%
% 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 = 'baseline_V1_HHlender';
M_.dynare_version = '4.5.5';
oo_.dynare_version = '4.5.5';
options_.dynare_version = '4.5.5';
%
% Some global variables initialization
%
global_initialization;
diary off;
diary('baseline_V1_HHlender.log');
M_.exo_names = 'eps_a';
M_.exo_names_tex = 'eps\_a';
M_.exo_names_long = 'eps_a';
M_.exo_names = char(M_.exo_names, 'eps_ag');
M_.exo_names_tex = char(M_.exo_names_tex, 'eps\_ag');
M_.exo_names_long = char(M_.exo_names_long, 'eps_ag');
M_.exo_names = char(M_.exo_names, 'eps_rf');
M_.exo_names_tex = char(M_.exo_names_tex, 'eps\_rf');
M_.exo_names_long = char(M_.exo_names_long, 'eps_rf');
M_.exo_names = char(M_.exo_names, 'eps_mu');
M_.exo_names_tex = char(M_.exo_names_tex, 'eps\_mu');
M_.exo_names_long = char(M_.exo_names_long, 'eps_mu');
M_.exo_names = char(M_.exo_names, 'eps_lg');
M_.exo_names_tex = char(M_.exo_names_tex, 'eps\_lg');
M_.exo_names_long = char(M_.exo_names_long, 'eps_lg');
M_.exo_names = char(M_.exo_names, 'eps_g');
M_.exo_names_tex = char(M_.exo_names_tex, 'eps\_g');
M_.exo_names_long = char(M_.exo_names_long, 'eps_g');
M_.exo_names = char(M_.exo_names, 'eps_ig');
M_.exo_names_tex = char(M_.exo_names_tex, 'eps\_ig');
M_.exo_names_long = char(M_.exo_names_long, 'eps_ig');
M_.exo_names = char(M_.exo_names, 'eps_tauk');
M_.exo_names_tex = char(M_.exo_names_tex, 'eps\_tauk');
M_.exo_names_long = char(M_.exo_names_long, 'eps_tauk');
M_.exo_names = char(M_.exo_names, 'eps_taun');
M_.exo_names_tex = char(M_.exo_names_tex, 'eps\_taun');
M_.exo_names_long = char(M_.exo_names_long, 'eps_taun');
M_.exo_names = char(M_.exo_names, 'eps_tauc');
M_.exo_names_tex = char(M_.exo_names_tex, 'eps\_tauc');
M_.exo_names_long = char(M_.exo_names_long, 'eps_tauc');
M_.exo_names = char(M_.exo_names, 'eps_lst');
M_.exo_names_tex = char(M_.exo_names_tex, 'eps\_lst');
M_.exo_names_long = char(M_.exo_names_long, 'eps_lst');
M_.endo_names = 'c';
M_.endo_names_tex = 'c';
M_.endo_names_long = 'c';
M_.endo_names = char(M_.endo_names, 'lp');
M_.endo_names_tex = char(M_.endo_names_tex, 'lp');
M_.endo_names_long = char(M_.endo_names_long, 'lp');
M_.endo_names = char(M_.endo_names, 'it');
M_.endo_names_tex = char(M_.endo_names_tex, 'it');
M_.endo_names_long = char(M_.endo_names_long, 'it');
M_.endo_names = char(M_.endo_names, 'int');
M_.endo_names_tex = char(M_.endo_names_tex, 'int');
M_.endo_names_long = char(M_.endo_names_long, 'int');
M_.endo_names = char(M_.endo_names, 'kt');
M_.endo_names_tex = char(M_.endo_names_tex, 'kt');
M_.endo_names_long = char(M_.endo_names_long, 'kt');
M_.endo_names = char(M_.endo_names, 'knt');
M_.endo_names_tex = char(M_.endo_names_tex, 'knt');
M_.endo_names_long = char(M_.endo_names_long, 'knt');
M_.endo_names = char(M_.endo_names, 'inve');
M_.endo_names_tex = char(M_.endo_names_tex, 'inve');
M_.endo_names_long = char(M_.endo_names_long, 'inve');
M_.endo_names = char(M_.endo_names, 'lambda1');
M_.endo_names_tex = char(M_.endo_names_tex, 'lambda1');
M_.endo_names_long = char(M_.endo_names_long, 'lambda1');
M_.endo_names = char(M_.endo_names, 'lambda2');
M_.endo_names_tex = char(M_.endo_names_tex, 'lambda2');
M_.endo_names_long = char(M_.endo_names_long, 'lambda2');
M_.endo_names = char(M_.endo_names, 'lambda3');
M_.endo_names_tex = char(M_.endo_names_tex, 'lambda3');
M_.endo_names_long = char(M_.endo_names_long, 'lambda3');
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, 'rstar');
M_.endo_names_tex = char(M_.endo_names_tex, 'rstar');
M_.endo_names_long = char(M_.endo_names_long, 'rstar');
M_.endo_names = char(M_.endo_names, 'ltt');
M_.endo_names_tex = char(M_.endo_names_tex, 'ltt');
M_.endo_names_long = char(M_.endo_names_long, 'ltt');
M_.endo_names = char(M_.endo_names, 'lnt');
M_.endo_names_tex = char(M_.endo_names_tex, 'lnt');
M_.endo_names_long = char(M_.endo_names_long, 'lnt');
M_.endo_names = char(M_.endo_names, 'rt');
M_.endo_names_tex = char(M_.endo_names_tex, 'rt');
M_.endo_names_long = char(M_.endo_names_long, 'rt');
M_.endo_names = char(M_.endo_names, 'rnt');
M_.endo_names_tex = char(M_.endo_names_tex, 'rnt');
M_.endo_names_long = char(M_.endo_names_long, 'rnt');
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, 'ynt');
M_.endo_names_tex = char(M_.endo_names_tex, 'ynt');
M_.endo_names_long = char(M_.endo_names_long, 'ynt');
M_.endo_names = char(M_.endo_names, 'yt');
M_.endo_names_tex = char(M_.endo_names_tex, 'yt');
M_.endo_names_long = char(M_.endo_names_long, 'yt');
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, 'et');
M_.endo_names_tex = char(M_.endo_names_tex, 'et');
M_.endo_names_long = char(M_.endo_names_long, 'et');
M_.endo_names = char(M_.endo_names, 'pn');
M_.endo_names_tex = char(M_.endo_names_tex, 'pn');
M_.endo_names_long = char(M_.endo_names_long, 'pn');
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, 'infnt');
M_.endo_names_tex = char(M_.endo_names_tex, 'infnt');
M_.endo_names_long = char(M_.endo_names_long, 'infnt');
M_.endo_names = char(M_.endo_names, 'mc');
M_.endo_names_tex = char(M_.endo_names_tex, 'mc');
M_.endo_names_long = char(M_.endo_names_long, 'mc');
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, 'infntopt');
M_.endo_names_tex = char(M_.endo_names_tex, 'infntopt');
M_.endo_names_long = char(M_.endo_names_long, 'infntopt');
M_.endo_names = char(M_.endo_names, 'x1');
M_.endo_names_tex = char(M_.endo_names_tex, 'x1');
M_.endo_names_long = char(M_.endo_names_long, 'x1');
M_.endo_names = char(M_.endo_names, 'x2');
M_.endo_names_tex = char(M_.endo_names_tex, 'x2');
M_.endo_names_long = char(M_.endo_names_long, 'x2');
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, 'kg');
M_.endo_names_tex = char(M_.endo_names_tex, 'kg');
M_.endo_names_long = char(M_.endo_names_long, 'kg');
M_.endo_names = char(M_.endo_names, 'zg');
M_.endo_names_tex = char(M_.endo_names_tex, 'zg');
M_.endo_names_long = char(M_.endo_names_long, 'zg');
M_.endo_names = char(M_.endo_names, 'tauc');
M_.endo_names_tex = char(M_.endo_names_tex, 'tauc');
M_.endo_names_long = char(M_.endo_names_long, 'tauc');
M_.endo_names = char(M_.endo_names, 'taun');
M_.endo_names_tex = char(M_.endo_names_tex, 'taun');
M_.endo_names_long = char(M_.endo_names_long, 'taun');
M_.endo_names = char(M_.endo_names, 'tauk');
M_.endo_names_tex = char(M_.endo_names_tex, 'tauk');
M_.endo_names_long = char(M_.endo_names_long, 'tauk');
M_.endo_names = char(M_.endo_names, 'bstar');
M_.endo_names_tex = char(M_.endo_names_tex, 'bstar');
M_.endo_names_long = char(M_.endo_names_long, 'bstar');
M_.endo_names = char(M_.endo_names, 'bgstar');
M_.endo_names_tex = char(M_.endo_names_tex, 'bgstar');
M_.endo_names_long = char(M_.endo_names_long, 'bgstar');
M_.endo_names = char(M_.endo_names, 'b');
M_.endo_names_tex = char(M_.endo_names_tex, 'b');
M_.endo_names_long = char(M_.endo_names_long, 'b');
M_.endo_names = char(M_.endo_names, 'bg');
M_.endo_names_tex = char(M_.endo_names_tex, 'bg');
M_.endo_names_long = char(M_.endo_names_long, 'bg');
M_.endo_names = char(M_.endo_names, 'wg');
M_.endo_names_tex = char(M_.endo_names_tex, 'wg');
M_.endo_names_long = char(M_.endo_names_long, 'wg');
M_.endo_names = char(M_.endo_names, 'by');
M_.endo_names_tex = char(M_.endo_names_tex, 'by');
M_.endo_names_long = char(M_.endo_names_long, 'by');
M_.endo_names = char(M_.endo_names, 'lst');
M_.endo_names_tex = char(M_.endo_names_tex, 'lst');
M_.endo_names_long = char(M_.endo_names_long, 'lst');
M_.endo_names = char(M_.endo_names, 'tr');
M_.endo_names_tex = char(M_.endo_names_tex, 'tr');
M_.endo_names_long = char(M_.endo_names_long, 'tr');
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, 'ag');
M_.endo_names_tex = char(M_.endo_names_tex, 'ag');
M_.endo_names_long = char(M_.endo_names_long, 'ag');
M_.endo_names = char(M_.endo_names, 'rf');
M_.endo_names_tex = char(M_.endo_names_tex, 'rf');
M_.endo_names_long = char(M_.endo_names_long, 'rf');
M_.endo_names = char(M_.endo_names, 'mu');
M_.endo_names_tex = char(M_.endo_names_tex, 'mu');
M_.endo_names_long = char(M_.endo_names_long, 'mu');
M_.endo_names = char(M_.endo_names, 'lg');
M_.endo_names_tex = char(M_.endo_names_tex, 'lg');
M_.endo_names_long = char(M_.endo_names_long, 'lg');
M_.endo_names = char(M_.endo_names, 'g');
M_.endo_names_tex = char(M_.endo_names_tex, 'g');
M_.endo_names_long = char(M_.endo_names_long, 'g');
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_partitions = struct();
M_.param_names = 'betta';
M_.param_names_tex = 'betta';
M_.param_names_long = 'betta';
M_.param_names = char(M_.param_names, 'alphat');
M_.param_names_tex = char(M_.param_names_tex, 'alphat');
M_.param_names_long = char(M_.param_names_long, 'alphat');
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, 'alphant');
M_.param_names_tex = char(M_.param_names_tex, 'alphant');
M_.param_names_long = char(M_.param_names_long, 'alphant');
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, 'deltaz');
M_.param_names_tex = char(M_.param_names_tex, 'deltaz');
M_.param_names_long = char(M_.param_names_long, 'deltaz');
M_.param_names = char(M_.param_names, 'eps');
M_.param_names_tex = char(M_.param_names_tex, 'eps');
M_.param_names_long = char(M_.param_names_long, 'eps');
M_.param_names = char(M_.param_names, 'epsp');
M_.param_names_tex = char(M_.param_names_tex, 'epsp');
M_.param_names_long = char(M_.param_names_long, 'epsp');
M_.param_names = char(M_.param_names, 'espw');
M_.param_names_tex = char(M_.param_names_tex, 'espw');
M_.param_names_long = char(M_.param_names_long, 'espw');
M_.param_names = char(M_.param_names, 'psil');
M_.param_names_tex = char(M_.param_names_tex, 'psil');
M_.param_names_long = char(M_.param_names_long, 'psil');
M_.param_names = char(M_.param_names, 'nut');
M_.param_names_tex = char(M_.param_names_tex, 'nut');
M_.param_names_long = char(M_.param_names_long, 'nut');
M_.param_names = char(M_.param_names, 'nunt');
M_.param_names_tex = char(M_.param_names_tex, 'nunt');
M_.param_names_long = char(M_.param_names_long, 'nunt');
M_.param_names = char(M_.param_names, 'nub');
M_.param_names_tex = char(M_.param_names_tex, 'nub');
M_.param_names_long = char(M_.param_names_long, 'nub');
M_.param_names = char(M_.param_names, 'phiw');
M_.param_names_tex = char(M_.param_names_tex, 'phiw');
M_.param_names_long = char(M_.param_names_long, 'phiw');
M_.param_names = char(M_.param_names, 'ac');
M_.param_names_tex = char(M_.param_names_tex, 'ac');
M_.param_names_long = char(M_.param_names_long, 'ac');
M_.param_names = char(M_.param_names, 'alphag');
M_.param_names_tex = char(M_.param_names_tex, 'alphag');
M_.param_names_long = char(M_.param_names_long, 'alphag');
M_.param_names = char(M_.param_names, 'psif');
M_.param_names_tex = char(M_.param_names_tex, 'psif');
M_.param_names_long = char(M_.param_names_long, 'psif');
M_.param_names = char(M_.param_names, 'alphaz');
M_.param_names_tex = char(M_.param_names_tex, 'alphaz');
M_.param_names_long = char(M_.param_names_long, 'alphaz');
M_.param_names = char(M_.param_names, 'alphas');
M_.param_names_tex = char(M_.param_names_tex, 'alphas');
M_.param_names_long = char(M_.param_names_long, 'alphas');
M_.param_names = char(M_.param_names, 'zeta');
M_.param_names_tex = char(M_.param_names_tex, 'zeta');
M_.param_names_long = char(M_.param_names_long, 'zeta');
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, 'phib');
M_.param_names_tex = char(M_.param_names_tex, 'phib');
M_.param_names_long = char(M_.param_names_long, 'phib');
M_.param_names = char(M_.param_names, 'chil');
M_.param_names_tex = char(M_.param_names_tex, 'chil');
M_.param_names_long = char(M_.param_names_long, 'chil');
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, 'rho_ag');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_ag');
M_.param_names_long = char(M_.param_names_long, 'rho_ag');
M_.param_names = char(M_.param_names, 'rho_g');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_g');
M_.param_names_long = char(M_.param_names_long, 'rho_g');
M_.param_names = char(M_.param_names, 'rho_ig');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_ig');
M_.param_names_long = char(M_.param_names_long, 'rho_ig');
M_.param_names = char(M_.param_names, 'rho_lg');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_lg');
M_.param_names_long = char(M_.param_names_long, 'rho_lg');
M_.param_names = char(M_.param_names, 'rho_rf');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_rf');
M_.param_names_long = char(M_.param_names_long, 'rho_rf');
M_.param_names = char(M_.param_names, 'rho_r');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_r');
M_.param_names_long = char(M_.param_names_long, 'rho_r');
M_.param_names = char(M_.param_names, 'rho_mu');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_mu');
M_.param_names_long = char(M_.param_names_long, 'rho_mu');
M_.param_names = char(M_.param_names, 'rho_infstar');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_infstar');
M_.param_names_long = char(M_.param_names_long, 'rho_infstar');
M_.param_names = char(M_.param_names, 'rho_tauk');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_tauk');
M_.param_names_long = char(M_.param_names_long, 'rho_tauk');
M_.param_names = char(M_.param_names, 'rho_taun');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_taun');
M_.param_names_long = char(M_.param_names_long, 'rho_taun');
M_.param_names = char(M_.param_names, 'rho_tauc');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_tauc');
M_.param_names_long = char(M_.param_names_long, 'rho_tauc');
M_.param_names = char(M_.param_names, 'rho_lst');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_lst');
M_.param_names_long = char(M_.param_names_long, 'rho_lst');
M_.param_names = char(M_.param_names, 'gamma_c');
M_.param_names_tex = char(M_.param_names_tex, 'gamma\_c');
M_.param_names_long = char(M_.param_names_long, 'gamma_c');
M_.param_names = char(M_.param_names, 'gamma_n');
M_.param_names_tex = char(M_.param_names_tex, 'gamma\_n');
M_.param_names_long = char(M_.param_names_long, 'gamma_n');
M_.param_names = char(M_.param_names, 'gamma_k');
M_.param_names_tex = char(M_.param_names_tex, 'gamma\_k');
M_.param_names_long = char(M_.param_names_long, 'gamma_k');
M_.param_names = char(M_.param_names, 'gamma_t');
M_.param_names_tex = char(M_.param_names_tex, 'gamma\_t');
M_.param_names_long = char(M_.param_names_long, 'gamma_t');
M_.param_names = char(M_.param_names, 'aratio');
M_.param_names_tex = char(M_.param_names_tex, 'aratio');
M_.param_names_long = char(M_.param_names_long, 'aratio');
M_.param_names = char(M_.param_names, 'bratio');
M_.param_names_tex = char(M_.param_names_tex, 'bratio');
M_.param_names_long = char(M_.param_names_long, 'bratio');
M_.param_names = char(M_.param_names, 'a_bar');
M_.param_names_tex = char(M_.param_names_tex, 'a\_bar');
M_.param_names_long = char(M_.param_names_long, 'a_bar');
M_.param_names = char(M_.param_names, 'ag_bar');
M_.param_names_tex = char(M_.param_names_tex, 'ag\_bar');
M_.param_names_long = char(M_.param_names_long, 'ag_bar');
M_.param_names = char(M_.param_names, 'y_bar');
M_.param_names_tex = char(M_.param_names_tex, 'y\_bar');
M_.param_names_long = char(M_.param_names_long, 'y_bar');
M_.param_names = char(M_.param_names, 'g_bar');
M_.param_names_tex = char(M_.param_names_tex, 'g\_bar');
M_.param_names_long = char(M_.param_names_long, 'g_bar');
M_.param_names = char(M_.param_names, 'c_bar');
M_.param_names_tex = char(M_.param_names_tex, 'c\_bar');
M_.param_names_long = char(M_.param_names_long, 'c_bar');
M_.param_names = char(M_.param_names, 'i_bar');
M_.param_names_tex = char(M_.param_names_tex, 'i\_bar');
M_.param_names_long = char(M_.param_names_long, 'i_bar');
M_.param_names = char(M_.param_names, 'k_bar');
M_.param_names_tex = char(M_.param_names_tex, 'k\_bar');
M_.param_names_long = char(M_.param_names_long, 'k_bar');
M_.param_names = char(M_.param_names, 'rnt_bar');
M_.param_names_tex = char(M_.param_names_tex, 'rnt\_bar');
M_.param_names_long = char(M_.param_names_long, 'rnt_bar');
M_.param_names = char(M_.param_names, 'rt_bar');
M_.param_names_tex = char(M_.param_names_tex, 'rt\_bar');
M_.param_names_long = char(M_.param_names_long, 'rt_bar');
M_.param_names = char(M_.param_names, 'et_bar');
M_.param_names_tex = char(M_.param_names_tex, 'et\_bar');
M_.param_names_long = char(M_.param_names_long, 'et_bar');
M_.param_names = char(M_.param_names, 'pn_bar');
M_.param_names_tex = char(M_.param_names_tex, 'pn\_bar');
M_.param_names_long = char(M_.param_names_long, 'pn_bar');
M_.param_names = char(M_.param_names, 'infstar_bar');
M_.param_names_tex = char(M_.param_names_tex, 'infstar\_bar');
M_.param_names_long = char(M_.param_names_long, 'infstar_bar');
M_.param_names = char(M_.param_names, 'infnt_bar');
M_.param_names_tex = char(M_.param_names_tex, 'infnt\_bar');
M_.param_names_long = char(M_.param_names_long, 'infnt_bar');
M_.param_names = char(M_.param_names, 'infl_bar');
M_.param_names_tex = char(M_.param_names_tex, 'infl\_bar');
M_.param_names_long = char(M_.param_names_long, 'infl_bar');
M_.param_names = char(M_.param_names, 'rf_bar');
M_.param_names_tex = char(M_.param_names_tex, 'rf\_bar');
M_.param_names_long = char(M_.param_names_long, 'rf_bar');
M_.param_names = char(M_.param_names, 'r_bar');
M_.param_names_tex = char(M_.param_names_tex, 'r\_bar');
M_.param_names_long = char(M_.param_names_long, 'r_bar');
M_.param_names = char(M_.param_names, 'mu_bar');
M_.param_names_tex = char(M_.param_names_tex, 'mu\_bar');
M_.param_names_long = char(M_.param_names_long, 'mu_bar');
M_.param_names = char(M_.param_names, 'mc_bar');
M_.param_names_tex = char(M_.param_names_tex, 'mc\_bar');
M_.param_names_long = char(M_.param_names_long, 'mc_bar');
M_.param_names = char(M_.param_names, 'kt_bar');
M_.param_names_tex = char(M_.param_names_tex, 'kt\_bar');
M_.param_names_long = char(M_.param_names_long, 'kt_bar');
M_.param_names = char(M_.param_names, 'it_bar');
M_.param_names_tex = char(M_.param_names_tex, 'it\_bar');
M_.param_names_long = char(M_.param_names_long, 'it_bar');
M_.param_names = char(M_.param_names, 'knt_bar');
M_.param_names_tex = char(M_.param_names_tex, 'knt\_bar');
M_.param_names_long = char(M_.param_names_long, 'knt_bar');
M_.param_names = char(M_.param_names, 'int_bar');
M_.param_names_tex = char(M_.param_names_tex, 'int\_bar');
M_.param_names_long = char(M_.param_names_long, 'int_bar');
M_.param_names = char(M_.param_names, 'v_bar');
M_.param_names_tex = char(M_.param_names_tex, 'v\_bar');
M_.param_names_long = char(M_.param_names_long, 'v_bar');
M_.param_names = char(M_.param_names, 'x1_bar');
M_.param_names_tex = char(M_.param_names_tex, 'x1\_bar');
M_.param_names_long = char(M_.param_names_long, 'x1_bar');
M_.param_names = char(M_.param_names, 'x2_var');
M_.param_names_tex = char(M_.param_names_tex, 'x2\_var');
M_.param_names_long = char(M_.param_names_long, 'x2_var');
M_.param_names = char(M_.param_names, 'infntopt_bar');
M_.param_names_tex = char(M_.param_names_tex, 'infntopt\_bar');
M_.param_names_long = char(M_.param_names_long, 'infntopt_bar');
M_.param_names = char(M_.param_names, 'lg_bar');
M_.param_names_tex = char(M_.param_names_tex, 'lg\_bar');
M_.param_names_long = char(M_.param_names_long, 'lg_bar');
M_.param_names = char(M_.param_names, 'ltt_bar');
M_.param_names_tex = char(M_.param_names_tex, 'ltt\_bar');
M_.param_names_long = char(M_.param_names_long, 'ltt_bar');
M_.param_names = char(M_.param_names, 'lnt_bar');
M_.param_names_tex = char(M_.param_names_tex, 'lnt\_bar');
M_.param_names_long = char(M_.param_names_long, 'lnt_bar');
M_.param_names = char(M_.param_names, 'ynt_bar');
M_.param_names_tex = char(M_.param_names_tex, 'ynt\_bar');
M_.param_names_long = char(M_.param_names_long, 'ynt_bar');
M_.param_names = char(M_.param_names, 'yt_bar');
M_.param_names_tex = char(M_.param_names_tex, 'yt\_bar');
M_.param_names_long = char(M_.param_names_long, 'yt_bar');
M_.param_names = char(M_.param_names, 'lp_bar');
M_.param_names_tex = char(M_.param_names_tex, 'lp\_bar');
M_.param_names_long = char(M_.param_names_long, 'lp_bar');
M_.param_names = char(M_.param_names, 'w_bar');
M_.param_names_tex = char(M_.param_names_tex, 'w\_bar');
M_.param_names_long = char(M_.param_names_long, 'w_bar');
M_.param_names = char(M_.param_names, 'wg_bar');
M_.param_names_tex = char(M_.param_names_tex, 'wg\_bar');
M_.param_names_long = char(M_.param_names_long, 'wg_bar');
M_.param_names = char(M_.param_names, 'lambda1_bar');
M_.param_names_tex = char(M_.param_names_tex, 'lambda1\_bar');
M_.param_names_long = char(M_.param_names_long, 'lambda1_bar');
M_.param_names = char(M_.param_names, 'lambda2_bar');
M_.param_names_tex = char(M_.param_names_tex, 'lambda2\_bar');
M_.param_names_long = char(M_.param_names_long, 'lambda2_bar');
M_.param_names = char(M_.param_names, 'lambda3_bar');
M_.param_names_tex = char(M_.param_names_tex, 'lambda3\_bar');
M_.param_names_long = char(M_.param_names_long, 'lambda3_bar');
M_.param_names = char(M_.param_names, 'p_bar');
M_.param_names_tex = char(M_.param_names_tex, 'p\_bar');
M_.param_names_long = char(M_.param_names_long, 'p_bar');
M_.param_names = char(M_.param_names, 'zg_bar');
M_.param_names_tex = char(M_.param_names_tex, 'zg\_bar');
M_.param_names_long = char(M_.param_names_long, 'zg_bar');
M_.param_names = char(M_.param_names, 's_bar');
M_.param_names_tex = char(M_.param_names_tex, 's\_bar');
M_.param_names_long = char(M_.param_names_long, 's_bar');
M_.param_names = char(M_.param_names, 'b_bar');
M_.param_names_tex = char(M_.param_names_tex, 'b\_bar');
M_.param_names_long = char(M_.param_names_long, 'b_bar');
M_.param_names = char(M_.param_names, 'bg_bar');
M_.param_names_tex = char(M_.param_names_tex, 'bg\_bar');
M_.param_names_long = char(M_.param_names_long, 'bg_bar');
M_.param_names = char(M_.param_names, 'bgstar_bar');
M_.param_names_tex = char(M_.param_names_tex, 'bgstar\_bar');
M_.param_names_long = char(M_.param_names_long, 'bgstar_bar');
M_.param_names = char(M_.param_names, 'bstar_bar');
M_.param_names_tex = char(M_.param_names_tex, 'bstar\_bar');
M_.param_names_long = char(M_.param_names_long, 'bstar_bar');
M_.param_names = char(M_.param_names, 'lst_bar');
M_.param_names_tex = char(M_.param_names_tex, 'lst\_bar');
M_.param_names_long = char(M_.param_names_long, 'lst_bar');
M_.param_names = char(M_.param_names, 'tauc_bar');
M_.param_names_tex = char(M_.param_names_tex, 'tauc\_bar');
M_.param_names_long = char(M_.param_names_long, 'tauc_bar');
M_.param_names = char(M_.param_names, 'taun_bar');
M_.param_names_tex = char(M_.param_names_tex, 'taun\_bar');
M_.param_names_long = char(M_.param_names_long, 'taun_bar');
M_.param_names = char(M_.param_names, 'tauk_bar');
M_.param_names_tex = char(M_.param_names_tex, 'tauk\_bar');
M_.param_names_long = char(M_.param_names_long, 'tauk_bar');
M_.param_names = char(M_.param_names, 'by_bar');
M_.param_names_tex = char(M_.param_names_tex, 'by\_bar');
M_.param_names_long = char(M_.param_names_long, 'by_bar');
M_.param_names = char(M_.param_names, 'ig_bar');
M_.param_names_tex = char(M_.param_names_tex, 'ig\_bar');
M_.param_names_long = char(M_.param_names_long, 'ig_bar');
M_.param_names = char(M_.param_names, 'kg_bar');
M_.param_names_tex = char(M_.param_names_tex, 'kg\_bar');
M_.param_names_long = char(M_.param_names_long, 'kg_bar');
M_.param_names = char(M_.param_names, 'zy_bar');
M_.param_names_tex = char(M_.param_names_tex, 'zy\_bar');
M_.param_names_long = char(M_.param_names_long, 'zy_bar');
M_.param_names = char(M_.param_names, 'z_bar');
M_.param_names_tex = char(M_.param_names_tex, 'z\_bar');
M_.param_names_long = char(M_.param_names_long, 'z_bar');
M_.param_partitions = struct();
M_.exo_det_nbr = 0;
M_.exo_nbr = 11;
M_.endo_nbr = 50;
M_.param_nbr = 95;
M_.orig_endo_nbr = 50;
M_.aux_vars = [];
M_.Sigma_e = zeros(11, 11);
M_.Correlation_matrix = eye(11, 11);
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 = 1;
erase_compiled_function('baseline_V1_HHlender_static');
erase_compiled_function('baseline_V1_HHlender_dynamic');
M_.orig_eq_nbr = 50;
M_.eq_nbr = 50;
M_.ramsey_eq_nbr = 0;
M_.set_auxiliary_variables = exist(['./' M_.fname '_set_auxiliary_variables.m'], 'file') == 2;
M_.lead_lag_incidence = [
 0 23 0;
 0 24 0;
 1 25 73;
 2 26 74;
 0 27 75;
 0 28 76;
 0 29 0;
 0 30 77;
 0 31 78;
 0 32 79;
 3 33 0;
 4 34 0;
 0 35 0;
 0 36 0;
 0 37 80;
 0 38 81;
 0 39 0;
 0 40 0;
 0 41 0;
 0 42 0;
 5 43 82;
 6 44 0;
 0 45 83;
 0 46 84;
 0 47 0;
 7 48 0;
 0 49 0;
 0 50 85;
 0 51 86;
 0 52 0;
 0 53 87;
 0 54 0;
 8 55 0;
 9 56 0;
 10 57 88;
 11 58 0;
 12 59 0;
 0 60 0;
 13 61 0;
 14 62 0;
 0 63 0;
 15 64 0;
 0 65 0;
 16 66 0;
 17 67 0;
 18 68 0;
 19 69 0;
 20 70 0;
 21 71 0;
 22 72 0;]';
