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

clear all
clear global
tic;
global M_ oo_ options_ ys0_ ex0_ estimation_info
options_ = [];
M_.fname = 'MCF';
%
% Some global variables initialization
%
global_initialization;
diary off;
logname_ = 'MCF.log';
if exist(logname_, 'file')
    delete(logname_)
end
diary(logname_)
M_.exo_names = 'varepsilonm';
M_.exo_names_tex = 'varepsilonm';
M_.endo_names = 'lambda';
M_.endo_names_tex = 'lambda';
M_.endo_names = char(M_.endo_names, 'y');
M_.endo_names_tex = char(M_.endo_names_tex, 'y');
M_.endo_names = char(M_.endo_names, 'c');
M_.endo_names_tex = char(M_.endo_names_tex, 'c');
M_.endo_names = char(M_.endo_names, 'r');
M_.endo_names_tex = char(M_.endo_names_tex, 'r');
M_.endo_names = char(M_.endo_names, 'm');
M_.endo_names_tex = char(M_.endo_names_tex, 'm');
M_.endo_names = char(M_.endo_names, 'u');
M_.endo_names_tex = char(M_.endo_names_tex, 'u');
M_.endo_names = char(M_.endo_names, 'v');
M_.endo_names_tex = char(M_.endo_names_tex, 'v');
M_.endo_names = char(M_.endo_names, 'q');
M_.endo_names_tex = char(M_.endo_names_tex, 'q');
M_.endo_names = char(M_.endo_names, 's');
M_.endo_names_tex = char(M_.endo_names_tex, 's');
M_.endo_names = char(M_.endo_names, 'n');
M_.endo_names_tex = char(M_.endo_names_tex, 'n');
M_.endo_names = char(M_.endo_names, 'x');
M_.endo_names_tex = char(M_.endo_names_tex, 'x');
M_.endo_names = char(M_.endo_names, 'h');
M_.endo_names_tex = char(M_.endo_names_tex, 'h');
M_.endo_names = char(M_.endo_names, 'w');
M_.endo_names_tex = char(M_.endo_names_tex, 'w');
M_.endo_names = char(M_.endo_names, 'pi');
M_.endo_names_tex = char(M_.endo_names_tex, 'pi');
M_.endo_names = char(M_.endo_names, 'rn');
M_.endo_names_tex = char(M_.endo_names_tex, 'rn');
M_.endo_names = char(M_.endo_names, 'mpl');
M_.endo_names_tex = char(M_.endo_names_tex, 'mpl');
M_.endo_names = char(M_.endo_names, 'mrs');
M_.endo_names_tex = char(M_.endo_names_tex, 'mrs');
M_.endo_names = char(M_.endo_names, 'theta');
M_.endo_names_tex = char(M_.endo_names_tex, 'theta');
M_.param_names = 'e';
M_.param_names_tex = 'e';
M_.param_names = char(M_.param_names, 'phi');
M_.param_names_tex = char(M_.param_names_tex, 'phi');
M_.param_names = char(M_.param_names, 'beta');
M_.param_names_tex = char(M_.param_names_tex, 'beta');
M_.param_names = char(M_.param_names, 'rho');
M_.param_names_tex = char(M_.param_names_tex, 'rho');
M_.param_names = char(M_.param_names, 'kappa');
M_.param_names_tex = char(M_.param_names_tex, 'kappa');
M_.param_names = char(M_.param_names, 'eta');
M_.param_names_tex = char(M_.param_names_tex, 'eta');
M_.param_names = char(M_.param_names, 'sigma');
M_.param_names_tex = char(M_.param_names_tex, 'sigma');
M_.param_names = char(M_.param_names, 'rho_u');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_u');
M_.param_names = char(M_.param_names, 'rho_a');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_a');
M_.param_names = char(M_.param_names, 'varphi');
M_.param_names_tex = char(M_.param_names_tex, 'varphi');
M_.param_names = char(M_.param_names, 'alpha');
M_.param_names_tex = char(M_.param_names_tex, 'alpha');
M_.param_names = char(M_.param_names, 'rho_m');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_m');
M_.param_names = char(M_.param_names, 'gamma_pi');
M_.param_names_tex = char(M_.param_names_tex, 'gamma\_pi');
M_.param_names = char(M_.param_names, 'gamma_y');
M_.param_names_tex = char(M_.param_names_tex, 'gamma\_y');
M_.param_names = char(M_.param_names, 'F');
M_.param_names_tex = char(M_.param_names_tex, 'F');
M_.param_names = char(M_.param_names, 'S');
M_.param_names_tex = char(M_.param_names_tex, 'S');
M_.param_names = char(M_.param_names, 'b');
M_.param_names_tex = char(M_.param_names_tex, 'b');
M_.exo_det_nbr = 0;
M_.exo_nbr = 1;
M_.endo_nbr = 18;
M_.param_nbr = 17;
M_.orig_endo_nbr = 18;
M_.aux_vars = [];
M_.Sigma_e = zeros(1, 1);
M_.H = 0;
options_.linear = 1;
options_.block=0;
options_.bytecode=0;
options_.use_dll=0;
erase_compiled_function('MCF_static');
erase_compiled_function('MCF_dynamic');
M_.lead_lag_incidence = [
 0 5 23;
 0 6 0;
 1 7 24;
 2 8 0;
 3 9 0;
 0 10 0;
 0 11 0;
 0 12 25;
 0 13 0;
 4 14 0;
 0 15 26;
 0 16 27;
 0 17 28;
 0 18 29;
 0 19 0;
 0 20 0;
 0 21 0;
 0 22 0;]';
