%
% Status : main Dynare file 
%
% Warning : this file is generated automatically by Dynare
%           from model file (.mod)

clear all
tic;
global M_ oo_ options_ ys0_ ex0_ estimation_info
options_ = [];
M_.fname = 'rbcb_monprud_nobias';
%
% Some global variables initialization
%
global_initialization;
diary off;
diary('rbcb_monprud_nobias.log');
M_.exo_names = 'errz';
M_.exo_names_tex = 'errz';
M_.exo_names_long = 'errz';
M_.exo_names = char(M_.exo_names, 'errmon');
M_.exo_names_tex = char(M_.exo_names_tex, 'errmon');
M_.exo_names_long = char(M_.exo_names_long, 'errmon');
M_.exo_names = char(M_.exo_names, 'errlambda');
M_.exo_names_tex = char(M_.exo_names_tex, 'errlambda');
M_.exo_names_long = char(M_.exo_names_long, 'errlambda');
M_.exo_names = char(M_.exo_names, 'errbt');
M_.exo_names_tex = char(M_.exo_names_tex, 'errbt');
M_.exo_names_long = char(M_.exo_names_long, 'errbt');
M_.exo_names = char(M_.exo_names, 'errbetavar');
M_.exo_names_tex = char(M_.exo_names_tex, 'errbetavar');
M_.exo_names_long = char(M_.exo_names_long, 'errbetavar');
M_.exo_names = char(M_.exo_names, 'errcpref');
M_.exo_names_tex = char(M_.exo_names_tex, 'errcpref');
M_.exo_names_long = char(M_.exo_names_long, 'errcpref');
M_.exo_names = char(M_.exo_names, 'errmarkup');
M_.exo_names_tex = char(M_.exo_names_tex, 'errmarkup');
M_.exo_names_long = char(M_.exo_names_long, 'errmarkup');
M_.endo_names = 'k';
M_.endo_names_tex = 'k';
M_.endo_names_long = 'k';
M_.endo_names = char(M_.endo_names, 'y');
M_.endo_names_tex = char(M_.endo_names_tex, 'y');
M_.endo_names_long = char(M_.endo_names_long, 'y');
M_.endo_names = char(M_.endo_names, 'w');
M_.endo_names_tex = char(M_.endo_names_tex, 'w');
M_.endo_names_long = char(M_.endo_names_long, 'w');
M_.endo_names = char(M_.endo_names, 'rs');
M_.endo_names_tex = char(M_.endo_names_tex, 'rs');
M_.endo_names_long = char(M_.endo_names_long, 'rs');
M_.endo_names = char(M_.endo_names, 'lambdac');
M_.endo_names_tex = char(M_.endo_names_tex, 'lambdac');
M_.endo_names_long = char(M_.endo_names_long, 'lambdac');
M_.endo_names = char(M_.endo_names, 'l');
M_.endo_names_tex = char(M_.endo_names_tex, 'l');
M_.endo_names_long = char(M_.endo_names_long, 'l');
M_.endo_names = char(M_.endo_names, 'r');
M_.endo_names_tex = char(M_.endo_names_tex, 'r');
M_.endo_names_long = char(M_.endo_names_long, 'r');
M_.endo_names = char(M_.endo_names, 's');
M_.endo_names_tex = char(M_.endo_names_tex, 's');
M_.endo_names_long = char(M_.endo_names_long, 's');
M_.endo_names = char(M_.endo_names, 'n');
M_.endo_names_tex = char(M_.endo_names_tex, 'n');
M_.endo_names_long = char(M_.endo_names_long, 'n');
M_.endo_names = char(M_.endo_names, 'v');
M_.endo_names_tex = char(M_.endo_names_tex, 'v');
M_.endo_names_long = char(M_.endo_names_long, 'v');
M_.endo_names = char(M_.endo_names, 'eta');
M_.endo_names_tex = char(M_.endo_names_tex, 'eta');
M_.endo_names_long = char(M_.endo_names_long, 'eta');
M_.endo_names = char(M_.endo_names, 'in');
M_.endo_names_tex = char(M_.endo_names_tex, 'in');
M_.endo_names_long = char(M_.endo_names_long, 'in');
M_.endo_names = char(M_.endo_names, 'ig');
M_.endo_names_tex = char(M_.endo_names_tex, 'ig');
M_.endo_names_long = char(M_.endo_names_long, 'ig');
M_.endo_names = char(M_.endo_names, 'q');
M_.endo_names_tex = char(M_.endo_names_tex, 'q');
M_.endo_names_long = char(M_.endo_names_long, 'q');
M_.endo_names = char(M_.endo_names, 'c');
M_.endo_names_tex = char(M_.endo_names_tex, 'c');
M_.endo_names_long = char(M_.endo_names_long, 'c');
M_.endo_names = char(M_.endo_names, 'gc');
M_.endo_names_tex = char(M_.endo_names_tex, 'gc');
M_.endo_names_long = char(M_.endo_names_long, 'gc');
M_.endo_names = char(M_.endo_names, 'spread');
M_.endo_names_tex = char(M_.endo_names_tex, 'spread');
M_.endo_names_long = char(M_.endo_names_long, 'spread');
M_.endo_names = char(M_.endo_names, 'l2n');
M_.endo_names_tex = char(M_.endo_names_tex, 'l2n');
M_.endo_names_long = char(M_.endo_names_long, 'l2n');
M_.endo_names = char(M_.endo_names, 'rmon');
M_.endo_names_tex = char(M_.endo_names_tex, 'rmon');
M_.endo_names_long = char(M_.endo_names_long, 'rmon');
M_.endo_names = char(M_.endo_names, 'infl');
M_.endo_names_tex = char(M_.endo_names_tex, 'infl');
M_.endo_names_long = char(M_.endo_names_long, 'infl');
M_.endo_names = char(M_.endo_names, 'phi');
M_.endo_names_tex = char(M_.endo_names_tex, 'phi');
M_.endo_names_long = char(M_.endo_names_long, 'phi');
M_.endo_names = char(M_.endo_names, 'dphidp');
M_.endo_names_tex = char(M_.endo_names_tex, 'dphidp');
M_.endo_names_long = char(M_.endo_names_long, 'dphidp');
M_.endo_names = char(M_.endo_names, 'dphidlp');
M_.endo_names_tex = char(M_.endo_names_tex, 'dphidlp');
M_.endo_names_long = char(M_.endo_names_long, 'dphidlp');
M_.endo_names = char(M_.endo_names, 'sigma');
M_.endo_names_tex = char(M_.endo_names_tex, 'sigma');
M_.endo_names_long = char(M_.endo_names_long, 'sigma');
M_.endo_names = char(M_.endo_names, 'n2y');
M_.endo_names_tex = char(M_.endo_names_tex, 'n2y');
M_.endo_names_long = char(M_.endo_names_long, 'n2y');
M_.endo_names = char(M_.endo_names, 'Util1');
M_.endo_names_tex = char(M_.endo_names_tex, 'Util1');
M_.endo_names_long = char(M_.endo_names_long, 'Util1');
M_.endo_names = char(M_.endo_names, 'Util2');
M_.endo_names_tex = char(M_.endo_names_tex, 'Util2');
M_.endo_names_long = char(M_.endo_names_long, 'Util2');
M_.endo_names = char(M_.endo_names, 'bt');
M_.endo_names_tex = char(M_.endo_names_tex, 'bt');
M_.endo_names_long = char(M_.endo_names_long, 'bt');
M_.endo_names = char(M_.endo_names, 'z');
M_.endo_names_tex = char(M_.endo_names_tex, 'z');
M_.endo_names_long = char(M_.endo_names_long, 'z');
M_.endo_names = char(M_.endo_names, 'lambda');
M_.endo_names_tex = char(M_.endo_names_tex, 'lambda');
M_.endo_names_long = char(M_.endo_names_long, 'lambda');
M_.endo_names = char(M_.endo_names, 'betavar');
M_.endo_names_tex = char(M_.endo_names_tex, 'betavar');
M_.endo_names_long = char(M_.endo_names_long, 'betavar');
M_.endo_names = char(M_.endo_names, 'cpref');
M_.endo_names_tex = char(M_.endo_names_tex, 'cpref');
M_.endo_names_long = char(M_.endo_names_long, 'cpref');
M_.endo_names = char(M_.endo_names, 'markup');
M_.endo_names_tex = char(M_.endo_names_tex, 'markup');
M_.endo_names_long = char(M_.endo_names_long, 'markup');
M_.param_names = 'alphap';
M_.param_names_tex = 'alphap';
M_.param_names_long = 'alphap';
M_.param_names = char(M_.param_names, 'gammap');
M_.param_names_tex = char(M_.param_names_tex, 'gammap');
M_.param_names_long = char(M_.param_names_long, 'gammap');
M_.param_names = char(M_.param_names, 'chi0p');
M_.param_names_tex = char(M_.param_names_tex, 'chi0p');
M_.param_names_long = char(M_.param_names_long, 'chi0p');
M_.param_names = char(M_.param_names, 'chip');
M_.param_names_tex = char(M_.param_names_tex, 'chip');
M_.param_names_long = char(M_.param_names_long, 'chip');
M_.param_names = char(M_.param_names, 'betap');
M_.param_names_tex = char(M_.param_names_tex, 'betap');
M_.param_names_long = char(M_.param_names_long, 'betap');
M_.param_names = char(M_.param_names, 'deltap');
M_.param_names_tex = char(M_.param_names_tex, 'deltap');
M_.param_names_long = char(M_.param_names_long, 'deltap');
M_.param_names = char(M_.param_names, 'psip');
M_.param_names_tex = char(M_.param_names_tex, 'psip');
M_.param_names_long = char(M_.param_names_long, 'psip');
M_.param_names = char(M_.param_names, 'omegagyp');
M_.param_names_tex = char(M_.param_names_tex, 'omegagyp');
M_.param_names_long = char(M_.param_names_long, 'omegagyp');
M_.param_names = char(M_.param_names, 'lambdap');
M_.param_names_tex = char(M_.param_names_tex, 'lambdap');
M_.param_names_long = char(M_.param_names_long, 'lambdap');
M_.param_names = char(M_.param_names, 'thetap');
M_.param_names_tex = char(M_.param_names_tex, 'thetap');
M_.param_names_long = char(M_.param_names_long, 'thetap');
M_.param_names = char(M_.param_names, 'baromegap');
M_.param_names_tex = char(M_.param_names_tex, 'baromegap');
M_.param_names_long = char(M_.param_names_long, 'baromegap');
M_.param_names = char(M_.param_names, 's2n');
M_.param_names_tex = char(M_.param_names_tex, 's2n');
M_.param_names_long = char(M_.param_names_long, 's2n');
M_.param_names = char(M_.param_names, 'phip');
M_.param_names_tex = char(M_.param_names_tex, 'phip');
M_.param_names_long = char(M_.param_names_long, 'phip');
M_.param_names = char(M_.param_names, 'spreadp');
M_.param_names_tex = char(M_.param_names_tex, 'spreadp');
M_.param_names_long = char(M_.param_names_long, 'spreadp');
M_.param_names = char(M_.param_names, 'nup');
M_.param_names_tex = char(M_.param_names_tex, 'nup');
M_.param_names_long = char(M_.param_names_long, 'nup');
M_.param_names = char(M_.param_names, 'taup');
M_.param_names_tex = char(M_.param_names_tex, 'taup');
M_.param_names_long = char(M_.param_names_long, 'taup');
M_.param_names = char(M_.param_names, 'indexp');
M_.param_names_tex = char(M_.param_names_tex, 'indexp');
M_.param_names_long = char(M_.param_names_long, 'indexp');
M_.param_names = char(M_.param_names, 'mu_mpr');
M_.param_names_tex = char(M_.param_names_tex, 'mu\_mpr');
M_.param_names_long = char(M_.param_names_long, 'mu_mpr');
M_.param_names = char(M_.param_names, 'mu_cb');
M_.param_names_tex = char(M_.param_names_tex, 'mu\_cb');
M_.param_names_long = char(M_.param_names_long, 'mu_cb');
M_.param_names = char(M_.param_names, 'xip');
M_.param_names_tex = char(M_.param_names_tex, 'xip');
M_.param_names_long = char(M_.param_names_long, 'xip');
M_.param_names = char(M_.param_names, 'blife');
M_.param_names_tex = char(M_.param_names_tex, 'blife');
M_.param_names_long = char(M_.param_names_long, 'blife');
M_.param_names = char(M_.param_names, 'rs_rslp');
M_.param_names_tex = char(M_.param_names_tex, 'rs\_rslp');
M_.param_names_long = char(M_.param_names_long, 'rs_rslp');
M_.param_names = char(M_.param_names, 'rs_inflp');
M_.param_names_tex = char(M_.param_names_tex, 'rs\_inflp');
M_.param_names_long = char(M_.param_names_long, 'rs_inflp');
M_.param_names = char(M_.param_names, 'rs_yp');
M_.param_names_tex = char(M_.param_names_tex, 'rs\_yp');
M_.param_names_long = char(M_.param_names_long, 'rs_yp');
M_.param_names = char(M_.param_names, 'bt_btlp');
M_.param_names_tex = char(M_.param_names_tex, 'bt\_btlp');
M_.param_names_long = char(M_.param_names_long, 'bt_btlp');
M_.param_names = char(M_.param_names, 'bt_sp');
M_.param_names_tex = char(M_.param_names_tex, 'bt\_sp');
M_.param_names_long = char(M_.param_names_long, 'bt_sp');
M_.param_names = char(M_.param_names, 'rhozp');
M_.param_names_tex = char(M_.param_names_tex, 'rhozp');
M_.param_names_long = char(M_.param_names_long, 'rhozp');
M_.param_names = char(M_.param_names, 'sigmazp');
M_.param_names_tex = char(M_.param_names_tex, 'sigmazp');
M_.param_names_long = char(M_.param_names_long, 'sigmazp');
M_.param_names = char(M_.param_names, 'rholambdap');
M_.param_names_tex = char(M_.param_names_tex, 'rholambdap');
M_.param_names_long = char(M_.param_names_long, 'rholambdap');
M_.param_names = char(M_.param_names, 'sigmalambdap');
M_.param_names_tex = char(M_.param_names_tex, 'sigmalambdap');
M_.param_names_long = char(M_.param_names_long, 'sigmalambdap');
M_.param_names = char(M_.param_names, 'rhobtp');
M_.param_names_tex = char(M_.param_names_tex, 'rhobtp');
M_.param_names_long = char(M_.param_names_long, 'rhobtp');
M_.param_names = char(M_.param_names, 'sigmabtp');
M_.param_names_tex = char(M_.param_names_tex, 'sigmabtp');
M_.param_names_long = char(M_.param_names_long, 'sigmabtp');
M_.param_names = char(M_.param_names, 'rhobetavarp');
M_.param_names_tex = char(M_.param_names_tex, 'rhobetavarp');
M_.param_names_long = char(M_.param_names_long, 'rhobetavarp');
M_.param_names = char(M_.param_names, 'sigmabetavarp');
M_.param_names_tex = char(M_.param_names_tex, 'sigmabetavarp');
M_.param_names_long = char(M_.param_names_long, 'sigmabetavarp');
M_.param_names = char(M_.param_names, 'rhocprefp');
M_.param_names_tex = char(M_.param_names_tex, 'rhocprefp');
M_.param_names_long = char(M_.param_names_long, 'rhocprefp');
M_.param_names = char(M_.param_names, 'sigmacprefp');
M_.param_names_tex = char(M_.param_names_tex, 'sigmacprefp');
M_.param_names_long = char(M_.param_names_long, 'sigmacprefp');
M_.param_names = char(M_.param_names, 'l_ss');
M_.param_names_tex = char(M_.param_names_tex, 'l\_ss');
M_.param_names_long = char(M_.param_names_long, 'l_ss');
M_.param_names = char(M_.param_names, 'q_ss');
M_.param_names_tex = char(M_.param_names_tex, 'q\_ss');
M_.param_names_long = char(M_.param_names_long, 'q_ss');
M_.param_names = char(M_.param_names, 'r_ss');
M_.param_names_tex = char(M_.param_names_tex, 'r\_ss');
M_.param_names_long = char(M_.param_names_long, 'r_ss');
M_.param_names = char(M_.param_names, 'y_ss');
M_.param_names_tex = char(M_.param_names_tex, 'y\_ss');
M_.param_names_long = char(M_.param_names_long, 'y_ss');
M_.param_names = char(M_.param_names, 'w_ss');
M_.param_names_tex = char(M_.param_names_tex, 'w\_ss');
M_.param_names_long = char(M_.param_names_long, 'w_ss');
M_.param_names = char(M_.param_names, 'k_ss');
M_.param_names_tex = char(M_.param_names_tex, 'k\_ss');
M_.param_names_long = char(M_.param_names_long, 'k_ss');
M_.param_names = char(M_.param_names, 's_ss');
M_.param_names_tex = char(M_.param_names_tex, 's\_ss');
M_.param_names_long = char(M_.param_names_long, 's_ss');
M_.param_names = char(M_.param_names, 'n_ss');
M_.param_names_tex = char(M_.param_names_tex, 'n\_ss');
M_.param_names_long = char(M_.param_names_long, 'n_ss');
M_.param_names = char(M_.param_names, 'in_ss');
M_.param_names_tex = char(M_.param_names_tex, 'in\_ss');
M_.param_names_long = char(M_.param_names_long, 'in_ss');
M_.param_names = char(M_.param_names, 'ig_ss');
M_.param_names_tex = char(M_.param_names_tex, 'ig\_ss');
M_.param_names_long = char(M_.param_names_long, 'ig_ss');
M_.param_names = char(M_.param_names, 'gc_ss');
M_.param_names_tex = char(M_.param_names_tex, 'gc\_ss');
M_.param_names_long = char(M_.param_names_long, 'gc_ss');
M_.param_names = char(M_.param_names, 'infl_ss');
M_.param_names_tex = char(M_.param_names_tex, 'infl\_ss');
M_.param_names_long = char(M_.param_names_long, 'infl_ss');
M_.param_names = char(M_.param_names, 'sigmamarkup');
M_.param_names_tex = char(M_.param_names_tex, 'sigmamarkup');
M_.param_names_long = char(M_.param_names_long, 'sigmamarkup');
M_.param_names = char(M_.param_names, 'rhomarkup');
M_.param_names_tex = char(M_.param_names_tex, 'rhomarkup');
M_.param_names_long = char(M_.param_names_long, 'rhomarkup');
M_.exo_det_nbr = 0;
M_.exo_nbr = 7;
M_.endo_nbr = 33;
M_.param_nbr = 50;
M_.orig_endo_nbr = 33;
M_.aux_vars = [];
M_.Sigma_e = zeros(7, 7);
M_.Correlation_matrix = eye(7, 7);
M_.H = 0;
M_.Correlation_matrix_ME = 1;
M_.sigma_e_is_diagonal = 1;
options_.block=0;
options_.bytecode=0;
options_.use_dll=0;
erase_compiled_function('rbcb_monprud_nobias_static');
erase_compiled_function('rbcb_monprud_nobias_dynamic');
M_.lead_lag_incidence = [
 1 19 0;
 2 20 52;
 0 21 0;
 0 22 53;
 0 23 54;
 0 24 0;
 3 25 0;
 4 26 0;
 5 27 0;
 6 28 55;
 7 29 56;
 0 30 0;
 8 31 57;
 9 32 58;
 10 33 59;
 0 34 0;
 0 35 0;
 0 36 0;
 11 37 0;
 12 38 60;
 0 39 0;
 0 40 0;
 0 41 61;
 0 42 0;
 0 43 0;
 0 44 0;
 0 45 0;
 13 46 62;
 14 47 0;
 15 48 0;
 16 49 0;
 17 50 63;
 18 51 0;]';
