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

if isoctave || matlab_ver_less_than('8.6')
    clear all
else
    clearvars -global
    clear_persistent_variables(fileparts(which('dynare')), false)
end
tic0 = tic;
% Save empty dates and dseries objects in memory.
dates('initialize');
dseries('initialize');
% Define global variables.
global M_ options_ oo_ estim_params_ bayestopt_ dataset_ dataset_info estimation_info ys0_ ex0_
options_ = [];
M_.fname = 'rep_lishu';
M_.dynare_version = '4.5.4';
oo_.dynare_version = '4.5.4';
options_.dynare_version = '4.5.4';
%
% Some global variables initialization
%
global_initialization;
diary off;
diary('rep_lishu.log');
M_.exo_names = 'eps_mp';
M_.exo_names_tex = 'eps\_mp';
M_.exo_names_long = 'eps_mp';
M_.exo_names = char(M_.exo_names, 'e_a');
M_.exo_names_tex = char(M_.exo_names_tex, 'e\_a');
M_.exo_names_long = char(M_.exo_names_long, 'e_a');
M_.exo_names = char(M_.exo_names, 'e_wg');
M_.exo_names_tex = char(M_.exo_names_tex, 'e\_wg');
M_.exo_names_long = char(M_.exo_names_long, 'e_wg');
M_.endo_names = 'pie';
M_.endo_names_tex = 'pie';
M_.endo_names_long = 'pie';
M_.endo_names = char(M_.endo_names, 'r_kh');
M_.endo_names_tex = char(M_.endo_names_tex, 'r\_kh');
M_.endo_names_long = char(M_.endo_names_long, 'r_kh');
M_.endo_names = char(M_.endo_names, 'u_h');
M_.endo_names_tex = char(M_.endo_names_tex, 'u\_h');
M_.endo_names_long = char(M_.endo_names_long, 'u_h');
M_.endo_names = char(M_.endo_names, 'R_KH');
M_.endo_names_tex = char(M_.endo_names_tex, 'R\_KH');
M_.endo_names_long = char(M_.endo_names_long, 'R_KH');
M_.endo_names = char(M_.endo_names, 'omegarbar');
M_.endo_names_tex = char(M_.endo_names_tex, 'omegarbar');
M_.endo_names_long = char(M_.endo_names_long, 'omegarbar');
M_.endo_names = char(M_.endo_names, 'k_h');
M_.endo_names_tex = char(M_.endo_names_tex, 'k\_h');
M_.endo_names_long = char(M_.endo_names_long, 'k_h');
M_.endo_names = char(M_.endo_names, 'n_h');
M_.endo_names_tex = char(M_.endo_names_tex, 'n\_h');
M_.endo_names_long = char(M_.endo_names_long, 'n_h');
M_.endo_names = char(M_.endo_names, 'b_h');
M_.endo_names_tex = char(M_.endo_names_tex, 'b\_h');
M_.endo_names_long = char(M_.endo_names_long, 'b_h');
M_.endo_names = char(M_.endo_names, 'omegarbar_b');
M_.endo_names_tex = char(M_.endo_names_tex, 'omegarbar\_b');
M_.endo_names_long = char(M_.endo_names_long, 'omegarbar_b');
M_.endo_names = char(M_.endo_names, 'chi_sb');
M_.endo_names_tex = char(M_.endo_names_tex, 'chi\_sb');
M_.endo_names_long = char(M_.endo_names_long, 'chi_sb');
M_.endo_names = char(M_.endo_names, 'n_sb');
M_.endo_names_tex = char(M_.endo_names_tex, 'n\_sb');
M_.endo_names_long = char(M_.endo_names_long, 'n_sb');
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, 'kap_sb');
M_.endo_names_tex = char(M_.endo_names_tex, 'kap\_sb');
M_.endo_names_long = char(M_.endo_names_long, 'kap_sb');
M_.endo_names = char(M_.endo_names, 'P_efp');
M_.endo_names_tex = char(M_.endo_names_tex, 'P\_efp');
M_.endo_names_long = char(M_.endo_names_long, 'P_efp');
M_.endo_names = char(M_.endo_names, 'R_hsb');
M_.endo_names_tex = char(M_.endo_names_tex, 'R\_hsb');
M_.endo_names_long = char(M_.endo_names_long, 'R_hsb');
M_.endo_names = char(M_.endo_names, 'r_kl');
M_.endo_names_tex = char(M_.endo_names_tex, 'r\_kl');
M_.endo_names_long = char(M_.endo_names_long, 'r_kl');
M_.endo_names = char(M_.endo_names, 'R_KL');
M_.endo_names_tex = char(M_.endo_names_tex, 'R\_KL');
M_.endo_names_long = char(M_.endo_names_long, 'R_KL');
M_.endo_names = char(M_.endo_names, 'u_l');
M_.endo_names_tex = char(M_.endo_names_tex, 'u\_l');
M_.endo_names_long = char(M_.endo_names_long, 'u_l');
M_.endo_names = char(M_.endo_names, 'k_l');
M_.endo_names_tex = char(M_.endo_names_tex, 'k\_l');
M_.endo_names_long = char(M_.endo_names_long, 'k_l');
M_.endo_names = char(M_.endo_names, 'n_l');
M_.endo_names_tex = char(M_.endo_names_tex, 'n\_l');
M_.endo_names_long = char(M_.endo_names_long, 'n_l');
M_.endo_names = char(M_.endo_names, 'b_l');
M_.endo_names_tex = char(M_.endo_names_tex, 'b\_l');
M_.endo_names_long = char(M_.endo_names_long, 'b_l');
M_.endo_names = char(M_.endo_names, 'r_bl');
M_.endo_names_tex = char(M_.endo_names_tex, 'r\_bl');
M_.endo_names_long = char(M_.endo_names_long, 'r_bl');
M_.endo_names = char(M_.endo_names, 'r_e');
M_.endo_names_tex = char(M_.endo_names_tex, 'r\_e');
M_.endo_names_long = char(M_.endo_names_long, 'r_e');
M_.endo_names = char(M_.endo_names, 'profit_b');
M_.endo_names_tex = char(M_.endo_names_tex, 'profit\_b');
M_.endo_names_long = char(M_.endo_names_long, 'profit_b');
M_.endo_names = char(M_.endo_names, 'chi_rb');
M_.endo_names_tex = char(M_.endo_names_tex, 'chi\_rb');
M_.endo_names_long = char(M_.endo_names_long, 'chi_rb');
M_.endo_names = char(M_.endo_names, 'eps_lop');
M_.endo_names_tex = char(M_.endo_names_tex, 'eps\_lop');
M_.endo_names_long = char(M_.endo_names_long, 'eps_lop');
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, 'r_l');
M_.endo_names_tex = char(M_.endo_names_tex, 'r\_l');
M_.endo_names_long = char(M_.endo_names_long, 'r_l');
M_.endo_names = char(M_.endo_names, 'a');
M_.endo_names_tex = char(M_.endo_names_tex, 'a');
M_.endo_names_long = char(M_.endo_names_long, 'a');
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, 'k');
M_.endo_names_tex = char(M_.endo_names_tex, 'k');
M_.endo_names_long = char(M_.endo_names_long, 'k');
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, '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, 'i');
M_.endo_names_tex = char(M_.endo_names_tex, 'i');
M_.endo_names_long = char(M_.endo_names_long, 'i');
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, 'Fp');
M_.endo_names_tex = char(M_.endo_names_tex, 'Fp');
M_.endo_names_long = char(M_.endo_names_long, 'Fp');
M_.endo_names = char(M_.endo_names, 'Kp');
M_.endo_names_tex = char(M_.endo_names_tex, 'Kp');
M_.endo_names_long = char(M_.endo_names_long, 'Kp');
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, '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, 'D');
M_.endo_names_tex = char(M_.endo_names_tex, 'D');
M_.endo_names_long = char(M_.endo_names_long, 'D');
M_.endo_names = char(M_.endo_names, 'h');
M_.endo_names_tex = char(M_.endo_names_tex, 'h');
M_.endo_names_long = char(M_.endo_names_long, 'h');
M_.endo_names = char(M_.endo_names, 'r_esb');
M_.endo_names_tex = char(M_.endo_names_tex, 'r\_esb');
M_.endo_names_long = char(M_.endo_names_long, 'r_esb');
M_.endo_names = char(M_.endo_names, 'Fw');
M_.endo_names_tex = char(M_.endo_names_tex, 'Fw');
M_.endo_names_long = char(M_.endo_names_long, 'Fw');
M_.endo_names = char(M_.endo_names, 'Kw');
M_.endo_names_tex = char(M_.endo_names_tex, 'Kw');
M_.endo_names_long = char(M_.endo_names_long, 'Kw');
M_.endo_names = char(M_.endo_names, 'R_p');
M_.endo_names_tex = char(M_.endo_names_tex, 'R\_p');
M_.endo_names_long = char(M_.endo_names_long, 'R_p');
M_.endo_names = char(M_.endo_names, 'L_cb');
M_.endo_names_tex = char(M_.endo_names_tex, 'L\_cb');
M_.endo_names_long = char(M_.endo_names_long, 'L_cb');
M_.endo_names = char(M_.endo_names, 'r_lcb');
M_.endo_names_tex = char(M_.endo_names_tex, 'r\_lcb');
M_.endo_names_long = char(M_.endo_names_long, 'r_lcb');
M_.endo_names = char(M_.endo_names, 'r_dcb');
M_.endo_names_tex = char(M_.endo_names_tex, 'r\_dcb');
M_.endo_names_long = char(M_.endo_names_long, 'r_dcb');
M_.endo_names = char(M_.endo_names, 'n_ag');
M_.endo_names_tex = char(M_.endo_names_tex, 'n\_ag');
M_.endo_names_long = char(M_.endo_names_long, 'n_ag');
M_.endo_names = char(M_.endo_names, 'lev_h');
M_.endo_names_tex = char(M_.endo_names_tex, 'lev\_h');
M_.endo_names_long = char(M_.endo_names_long, 'lev_h');
M_.endo_names = char(M_.endo_names, 'lev_l');
M_.endo_names_tex = char(M_.endo_names_tex, 'lev\_l');
M_.endo_names_long = char(M_.endo_names_long, 'lev_l');
M_.endo_names = char(M_.endo_names, 'soc_lev');
M_.endo_names_tex = char(M_.endo_names_tex, 'soc\_lev');
M_.endo_names_long = char(M_.endo_names_long, 'soc_lev');
M_.endo_names = char(M_.endo_names, 'TSF');
M_.endo_names_tex = char(M_.endo_names_tex, 'TSF');
M_.endo_names_long = char(M_.endo_names_long, 'TSF');
M_.endo_names = char(M_.endo_names, 'AUX_ENDO_LAG_33_1');
M_.endo_names_tex = char(M_.endo_names_tex, 'AUX\_ENDO\_LAG\_33\_1');
M_.endo_names_long = char(M_.endo_names_long, 'AUX_ENDO_LAG_33_1');
M_.endo_partitions = struct();
M_.param_names = 'beta';
M_.param_names_tex = 'beta';
M_.param_names_long = 'beta';
M_.param_names = char(M_.param_names, 'r_kh_ss');
M_.param_names_tex = char(M_.param_names_tex, 'r\_kh\_ss');
M_.param_names_long = char(M_.param_names_long, 'r_kh_ss');
M_.param_names = char(M_.param_names, 'delta_h');
M_.param_names_tex = char(M_.param_names_tex, 'delta\_h');
M_.param_names_long = char(M_.param_names_long, 'delta_h');
M_.param_names = char(M_.param_names, 'sigmaa_h');
M_.param_names_tex = char(M_.param_names_tex, 'sigmaa\_h');
M_.param_names_long = char(M_.param_names_long, 'sigmaa_h');
M_.param_names = char(M_.param_names, 'miu');
M_.param_names_tex = char(M_.param_names_tex, 'miu');
M_.param_names_long = char(M_.param_names_long, 'miu');
M_.param_names = char(M_.param_names, 'gamma_h');
M_.param_names_tex = char(M_.param_names_tex, 'gamma\_h');
M_.param_names_long = char(M_.param_names_long, 'gamma_h');
M_.param_names = char(M_.param_names, 'sigma_ss');
M_.param_names_tex = char(M_.param_names_tex, 'sigma\_ss');
M_.param_names_long = char(M_.param_names_long, 'sigma_ss');
M_.param_names = char(M_.param_names, 'w_h');
M_.param_names_tex = char(M_.param_names_tex, 'w\_h');
M_.param_names_long = char(M_.param_names_long, 'w_h');
M_.param_names = char(M_.param_names, 'rho_chi_sb');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_chi\_sb');
M_.param_names_long = char(M_.param_names_long, 'rho_chi_sb');
M_.param_names = char(M_.param_names, 'alpha_sb');
M_.param_names_tex = char(M_.param_names_tex, 'alpha\_sb');
M_.param_names_long = char(M_.param_names_long, 'alpha_sb');
M_.param_names = char(M_.param_names, 'delta_sb');
M_.param_names_tex = char(M_.param_names_tex, 'delta\_sb');
M_.param_names_long = char(M_.param_names_long, 'delta_sb');
M_.param_names = char(M_.param_names, 'w_sb');
M_.param_names_tex = char(M_.param_names_tex, 'w\_sb');
M_.param_names_long = char(M_.param_names_long, 'w_sb');
M_.param_names = char(M_.param_names, 'sigma');
M_.param_names_tex = char(M_.param_names_tex, 'sigma');
M_.param_names_long = char(M_.param_names_long, 'sigma');
M_.param_names = char(M_.param_names, 'r_kl_ss');
M_.param_names_tex = char(M_.param_names_tex, 'r\_kl\_ss');
M_.param_names_long = char(M_.param_names_long, 'r_kl_ss');
M_.param_names = char(M_.param_names, 'sigmaa_l');
M_.param_names_tex = char(M_.param_names_tex, 'sigmaa\_l');
M_.param_names_long = char(M_.param_names_long, 'sigmaa_l');
M_.param_names = char(M_.param_names, 'delta_ll');
M_.param_names_tex = char(M_.param_names_tex, 'delta\_ll');
M_.param_names_long = char(M_.param_names_long, 'delta_ll');
M_.param_names = char(M_.param_names, 'gamma_l');
M_.param_names_tex = char(M_.param_names_tex, 'gamma\_l');
M_.param_names_long = char(M_.param_names_long, 'gamma_l');
M_.param_names = char(M_.param_names, 'w_l');
M_.param_names_tex = char(M_.param_names_tex, 'w\_l');
M_.param_names_long = char(M_.param_names_long, 'w_l');
M_.param_names = char(M_.param_names, 'v');
M_.param_names_tex = char(M_.param_names_tex, 'v');
M_.param_names_long = char(M_.param_names_long, 'v');
M_.param_names = char(M_.param_names, 'epsilon_rl');
M_.param_names_tex = char(M_.param_names_tex, 'epsilon\_rl');
M_.param_names_long = char(M_.param_names_long, 'epsilon_rl');
M_.param_names = char(M_.param_names, 'eta_cost');
M_.param_names_tex = char(M_.param_names_tex, 'eta\_cost');
M_.param_names_long = char(M_.param_names_long, 'eta_cost');
M_.param_names = char(M_.param_names, 'rho_chi_rb');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_chi\_rb');
M_.param_names_long = char(M_.param_names_long, 'rho_chi_rb');
M_.param_names = char(M_.param_names, 'alpha_rb');
M_.param_names_tex = char(M_.param_names_tex, 'alpha\_rb');
M_.param_names_long = char(M_.param_names_long, 'alpha_rb');
M_.param_names = char(M_.param_names, 'kap_l');
M_.param_names_tex = char(M_.param_names_tex, 'kap\_l');
M_.param_names_long = char(M_.param_names_long, 'kap_l');
M_.param_names = char(M_.param_names, 'epsilon_d');
M_.param_names_tex = char(M_.param_names_tex, 'epsilon\_d');
M_.param_names_long = char(M_.param_names_long, 'epsilon_d');
M_.param_names = char(M_.param_names, 'kap_w');
M_.param_names_tex = char(M_.param_names_tex, 'kap\_w');
M_.param_names_long = char(M_.param_names_long, 'kap_w');
M_.param_names = char(M_.param_names, 'rho_a');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_a');
M_.param_names_long = char(M_.param_names_long, 'rho_a');
M_.param_names = char(M_.param_names, 'alpha');
M_.param_names_tex = char(M_.param_names_tex, 'alpha');
M_.param_names_long = char(M_.param_names_long, 'alpha');
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, 'rho');
M_.param_names_tex = char(M_.param_names_tex, 'rho');
M_.param_names_long = char(M_.param_names_long, 'rho');
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, 'pietarget');
M_.param_names_tex = char(M_.param_names_tex, 'pietarget');
M_.param_names_long = char(M_.param_names_long, 'pietarget');
M_.param_names = char(M_.param_names, 'lambdaf');
M_.param_names_tex = char(M_.param_names_tex, 'lambdaf');
M_.param_names_long = char(M_.param_names_long, 'lambdaf');
M_.param_names = char(M_.param_names, 'iota');
M_.param_names_tex = char(M_.param_names_tex, 'iota');
M_.param_names_long = char(M_.param_names_long, 'iota');
M_.param_names = char(M_.param_names, 'kappa_i');
M_.param_names_tex = char(M_.param_names_tex, 'kappa\_i');
M_.param_names_long = char(M_.param_names_long, 'kappa_i');
M_.param_names = char(M_.param_names, 'delta');
M_.param_names_tex = char(M_.param_names_tex, 'delta');
M_.param_names_long = char(M_.param_names_long, 'delta');
M_.param_names = char(M_.param_names, 'phi_l');
M_.param_names_tex = char(M_.param_names_tex, 'phi\_l');
M_.param_names_long = char(M_.param_names_long, 'phi_l');
M_.param_names = char(M_.param_names, 'delta_l');
M_.param_names_tex = char(M_.param_names_tex, 'delta\_l');
M_.param_names_long = char(M_.param_names_long, 'delta_l');
M_.param_names = char(M_.param_names, 'b');
M_.param_names_tex = char(M_.param_names_tex, 'b');
M_.param_names_long = char(M_.param_names_long, 'b');
M_.param_names = char(M_.param_names, 'sigma_c');
M_.param_names_tex = char(M_.param_names_tex, 'sigma\_c');
M_.param_names_long = char(M_.param_names_long, 'sigma_c');
M_.param_names = char(M_.param_names, 'xiw');
M_.param_names_tex = char(M_.param_names_tex, 'xiw');
M_.param_names_long = char(M_.param_names_long, 'xiw');
M_.param_names = char(M_.param_names, 'lambdaw');
M_.param_names_tex = char(M_.param_names_tex, 'lambdaw');
M_.param_names_long = char(M_.param_names_long, 'lambdaw');
M_.param_names = char(M_.param_names, 'iotaw');
M_.param_names_tex = char(M_.param_names_tex, 'iotaw');
M_.param_names_long = char(M_.param_names_long, 'iotaw');
M_.param_names = char(M_.param_names, 'etag');
M_.param_names_tex = char(M_.param_names_tex, 'etag');
M_.param_names_long = char(M_.param_names_long, 'etag');
M_.param_names = char(M_.param_names, 'rho_rd');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_rd');
M_.param_names_long = char(M_.param_names_long, 'rho_rd');
M_.param_names = char(M_.param_names, 'phi_rpie');
M_.param_names_tex = char(M_.param_names_tex, 'phi\_rpie');
M_.param_names_long = char(M_.param_names_long, 'phi_rpie');
M_.param_names = char(M_.param_names, 'phi_ry');
M_.param_names_tex = char(M_.param_names_tex, 'phi\_ry');
M_.param_names_long = char(M_.param_names_long, 'phi_ry');
M_.param_names = char(M_.param_names, 'R_pss');
M_.param_names_tex = char(M_.param_names_tex, 'R\_pss');
M_.param_names_long = char(M_.param_names_long, 'R_pss');
M_.param_names = char(M_.param_names, 'c_d');
M_.param_names_tex = char(M_.param_names_tex, 'c\_d');
M_.param_names_long = char(M_.param_names_long, 'c_d');
M_.param_names = char(M_.param_names, 'c_l');
M_.param_names_tex = char(M_.param_names_tex, 'c\_l');
M_.param_names_long = char(M_.param_names_long, 'c_l');
M_.param_names = char(M_.param_names, 'phi_lcb');
M_.param_names_tex = char(M_.param_names_tex, 'phi\_lcb');
M_.param_names_long = char(M_.param_names_long, 'phi_lcb');
M_.param_names = char(M_.param_names, 'phi_ll');
M_.param_names_tex = char(M_.param_names_tex, 'phi\_ll');
M_.param_names_long = char(M_.param_names_long, 'phi_ll');
M_.param_names = char(M_.param_names, 'phi_ly');
M_.param_names_tex = char(M_.param_names_tex, 'phi\_ly');
M_.param_names_long = char(M_.param_names_long, 'phi_ly');
M_.param_names = char(M_.param_names, 'phi_lpie');
M_.param_names_tex = char(M_.param_names_tex, 'phi\_lpie');
M_.param_names_long = char(M_.param_names_long, 'phi_lpie');
M_.param_partitions = struct();
M_.exo_det_nbr = 0;
M_.exo_nbr = 3;
M_.endo_nbr = 54;
M_.param_nbr = 54;
M_.orig_endo_nbr = 53;
M_.aux_vars(1).endo_index = 54;
M_.aux_vars(1).type = 1;
M_.aux_vars(1).orig_index = 34;
M_.aux_vars(1).orig_lead_lag = -1;
M_.Sigma_e = zeros(3, 3);
M_.Correlation_matrix = eye(3, 3);
M_.H = 0;
M_.Correlation_matrix_ME = 1;
M_.sigma_e_is_diagonal = 1;
M_.det_shocks = [];
options_.block=0;
options_.bytecode=0;
options_.use_dll=0;
M_.hessian_eq_zero = 0;
erase_compiled_function('rep_lishu_static');
erase_compiled_function('rep_lishu_dynamic');
M_.orig_eq_nbr = 53;
M_.eq_nbr = 54;
M_.ramsey_eq_nbr = 0;
M_.set_auxiliary_variables = exist(['./' M_.fname '_set_auxiliary_variables.m'], 'file') == 2;
M_.lead_lag_incidence = [
 1 20 74;
 0 21 0;
 0 22 0;
 0 23 75;
 0 24 76;
 2 25 0;
 3 26 0;
 0 27 0;
 0 28 0;
 4 29 0;
 5 30 0;
 6 31 0;
 0 32 0;
 0 33 0;
 0 34 0;
 0 35 0;
 0 36 0;
 0 37 0;
 7 38 0;
 8 39 0;
 0 40 0;
 9 41 0;
 0 42 77;
 0 43 0;
 10 44 0;
 0 45 0;
 0 46 0;
 0 47 0;
 11 48 0;
 0 49 0;
 12 50 0;
 13 51 78;
 0 52 0;
 14 53 79;
 15 54 80;
 0 55 81;
 0 56 82;
 16 57 83;
 0 58 84;
 0 59 0;
 0 60 0;
 0 61 0;
 0 62 85;
 0 63 86;
 17 64 0;
 18 65 0;
 0 66 0;
 0 67 0;
 0 68 0;
 0 69 0;
 0 70 0;
 0 71 0;
 0 72 0;
 19 73 0;]';