M_.nstatic = 8;
M_.nfwrd   = 6;
M_.npred   = 3;
M_.nboth   = 1;
M_.equations_tags = {
};
M_.exo_names_orig_ord = [1:1];
M_.maximum_lag = 1;
M_.maximum_lead = 1;
M_.maximum_endo_lag = 1;
M_.maximum_endo_lead = 1;
oo_.steady_state = zeros(18, 1);
M_.maximum_exo_lag = 0;
M_.maximum_exo_lead = 0;
oo_.exo_steady_state = zeros(1, 1);
M_.params = NaN(17, 1);
M_.NNZDerivatives = zeros(3, 1);
M_.NNZDerivatives(1) = 63;
M_.NNZDerivatives(2) = 3;
M_.NNZDerivatives(3) = -1;
M_.params( 1 ) = 0.7;
e = M_.params( 1 );
M_.params( 2 ) = 10;
phi = M_.params( 2 );
M_.params( 3 ) = 0.99;
beta = M_.params( 3 );
M_.params( 4 ) = 0.08;
rho = M_.params( 4 );
M_.params( 5 ) = 0.16;
kappa = M_.params( 5 );
M_.params( 6 ) = 0.5;
eta = M_.params( 6 );
M_.params( 7 ) = 0.5;
sigma = M_.params( 7 );
M_.params( 8 ) = 0.5;
rho_u = M_.params( 8 );
M_.params( 9 ) = 0.5;
rho_a = M_.params( 9 );
M_.params( 10 ) = 0.75;
varphi = M_.params( 10 );
M_.params( 11 ) = 0.99;
alpha = M_.params( 11 );
M_.params( 12 ) = 0.85;
rho_m = M_.params( 12 );
M_.params( 13 ) = 1.5;
gamma_pi = M_.params( 13 );
M_.params( 14 ) = 0.5;
gamma_y = M_.params( 14 );
M_.params( 15 ) = 0.844175;
F = M_.params( 15 );
M_.params( 16 ) = 3.6455;
S = M_.params( 16 );
M_.params( 17 ) = 0.62153;
b = M_.params( 17 );
%
% INITVAL instructions
%
options_.initval_file = 0;
oo_.steady_state( 2 ) = 1;
oo_.steady_state( 3 ) = 1;
oo_.steady_state( 12 ) = 1;
oo_.steady_state( 1 ) = 1;
oo_.steady_state( 5 ) = 0.0349;
oo_.steady_state( 8 ) = 0.5626;
oo_.steady_state( 9 ) = 0.3878;
oo_.steady_state( 18 ) = 0.9489;
oo_.steady_state( 10 ) = 0.9891;
oo_.steady_state( 6 ) = 0.09;
oo_.steady_state( 7 ) = 0.0854;
oo_.steady_state( 16 ) = 1.08794;
oo_.steady_state( 17 ) = 0.3333333333333333;
oo_.steady_state( 11 ) = 0.3064;
oo_.steady_state( 13 ) = 1.37066;
if M_.exo_nbr > 0;
	oo_.exo_simul = [ones(M_.maximum_lag,1)*oo_.exo_steady_state'];
end;
if M_.exo_det_nbr > 0;
	oo_.exo_det_simul = [ones(M_.maximum_lag,1)*oo_.exo_det_steady_state'];
end;
oo_.dr.eigval = check(M_,options_,oo_);
%
% SHOCKS instructions
%
make_ex_;
M_.exo_det_length = 0;
M_.Sigma_e(1, 1) = (-0.025);
M_.sigma_e_is_diagonal = 1;
options_.periods = 2100;
var_list_=[];
info = stoch_simul(var_list_);
save('MCF_results.mat', 'oo_', 'M_', 'options_');


disp(['Total computing time : ' dynsec2hms(toc) ]);
diary off