%set_param_value('A_ss',0.5);

dynare Basu_basic_EZ.mod;
%%
options_.k_order_solver=1;

SimLength=20000;

set_param_value('beta',0.99);
rng(1)
shocks=randn(SimLength,2);
[oo_.dr,info,M_,options_,oo_] = resol(0,M_,options_,oo_); %run model solution in Dynare


SimData = simult_(oo_.steady_state,oo_.dr,shocks,3);



%% Construct Moments
endo_names=cellstr(M_.endo_names);
%names={'Y_t','C_t','I_t','N_t','prod_t','w_t','R_t','A_t'};

means=mean(SimData(:,ceil(SimLength/2):end),2);

sim_mom=[means(strcmp('r',endo_names))]