%
% 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 = 'JamTest1';
%
% Some global variables initialization
%
global_initialization;
diary off;
diary('JamTest1.log');
M_.exo_names = 'rer';
M_.exo_names_tex = 'rer';
M_.exo_names_long = 'rer';
M_.exo_names = char(M_.exo_names, 'pr_H');
M_.exo_names_tex = char(M_.exo_names_tex, 'pr\_H');
M_.exo_names_long = char(M_.exo_names_long, 'pr_H');
M_.exo_names = char(M_.exo_names, 'pr_O');
M_.exo_names_tex = char(M_.exo_names_tex, 'pr\_O');
M_.exo_names_long = char(M_.exo_names_long, 'pr_O');
M_.exo_names = char(M_.exo_names, 'i');
M_.exo_names_tex = char(M_.exo_names_tex, 'i');
M_.exo_names_long = char(M_.exo_names_long, 'i');
M_.exo_names = char(M_.exo_names, 'ppi');
M_.exo_names_tex = char(M_.exo_names_tex, 'ppi');
M_.exo_names_long = char(M_.exo_names_long, 'ppi');
M_.endo_names = 'c_H';
M_.endo_names_tex = 'c\_H';
M_.endo_names_long = 'c_H';
M_.endo_names = char(M_.endo_names, 'c_F');
M_.endo_names_tex = char(M_.endo_names_tex, 'c\_F');
M_.endo_names_long = char(M_.endo_names_long, 'c_F');
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, 'o_C');
M_.endo_names_tex = char(M_.endo_names_tex, 'o\_C');
M_.endo_names_long = char(M_.endo_names_long, 'o_C');
M_.param_names = 'gamma';
M_.param_names_tex = 'gamma';
M_.param_names_long = 'gamma';
M_.param_names = char(M_.param_names, 'theta');
M_.param_names_tex = char(M_.param_names_tex, 'theta');
M_.param_names_long = char(M_.param_names_long, 'theta');
M_.param_names = char(M_.param_names, 'eta');
M_.param_names_tex = char(M_.param_names_tex, 'eta');
M_.param_names_long = char(M_.param_names_long, 'eta');
M_.param_names = char(M_.param_names, 'h');
M_.param_names_tex = char(M_.param_names_tex, 'h');
M_.param_names_long = char(M_.param_names_long, 'h');
M_.exo_det_nbr = 0;
M_.exo_nbr = 5;
M_.endo_nbr = 4;
M_.param_nbr = 4;
M_.orig_endo_nbr = 4;
M_.aux_vars = [];
M_.Sigma_e = zeros(5, 5);
M_.Correlation_matrix = eye(5, 5);
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('JamTest1_static');
erase_compiled_function('JamTest1_dynamic');
M_.lead_lag_incidence = [
 0 2 0;
 0 3 0;
 1 4 6;
 0 5 0;]';
M_.nstatic = 3;
M_.nfwrd   = 0;
M_.npred   = 0;
M_.nboth   = 1;
M_.nsfwrd   = 1;
M_.nspred   = 1;
M_.ndynamic   = 1;
M_.equations_tags = {
};
M_.static_and_dynamic_models_differ = 0;
M_.exo_names_orig_ord = [1:5];
M_.maximum_lag = 1;
M_.maximum_lead = 1;
M_.maximum_endo_lag = 1;
M_.maximum_endo_lead = 1;
oo_.steady_state = zeros(4, 1);
M_.maximum_exo_lag = 0;
M_.maximum_exo_lead = 1;
oo_.exo_steady_state = zeros(5, 1);
M_.params = NaN(4, 1);
M_.NNZDerivatives = zeros(3, 1);
M_.NNZDerivatives(1) = 16;
M_.NNZDerivatives(2) = 0;
M_.NNZDerivatives(3) = -1;
M_.params( 1 ) = .4;
gamma = M_.params( 1 );
M_.params( 2 ) = 1;
theta = M_.params( 2 );
M_.params( 3 ) = .15;
eta = M_.params( 3 );
M_.params( 4 ) = .5;
h = M_.params( 4 );
end;
%
% INITVAL instructions
%
options_.initval_file = 0;
oo_.steady_state( 1 ) = 0;
oo_.steady_state( 2 ) = 0;
oo_.steady_state( 3 ) = 0;
oo_.steady_state( 4 ) = 0;
oo_.exo_steady_state( 1 ) = 0;
oo_.exo_steady_state( 2 ) = 0;
oo_.exo_steady_state( 3 ) = 0;
oo_.exo_steady_state( 4 ) = 0;
oo_.exo_steady_state( 5 ) = 0;
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;
%
% SHOCKS instructions
%
make_ex_;
set_shocks(0,1:9, 1, 0.1);
set_shocks(0,1:9, 2, 0.1);
set_shocks(0,1:9, 3, 0.1);
set_shocks(0,1:9, 4, 0.1);
set_shocks(0,1:9, 5, 0.1);
M_.exo_det_length = 0;
steady;
oo_.dr.eigval = check(M_,options_,oo_);
options_.periods = 2100;
simul();
save('JamTest1_results.mat', 'oo_', 'M_', 'options_');
if exist('estim_params_', 'var') == 1
  save('JamTest1_results.mat', 'estim_params_', '-append');
end
if exist('bayestopt_', 'var') == 1
  save('JamTest1_results.mat', 'bayestopt_', '-append');
end
if exist('dataset_', 'var') == 1
  save('JamTest1_results.mat', 'dataset_', '-append');
end
if exist('estimation_info', 'var') == 1
  save('JamTest1_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