%% Specify the model name
model = 'model3';
num_p = '[0 1]';
den_p = '[0.1 1]';
set_param('model3/Transfer', 'Numerator', num_p);
set_param('model3/Transfer', 'Denominator',  den_p);

%% Specify the analysis I/Os
% Get the analysis I/Os from the model
io = getlinio(model);

%% Specify the operating point
% Use the model initial condition
op = operpoint(model);


%% Linearize the model
sys = linearize(model,io,op);

%% Plot the resulting linearization
step(sys)