%
% 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 = 'BGG_financial_3_p45_phi2';
%
% Some global variables initialization
%
global_initialization;
diary off;
diary('BGG_financial_3_p45_phi2.log');
M_.exo_names = 'epsA';
M_.exo_names_tex = 'epsA';
M_.exo_names_long = 'epsA';
M_.exo_names = char(M_.exo_names, 'epsG');
M_.exo_names_tex = char(M_.exo_names_tex, 'epsG');
M_.exo_names_long = char(M_.exo_names_long, 'epsG');
M_.exo_names = char(M_.exo_names, 'epsMCS');
M_.exo_names_tex = char(M_.exo_names_tex, 'epsMCS');
M_.exo_names_long = char(M_.exo_names_long, 'epsMCS');
M_.exo_names = char(M_.exo_names, 'epsMPS');
M_.exo_names_tex = char(M_.exo_names_tex, 'epsMPS');
M_.exo_names_long = char(M_.exo_names_long, 'epsMPS');
M_.exo_names = char(M_.exo_names, 'epsIS');
M_.exo_names_tex = char(M_.exo_names_tex, 'epsIS');
M_.exo_names_long = char(M_.exo_names_long, 'epsIS');
M_.exo_names = char(M_.exo_names, 'epsMRSS');
M_.exo_names_tex = char(M_.exo_names_tex, 'epsMRSS');
M_.exo_names_long = char(M_.exo_names_long, 'epsMRSS');
M_.exo_names = char(M_.exo_names, 'epsRPS');
M_.exo_names_tex = char(M_.exo_names_tex, 'epsRPS');
M_.exo_names_long = char(M_.exo_names_long, 'epsRPS');
M_.endo_names = 'V';
M_.endo_names_tex = 'V';
M_.endo_names_long = 'V';
M_.endo_names = char(M_.endo_names, 'Vmodified');
M_.endo_names_tex = char(M_.endo_names_tex, 'Vmodified');
M_.endo_names_long = char(M_.endo_names_long, 'Vmodified');
M_.endo_names = char(M_.endo_names, 'U');
M_.endo_names_tex = char(M_.endo_names_tex, 'U');
M_.endo_names_long = char(M_.endo_names_long, 'U');
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, '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, '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, '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, 'X');
M_.endo_names_tex = char(M_.endo_names_tex, 'X');
M_.endo_names_long = char(M_.endo_names_long, 'X');
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, 'IS');
M_.endo_names_tex = char(M_.endo_names_tex, 'IS');
M_.endo_names_long = char(M_.endo_names_long, 'IS');
M_.endo_names = char(M_.endo_names, 'Sdash');
M_.endo_names_tex = char(M_.endo_names_tex, 'Sdash');
M_.endo_names_long = char(M_.endo_names_long, 'Sdash');
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, 'Wh');
M_.endo_names_tex = char(M_.endo_names_tex, 'Wh');
M_.endo_names_long = char(M_.endo_names_long, 'Wh');
M_.endo_names = char(M_.endo_names, 'rK');
M_.endo_names_tex = char(M_.endo_names_tex, 'rK');
M_.endo_names_long = char(M_.endo_names_long, 'rK');
M_.endo_names = char(M_.endo_names, 'RPS');
M_.endo_names_tex = char(M_.endo_names_tex, 'RPS');
M_.endo_names_long = char(M_.endo_names_long, 'RPS');
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, '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, 'Rn');
M_.endo_names_tex = char(M_.endo_names_tex, 'Rn');
M_.endo_names_long = char(M_.endo_names_long, 'Rn');
M_.endo_names = char(M_.endo_names, 'PIE');
M_.endo_names_tex = char(M_.endo_names_tex, 'PIE');
M_.endo_names_long = char(M_.endo_names_long, 'PIE');
M_.endo_names = char(M_.endo_names, 'PIEw');
M_.endo_names_tex = char(M_.endo_names_tex, 'PIEw');
M_.endo_names_long = char(M_.endo_names_long, 'PIEw');
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, 'Jw');
M_.endo_names_tex = char(M_.endo_names_tex, 'Jw');
M_.endo_names_long = char(M_.endo_names_long, 'Jw');
M_.endo_names = char(M_.endo_names, 'Hd');
M_.endo_names_tex = char(M_.endo_names_tex, 'Hd');
M_.endo_names_long = char(M_.endo_names_long, 'Hd');
M_.endo_names = char(M_.endo_names, 'MRSS');
M_.endo_names_tex = char(M_.endo_names_tex, 'MRSS');
M_.endo_names_long = char(M_.endo_names_long, 'MRSS');
M_.endo_names = char(M_.endo_names, 'JJw');
M_.endo_names_tex = char(M_.endo_names_tex, 'JJw');
M_.endo_names_long = char(M_.endo_names_long, 'JJw');
M_.endo_names = char(M_.endo_names, 'WO');
M_.endo_names_tex = char(M_.endo_names_tex, 'WO');
M_.endo_names_long = char(M_.endo_names_long, 'WO');
M_.endo_names = char(M_.endo_names, 'Deltaw');
M_.endo_names_tex = char(M_.endo_names_tex, 'Deltaw');
M_.endo_names_long = char(M_.endo_names_long, 'Deltaw');
M_.endo_names = char(M_.endo_names, 'YW');
M_.endo_names_tex = char(M_.endo_names_tex, 'YW');
M_.endo_names_long = char(M_.endo_names_long, 'YW');
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, 'PWP');
M_.endo_names_tex = char(M_.endo_names_tex, 'PWP');
M_.endo_names_long = char(M_.endo_names_long, 'PWP');
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, 'Jp');
M_.endo_names_tex = char(M_.endo_names_tex, 'Jp');
M_.endo_names_long = char(M_.endo_names_long, 'Jp');
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, 'MCS');
M_.endo_names_tex = char(M_.endo_names_tex, 'MCS');
M_.endo_names_long = char(M_.endo_names_long, 'MCS');
M_.endo_names = char(M_.endo_names, 'JJp');
M_.endo_names_tex = char(M_.endo_names_tex, 'JJp');
M_.endo_names_long = char(M_.endo_names_long, 'JJp');
M_.endo_names = char(M_.endo_names, 'PO');
M_.endo_names_tex = char(M_.endo_names_tex, 'PO');
M_.endo_names_long = char(M_.endo_names_long, 'PO');
M_.endo_names = char(M_.endo_names, 'Deltap');
M_.endo_names_tex = char(M_.endo_names_tex, 'Deltap');
M_.endo_names_long = char(M_.endo_names_long, 'Deltap');
M_.endo_names = char(M_.endo_names, 'MPS');
M_.endo_names_tex = char(M_.endo_names_tex, 'MPS');
M_.endo_names_long = char(M_.endo_names_long, 'MPS');
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, 'RK');
M_.endo_names_tex = char(M_.endo_names_tex, 'RK');
M_.endo_names_long = char(M_.endo_names_long, 'RK');
M_.endo_names = char(M_.endo_names, 'spread');
M_.endo_names_tex = char(M_.endo_names_tex, 'spread');
M_.endo_names_long = char(M_.endo_names_long, 'spread');
M_.endo_names = char(M_.endo_names, 'CEquiv');
M_.endo_names_tex = char(M_.endo_names_tex, 'CEquiv');
M_.endo_names_long = char(M_.endo_names_long, 'CEquiv');
M_.endo_names = char(M_.endo_names, 'rhoofpsi');
M_.endo_names_tex = char(M_.endo_names_tex, 'rhoofpsi');
M_.endo_names_long = char(M_.endo_names_long, 'rhoofpsi');
M_.endo_names = char(M_.endo_names, 'psiFF');
M_.endo_names_tex = char(M_.endo_names_tex, 'psiFF');
M_.endo_names_long = char(M_.endo_names_long, 'psiFF');
M_.endo_names = char(M_.endo_names, 'NE');
M_.endo_names_tex = char(M_.endo_names_tex, 'NE');
M_.endo_names_long = char(M_.endo_names_long, 'NE');
M_.endo_names = char(M_.endo_names, 'CE');
M_.endo_names_tex = char(M_.endo_names_tex, 'CE');
M_.endo_names_long = char(M_.endo_names_long, 'CE');
M_.endo_names = char(M_.endo_names, 'L');
M_.endo_names_tex = char(M_.endo_names_tex, 'L');
M_.endo_names_long = char(M_.endo_names_long, 'L');
M_.endo_names = char(M_.endo_names, 'Rl');
M_.endo_names_tex = char(M_.endo_names_tex, 'Rl');
M_.endo_names_long = char(M_.endo_names_long, 'Rl');
M_.endo_names = char(M_.endo_names, 'p');
M_.endo_names_tex = char(M_.endo_names_tex, 'p');
M_.endo_names_long = char(M_.endo_names_long, 'p');
M_.endo_names = char(M_.endo_names, 'Gofpsi');
M_.endo_names_tex = char(M_.endo_names_tex, 'Gofpsi');
M_.endo_names_long = char(M_.endo_names_long, 'Gofpsi');
M_.endo_names = char(M_.endo_names, 'Gammaofpsi');
M_.endo_names_tex = char(M_.endo_names_tex, 'Gammaofpsi');
M_.endo_names_long = char(M_.endo_names_long, 'Gammaofpsi');
M_.endo_names = char(M_.endo_names, 'Gammadashofpsi');
M_.endo_names_tex = char(M_.endo_names_tex, 'Gammadashofpsi');
M_.endo_names_long = char(M_.endo_names_long, 'Gammadashofpsi');
M_.endo_names = char(M_.endo_names, 'Gdashofpsi');
M_.endo_names_tex = char(M_.endo_names_tex, 'Gdashofpsi');
M_.endo_names_long = char(M_.endo_names_long, 'Gdashofpsi');
M_.endo_names = char(M_.endo_names, 'sigmaE');
M_.endo_names_tex = char(M_.endo_names_tex, 'sigmaE');
M_.endo_names_long = char(M_.endo_names_long, 'sigmaE');
M_.endo_names = char(M_.endo_names, 'xiE');
M_.endo_names_tex = char(M_.endo_names_tex, 'xiE');
M_.endo_names_long = char(M_.endo_names_long, 'xiE');
M_.endo_names = char(M_.endo_names, 'muFF');
M_.endo_names_tex = char(M_.endo_names_tex, 'muFF');
M_.endo_names_long = char(M_.endo_names_long, 'muFF');
M_.endo_names = char(M_.endo_names, 'sigmapsi');
M_.endo_names_tex = char(M_.endo_names_tex, 'sigmapsi');
M_.endo_names_long = char(M_.endo_names_long, 'sigmapsi');
M_.endo_names = char(M_.endo_names, 'dyobs');
M_.endo_names_tex = char(M_.endo_names_tex, 'dyobs');
M_.endo_names_long = char(M_.endo_names_long, 'dyobs');
M_.endo_names = char(M_.endo_names, 'dcobs');
M_.endo_names_tex = char(M_.endo_names_tex, 'dcobs');
M_.endo_names_long = char(M_.endo_names_long, 'dcobs');
M_.endo_names = char(M_.endo_names, 'pinfobs');
M_.endo_names_tex = char(M_.endo_names_tex, 'pinfobs');
M_.endo_names_long = char(M_.endo_names_long, 'pinfobs');
M_.endo_names = char(M_.endo_names, 'robs');
M_.endo_names_tex = char(M_.endo_names_tex, 'robs');
M_.endo_names_long = char(M_.endo_names_long, 'robs');
M_.endo_names = char(M_.endo_names, 'labobs');
M_.endo_names_tex = char(M_.endo_names_tex, 'labobs');
M_.endo_names_long = char(M_.endo_names_long, 'labobs');
M_.endo_names = char(M_.endo_names, 'spreadobs');
M_.endo_names_tex = char(M_.endo_names_tex, 'spreadobs');
M_.endo_names_long = char(M_.endo_names_long, 'spreadobs');
M_.endo_names = char(M_.endo_names, 'dinvobs');
M_.endo_names_tex = char(M_.endo_names_tex, 'dinvobs');
M_.endo_names_long = char(M_.endo_names_long, 'dinvobs');
M_.endo_names = char(M_.endo_names, 'dwobs');
M_.endo_names_tex = char(M_.endo_names_tex, 'dwobs');
M_.endo_names_long = char(M_.endo_names_long, 'dwobs');
M_.endo_names = char(M_.endo_names, 'UH');
M_.endo_names_tex = char(M_.endo_names_tex, 'UH');
M_.endo_names_long = char(M_.endo_names_long, 'UH');
M_.param_names = 'gy';
M_.param_names_tex = 'gy';
M_.param_names_long = 'gy';
M_.param_names = char(M_.param_names, 'alp');
M_.param_names_tex = char(M_.param_names_tex, 'alp');
M_.param_names_long = char(M_.param_names_long, 'alp');
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, 'zzetap');
M_.param_names_tex = char(M_.param_names_tex, 'zzetap');
M_.param_names_long = char(M_.param_names_long, 'zzetap');
M_.param_names = char(M_.param_names, 'zzetaw');
M_.param_names_tex = char(M_.param_names_tex, 'zzetaw');
M_.param_names_long = char(M_.param_names_long, 'zzetaw');
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, 'chi');
M_.param_names_tex = char(M_.param_names_tex, 'chi');
M_.param_names_long = char(M_.param_names_long, 'chi');
M_.param_names = char(M_.param_names, 'phiX');
M_.param_names_tex = char(M_.param_names_tex, 'phiX');
M_.param_names_long = char(M_.param_names_long, 'phiX');
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, '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, 'gammap');
M_.param_names_tex = char(M_.param_names_tex, 'gammap');
M_.param_names_long = char(M_.param_names_long, 'gammap');
M_.param_names = char(M_.param_names, 'gammaw');
M_.param_names_tex = char(M_.param_names_tex, 'gammaw');
M_.param_names_long = char(M_.param_names_long, 'gammaw');
M_.param_names = char(M_.param_names, 'PIE_bar');
M_.param_names_tex = char(M_.param_names_tex, 'PIE\_bar');
M_.param_names_long = char(M_.param_names_long, 'PIE_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, 'Rn_bar');
M_.param_names_tex = char(M_.param_names_tex, 'Rn\_bar');
M_.param_names_long = char(M_.param_names_long, 'Rn_bar');
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, 'theta_pie');
M_.param_names_tex = char(M_.param_names_tex, 'theta\_pie');
M_.param_names_long = char(M_.param_names_long, 'theta_pie');
M_.param_names = char(M_.param_names, 'theta_y');
M_.param_names_tex = char(M_.param_names_tex, 'theta\_y');
M_.param_names_long = char(M_.param_names_long, 'theta_y');
M_.param_names = char(M_.param_names, 'theta_dy');
M_.param_names_tex = char(M_.param_names_tex, 'theta\_dy');
M_.param_names_long = char(M_.param_names_long, 'theta_dy');
M_.param_names = char(M_.param_names, 'rhoA');
M_.param_names_tex = char(M_.param_names_tex, 'rhoA');
M_.param_names_long = char(M_.param_names_long, 'rhoA');
M_.param_names = char(M_.param_names, 'rhoG');
M_.param_names_tex = char(M_.param_names_tex, 'rhoG');
M_.param_names_long = char(M_.param_names_long, 'rhoG');
M_.param_names = char(M_.param_names, 'rhoMCS');
M_.param_names_tex = char(M_.param_names_tex, 'rhoMCS');
M_.param_names_long = char(M_.param_names_long, 'rhoMCS');
M_.param_names = char(M_.param_names, 'rhoMRSS');
M_.param_names_tex = char(M_.param_names_tex, 'rhoMRSS');
M_.param_names_long = char(M_.param_names_long, 'rhoMRSS');
M_.param_names = char(M_.param_names, 'rhoMPS');
M_.param_names_tex = char(M_.param_names_tex, 'rhoMPS');
M_.param_names_long = char(M_.param_names_long, 'rhoMPS');
M_.param_names = char(M_.param_names, 'rhoRPS');
M_.param_names_tex = char(M_.param_names_tex, 'rhoRPS');
M_.param_names_long = char(M_.param_names_long, 'rhoRPS');
M_.param_names = char(M_.param_names, 'rhoIS');
M_.param_names_tex = char(M_.param_names_tex, 'rhoIS');
M_.param_names_long = char(M_.param_names_long, 'rhoIS');
M_.param_names = char(M_.param_names, 'ptarget');
M_.param_names_tex = char(M_.param_names_tex, 'ptarget');
M_.param_names_long = char(M_.param_names_long, 'ptarget');
M_.param_names = char(M_.param_names, 'rhotarget');
M_.param_names_tex = char(M_.param_names_tex, 'rhotarget');
M_.param_names_long = char(M_.param_names_long, 'rhotarget');
M_.param_names = char(M_.param_names, 'phitarget');
M_.param_names_tex = char(M_.param_names_tex, 'phitarget');
M_.param_names_long = char(M_.param_names_long, 'phitarget');
M_.param_names = char(M_.param_names, 'CEtarget');
M_.param_names_tex = char(M_.param_names_tex, 'CEtarget');
M_.param_names_long = char(M_.param_names_long, 'CEtarget');
M_.param_names = char(M_.param_names, 'sigma_s');
M_.param_names_tex = char(M_.param_names_tex, 'sigma\_s');
M_.param_names_long = char(M_.param_names_long, 'sigma_s');
M_.param_names = char(M_.param_names, 'wr');
M_.param_names_tex = char(M_.param_names_tex, 'wr');
M_.param_names_long = char(M_.param_names_long, 'wr');
M_.param_names = char(M_.param_names, 'psi');
M_.param_names_tex = char(M_.param_names_tex, 'psi');
M_.param_names_long = char(M_.param_names_long, 'psi');
M_.exo_det_nbr = 0;
M_.exo_nbr = 7;
M_.endo_nbr = 66;
M_.param_nbr = 33;
M_.orig_endo_nbr = 66;
M_.aux_vars = [];
options_.varobs = [];
options_.varobs = 'dyobs';
options_.varobs = char(options_.varobs, 'dcobs');
options_.varobs = char(options_.varobs, 'dinvobs');
options_.varobs = char(options_.varobs, 'dwobs');
options_.varobs = char(options_.varobs, 'robs');
options_.varobs = char(options_.varobs, 'pinfobs');
options_.varobs = char(options_.varobs, 'labobs');
options_.varobs_id = [ 58 59 64 65 61 60 62  ];
M_.Sigma_e = zeros(7, 7);
M_.Correlation_matrix = eye(7, 7);
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('BGG_financial_3_p45_phi2_static');
erase_compiled_function('BGG_financial_3_p45_phi2_dynamic');
M_.lead_lag_incidence = [
 0 20 85;
 0 21 86;
 0 22 87;
 1 23 88;
 0 24 0;
 2 25 0;
 0 26 0;
 0 27 89;
 3 28 0;
 4 29 90;
 0 30 91;
 0 31 92;
 0 32 0;
 0 33 0;
 5 34 0;
 0 35 93;
 6 36 94;
 7 37 0;
 8 38 95;
 0 39 96;
 9 40 0;
 0 41 97;
 0 42 0;
 10 43 0;
 0 44 98;
 0 45 0;
 11 46 0;
 0 47 0;
 12 48 0;
 0 49 0;
 0 50 0;
 0 51 99;
 13 52 0;
 14 53 0;
 0 54 100;
 0 55 0;
 15 56 0;
 16 57 0;
 17 58 0;
 0 59 101;
 0 60 0;
 0 61 102;
 0 62 103;
 0 63 0;
 0 64 0;
 0 65 0;
 18 66 0;
 19 0 0;
 0 67 0;
 0 68 104;
 0 69 105;
 0 70 0;
 0 71 0;
 0 72 0;
 0 73 0;
 0 74 0;
 0 75 0;
 0 76 0;
 0 77 0;
 0 78 0;
 0 79 0;
 0 80 0;
 0 81 0;
 0 82 0;
 0 83 0;
 0 84 0;]';
