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

tic;
global M_ oo_ options_ ys0_ ex0_ estimation_info
options_ = [];
M_.fname = 'usmodel';
%
% Some global variables initialization
%
global_initialization;
diary off;
M_.exo_names = 'ea';
M_.exo_names_tex = 'ea';
M_.exo_names_long = 'ea';
M_.exo_names = char(M_.exo_names, 'eb');
M_.exo_names_tex = char(M_.exo_names_tex, 'eb');
M_.exo_names_long = char(M_.exo_names_long, 'eb');
M_.exo_names = char(M_.exo_names, 'eg');
M_.exo_names_tex = char(M_.exo_names_tex, 'eg');
M_.exo_names_long = char(M_.exo_names_long, 'eg');
M_.exo_names = char(M_.exo_names, 'eqs');
M_.exo_names_tex = char(M_.exo_names_tex, 'eqs');
M_.exo_names_long = char(M_.exo_names_long, 'eqs');
M_.exo_names = char(M_.exo_names, 'em');
M_.exo_names_tex = char(M_.exo_names_tex, 'em');
M_.exo_names_long = char(M_.exo_names_long, 'em');
M_.exo_names = char(M_.exo_names, 'epinf');
M_.exo_names_tex = char(M_.exo_names_tex, 'epinf');
M_.exo_names_long = char(M_.exo_names_long, 'epinf');
M_.exo_names = char(M_.exo_names, 'ew');
M_.exo_names_tex = char(M_.exo_names_tex, 'ew');
M_.exo_names_long = char(M_.exo_names_long, 'ew');
M_.endo_names = 'labobs';
M_.endo_names_tex = 'labobs';
M_.endo_names_long = 'labobs';
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, '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, 'dy');
M_.endo_names_tex = char(M_.endo_names_tex, 'dy');
M_.endo_names_long = char(M_.endo_names_long, 'dy');
M_.endo_names = char(M_.endo_names, 'dc');
M_.endo_names_tex = char(M_.endo_names_tex, 'dc');
M_.endo_names_long = char(M_.endo_names_long, 'dc');
M_.endo_names = char(M_.endo_names, 'dinve');
M_.endo_names_tex = char(M_.endo_names_tex, 'dinve');
M_.endo_names_long = char(M_.endo_names_long, 'dinve');
M_.endo_names = char(M_.endo_names, 'dw');
M_.endo_names_tex = char(M_.endo_names_tex, 'dw');
M_.endo_names_long = char(M_.endo_names_long, 'dw');
M_.endo_names = char(M_.endo_names, 'ewma');
M_.endo_names_tex = char(M_.endo_names_tex, 'ewma');
M_.endo_names_long = char(M_.endo_names_long, 'ewma');
M_.endo_names = char(M_.endo_names, 'epinfma');
M_.endo_names_tex = char(M_.endo_names_tex, 'epinfma');
M_.endo_names_long = char(M_.endo_names_long, 'epinfma');
M_.endo_names = char(M_.endo_names, 'zcapf');
M_.endo_names_tex = char(M_.endo_names_tex, 'zcapf');
M_.endo_names_long = char(M_.endo_names_long, 'zcapf');
M_.endo_names = char(M_.endo_names, 'rkf');
M_.endo_names_tex = char(M_.endo_names_tex, 'rkf');
M_.endo_names_long = char(M_.endo_names_long, 'rkf');
M_.endo_names = char(M_.endo_names, 'kf');
M_.endo_names_tex = char(M_.endo_names_tex, 'kf');
M_.endo_names_long = char(M_.endo_names_long, 'kf');
M_.endo_names = char(M_.endo_names, 'pkf');
M_.endo_names_tex = char(M_.endo_names_tex, 'pkf');
M_.endo_names_long = char(M_.endo_names_long, 'pkf');
M_.endo_names = char(M_.endo_names, 'cf');
M_.endo_names_tex = char(M_.endo_names_tex, 'cf');
M_.endo_names_long = char(M_.endo_names_long, 'cf');
M_.endo_names = char(M_.endo_names, 'invef');
M_.endo_names_tex = char(M_.endo_names_tex, 'invef');
M_.endo_names_long = char(M_.endo_names_long, 'invef');
M_.endo_names = char(M_.endo_names, 'yf');
M_.endo_names_tex = char(M_.endo_names_tex, 'yf');
M_.endo_names_long = char(M_.endo_names_long, 'yf');
M_.endo_names = char(M_.endo_names, 'labf');
M_.endo_names_tex = char(M_.endo_names_tex, 'labf');
M_.endo_names_long = char(M_.endo_names_long, 'labf');
M_.endo_names = char(M_.endo_names, 'wf');
M_.endo_names_tex = char(M_.endo_names_tex, 'wf');
M_.endo_names_long = char(M_.endo_names_long, 'wf');
M_.endo_names = char(M_.endo_names, 'rrf');
M_.endo_names_tex = char(M_.endo_names_tex, 'rrf');
M_.endo_names_long = char(M_.endo_names_long, 'rrf');
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, 'zcap');
M_.endo_names_tex = char(M_.endo_names_tex, 'zcap');
M_.endo_names_long = char(M_.endo_names_long, 'zcap');
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, '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, 'pk');
M_.endo_names_tex = char(M_.endo_names_tex, 'pk');
M_.endo_names_long = char(M_.endo_names_long, 'pk');
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, '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, '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, 'lab');
M_.endo_names_tex = char(M_.endo_names_tex, 'lab');
M_.endo_names_long = char(M_.endo_names_long, 'lab');
M_.endo_names = char(M_.endo_names, 'pinf');
M_.endo_names_tex = char(M_.endo_names_tex, 'pinf');
M_.endo_names_long = char(M_.endo_names_long, 'pinf');
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, '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, '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, '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, '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, 'qs');
M_.endo_names_tex = char(M_.endo_names_tex, 'qs');
M_.endo_names_long = char(M_.endo_names_long, 'qs');
M_.endo_names = char(M_.endo_names, 'ms');
M_.endo_names_tex = char(M_.endo_names_tex, 'ms');
M_.endo_names_long = char(M_.endo_names_long, 'ms');
M_.endo_names = char(M_.endo_names, 'spinf');
M_.endo_names_tex = char(M_.endo_names_tex, 'spinf');
M_.endo_names_long = char(M_.endo_names_long, 'spinf');
M_.endo_names = char(M_.endo_names, 'sw');
M_.endo_names_tex = char(M_.endo_names_tex, 'sw');
M_.endo_names_long = char(M_.endo_names_long, 'sw');
M_.endo_names = char(M_.endo_names, 'kpf');
M_.endo_names_tex = char(M_.endo_names_tex, 'kpf');
M_.endo_names_long = char(M_.endo_names_long, 'kpf');
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, 'rnot');
M_.endo_names_tex = char(M_.endo_names_tex, 'rnot');
M_.endo_names_long = char(M_.endo_names_long, 'rnot');
M_.param_names = 'curvw';
M_.param_names_tex = 'curvw';
M_.param_names_long = 'curvw';
M_.param_names = char(M_.param_names, 'cgy');
M_.param_names_tex = char(M_.param_names_tex, 'cgy');
M_.param_names_long = char(M_.param_names_long, 'cgy');
M_.param_names = char(M_.param_names, 'curvp');
M_.param_names_tex = char(M_.param_names_tex, 'curvp');
M_.param_names_long = char(M_.param_names_long, 'curvp');
M_.param_names = char(M_.param_names, 'constelab');
M_.param_names_tex = char(M_.param_names_tex, 'constelab');
M_.param_names_long = char(M_.param_names_long, 'constelab');
M_.param_names = char(M_.param_names, 'constepinf');
M_.param_names_tex = char(M_.param_names_tex, 'constepinf');
M_.param_names_long = char(M_.param_names_long, 'constepinf');
M_.param_names = char(M_.param_names, 'constebeta');
M_.param_names_tex = char(M_.param_names_tex, 'constebeta');
M_.param_names_long = char(M_.param_names_long, 'constebeta');
M_.param_names = char(M_.param_names, 'cmaw');
M_.param_names_tex = char(M_.param_names_tex, 'cmaw');
M_.param_names_long = char(M_.param_names_long, 'cmaw');
M_.param_names = char(M_.param_names, 'cmap');
M_.param_names_tex = char(M_.param_names_tex, 'cmap');
M_.param_names_long = char(M_.param_names_long, 'cmap');
M_.param_names = char(M_.param_names, 'calfa');
M_.param_names_tex = char(M_.param_names_tex, 'calfa');
M_.param_names_long = char(M_.param_names_long, 'calfa');
M_.param_names = char(M_.param_names, 'czcap');
M_.param_names_tex = char(M_.param_names_tex, 'czcap');
M_.param_names_long = char(M_.param_names_long, 'czcap');
M_.param_names = char(M_.param_names, 'cbeta');
M_.param_names_tex = char(M_.param_names_tex, 'cbeta');
M_.param_names_long = char(M_.param_names_long, 'cbeta');
M_.param_names = char(M_.param_names, 'csadjcost');
M_.param_names_tex = char(M_.param_names_tex, 'csadjcost');
M_.param_names_long = char(M_.param_names_long, 'csadjcost');
M_.param_names = char(M_.param_names, 'ctou');
M_.param_names_tex = char(M_.param_names_tex, 'ctou');
M_.param_names_long = char(M_.param_names_long, 'ctou');
M_.param_names = char(M_.param_names, 'csigma');
M_.param_names_tex = char(M_.param_names_tex, 'csigma');
M_.param_names_long = char(M_.param_names_long, 'csigma');
M_.param_names = char(M_.param_names, 'chabb');
M_.param_names_tex = char(M_.param_names_tex, 'chabb');
M_.param_names_long = char(M_.param_names_long, 'chabb');
M_.param_names = char(M_.param_names, 'cfc');
M_.param_names_tex = char(M_.param_names_tex, 'cfc');
M_.param_names_long = char(M_.param_names_long, 'cfc');
M_.param_names = char(M_.param_names, 'cindw');
M_.param_names_tex = char(M_.param_names_tex, 'cindw');
M_.param_names_long = char(M_.param_names_long, 'cindw');
M_.param_names = char(M_.param_names, 'cprobw');
M_.param_names_tex = char(M_.param_names_tex, 'cprobw');
M_.param_names_long = char(M_.param_names_long, 'cprobw');
M_.param_names = char(M_.param_names, 'cindp');
M_.param_names_tex = char(M_.param_names_tex, 'cindp');
M_.param_names_long = char(M_.param_names_long, 'cindp');
M_.param_names = char(M_.param_names, 'cprobp');
M_.param_names_tex = char(M_.param_names_tex, 'cprobp');
M_.param_names_long = char(M_.param_names_long, 'cprobp');
M_.param_names = char(M_.param_names, 'csigl');
M_.param_names_tex = char(M_.param_names_tex, 'csigl');
M_.param_names_long = char(M_.param_names_long, 'csigl');
M_.param_names = char(M_.param_names, 'clandaw');
M_.param_names_tex = char(M_.param_names_tex, 'clandaw');
M_.param_names_long = char(M_.param_names_long, 'clandaw');
M_.param_names = char(M_.param_names, 'crpi');
M_.param_names_tex = char(M_.param_names_tex, 'crpi');
M_.param_names_long = char(M_.param_names_long, 'crpi');
M_.param_names = char(M_.param_names, 'crdy');
M_.param_names_tex = char(M_.param_names_tex, 'crdy');
M_.param_names_long = char(M_.param_names_long, 'crdy');
M_.param_names = char(M_.param_names, 'cry');
M_.param_names_tex = char(M_.param_names_tex, 'cry');
M_.param_names_long = char(M_.param_names_long, 'cry');
M_.param_names = char(M_.param_names, 'crr');
M_.param_names_tex = char(M_.param_names_tex, 'crr');
M_.param_names_long = char(M_.param_names_long, 'crr');
M_.param_names = char(M_.param_names, 'crhoa');
M_.param_names_tex = char(M_.param_names_tex, 'crhoa');
M_.param_names_long = char(M_.param_names_long, 'crhoa');
M_.param_names = char(M_.param_names, 'crhoas');
M_.param_names_tex = char(M_.param_names_tex, 'crhoas');
M_.param_names_long = char(M_.param_names_long, 'crhoas');
M_.param_names = char(M_.param_names, 'crhob');
M_.param_names_tex = char(M_.param_names_tex, 'crhob');
M_.param_names_long = char(M_.param_names_long, 'crhob');
M_.param_names = char(M_.param_names, 'crhog');
M_.param_names_tex = char(M_.param_names_tex, 'crhog');
M_.param_names_long = char(M_.param_names_long, 'crhog');
M_.param_names = char(M_.param_names, 'crhols');
M_.param_names_tex = char(M_.param_names_tex, 'crhols');
M_.param_names_long = char(M_.param_names_long, 'crhols');
M_.param_names = char(M_.param_names, 'crhoqs');
M_.param_names_tex = char(M_.param_names_tex, 'crhoqs');
M_.param_names_long = char(M_.param_names_long, 'crhoqs');
M_.param_names = char(M_.param_names, 'crhoms');
M_.param_names_tex = char(M_.param_names_tex, 'crhoms');
M_.param_names_long = char(M_.param_names_long, 'crhoms');
M_.param_names = char(M_.param_names, 'crhopinf');
M_.param_names_tex = char(M_.param_names_tex, 'crhopinf');
M_.param_names_long = char(M_.param_names_long, 'crhopinf');
M_.param_names = char(M_.param_names, 'crhow');
M_.param_names_tex = char(M_.param_names_tex, 'crhow');
M_.param_names_long = char(M_.param_names_long, 'crhow');
M_.param_names = char(M_.param_names, 'ctrend');
M_.param_names_tex = char(M_.param_names_tex, 'ctrend');
M_.param_names_long = char(M_.param_names_long, 'ctrend');
M_.param_names = char(M_.param_names, 'conster');
M_.param_names_tex = char(M_.param_names_tex, 'conster');
M_.param_names_long = char(M_.param_names_long, 'conster');
M_.param_names = char(M_.param_names, 'cg');
M_.param_names_tex = char(M_.param_names_tex, 'cg');
M_.param_names_long = char(M_.param_names_long, 'cg');
M_.param_names = char(M_.param_names, 'cgamma');
M_.param_names_tex = char(M_.param_names_tex, 'cgamma');
M_.param_names_long = char(M_.param_names_long, 'cgamma');
M_.param_names = char(M_.param_names, 'clandap');
M_.param_names_tex = char(M_.param_names_tex, 'clandap');
M_.param_names_long = char(M_.param_names_long, 'clandap');
M_.param_names = char(M_.param_names, 'cbetabar');
M_.param_names_tex = char(M_.param_names_tex, 'cbetabar');
M_.param_names_long = char(M_.param_names_long, 'cbetabar');
M_.param_names = char(M_.param_names, 'cr');
M_.param_names_tex = char(M_.param_names_tex, 'cr');
M_.param_names_long = char(M_.param_names_long, 'cr');
M_.param_names = char(M_.param_names, 'cpie');
M_.param_names_tex = char(M_.param_names_tex, 'cpie');
M_.param_names_long = char(M_.param_names_long, 'cpie');
M_.param_names = char(M_.param_names, 'crk');
M_.param_names_tex = char(M_.param_names_tex, 'crk');
M_.param_names_long = char(M_.param_names_long, 'crk');
M_.param_names = char(M_.param_names, 'cw');
M_.param_names_tex = char(M_.param_names_tex, 'cw');
M_.param_names_long = char(M_.param_names_long, 'cw');
M_.param_names = char(M_.param_names, 'cikbar');
M_.param_names_tex = char(M_.param_names_tex, 'cikbar');
M_.param_names_long = char(M_.param_names_long, 'cikbar');
M_.param_names = char(M_.param_names, 'cik');
M_.param_names_tex = char(M_.param_names_tex, 'cik');
M_.param_names_long = char(M_.param_names_long, 'cik');
M_.param_names = char(M_.param_names, 'clk');
M_.param_names_tex = char(M_.param_names_tex, 'clk');
M_.param_names_long = char(M_.param_names_long, 'clk');
M_.param_names = char(M_.param_names, 'cky');
M_.param_names_tex = char(M_.param_names_tex, 'cky');
M_.param_names_long = char(M_.param_names_long, 'cky');
M_.param_names = char(M_.param_names, 'ciy');
M_.param_names_tex = char(M_.param_names_tex, 'ciy');
M_.param_names_long = char(M_.param_names_long, 'ciy');
M_.param_names = char(M_.param_names, 'ccy');
M_.param_names_tex = char(M_.param_names_tex, 'ccy');
M_.param_names_long = char(M_.param_names_long, 'ccy');
M_.param_names = char(M_.param_names, 'crkky');
M_.param_names_tex = char(M_.param_names_tex, 'crkky');
M_.param_names_long = char(M_.param_names_long, 'crkky');
M_.param_names = char(M_.param_names, 'cwhlc');
M_.param_names_tex = char(M_.param_names_tex, 'cwhlc');
M_.param_names_long = char(M_.param_names_long, 'cwhlc');
M_.param_names = char(M_.param_names, 'cwly');
M_.param_names_tex = char(M_.param_names_tex, 'cwly');
M_.param_names_long = char(M_.param_names_long, 'cwly');
M_.exo_det_nbr = 0;
M_.exo_nbr = 7;
M_.endo_nbr = 41;
M_.param_nbr = 54;
M_.orig_endo_nbr = 41;
M_.aux_vars = [];
M_.Sigma_e = zeros(7, 7);
M_.Correlation_matrix = eye(7, 7);
M_.H = 0;
M_.Correlation_matrix_ME = 1;
options_.linear = 1;
options_.block=0;
options_.bytecode=0;
options_.use_dll=0;
erase_compiled_function('usmodel_static');
erase_compiled_function('usmodel_dynamic');
M_.lead_lag_incidence = [
 0 21 0;
 0 22 0;
 0 23 0;
 0 24 0;
 0 25 0;
 0 26 0;
 0 27 0;
 1 28 0;
 2 29 0;
 0 30 0;
 0 31 62;
 0 32 0;
 0 33 63;
 3 34 64;
 4 35 65;
 5 36 0;
 0 37 66;
 0 38 0;
 0 39 0;
 0 40 0;
 0 41 0;
 0 42 67;
 0 43 0;
 0 44 68;
 6 45 69;
 7 46 70;
 8 47 0;
 0 48 71;
 9 49 72;
 10 50 73;
 11 51 0;
 12 52 0;
 13 53 0;
 14 54 0;
 15 55 0;
 16 56 0;
 17 57 0;
 18 58 0;
 19 59 0;
 20 60 0;
 0 61 0;]';
