
%--------------------------------------------------------------------------
% Simulate Model                            
%--------------------------------------------------------------------------

dynare benchmark

% Change parameter value
init_val = 2;
M_.params(strcmp('nu_p',M_.param_names)) = init_val;

% Compute steady state with new value
exo = zeros(M_.exo_nbr,1); ys = zeros(M_.endo_nbr,1);
oo_.dr.ys = benchmark_steadystate(ys,exo,M_,options_);

Steady = exp(oo_.dr.ys);