M_.nstatic = 30;
M_.nfwrd   = 17;
M_.npred   = 15;
M_.nboth   = 4;
M_.nsfwrd   = 21;
M_.nspred   = 19;
M_.ndynamic   = 36;
M_.equations_tags = {
  1 , 'name' , '(226) Consumption equivalent variation' ;
  2 , 'name' , '(185) Intertemporal welfare' ;
  3 , 'name' , '(185b) Modified intertemporal welfare' ;
  4 , 'name' , '(186) Flow utility' ;
  5 , 'name' , '(187) Law of motion for capital' ;
  6 , 'name' , '(188) Investment' ;
  7 , 'name' , '(189) Capital adjustment cost' ;
  8 , 'name' , '(190) Capital adjustment cost' ;
  9 , 'name' , '(191) MUC' ;
  10 , 'name' , '(192) Household foc for hours' ;
  11 , 'name' , '(193) Euler equation' ;
  12 , 'name' , '(195) Household foc for investment' ;
  13 , 'name' , '(197) Fischer equation' ;
  15 , 'name' , '(198) Wage inflation' ;
  16 , 'name' , '(199) Wage setting condition I' ;
  17 , 'name' , '(200) Wage setting condition II' ;
  18 , 'name' , '(201) Wage setting condition III' ;
  19 , 'name' , '(202) Law of motion for wage level' ;
  20 , 'name' , '(203) Law of motion for wage dispersion' ;
  21 , 'name' , '(204) Production function' ;
  22 , 'name' , '(205) Firm for for labour' ;
  23 , 'name' , '(206) Firm foc: return on capital' ;
  24 , 'name' , '(207) Marginal cost' ;
  25 , 'name' , '(208) Price setting condition' ;
  26 , 'name' , '(209) Price setting condition' ;
  27 , 'name' , '(210) Price setting condition' ;
  28 , 'name' , '(211) Law of motion for price level' ;
  29 , 'name' , '(212) Law of motion for price dispersion' ;
  30 , 'name' , '(213) Taylor rule' ;
  31 , 'name' , '(215) Market clearing for labour)' ;
  32 , 'name' , '(216) Market clearing for goods' ;
  33 , 'name' , '(217) Gross real return on capital' ;
  34 , 'name' , '(218)' ;
  35 , 'name' , '(439)' ;
  36 , 'name' , '(440)' ;
  37 , 'name' , '(441)' ;
  38 , 'name' , '(442)' ;
  39 , 'name' , '(443) Resource constraint' ;
  40 , 'name' , '(444)' ;
  41 , 'name' , '(445)' ;
  42 , 'name' , '(446)' ;
  43 , 'name' , '(446)' ;
  44 , 'name' , '(451)' ;
  45 , 'name' , '(452)' ;
  46 , 'name' , '(460)' ;
  47 , 'name' , '(454)' ;
  48 , 'name' , '(219) Law of motion for productivity' ;
  49 , 'name' , '(220) Law of motion for government expenditure' ;
  50 , 'name' , '(221) Law of motion for price markup shock' ;
  51 , 'name' , '(232) Law of motion for wage markup shock' ;
  52 , 'name' , '(233) Law of motion for capital quality shock' ;
  53 , 'name' , '(234) Law of motion for monetary policy shock' ;
  54 , 'name' , '(235) Law of motion for risk premium shock' ;
  55 , 'name' , '(264) Measurement equation for output growth' ;
  56 , 'name' , '(265) Measurement equation for consumption growth' ;
  57 , 'name' , '(266) Measurement equation for investment growth' ;
  58 , 'name' , '(267) Measurement equation for real wage growth' ;
  59 , 'name' , '(268) Measurement equation for hours' ;
  60 , 'name' , '(269) Measurement equation for policy rate' ;
  61 , 'name' , '(270) Measurement equation for inflation' ;
  62 , 'name' , '(271) Measurement equation for spread' ;
};
M_.static_and_dynamic_models_differ = 0;
M_.exo_names_orig_ord = [1:7];
M_.maximum_lag = 1;
M_.maximum_lead = 1;
M_.maximum_endo_lag = 1;
M_.maximum_endo_lead = 1;
oo_.steady_state = zeros(66, 1);
M_.maximum_exo_lag = 0;
M_.maximum_exo_lead = 0;
oo_.exo_steady_state = zeros(7, 1);
M_.params = NaN(33, 1);
M_.NNZDerivatives = zeros(3, 1);
M_.NNZDerivatives(1) = 261;
M_.NNZDerivatives(2) = 524;
M_.NNZDerivatives(3) = -1;
M_.params( 1 ) = 0.18;
gy = M_.params( 1 );
M_.params( 2 ) = 0.67;
alp = M_.params( 2 );
M_.params( 3 ) = 0.0250;
delta = M_.params( 3 );
M_.params( 4 ) = 7.0;
zzetap = M_.params( 4 );
M_.params( 5 ) = 3.0;
zzetaw = M_.params( 5 );
M_.params( 13 ) = 0.8662984461054001;
PIE_bar = M_.params( 13 );
M_.params( 14 ) = 0.3550693291143;
g_bar = M_.params( 14 );
M_.params( 15 ) = 1.3318469551282;
Rn_bar = M_.params( 15 );
M_.params( 27 ) = 0.03;
ptarget = M_.params( 27 );
M_.params( 28 ) = 1.006901;
rhotarget = M_.params( 28 );
M_.params( 29 ) = 2.0;
phitarget = M_.params( 29 );
M_.params( 30 ) = 0.075;
CEtarget = M_.params( 30 );
M_.params( 6 ) = 1.5;
sigma_c = M_.params( 6 );
M_.params( 7 ) = 0.7;
chi = M_.params( 7 );
M_.params( 8 ) = 2.0;
phiX = M_.params( 8 );
M_.params( 9 ) = 0.7;
xip = M_.params( 9 );
M_.params( 10 ) = 0.7;
xiw = M_.params( 10 );
M_.params( 11 ) = 0.5;
gammap = M_.params( 11 );
M_.params( 12 ) = 0.5;
gammaw = M_.params( 12 );
M_.params( 33 ) = 2;
psi = M_.params( 33 );
M_.params( 16 ) = 0.7;
rho_r = M_.params( 16 );
M_.params( 17 ) = 1.5;
theta_pie = M_.params( 17 );
M_.params( 18 ) = 0.3;
theta_y = M_.params( 18 );
M_.params( 19 ) = 0.3;
theta_dy = M_.params( 19 );
M_.params( 20 ) = 0.7;
rhoA = M_.params( 20 );
M_.params( 21 ) = 0.7;
rhoG = M_.params( 21 );
M_.params( 22 ) = 0.7;
rhoMCS = M_.params( 22 );
M_.params( 23 ) = 0.7;
rhoMRSS = M_.params( 23 );
M_.params( 24 ) = 0.7;
rhoMPS = M_.params( 24 );
M_.params( 25 ) = 0.7;
rhoRPS = M_.params( 25 );
M_.params( 26 ) = 0.7;
rhoIS = M_.params( 26 );
M_.params( 31 ) = 0.01;
sigma_s = M_.params( 31 );
M_.params( 32 ) = 0.0;
wr = M_.params( 32 );
%
% SHOCKS instructions
%
make_ex_;
M_.exo_det_length = 0;
M_.Sigma_e(1, 1) = (M_.params(31))^2;
M_.Sigma_e(2, 2) = (M_.params(31))^2;
M_.Sigma_e(3, 3) = (M_.params(31))^2;
M_.Sigma_e(4, 4) = (M_.params(31))^2;
M_.Sigma_e(5, 5) = (M_.params(31))^2;
M_.Sigma_e(6, 6) = (M_.params(31))^2;
M_.Sigma_e(7, 7) = (M_.params(31))^2;
steady;
oo_.dr.eigval = check(M_,options_,oo_);
global estim_params_
estim_params_.var_exo = [];
estim_params_.var_endo = [];
estim_params_.corrx = [];
estim_params_.corrn = [];
estim_params_.param_vals = [];
estim_params_.param_vals = [estim_params_.param_vals; 6, NaN, (-Inf), Inf, 3, 1.50, 0.375, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 33, NaN, 0, NaN, 3, 2.00, 0.75, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 7, NaN, (-Inf), Inf, 1, 0.50, 0.10, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 8, NaN, (-Inf), Inf, 3, 2.00, 0.75, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 9, NaN, (-Inf), Inf, 1, 0.50, 0.10, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 10, NaN, (-Inf), Inf, 1, 0.50, 0.10, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 11, NaN, (-Inf), Inf, 1, 0.50, 0.10, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 12, NaN, (-Inf), Inf, 1, 0.50, 0.10, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 16, NaN, (-Inf), Inf, 1, 0.75, 0.10, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 17, NaN, (-Inf), Inf, 3, 1.50, 0.25, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 18, NaN, (-Inf), Inf, 3, 0.12, 0.05, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 19, NaN, (-Inf), Inf, 3, 0.12, 0.05, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 20, NaN, (-Inf), Inf, 1, 0.50, 0.20, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 21, NaN, (-Inf), Inf, 1, 0.50, 0.20, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 22, NaN, (-Inf), Inf, 1, 0.50, 0.20, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 23, NaN, (-Inf), Inf, 1, 0.50, 0.20, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 24, NaN, (-Inf), Inf, 1, 0.50, 0.20, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 25, NaN, (-Inf), Inf, 1, 0.50, 0.20, NaN, NaN, NaN ];
estim_params_.param_vals = [estim_params_.param_vals; 26, NaN, (-Inf), Inf, 1, 0.50, 0.20, NaN, NaN, NaN ];
estim_params_.var_exo = [estim_params_.var_exo; 1, NaN, (-Inf), Inf, 4, 0.001, 0.020, NaN, NaN, NaN ];
estim_params_.var_exo = [estim_params_.var_exo; 2, NaN, (-Inf), Inf, 4, 0.001, 0.020, NaN, NaN, NaN ];
estim_params_.var_exo = [estim_params_.var_exo; 3, NaN, (-Inf), Inf, 4, 0.001, 0.020, NaN, NaN, NaN ];
estim_params_.var_exo = [estim_params_.var_exo; 6, NaN, (-Inf), Inf, 4, 0.001, 0.020, NaN, NaN, NaN ];
estim_params_.var_exo = [estim_params_.var_exo; 4, NaN, (-Inf), Inf, 4, 0.001, 0.020, NaN, NaN, NaN ];
estim_params_.var_exo = [estim_params_.var_exo; 7, NaN, (-Inf), Inf, 4, 0.001, 0.020, NaN, NaN, NaN ];
estim_params_.var_exo = [estim_params_.var_exo; 5, NaN, (-Inf), Inf, 4, 0.001, 0.020, NaN, NaN, NaN ];
options_.first_obs = 75;
options_.lik_init = 1;
options_.mh_drop = 0.25;
options_.mh_jscale = 0.4;
options_.mh_nblck = 2;
options_.mh_replic = 500000;
options_.mode_compute = 6;
options_.plot_priors = 0;
options_.prefilter = 0;
options_.presample = 0;
options_.datafile = 'usdata_47q1_17q4_unscaled';
options_.nobs = 208;
options_.order = 1;
var_list_=[];
dynare_estimation(var_list_);
save('BGG_financial_3_p45_phi2_results.mat', 'oo_', 'M_', 'options_');
if exist('estim_params_', 'var') == 1
  save('BGG_financial_3_p45_phi2_results.mat', 'estim_params_', '-append');
end
if exist('bayestopt_', 'var') == 1
  save('BGG_financial_3_p45_phi2_results.mat', 'bayestopt_', '-append');
end
if exist('dataset_', 'var') == 1
  save('BGG_financial_3_p45_phi2_results.mat', 'dataset_', '-append');
end
if exist('estimation_info', 'var') == 1
  save('BGG_financial_3_p45_phi2_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