M_.nstatic = 27;
M_.nfwrd   = 8;
M_.npred   = 14;
M_.nboth   = 5;
M_.nsfwrd   = 13;
M_.nspred   = 19;
M_.ndynamic   = 27;
M_.equations_tags = {
};
M_.static_and_dynamic_models_differ = 0;
M_.exo_names_orig_ord = [1:3];
M_.maximum_lag = 1;
M_.maximum_lead = 1;
M_.maximum_endo_lag = 1;
M_.maximum_endo_lead = 1;
oo_.steady_state = zeros(54, 1);
M_.maximum_exo_lag = 0;
M_.maximum_exo_lead = 0;
oo_.exo_steady_state = zeros(3, 1);
M_.params = NaN(54, 1);
M_.NNZDerivatives = [249; 527; -1];
*****************************
M_.params( 1 ) = 0.9925;
beta = M_.params( 1 );
M_.params( 2 ) = 0.07;
r_kh_ss = M_.params( 2 );
M_.params( 3 ) = 0.03;
delta_h = M_.params( 3 );
M_.params( 4 ) = 18.9;
sigmaa_h = M_.params( 4 );
M_.params( 7 ) = 0.65;
sigma_ss = M_.params( 7 );
M_.params( 5 ) = 0.3;
miu = M_.params( 5 );
M_.params( 6 ) = 0.05;
gamma_h = M_.params( 6 );
M_.params( 8 ) = 0.14;
w_h = M_.params( 8 );
M_.params( 9 ) = 0.8;
rho_chi_sb = M_.params( 9 );
M_.params( 10 ) = 0.9;
alpha_sb = M_.params( 10 );
M_.params( 11 ) = 0.2;
delta_sb = M_.params( 11 );
M_.params( 12 ) = 0.116282119029118;
w_sb = M_.params( 12 );
M_.params( 13 ) = 0.8246211251235321;
sigma = M_.params( 13 );
M_.params( 15 ) = 18.9;
sigmaa_l = M_.params( 15 );
M_.params( 16 ) = 0.025;
delta_ll = M_.params( 16 );
M_.params( 17 ) = 0.03;
gamma_l = M_.params( 17 );
M_.params( 14 ) = 0.045;
r_kl_ss = M_.params( 14 );
M_.params( 18 ) = 1.12;
w_l = M_.params( 18 );
M_.params( 19 ) = 0.2;
v = M_.params( 19 );
M_.params( 20 ) = 389;
epsilon_rl = M_.params( 20 );
M_.params( 21 ) = 0.6;
eta_cost = M_.params( 21 );
M_.params( 9 ) = 0.9;
rho_chi_sb = M_.params( 9 );
M_.params( 23 ) = 40;
alpha_rb = M_.params( 23 );
M_.params( 24 ) = M_.params(20);
kap_l = M_.params( 24 );
M_.params( 25 ) = (-222);
epsilon_d = M_.params( 25 );
M_.params( 26 ) = 0.4;
kap_w = M_.params( 26 );
M_.params( 27 ) = 0.95;
rho_a = M_.params( 27 );
M_.params( 28 ) = 0.6;
alpha = M_.params( 28 );
M_.params( 29 ) = 0.16;
eta = M_.params( 29 );
M_.params( 30 ) = 0.71;
rho = M_.params( 30 );
M_.params( 33 ) = 1.199129448472893;
lambdaf = M_.params( 33 );
M_.params( 31 ) = 0.8;
xip = M_.params( 31 );
M_.params( 32 ) = 1;
pietarget = M_.params( 32 );
M_.params( 34 ) = 0.16;
iota = M_.params( 34 );
M_.params( 35 ) = 2.5;
kappa_i = M_.params( 35 );
M_.params( 36 ) = 0.03;
delta = M_.params( 36 );
M_.params( 37 ) = 14.129780032855374;
phi_l = M_.params( 37 );
M_.params( 38 ) = 1;
delta_l = M_.params( 38 );
M_.params( 39 ) = 0.8;
b = M_.params( 39 );
M_.params( 40 ) = 1;
sigma_c = M_.params( 40 );
M_.params( 41 ) = 0.75;
xiw = M_.params( 41 );
M_.params( 42 ) = 1.05;
lambdaw = M_.params( 42 );
M_.params( 43 ) = 0.29;
iotaw = M_.params( 43 );
M_.params( 44 ) = 0.135;
etag = M_.params( 44 );
M_.params( 45 ) = 0.8;
rho_rd = M_.params( 45 );
M_.params( 48 ) = 1/M_.params(1)-1;
R_pss = M_.params( 48 );
M_.params( 46 ) = 1.2;
phi_rpie = M_.params( 46 );
M_.params( 47 ) = 0.4;
phi_ry = M_.params( 47 );
M_.params( 49 ) = 0.0001;
c_d = M_.params( 49 );
M_.params( 50 ) = 0.000158;
c_l = M_.params( 50 );
M_.params( 51 ) = 0.8;
phi_lcb = M_.params( 51 );
M_.params( 52 ) = 0.3;
phi_ll = M_.params( 52 );
M_.params( 53 ) = 0.1;
phi_ly = M_.params( 53 );
M_.params( 54 ) = 1.8;
phi_lpie = M_.params( 54 );
%
% INITVAL instructions
%
options_.initval_file = 0;
oo_.steady_state( 23 ) = 0.015;
oo_.steady_state( 2 ) = M_.params(2);
oo_.steady_state( 3 ) = 1;
oo_.steady_state( 35 ) = 1;
oo_.steady_state( 4 ) = oo_.steady_state(2)-M_.params(3);
oo_.steady_state( 5 ) = 0.230648545301527;
oo_.steady_state( 6 ) = 27.409417048714833;
oo_.steady_state( 7 ) = 21.008063698336578;
oo_.steady_state( 8 ) = oo_.steady_state(35)*oo_.steady_state(6)-oo_.steady_state(7);
oo_.steady_state( 9 ) = oo_.steady_state(5);
oo_.steady_state( 10 ) = 0;
oo_.steady_state( 11 ) = 0.581410883718849;
oo_.steady_state( 12 ) = 0.2;
oo_.steady_state( 13 ) = 0.03;
oo_.steady_state( 42 ) = 1;
oo_.steady_state( 16 ) = M_.params(14);
oo_.steady_state( 18 ) = 1;
oo_.steady_state( 17 ) = oo_.steady_state(16)-M_.params(16);
oo_.steady_state( 22 ) = 1+oo_.steady_state(17)-1/M_.params(1);
oo_.steady_state( 19 ) = (oo_.steady_state(2)/oo_.steady_state(16))^(1/(1-M_.params(30)))*oo_.steady_state(6)*oo_.steady_state(3)/oo_.steady_state(18);
oo_.steady_state( 20 ) = (oo_.steady_state(19)*(1-M_.params(17))*(oo_.steady_state(17)-oo_.steady_state(22))+M_.params(18))/(1-(1-M_.params(17))*(1+oo_.steady_state(22)));
oo_.steady_state( 21 ) = oo_.steady_state(35)*oo_.steady_state(19)-oo_.steady_state(20);
oo_.steady_state( 40 ) = oo_.steady_state(21)*M_.params(19)/(1-M_.params(19));
oo_.steady_state( 24 ) = oo_.steady_state(21)*oo_.steady_state(22)*(1-M_.params(29))-(1-M_.params(29))*M_.params(48)*oo_.steady_state(40);
oo_.steady_state( 25 ) = 0;
oo_.steady_state( 26 ) = M_.params(20);
oo_.steady_state( 27 ) = M_.params(48);
oo_.steady_state( 28 ) = M_.params(48);
oo_.steady_state( 31 ) = (M_.params(29)*(oo_.steady_state(6)*oo_.steady_state(3))^M_.params(30)+(1-M_.params(29))*(oo_.steady_state(18)*oo_.steady_state(19))^M_.params(30))^(1/M_.params(30));
oo_.steady_state( 34 ) = M_.params(36)*(oo_.steady_state(6)*oo_.steady_state(3)*M_.params(29)+(1-M_.params(29))*oo_.steady_state(18)*oo_.steady_state(19));
oo_.steady_state( 33 ) = 1/M_.params(33);
oo_.steady_state( 41 ) = 0.292139066660454;
oo_.steady_state( 32 ) = 22.849056125573080;
oo_.steady_state( 30 ) = oo_.steady_state(31)^M_.params(28)*oo_.steady_state(41)^(1-M_.params(28));
oo_.steady_state( 38 ) = (1-M_.params(44))*oo_.steady_state(30)-oo_.steady_state(34)-oo_.steady_state(6)*oo_.steady_state(35)*M_.params(29)*M_.params(5)*normcdf((log(oo_.steady_state(5))+M_.params(13)^2/2)/M_.params(13)-M_.params(13),0,1)*(1+oo_.steady_state(4))-(1-M_.params(21))*oo_.steady_state(24);
oo_.steady_state( 39 ) = (1-M_.params(1)*M_.params(39))/((1-M_.params(39))*oo_.steady_state(38));
oo_.steady_state( 36 ) = oo_.steady_state(30)*oo_.steady_state(39)/(1-M_.params(1)*M_.params(31));
oo_.steady_state( 37 ) = oo_.steady_state(36);
oo_.steady_state( 1 ) = M_.params(32);
oo_.steady_state( 43 ) = oo_.steady_state(41)*oo_.steady_state(39)/M_.params(42)/(1-M_.params(1)*M_.params(41));
oo_.steady_state( 44 ) = oo_.steady_state(41)^(1+M_.params(38))/(1-M_.params(1)*M_.params(41));
oo_.steady_state( 45 ) = M_.params(48);
oo_.steady_state( 46 ) = 0.3*oo_.steady_state(30);
oo_.steady_state( 47 ) = 0.06;
oo_.steady_state( 48 ) = 0.027;
oo_.steady_state( 49 ) = oo_.steady_state(7)*M_.params(29)+oo_.steady_state(20)*(1-M_.params(29));
oo_.steady_state( 50 ) = oo_.steady_state(35)*oo_.steady_state(6)/oo_.steady_state(7);
oo_.steady_state( 51 ) = oo_.steady_state(35)*oo_.steady_state(19)/oo_.steady_state(20);
oo_.steady_state( 52 ) = M_.params(29)*oo_.steady_state(50)+(1-M_.params(29))*oo_.steady_state(51);
oo_.steady_state( 53 ) = M_.params(29)*oo_.steady_state(8)+oo_.steady_state(21)*(1-M_.params(29));
oo_.steady_state(54)=oo_.steady_state(34);
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
%
% ENDVAL instructions
%
ys0_= oo_.steady_state;
ex0_ = oo_.exo_steady_state;
oo_.steady_state( 47 ) = 0.06;
oo_.steady_state( 48 ) = 0.027;
oo_.steady_state(54)=oo_.steady_state(34);
steady;
oo_.dr.eigval = check(M_,options_,oo_);
model_diagnostics(M_,options_,oo_);
save('rep_lishu_results.mat', 'oo_', 'M_', 'options_');
if exist('estim_params_', 'var') == 1
  save('rep_lishu_results.mat', 'estim_params_', '-append');
end
if exist('bayestopt_', 'var') == 1
  save('rep_lishu_results.mat', 'bayestopt_', '-append');
end
if exist('dataset_', 'var') == 1
  save('rep_lishu_results.mat', 'dataset_', '-append');
end
if exist('estimation_info', 'var') == 1
  save('rep_lishu_results.mat', 'estimation_info', '-append');
end
if exist('dataset_info', 'var') == 1
  save('rep_lishu_results.mat', 'dataset_info', '-append');
end
if exist('oo_recursive_', 'var') == 1
  save('rep_lishu_results.mat', 'oo_recursive_', '-append');
end


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