M_.nstatic = 12;
M_.nfwrd   = 3;
M_.npred   = 9;
M_.nboth   = 9;
M_.nsfwrd   = 12;
M_.nspred   = 18;
M_.ndynamic   = 21;
M_.equations_tags = {
};
M_.static_and_dynamic_models_differ = 0;
M_.exo_names_orig_ord = [1:7];
M_.maximum_lag = 1;
M_.maximum_lead = 1;
M_.maximum_endo_lag = 1;
M_.maximum_endo_lead = 1;
oo_.steady_state = zeros(33, 1);
M_.maximum_exo_lag = 0;
M_.maximum_exo_lead = 0;
oo_.exo_steady_state = zeros(7, 1);
M_.params = NaN(50, 1);
M_.NNZDerivatives = zeros(3, 1);
M_.NNZDerivatives(1) = 151;
M_.NNZDerivatives(2) = 374;
M_.NNZDerivatives(3) = -1;
%
% SHOCKS instructions
%
make_ex_;
M_.exo_det_length = 0;
M_.Sigma_e(1, 1) = (1)^2;
M_.Sigma_e(7, 7) = (1)^2;
options_.irf = 0;
options_.nomoments = 1;
options_.order = 2;
options_.periods = 0;
var_list_=[];
info = stoch_simul(var_list_);
save('rbcb_monprud_nobias_results.mat', 'oo_', 'M_', 'options_');
if exist('estim_params_', 'var') == 1
  save('rbcb_monprud_nobias_results.mat', 'estim_params_', '-append');
end
if exist('bayestopt_', 'var') == 1
  save('rbcb_monprud_nobias_results.mat', 'bayestopt_', '-append');
end
if exist('dataset_', 'var') == 1
  save('rbcb_monprud_nobias_results.mat', 'dataset_', '-append');
end
if exist('estimation_info', 'var') == 1
  save('rbcb_monprud_nobias_results.mat', 'estimation_info', '-append');
end


disp(['Total computing time : ' dynsec2hms(toc) ]);
if ~isempty(lastwarn)
  disp('Note: warning(s) encountered in MATLAB/Octave code')
end
diary off