M_.nstatic = 16;
M_.nfwrd   = 12;
M_.npred   = 18;
M_.nboth   = 4;
M_.nsfwrd   = 16;
M_.nspred   = 22;
M_.ndynamic   = 34;
M_.equations_tags = {
};
M_.static_and_dynamic_models_differ = 0;
M_.exo_names_orig_ord = [1:11];
M_.maximum_lag = 1;
M_.maximum_lead = 1;
M_.maximum_endo_lag = 1;
M_.maximum_endo_lead = 1;
oo_.steady_state = zeros(50, 1);
M_.maximum_exo_lag = 0;
M_.maximum_exo_lead = 0;
oo_.exo_steady_state = zeros(11, 1);
M_.params = NaN(95, 1);
M_.NNZDerivatives = [243; -1; -1];
dynare baseline_V1_HHlender
trace_lots = 0; 
close all;
clc; 
M_.params( 1 ) = 0.99;
betta = M_.params( 1 );
M_.params( 3 ) = 0.2;
theta = M_.params( 3 );
M_.params( 5 ) = 0.01425;
delta = M_.params( 5 );
M_.params( 6 ) = 0.01425;
deltaz = M_.params( 6 );
M_.params( 10 ) = 1;
psil = M_.params( 10 );
M_.params( 11 ) = 8.6202;
nut = M_.params( 11 );
M_.params( 12 ) = 9.7666;
nunt = M_.params( 12 );
M_.params( 13 ) = 0.007;
nub = M_.params( 13 );
M_.params( 14 ) = 1.2;
phiw = M_.params( 14 );
M_.params( 16 ) = 0.12;
alphag = M_.params( 16 );
M_.params( 17 ) = 0.01;
psif = M_.params( 17 );
M_.params( 20 ) = 0.4;
zeta = M_.params( 20 );
M_.params( 21 ) = 0.6225;
phip = M_.params( 21 );
M_.params( 22 ) = 0.58;
phib = M_.params( 22 );
chis = 0.5; 
M_.params( 8 ) = 10;
epsp = M_.params( 8 );
M_.params( 24 ) = 0.8382;
rho_a = M_.params( 24 );
M_.params( 25 ) = 0.5414;
rho_ag = M_.params( 25 );
M_.params( 26 ) = 0.6839;
rho_g = M_.params( 26 );
M_.params( 27 ) = 0.7246;
rho_ig = M_.params( 27 );
M_.params( 28 ) = 0.6266;
rho_lg = M_.params( 28 );
M_.params( 29 ) = 0.9095;
rho_rf = M_.params( 29 );
M_.params( 30 ) = 0.91;
rho_r = M_.params( 30 );
M_.params( 31 ) = 0.5;
rho_mu = M_.params( 31 );
M_.params( 32 ) = 0.5;
rho_infstar = M_.params( 32 );
M_.params( 39 ) = 0.2;
gamma_k = M_.params( 39 );
M_.params( 38 ) = 0.2;
gamma_n = M_.params( 38 );
M_.params( 37 ) = 0.2;
gamma_c = M_.params( 37 );
M_.params( 40 ) = 0.2;
gamma_t = M_.params( 40 );
M_.params( 33 ) = 0.5;
rho_tauk = M_.params( 33 );
M_.params( 34 ) = 0.5;
rho_taun = M_.params( 34 );
M_.params( 35 ) = 0.5;
rho_tauc = M_.params( 35 );
M_.params( 36 ) = 0.5;
rho_lst = M_.params( 36 );
iy_bar = 0.21;  
cy_bar = 0.6091; 
igy_bar = 0.035;  
M_.params( 91 ) = 2.8;
by_bar = M_.params( 91 );
M_.params( 94 ) = 0.2;
zy_bar = M_.params( 94 );
M_.params( 56 ) = 1;
infl_bar = M_.params( 56 );
inft_bar = infl_bar;
M_.params( 55 ) = M_.params(56);
infnt_bar = M_.params( 55 );
M_.params( 58 ) = M_.params(56)*1/M_.params(1);
r_bar = M_.params( 58 );
M_.params( 54 ) = 1.0028;
infstar_bar = M_.params( 54 );
rstar_bar = (1/betta)*infstar_bar; 
M_.params( 57 ) = rstar_bar;
rf_bar = M_.params( 57 );
M_.params( 8 ) = 10;
epsp = M_.params( 8 );
M_.params( 60 ) = (M_.params(8)-1)/M_.params(8);
mc_bar = M_.params( 60 );
M_.params( 90 ) = 0.245;
tauk_bar = M_.params( 90 );
M_.params( 89 ) = 0.349;
taun_bar = M_.params( 89 );
M_.params( 88 ) = 0.179;
tauc_bar = M_.params( 88 );
M_.params( 43 ) = 1;
a_bar = M_.params( 43 );
M_.params( 65 ) = 1;
v_bar = M_.params( 65 );
M_.params( 59 ) = 1;
mu_bar = M_.params( 59 );
ltotal_bar = 1/3; 
M_.params( 69 ) = 0.2*ltotal_bar;
lg_bar = M_.params( 69 );
M_.params( 74 ) = ltotal_bar*0.8;
lp_bar = M_.params( 74 );
lttlnt = 0.5; 
M_.params( 70 ) = lttlnt/(1+lttlnt)*M_.params(74);
ltt_bar = M_.params( 70 );
M_.params( 71 ) = M_.params(74)-M_.params(70);
lnt_bar = M_.params( 71 );
M_.params( 2 ) = 0.39;
alphat = M_.params( 2 );
M_.params( 4 ) = 0.44;
alphant = M_.params( 4 );
ktknt = (1-alphant)/(1-alphat)*alphat/alphant*lttlnt;
M_.params( 50 ) = (1/M_.params(1)-(1-M_.params(5))-M_.params(5)*M_.params(90))/(1-M_.params(90));
rnt_bar = M_.params( 50 );
M_.params( 51 ) = M_.params(50);
rt_bar = M_.params( 51 );
itint = ktknt;
ity_bar = iy_bar*ktknt/(1+ktknt);
inty_bar = iy_bar-ity_bar;
eytpnynt = lttlnt*mc_bar*(1-alphant)/(1-alphat); 
eyty = eytpnynt/(1+eytpnynt) 
pnynty = 1-eyty; 
M_.params( 41 ) = eyty;
aratio = M_.params( 41 );
M_.params( 42 ) = pnynty;
bratio = M_.params( 42 );
M_.params( 51 ) = M_.params(41)*M_.params(2)/(ity_bar/M_.params(5));
rt_bar = M_.params( 51 );
M_.params( 45 ) = 1;
y_bar = M_.params( 45 );
M_.params( 92 ) = igy_bar*M_.params(45);
ig_bar = M_.params( 92 );
M_.params( 47 ) = M_.params(45)*cy_bar;
c_bar = M_.params( 47 );
M_.params( 95 ) = M_.params(45)*M_.params(94);
z_bar = M_.params( 95 );
M_.params( 92 ) = igy_bar*M_.params(45);
ig_bar = M_.params( 92 );
M_.params( 93 ) = M_.params(92)/M_.params(6);
kg_bar = M_.params( 93 );
M_.params( 48 ) = M_.params(45)*iy_bar;
i_bar = M_.params( 48 );
M_.params( 49 ) = M_.params(48)/M_.params(5);
k_bar = M_.params( 49 );
M_.params( 62 ) = ity_bar*M_.params(45);
it_bar = M_.params( 62 );
M_.params( 64 ) = M_.params(45)*inty_bar;
int_bar = M_.params( 64 );
M_.params( 63 ) = M_.params(64)/M_.params(5);
knt_bar = M_.params( 63 );
M_.params( 61 ) = M_.params(62)/M_.params(5);
kt_bar = M_.params( 61 );
ytynt = (kt_bar^alphat/knt_bar^alphant)*(ltt_bar^(1-alphat)/lnt_bar^(1-alphant));
yty = ytynt/(1+ytynt);
ynty = 1-yty;
M_.params( 52 ) = M_.params(41)/yty;
et_bar = M_.params( 52 );
M_.params( 53 ) = M_.params(42)/ynty;
pn_bar = M_.params( 53 );
M_.params( 73 ) = M_.params(41)*M_.params(45)/M_.params(52);
yt_bar = M_.params( 73 );
M_.params( 72 ) = M_.params(45)*M_.params(42)/M_.params(53);
ynt_bar = M_.params( 72 );
M_.params( 15 ) = 0.33;
ac = M_.params( 15 );
M_.params( 7 ) = 0.99;
eps = M_.params( 7 );
gy_bar = (ynt_bar/y_bar)/((1-ac)*pn_bar^(-eps))-cy_bar-iy_bar-igy_bar; 
M_.params( 46 ) = M_.params(45)*gy_bar;
g_bar = M_.params( 46 );
(ac*et_bar^(1-eps)+(1-ac)*pn_bar^(1-eps))^(1/(1-eps)); 
M_.params( 44 ) = 0.1;
ag_bar = M_.params( 44 );
M_.params( 18 ) = 0.6;
alphaz = M_.params( 18 );
M_.params( 19 ) = 0.3;
alphas = M_.params( 19 );
M_.params( 81 ) = M_.params(44)*M_.params(93)^M_.params(18)*(0.5*M_.params(69))^(1-M_.params(18));
zg_bar = M_.params( 81 );
M_.params( 82 ) = M_.params(44)*M_.params(46)^M_.params(19)*(0.5*M_.params(69))^(1-M_.params(19));
s_bar = M_.params( 82 );
M_.params( 75 ) = M_.params(52)*(1-M_.params(2))*M_.params(73)/M_.params(70);
w_bar = M_.params( 75 );
M_.params( 76 ) = M_.params(14)*M_.params(75);
wg_bar = M_.params( 76 );
M_.params( 77 ) = 1/((M_.params(47)+M_.params(3)*M_.params(82))*(1+M_.params(88)));
lambda1_bar = M_.params( 77 );
M_.params( 78 ) = M_.params(77);
lambda2_bar = M_.params( 78 );
M_.params( 79 ) = M_.params(77);
lambda3_bar = M_.params( 79 );
M_.params( 23 ) = M_.params(75)*M_.params(77)*(1-M_.params(89))/(M_.params(74)+M_.params(69))^M_.params(10);
chil = M_.params( 23 );
M_.params( 83 ) = M_.params(45)*M_.params(91);
b_bar = M_.params( 83 );
M_.params( 84 ) = 0.58*M_.params(83);
bg_bar = M_.params( 84 );
M_.params( 85 ) = (M_.params(83)-M_.params(84))/M_.params(52);
bgstar_bar = M_.params( 85 );
M_.params( 87 ) = (-(M_.params(84)+M_.params(47)*M_.params(88)+M_.params(89)*(M_.params(74)*M_.params(75)+M_.params(69)*M_.params(76))+M_.params(90)*(M_.params(51)*M_.params(61)+M_.params(50)*M_.params(63))+M_.params(52)*M_.params(85)-M_.params(92)-M_.params(69)*M_.params(76)-M_.params(46)-M_.params(84)*M_.params(58)/M_.params(56)-M_.params(85)*rstar_bar*M_.params(52)/M_.params(54)-M_.params(5)*M_.params(90)*(M_.params(61)+M_.params(63))-M_.params(95)));
lst_bar = M_.params( 87 );
M_.params( 86 ) = (M_.params(45)+M_.params(52)*M_.params(85)-M_.params(46)-M_.params(92)-M_.params(48)-M_.params(47)-M_.params(85)*rstar_bar*M_.params(52)/M_.params(54))/(M_.params(52)*(1-rstar_bar/M_.params(54)));
bstar_bar = M_.params( 86 );
inve_bar = it_bar + int_bar;
tr_bar =  tauc_bar*c_bar + taun_bar*(w_bar*lp_bar + wg_bar*lg_bar) + tauk_bar*(rt_bar*kt_bar + rnt_bar*knt_bar) + lst_bar - tauk_bar*delta*(kt_bar+knt_bar); 
M_.params( 66 ) = 1/(M_.params(47)+M_.params(3)*M_.params(82))*M_.params(60)*M_.params(72)/(1-M_.params(1)*M_.params(21));
x1_bar = M_.params( 66 );
x2_bar = (1/(c_bar+theta*s_bar)*ynt_bar)/(1-phip*betta);
M_.params( 68 ) = 1;
infntopt_bar = M_.params( 68 );
%
% SHOCKS instructions
%
M_.exo_det_length = 0;
M_.Sigma_e(1, 1) = (1)^2;
M_.Sigma_e(2, 2) = (1)^2;
M_.Sigma_e(3, 3) = (1)^2;
M_.Sigma_e(4, 4) = (1)^2;
M_.Sigma_e(5, 5) = (1)^2;
M_.Sigma_e(6, 6) = (1)^2;
M_.Sigma_e(7, 7) = (1)^2;
M_.Sigma_e(8, 8) = (1)^2;
M_.Sigma_e(9, 9) = (1)^2;
M_.Sigma_e(10, 10) = (1)^2;
M_.Sigma_e(11, 11) = (1)^2;
%
% INITVAL instructions
%
options_.initval_file = 0;
oo_.steady_state( 1 ) = log(M_.params(47));
oo_.steady_state( 2 ) = log(M_.params(74));
oo_.steady_state( 3 ) = log(M_.params(62));
oo_.steady_state( 4 ) = log(M_.params(64));
oo_.steady_state( 8 ) = log(M_.params(77));
oo_.steady_state( 9 ) = log(M_.params(78));
oo_.steady_state( 10 ) = log(M_.params(79));
oo_.steady_state( 13 ) = log(M_.params(70));
oo_.steady_state( 14 ) = log(M_.params(71));
oo_.steady_state( 15 ) = log(M_.params(51));
oo_.steady_state( 16 ) = log(M_.params(50));
oo_.steady_state( 17 ) = log(M_.params(75));
oo_.steady_state( 18 ) = log(M_.params(72));
oo_.steady_state( 19 ) = log(M_.params(73));
oo_.steady_state( 20 ) = log(M_.params(45));
oo_.steady_state( 21 ) = log(M_.params(52));
oo_.steady_state( 22 ) = log(M_.params(53));
oo_.steady_state( 23 ) = log(M_.params(56));
oo_.steady_state( 24 ) = log(M_.params(55));
oo_.steady_state( 25 ) = log(M_.params(60));
oo_.steady_state( 32 ) = log(M_.params(81));
oo_.steady_state( 40 ) = log(M_.params(76));
oo_.steady_state( 30 ) = log(M_.params(82));
oo_.steady_state( 7 ) = log(inve_bar);
oo_.steady_state( 33 ) = log(M_.params(88));
oo_.steady_state( 34 ) = log(M_.params(89));
oo_.steady_state( 35 ) = log(M_.params(90));
oo_.steady_state( 36 ) = log(M_.params(86));
oo_.steady_state( 37 ) = log(M_.params(85));
oo_.steady_state( 12 ) = log(rstar_bar);
oo_.steady_state( 38 ) = log(M_.params(83));
oo_.steady_state( 41 ) = log(M_.params(83)/M_.params(45));
oo_.steady_state( 39 ) = log(M_.params(84));
oo_.steady_state( 44 ) = log(M_.params(43));
oo_.steady_state( 45 ) = log(M_.params(44));
oo_.steady_state( 11 ) = log(M_.params(58));
oo_.steady_state( 46 ) = log(M_.params(57));
oo_.steady_state( 47 ) = log(M_.params(59));
oo_.steady_state( 48 ) = log(M_.params(69));
oo_.steady_state( 49 ) = log(M_.params(46));
oo_.steady_state( 50 ) = log(M_.params(92));
oo_.steady_state( 42 ) = log(M_.params(87));
oo_.steady_state( 31 ) = log(M_.params(93));
oo_.steady_state( 5 ) = log(M_.params(61));
oo_.steady_state( 6 ) = log(M_.params(63));
oo_.steady_state( 42 ) = log(M_.params(87));
oo_.steady_state( 43 ) = log(tr_bar);
oo_.steady_state( 26 ) = log(M_.params(65));
oo_.steady_state( 28 ) = log(M_.params(66));
oo_.steady_state( 29 ) = log(x2_bar);
oo_.steady_state( 27 ) = log(M_.params(68));
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
options_.irf = 21;
options_.nograph = 1;
options_.order = 1;
var_list_ = char();
info = stoch_simul(var_list_);
oo_baseline = oo_;
M_baseline = M_;
options_baseline = options_;
shocks_baseline = {'eps_a'; 'eps_mu'; 'eps_lg'; 'eps_g'; 'eps_ig'};
VarsExo = {'a'; 'mu'; 'lg'; 'g'; 'ig'};
NameExo = {'TFP'; 'Risk premium'; 'public employment'; 'government consumption', 'government investment'};
nshocks = length(shocks_baseline);
shock_scale_vec = [1;1;1;1;1];
VarsEndo1 = {'c'; 'lp'; 'inve'; 'bg'; 'bstar'; 'bgstar'; 'y'; 'yt'; 'ynt'; 'ltt'; 'lnt'; ...
'w'; 'it'; 'int'; 'rt'; 'rnt'; 'et'; 'pn'; 'mc'; 'inf'; 'infl'; ...
'ig'; 'zg'; 'b'; 'bg'; 'by'; 'rstar'; 'r'; 'tauc'; 'taun'; 'tauk'};
NameEndo1 = {'consumption'; 'private labor'; 'investment'; 'HH govt bond holdings'; 'HH intl lending'; 'government international borrowing'; ...
'output'; 'tradables'; 'non-tradables'; 'labor: NT' ; 'labor: T'; ...
'wage'; 'T inv'; 'NT inv'; 'capital rental: T'; 'capital rental: NT'; 'real exchange rate'; 'NT real price'; 'mc'; 'NT inflation'; 'inflation'; ...
'public investment'; 'infrastructure'; ...
'total govt debt'; 'govt intl borrowing'; 'debt ratio'; 'int rate on foreign debt'; 'domestic int rate'; 'consumption tax rate'; 'labor tax rate'; 'capital tax rate'}; 
nvars1 = length(VarsEndo1); 
lirf = options_.irf;
for jj = 1:nshocks;
m_scale = shock_scale_vec(jj); 
figure('Name',char(NameExo{jj}),'NumberTitle','off');     
for ii = 1:nvars1
var1 = [char(VarsEndo1{ii}),'_',char(shocks_baseline{jj})];
temp1 = eval(['oo_baseline.irfs.',var1]) *100;
subplot(3, 4, ii);
plot(0:lirf-1, m_scale*temp1, '-b', 'LineWidth', 1.5);
title(char(NameEndo1{ii}));         
end
end
save('baseline_V1_HHlender_results.mat', 'oo_', 'M_', 'options_');
if exist('estim_params_', 'var') == 1
  save('baseline_V1_HHlender_results.mat', 'estim_params_', '-append');
end
if exist('bayestopt_', 'var') == 1
  save('baseline_V1_HHlender_results.mat', 'bayestopt_', '-append');
end
if exist('dataset_', 'var') == 1
  save('baseline_V1_HHlender_results.mat', 'dataset_', '-append');
end
if exist('estimation_info', 'var') == 1
  save('baseline_V1_HHlender_results.mat', 'estimation_info', '-append');
end
if exist('dataset_info', 'var') == 1
  save('baseline_V1_HHlender_results.mat', 'dataset_info', '-append');
end
if exist('oo_recursive_', 'var') == 1
  save('baseline_V1_HHlender_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
