% % 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_ options_ = []; M_.fname = 'Jermann_Quadrini_FM'; % % Some global variables initialization % global_initialization; diary off; logname_ = 'Jermann_Quadrini_FM.log'; if exist(logname_, 'file') delete(logname_) end diary(logname_) M_.exo_names = 'eps_z'; M_.exo_names_tex = 'eps\_z'; M_.exo_names = char(M_.exo_names, 'eps_xi'); M_.exo_names_tex = char(M_.exo_names_tex, 'eps\_xi'); M_.endo_names = 'z'; M_.endo_names_tex = 'z'; M_.endo_names = char(M_.endo_names, 'xi'); M_.endo_names_tex = char(M_.endo_names_tex, 'xi'); M_.endo_names = char(M_.endo_names, 'k'); M_.endo_names_tex = char(M_.endo_names_tex, 'k'); M_.endo_names = char(M_.endo_names, 'b'); M_.endo_names_tex = char(M_.endo_names_tex, 'b'); 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, 'w'); M_.endo_names_tex = char(M_.endo_names_tex, 'w'); 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, 'd'); M_.endo_names_tex = char(M_.endo_names_tex, 'd'); 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, 'm_tilde'); M_.endo_names_tex = char(M_.endo_names_tex, 'm\_tilde'); 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, 'mu'); M_.endo_names_tex = char(M_.endo_names_tex, 'mu'); M_.endo_names = char(M_.endo_names, 'phi_d'); M_.endo_names_tex = char(M_.endo_names_tex, 'phi\_d'); M_.param_names = 'beta'; M_.param_names_tex = 'beta'; M_.param_names = char(M_.param_names, 'tau'); M_.param_names_tex = char(M_.param_names_tex, 'tau'); M_.param_names = char(M_.param_names, 'd_bar'); M_.param_names_tex = char(M_.param_names_tex, 'd\_bar'); 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, 'theta'); M_.param_names_tex = char(M_.param_names_tex, 'theta'); M_.param_names = char(M_.param_names, 'delta'); M_.param_names_tex = char(M_.param_names_tex, 'delta'); M_.param_names = char(M_.param_names, 'xi_bar'); M_.param_names_tex = char(M_.param_names_tex, 'xi\_bar'); 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, 'rho_z'); M_.param_names_tex = char(M_.param_names_tex, 'rho\_z'); M_.param_names = char(M_.param_names, 'sigma_z'); M_.param_names_tex = char(M_.param_names_tex, 'sigma\_z'); M_.param_names = char(M_.param_names, 'rho_xi'); M_.param_names_tex = char(M_.param_names_tex, 'rho\_xi'); M_.param_names = char(M_.param_names, 'sigma_xi'); M_.param_names_tex = char(M_.param_names_tex, 'sigma\_xi'); M_.exo_det_nbr = 0; M_.exo_nbr = 2; M_.endo_nbr = 13; M_.param_nbr = 12; M_.orig_endo_nbr = 13; M_.aux_vars = []; M_.Sigma_e = zeros(2, 2); M_.H = 0; options_.block=0; options_.bytecode=0; options_.use_dll=0; erase_compiled_function('Jermann_Quadrini_FM_dynamic'); M_.lead_lag_incidence = [ 1 5 18; 2 6 0; 3 7 0; 4 8 0; 0 9 19; 0 10 0; 0 11 20; 0 12 0; 0 13 0; 0 14 0; 0 15 0; 0 16 21; 0 17 22;]'; M_.equations_tags = { }; M_.exo_names_orig_ord = [1:2]; M_.maximum_lag = 1; M_.maximum_lead = 1; M_.maximum_endo_lag = 1; M_.maximum_endo_lead = 1; oo_.steady_state = zeros(13, 1); M_.maximum_exo_lag = 0; M_.maximum_exo_lead = 0; oo_.exo_steady_state = zeros(2, 1); M_.params = NaN(12, 1); M_.NNZDerivatives = zeros(3, 1); M_.NNZDerivatives(1) = 64; M_.NNZDerivatives(2) = -1; M_.NNZDerivatives(3) = -1; M_.params( 1 ) = 0.9825; beta = M_.params( 1 ); M_.params( 2 ) = 0.3500; tau = M_.params( 2 ); M_.params( 3 ) = 0.1100; d_bar = M_.params( 3 ); M_.params( 4 ) = 1.8834; alpha = M_.params( 4 ); M_.params( 5 ) = 0.3600; theta = M_.params( 5 ); M_.params( 6 ) = 0.0250; delta = M_.params( 6 ); M_.params( 7 ) = 0.1634; xi_bar = M_.params( 7 ); M_.params( 8 ) = 0.1460; kappa = M_.params( 8 ); M_.params( 9 ) = 0.9457; rho_z = M_.params( 9 ); M_.params( 10 ) = 0.0045; sigma_z = M_.params( 10 ); M_.params( 11 ) = 0.9703; rho_xi = M_.params( 11 ); M_.params( 12 ) = 0.0098; sigma_xi = M_.params( 12 ); % % INITVAL instructions % options_.initval_file = 0; oo_.steady_state( 1 ) = 0; oo_.steady_state( 2 ) = 0; oo_.steady_state( 3 ) = 2.302585092994046; oo_.steady_state( 4 ) = 1.386294361119891; oo_.steady_state( 5 ) = (-1.203972804325936); oo_.steady_state( 6 ) = 0.6931471805599453; oo_.steady_state( 7 ) = (-0.2231435513142097); oo_.steady_state( 8 ) = (-2.302585092994046); oo_.steady_state( 9 ) = 0; oo_.steady_state( 10 ) = log(M_.params(1)); oo_.steady_state( 11 ) = log(1/M_.params(1)); oo_.steady_state( 12 ) = (-3.218875824868201); oo_.steady_state( 13 ) = 1+2*M_.params(8)*((-2.302585092994046)-M_.params(3)); oo_.endo_simul=[oo_.steady_state*ones(1,M_.maximum_lag)]; 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_; M_.exo_det_length = 0; M_.Sigma_e(1, 1) = (1)^2; M_.Sigma_e(2, 2) = (1)^2; M_.sigma_e_is_diagonal = 1; steady; check; options_.irf = 30; options_.order = 1; var_list_=[]; var_list_ = 'b'; var_list_ = char(var_list_, 'n'); var_list_ = char(var_list_, 'y'); var_list_ = char(var_list_, 'd'); info = stoch_simul(var_list_); save('Jermann_Quadrini_FM_results.mat', 'oo_', 'M_', 'options_'); diary off disp(['Total computing time : ' dynsec2hms(toc) ]);