M_.nstatic = 15;
M_.nfwrd   = 6;
M_.npred   = 14;
M_.nboth   = 6;
M_.nsfwrd   = 12;
M_.nspred   = 20;
M_.ndynamic   = 26;
M_.equations_tags = {
};
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(41, 1);
M_.maximum_exo_lag = 0;
M_.maximum_exo_lead = 0;
oo_.exo_steady_state = zeros(7, 1);
M_.params = NaN(54, 1);
M_.NNZDerivatives = zeros(3, 1);
M_.NNZDerivatives(1) = 162;
M_.NNZDerivatives(2) = 0;
M_.NNZDerivatives(3) = -1;
M_.params( 13 ) = .025;
ctou = M_.params( 13 );
M_.params( 22 ) = 1.5;
clandaw = M_.params( 22 );
M_.params( 38 ) = 0.18;
cg = M_.params( 38 );
M_.params( 3 ) = 10;
curvp = M_.params( 3 );
M_.params( 1 ) = 10;
curvw = M_.params( 1 );
M_.params( 9 ) = .24;
calfa = M_.params( 9 );
M_.params( 39 ) = 1.004;
cgamma = M_.params( 39 );
M_.params( 11 ) = .9995;
cbeta = M_.params( 11 );
M_.params( 14 ) = 1.5;
csigma = M_.params( 14 );
M_.params( 43 ) = 1.005;
cpie = M_.params( 43 );
M_.params( 16 ) = 1.5;
cfc = M_.params( 16 );
M_.params( 2 ) = 0.51;
cgy = M_.params( 2 );
M_.params( 12 ) = 6.0144;
csadjcost = M_.params( 12 );
M_.params( 15 ) = 0.6361;
chabb = M_.params( 15 );
M_.params( 18 ) = 0.8087;
cprobw = M_.params( 18 );
M_.params( 21 ) = 1.9423;
csigl = M_.params( 21 );
M_.params( 20 ) = 0.6;
cprobp = M_.params( 20 );
M_.params( 17 ) = 0.3243;
cindw = M_.params( 17 );
M_.params( 19 ) = 0.47;
cindp = M_.params( 19 );
M_.params( 10 ) = 0.2696;
czcap = M_.params( 10 );
M_.params( 23 ) = 1.488;
crpi = M_.params( 23 );
M_.params( 26 ) = 0.8762;
crr = M_.params( 26 );
M_.params( 25 ) = 0.0593;
cry = M_.params( 25 );
M_.params( 24 ) = 0.2347;
crdy = M_.params( 24 );
M_.params( 27 ) = 0.9977;
crhoa = M_.params( 27 );
M_.params( 29 ) = 0.5799;
crhob = M_.params( 29 );
M_.params( 30 ) = 0.9957;
crhog = M_.params( 30 );
M_.params( 31 ) = 0.9928;
crhols = M_.params( 31 );
M_.params( 32 ) = 0.7165;
crhoqs = M_.params( 32 );
M_.params( 28 ) = 1;
crhoas = M_.params( 28 );
M_.params( 33 ) = 0;
crhoms = M_.params( 33 );
M_.params( 34 ) = 0;
crhopinf = M_.params( 34 );
M_.params( 35 ) = 0;
crhow = M_.params( 35 );
M_.params( 8 ) = 0;
cmap = M_.params( 8 );
M_.params( 7 ) = 0;
cmaw = M_.params( 7 );
M_.params( 27 ) = .9676;
crhoa = M_.params( 27 );
M_.params( 29 ) = .2703;
crhob = M_.params( 29 );
M_.params( 30 ) = .9930;
crhog = M_.params( 30 );
M_.params( 32 ) = .5724;
crhoqs = M_.params( 32 );
M_.params( 33 ) = .3;
crhoms = M_.params( 33 );
M_.params( 34 ) = .8692;
crhopinf = M_.params( 34 );
M_.params( 35 ) = .9546;
crhow = M_.params( 35 );
M_.params( 8 ) = .7652;
cmap = M_.params( 8 );
M_.params( 7 ) = .8936;
cmaw = M_.params( 7 );
M_.params( 12 ) = 6.3325;
csadjcost = M_.params( 12 );
M_.params( 14 ) = 1.2312;
csigma = M_.params( 14 );
M_.params( 15 ) = 0.7205;
chabb = M_.params( 15 );
M_.params( 18 ) = 0.7937;
cprobw = M_.params( 18 );
M_.params( 21 ) = 2.8401;
csigl = M_.params( 21 );
M_.params( 20 ) = 0.7813;
cprobp = M_.params( 20 );
M_.params( 17 ) = 0.4425;
cindw = M_.params( 17 );
M_.params( 19 ) = 0.3291;
cindp = M_.params( 19 );
M_.params( 10 ) = 0.2648;
czcap = M_.params( 10 );
M_.params( 16 ) = 1.4672;
cfc = M_.params( 16 );
M_.params( 23 ) = 1.7985;
crpi = M_.params( 23 );
M_.params( 26 ) = 0.8258;
crr = M_.params( 26 );
M_.params( 25 ) = 0.0893;
cry = M_.params( 25 );
M_.params( 24 ) = 0.2239;
crdy = M_.params( 24 );
M_.params( 5 ) = 0.7;
constepinf = M_.params( 5 );
M_.params( 6 ) = 0.7420;
constebeta = M_.params( 6 );
M_.params( 4 ) = 1.2918;
constelab = M_.params( 4 );
M_.params( 36 ) = 0.3982;
ctrend = M_.params( 36 );
M_.params( 2 ) = 0.05;
cgy = M_.params( 2 );
M_.params( 9 ) = 0.24;
calfa = M_.params( 9 );
M_.params( 40 ) = M_.params(16);
clandap = M_.params( 40 );
M_.params( 41 ) = M_.params(11)*M_.params(39)^(-M_.params(14));
cbetabar = M_.params( 41 );
M_.params( 42 ) = M_.params(43)/(M_.params(11)*M_.params(39)^(-M_.params(14)));
cr = M_.params( 42 );
M_.params( 44 ) = M_.params(11)^(-1)*M_.params(39)^M_.params(14)-(1-M_.params(13));
crk = M_.params( 44 );
M_.params( 45 ) = (M_.params(9)^M_.params(9)*(1-M_.params(9))^(1-M_.params(9))/(M_.params(40)*M_.params(44)^M_.params(9)))^(1/(1-M_.params(9)));
cw = M_.params( 45 );
M_.params( 46 ) = 1-(1-M_.params(13))/M_.params(39);
cikbar = M_.params( 46 );
M_.params( 47 ) = M_.params(39)*(1-(1-M_.params(13))/M_.params(39));
cik = M_.params( 47 );
M_.params( 48 ) = (1-M_.params(9))/M_.params(9)*M_.params(44)/M_.params(45);
clk = M_.params( 48 );
M_.params( 49 ) = M_.params(16)*M_.params(48)^(M_.params(9)-1);
cky = M_.params( 49 );
M_.params( 50 ) = M_.params(47)*M_.params(49);
ciy = M_.params( 50 );
M_.params( 51 ) = 1-M_.params(38)-M_.params(47)*M_.params(49);
ccy = M_.params( 51 );
M_.params( 52 ) = M_.params(44)*M_.params(49);
crkky = M_.params( 52 );
M_.params( 53 ) = M_.params(49)*M_.params(44)*(1-M_.params(9))*1/M_.params(22)/M_.params(9)/M_.params(51);
cwhlc = M_.params( 53 );
M_.params( 54 ) = 1-M_.params(44)*M_.params(49);
cwly = M_.params( 54 );
M_.params( 36 ) = (M_.params(39)-1)*100;
ctrend = M_.params( 36 );
M_.params( 37 ) = 100*(M_.params(42)-1);
conster = M_.params( 37 );
M_.params( 5 ) = 100*(M_.params(43)-1);
constepinf = M_.params( 5 );
M_.params( 4 ) = 0;
constelab = M_.params( 4 );
M_.params( 43 ) = 1+M_.params(5)/100;
cpie = M_.params( 43 );
M_.params( 39 ) = 1+M_.params(36)/100;
cgamma = M_.params( 39 );
M_.params( 11 ) = 1/(1+M_.params(6)/100);
cbeta = M_.params( 11 );
M_.params( 40 ) = M_.params(16);
clandap = M_.params( 40 );
M_.params( 41 ) = M_.params(11)*M_.params(39)^(-M_.params(14));
cbetabar = M_.params( 41 );
M_.params( 42 ) = M_.params(43)/(M_.params(11)*M_.params(39)^(-M_.params(14)));
cr = M_.params( 42 );
M_.params( 44 ) = M_.params(11)^(-1)*M_.params(39)^M_.params(14)-(1-M_.params(13));
crk = M_.params( 44 );
M_.params( 45 ) = (M_.params(9)^M_.params(9)*(1-M_.params(9))^(1-M_.params(9))/(M_.params(40)*M_.params(44)^M_.params(9)))^(1/(1-M_.params(9)));
cw = M_.params( 45 );
M_.params( 46 ) = 1-(1-M_.params(13))/M_.params(39);
cikbar = M_.params( 46 );
M_.params( 47 ) = M_.params(39)*(1-(1-M_.params(13))/M_.params(39));
cik = M_.params( 47 );
M_.params( 48 ) = (1-M_.params(9))/M_.params(9)*M_.params(44)/M_.params(45);
clk = M_.params( 48 );
M_.params( 49 ) = M_.params(16)*M_.params(48)^(M_.params(9)-1);
cky = M_.params( 49 );
M_.params( 50 ) = M_.params(47)*M_.params(49);
ciy = M_.params( 50 );
M_.params( 51 ) = 1-M_.params(38)-M_.params(47)*M_.params(49);
ccy = M_.params( 51 );
M_.params( 52 ) = M_.params(44)*M_.params(49);
crkky = M_.params( 52 );
M_.params( 53 ) = M_.params(49)*M_.params(44)*(1-M_.params(9))*1/M_.params(22)/M_.params(9)/M_.params(51);
cwhlc = M_.params( 53 );
M_.params( 54 ) = 1-M_.params(44)*M_.params(49);
cwly = M_.params( 54 );
M_.params( 37 ) = 100*(M_.params(42)-1);
conster = M_.params( 37 );
%
% SHOCKS instructions
%
make_ex_;
M_.exo_det_length = 0;
M_.Sigma_e(1, 1) = (0.4618)^2;
M_.Sigma_e(2, 2) = (1.8513)^2;
M_.Sigma_e(3, 3) = (0.6090)^2;
M_.Sigma_e(4, 4) = (0.6017)^2;
M_.Sigma_e(5, 5) = (0.2397)^2;
M_.Sigma_e(6, 6) = (0.1455)^2;
M_.Sigma_e(7, 7) = (0.2089)^2;
M_.sigma_e_is_diagonal = 1;
options_.irf = 0;
var_list_=[];
var_list_ = 'dy';
var_list_ = char(var_list_, 'pinfobs');
var_list_ = char(var_list_, 'robs');
info = stoch_simul(var_list_);
save('usmodel_results.mat', 'oo_', 'M_', 'options_');
if exist('estim_params_', 'var') == 1
  save('usmodel_results.mat', 'estim_params_', '-append');
end
if exist('bayestopt_', 'var') == 1
  save('usmodel_results.mat', 'bayestopt_', '-append');
end
if exist('dataset_', 'var') == 1
  save('usmodel_results.mat', 'dataset_', '-append');
end
if exist('estimation_info', 'var') == 1
  save('usmodel_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
