%
% 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 = 'LGK';
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('LGK.log');
M_.exo_names = 'ea';
M_.exo_names_tex = 'ea';
M_.exo_names_long = 'ea';
M_.exo_names = char(M_.exo_names, 'epsi');
M_.exo_names_tex = char(M_.exo_names_tex, 'epsi');
M_.exo_names_long = char(M_.exo_names_long, 'epsi');
M_.endo_names = 'y';
M_.endo_names_tex = 'y';
M_.endo_names_long = 'y';
M_.endo_names = char(M_.endo_names, 'k');
M_.endo_names_tex = char(M_.endo_names_tex, 'k');
M_.endo_names_long = char(M_.endo_names_long, 'k');
M_.endo_names = char(M_.endo_names, '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, 'z');
M_.endo_names_tex = char(M_.endo_names_tex, 'z');
M_.endo_names_long = char(M_.endo_names_long, 'z');
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, 'uc');
M_.endo_names_tex = char(M_.endo_names_tex, 'uc');
M_.endo_names_long = char(M_.endo_names_long, 'uc');
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, 'qi');
M_.endo_names_tex = char(M_.endo_names_tex, 'qi');
M_.endo_names_long = char(M_.endo_names_long, 'qi');
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, 'phin');
M_.endo_names_tex = char(M_.endo_names_tex, 'phin');
M_.endo_names_long = char(M_.endo_names_long, 'phin');
M_.endo_names = char(M_.endo_names, 'ni');
M_.endo_names_tex = char(M_.endo_names_tex, 'ni');
M_.endo_names_long = char(M_.endo_names_long, 'ni');
M_.endo_names = char(M_.endo_names, 'nn');
M_.endo_names_tex = char(M_.endo_names_tex, 'nn');
M_.endo_names_long = char(M_.endo_names_long, 'nn');
M_.endo_names = char(M_.endo_names, 'n');
M_.endo_names_tex = char(M_.endo_names_tex, 'n');
M_.endo_names_long = char(M_.endo_names_long, 'n');
M_.endo_names = char(M_.endo_names, 'd');
M_.endo_names_tex = char(M_.endo_names_tex, 'd');
M_.endo_names_long = char(M_.endo_names_long, 'd');
M_.endo_names = char(M_.endo_names, 'qbar');
M_.endo_names_tex = char(M_.endo_names_tex, 'qbar');
M_.endo_names_long = char(M_.endo_names_long, 'qbar');
M_.endo_names = char(M_.endo_names, 'phii');
M_.endo_names_tex = char(M_.endo_names_tex, 'phii');
M_.endo_names_long = char(M_.endo_names_long, 'phii');
M_.endo_names = char(M_.endo_names, 'nu');
M_.endo_names_tex = char(M_.endo_names_tex, 'nu');
M_.endo_names_long = char(M_.endo_names_long, 'nu');
M_.endo_names = char(M_.endo_names, 'omegai');
M_.endo_names_tex = char(M_.endo_names_tex, 'omegai');
M_.endo_names_long = char(M_.endo_names_long, 'omegai');
M_.endo_names = char(M_.endo_names, 'omegan');
M_.endo_names_tex = char(M_.endo_names_tex, 'omegan');
M_.endo_names_long = char(M_.endo_names_long, 'omegan');
M_.endo_names = char(M_.endo_names, 'mui');
M_.endo_names_tex = char(M_.endo_names_tex, 'mui');
M_.endo_names_long = char(M_.endo_names_long, 'mui');
M_.endo_names = char(M_.endo_names, 'qn');
M_.endo_names_tex = char(M_.endo_names_tex, 'qn');
M_.endo_names_long = char(M_.endo_names_long, 'qn');
M_.endo_names = char(M_.endo_names, 'rkii');
M_.endo_names_tex = char(M_.endo_names_tex, 'rkii');
M_.endo_names_long = char(M_.endo_names_long, 'rkii');
M_.endo_names = char(M_.endo_names, 'rkin');
M_.endo_names_tex = char(M_.endo_names_tex, 'rkin');
M_.endo_names_long = char(M_.endo_names_long, 'rkin');
M_.endo_names = char(M_.endo_names, 'rkni');
M_.endo_names_tex = char(M_.endo_names_tex, 'rkni');
M_.endo_names_long = char(M_.endo_names_long, 'rkni');
M_.endo_names = char(M_.endo_names, 'rknn');
M_.endo_names_tex = char(M_.endo_names_tex, 'rknn');
M_.endo_names_long = char(M_.endo_names_long, 'rknn');
M_.endo_names = char(M_.endo_names, 'wuc');
M_.endo_names_tex = char(M_.endo_names_tex, 'wuc');
M_.endo_names_long = char(M_.endo_names_long, 'wuc');
M_.endo_names = char(M_.endo_names, 'wc');
M_.endo_names_tex = char(M_.endo_names_tex, 'wc');
M_.endo_names_long = char(M_.endo_names_long, 'wc');
M_.endo_names = char(M_.endo_names, 'wi');
M_.endo_names_tex = char(M_.endo_names_tex, 'wi');
M_.endo_names_long = char(M_.endo_names_long, 'wi');
M_.endo_names = char(M_.endo_names, 'womegai');
M_.endo_names_tex = char(M_.endo_names_tex, 'womegai');
M_.endo_names_long = char(M_.endo_names_long, 'womegai');
M_.endo_names = char(M_.endo_names, 'womegan');
M_.endo_names_tex = char(M_.endo_names_tex, 'womegan');
M_.endo_names_long = char(M_.endo_names_long, 'womegan');
M_.endo_names = char(M_.endo_names, 'wrkii');
M_.endo_names_tex = char(M_.endo_names_tex, 'wrkii');
M_.endo_names_long = char(M_.endo_names_long, 'wrkii');
M_.endo_names = char(M_.endo_names, 'wrkin');
M_.endo_names_tex = char(M_.endo_names_tex, 'wrkin');
M_.endo_names_long = char(M_.endo_names_long, 'wrkin');
M_.endo_names = char(M_.endo_names, 'wrkni');
M_.endo_names_tex = char(M_.endo_names_tex, 'wrkni');
M_.endo_names_long = char(M_.endo_names_long, 'wrkni');
M_.endo_names = char(M_.endo_names, 'wrknn');
M_.endo_names_tex = char(M_.endo_names_tex, 'wrknn');
M_.endo_names_long = char(M_.endo_names_long, 'wrknn');
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, 'psi');
M_.endo_names_tex = char(M_.endo_names_tex, 'psi');
M_.endo_names_long = char(M_.endo_names_long, 'psi');
M_.endo_names = char(M_.endo_names, 'phic');
M_.endo_names_tex = char(M_.endo_names_tex, 'phic');
M_.endo_names_long = char(M_.endo_names_long, 'phic');
M_.endo_names = char(M_.endo_names, 'psig');
M_.endo_names_tex = char(M_.endo_names_tex, 'psig');
M_.endo_names_long = char(M_.endo_names_long, 'psig');
M_.endo_partitions = struct();
M_.param_names = 'alpha';
M_.param_names_tex = 'alpha';
M_.param_names_long = 'alpha';
M_.param_names = char(M_.param_names, 'gamma');
M_.param_names_tex = char(M_.param_names_tex, 'gamma');
M_.param_names_long = char(M_.param_names_long, 'gamma');
M_.param_names = char(M_.param_names, 'beta');
M_.param_names_tex = char(M_.param_names_tex, 'beta');
M_.param_names_long = char(M_.param_names_long, 'beta');
M_.param_names = char(M_.param_names, 'pii');
M_.param_names_tex = char(M_.param_names_tex, 'pii');
M_.param_names_long = char(M_.param_names_long, 'pii');
M_.param_names = char(M_.param_names, 'pin');
M_.param_names_tex = char(M_.param_names_tex, 'pin');
M_.param_names_long = char(M_.param_names_long, 'pin');
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, 'epsilon');
M_.param_names_tex = char(M_.param_names_tex, 'epsilon');
M_.param_names_long = char(M_.param_names_long, 'epsilon');
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, 'xi');
M_.param_names_tex = char(M_.param_names_tex, 'xi');
M_.param_names_long = char(M_.param_names_long, 'xi');
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, 'rhopsi');
M_.param_names_tex = char(M_.param_names_tex, 'rhopsi');
M_.param_names_long = char(M_.param_names_long, 'rhopsi');
M_.param_names = char(M_.param_names, 'sigmaa');
M_.param_names_tex = char(M_.param_names_tex, 'sigmaa');
M_.param_names_long = char(M_.param_names_long, 'sigmaa');
M_.param_names = char(M_.param_names, 'sigmapsi');
M_.param_names_tex = char(M_.param_names_tex, 'sigmapsi');
M_.param_names_long = char(M_.param_names_long, 'sigmapsi');
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, 'sigma');
M_.param_names_tex = char(M_.param_names_tex, 'sigma');
M_.param_names_long = char(M_.param_names_long, 'sigma');
M_.param_names = char(M_.param_names, 'kappa');
M_.param_names_tex = char(M_.param_names_tex, 'kappa');
M_.param_names_long = char(M_.param_names_long, 'kappa');
M_.param_names = char(M_.param_names, 'C');
M_.param_names_tex = char(M_.param_names_tex, 'C');
M_.param_names_long = char(M_.param_names_long, 'C');
M_.param_names = char(M_.param_names, 'Y');
M_.param_names_tex = char(M_.param_names_tex, 'Y');
M_.param_names_long = char(M_.param_names_long, 'Y');
M_.param_names = char(M_.param_names, 'Z');
M_.param_names_tex = char(M_.param_names_tex, 'Z');
M_.param_names_long = char(M_.param_names_long, 'Z');
M_.param_names = char(M_.param_names, 'I');
M_.param_names_tex = char(M_.param_names_tex, 'I');
M_.param_names_long = char(M_.param_names_long, 'I');
M_.param_names = char(M_.param_names, 'Ni');
M_.param_names_tex = char(M_.param_names_tex, 'Ni');
M_.param_names_long = char(M_.param_names_long, 'Ni');
M_.param_names = char(M_.param_names, 'Nn');
M_.param_names_tex = char(M_.param_names_tex, 'Nn');
M_.param_names_long = char(M_.param_names_long, 'Nn');
M_.param_names = char(M_.param_names, 'K');
M_.param_names_tex = char(M_.param_names_tex, 'K');
M_.param_names_long = char(M_.param_names_long, 'K');
M_.param_names = char(M_.param_names, 'Qn');
M_.param_names_tex = char(M_.param_names_tex, 'Qn');
M_.param_names_long = char(M_.param_names_long, 'Qn');
M_.param_names = char(M_.param_names, 'RKin');
M_.param_names_tex = char(M_.param_names_tex, 'RKin');
M_.param_names_long = char(M_.param_names_long, 'RKin');
M_.param_names = char(M_.param_names, 'RKii');
M_.param_names_tex = char(M_.param_names_tex, 'RKii');
M_.param_names_long = char(M_.param_names_long, 'RKii');
M_.param_names = char(M_.param_names, 'N');
M_.param_names_tex = char(M_.param_names_tex, 'N');
M_.param_names_long = char(M_.param_names_long, 'N');
M_.param_names = char(M_.param_names, 'D');
M_.param_names_tex = char(M_.param_names_tex, 'D');
M_.param_names_long = char(M_.param_names_long, 'D');
M_.param_names = char(M_.param_names, 'Qbar');
M_.param_names_tex = char(M_.param_names_tex, 'Qbar');
M_.param_names_long = char(M_.param_names_long, 'Qbar');
M_.param_names = char(M_.param_names, 'MUi');
M_.param_names_tex = char(M_.param_names_tex, 'MUi');
M_.param_names_long = char(M_.param_names_long, 'MUi');
M_.param_names = char(M_.param_names, 'PHIi');
M_.param_names_tex = char(M_.param_names_tex, 'PHIi');
M_.param_names_long = char(M_.param_names_long, 'PHIi');
M_.param_names = char(M_.param_names, 'PHIn');
M_.param_names_tex = char(M_.param_names_tex, 'PHIn');
M_.param_names_long = char(M_.param_names_long, 'PHIn');
M_.param_names = char(M_.param_names, 'OMEGAi');
M_.param_names_tex = char(M_.param_names_tex, 'OMEGAi');
M_.param_names_long = char(M_.param_names_long, 'OMEGAi');
M_.param_names = char(M_.param_names, 'OMEGAn');
M_.param_names_tex = char(M_.param_names_tex, 'OMEGAn');
M_.param_names_long = char(M_.param_names_long, 'OMEGAn');
M_.param_names = char(M_.param_names, 'NU');
M_.param_names_tex = char(M_.param_names_tex, 'NU');
M_.param_names_long = char(M_.param_names_long, 'NU');
M_.param_names = char(M_.param_names, 'RKni');
M_.param_names_tex = char(M_.param_names_tex, 'RKni');
M_.param_names_long = char(M_.param_names_long, 'RKni');
M_.param_names = char(M_.param_names, 'RKnn');
M_.param_names_tex = char(M_.param_names_tex, 'RKnn');
M_.param_names_long = char(M_.param_names_long, 'RKnn');
M_.param_names = char(M_.param_names, 'etac');
M_.param_names_tex = char(M_.param_names_tex, 'etac');
M_.param_names_long = char(M_.param_names_long, 'etac');
M_.param_names = char(M_.param_names, 'etauc');
M_.param_names_tex = char(M_.param_names_tex, 'etauc');
M_.param_names_long = char(M_.param_names_long, 'etauc');
M_.param_names = char(M_.param_names, 'etai');
M_.param_names_tex = char(M_.param_names_tex, 'etai');
M_.param_names_long = char(M_.param_names_long, 'etai');
M_.param_names = char(M_.param_names, 'etaomegai');
M_.param_names_tex = char(M_.param_names_tex, 'etaomegai');
M_.param_names_long = char(M_.param_names_long, 'etaomegai');
M_.param_names = char(M_.param_names, 'etaomegan');
M_.param_names_tex = char(M_.param_names_tex, 'etaomegan');
M_.param_names_long = char(M_.param_names_long, 'etaomegan');
M_.param_names = char(M_.param_names, 'etarkii');
M_.param_names_tex = char(M_.param_names_tex, 'etarkii');
M_.param_names_long = char(M_.param_names_long, 'etarkii');
M_.param_names = char(M_.param_names, 'etarkin');
M_.param_names_tex = char(M_.param_names_tex, 'etarkin');
M_.param_names_long = char(M_.param_names_long, 'etarkin');
M_.param_names = char(M_.param_names, 'etarkni');
M_.param_names_tex = char(M_.param_names_tex, 'etarkni');
M_.param_names_long = char(M_.param_names_long, 'etarkni');
M_.param_names = char(M_.param_names, 'etarknn');
M_.param_names_tex = char(M_.param_names_tex, 'etarknn');
M_.param_names_long = char(M_.param_names_long, 'etarknn');
M_.param_names = char(M_.param_names, 'nug');
M_.param_names_tex = char(M_.param_names_tex, 'nug');
M_.param_names_long = char(M_.param_names_long, 'nug');
M_.param_partitions = struct();
M_.exo_det_nbr = 0;
M_.exo_nbr = 2;
M_.endo_nbr = 38;
M_.param_nbr = 47;
M_.orig_endo_nbr = 38;
M_.aux_vars = [];
M_.Sigma_e = zeros(2, 2);
M_.Correlation_matrix = eye(2, 2);
M_.H = 0;
M_.Correlation_matrix_ME = 1;
M_.sigma_e_is_diagonal = 1;
M_.det_shocks = [];
options_.linear = 1;
options_.block=0;
options_.bytecode=0;
options_.use_dll=0;
M_.hessian_eq_zero = 1;
erase_compiled_function('LGK_static');
erase_compiled_function('LGK_dynamic');
M_.orig_eq_nbr = 38;
M_.eq_nbr = 38;
M_.ramsey_eq_nbr = 0;
M_.set_auxiliary_variables = exist(['./' M_.fname '_set_auxiliary_variables.m'], 'file') == 2;
M_.lead_lag_incidence = [
 0 7 0;
 1 8 41;
 2 9 42;
 0 10 43;
 0 11 0;
 0 12 44;
 0 13 0;
 0 14 45;
 3 15 46;
 0 16 0;
 0 17 0;
 0 18 0;
 0 19 0;
 4 20 0;
 0 21 0;
 0 22 0;
 0 23 0;
 0 24 47;
 0 25 48;
 0 26 0;
 0 27 49;
 0 0 50;
 0 0 51;
 0 0 52;
 0 0 53;
 0 28 0;
 0 29 0;
 0 30 0;
 0 31 0;
 0 32 0;
 0 33 0;
 0 34 0;
 0 35 0;
 0 36 0;
 5 37 0;
 6 38 54;
 0 39 0;
 0 40 0;]';
M_.nstatic = 22;
M_.nfwrd   = 10;
M_.npred   = 2;
M_.nboth   = 4;
M_.nsfwrd   = 14;
M_.nspred   = 6;
M_.ndynamic   = 16;
M_.equations_tags = {
};
M_.static_and_dynamic_models_differ = 0;
M_.exo_names_orig_ord = [1:2];
M_.maximum_lag = 1;
M_.maximum_lead = 1;
M_.maximum_endo_lag = 1;
M_.maximum_endo_lead = 1;
oo_.steady_state = zeros(38, 1);
M_.maximum_exo_lag = 0;
M_.maximum_exo_lead = 0;
oo_.exo_steady_state = zeros(2, 1);
M_.params = NaN(47, 1);
M_.NNZDerivatives = [133; 0; -1];
M_.params( 1 ) = 0.33;
alpha = M_.params( 1 );
M_.params( 2 ) = 0.5;
gamma = M_.params( 2 );
M_.params( 3 ) = 0.99;
beta = M_.params( 3 );
M_.params( 4 ) = 0.25;
pii = M_.params( 4 );
M_.params( 5 ) = 1-M_.params(4);
pin = M_.params( 5 );
M_.params( 6 ) = 5.013;
chi = M_.params( 6 );
M_.params( 7 ) = 0.1;
epsilon = M_.params( 7 );
M_.params( 8 ) = 0.1288;
theta = M_.params( 8 );
M_.params( 9 ) = 0.0019;
xi = M_.params( 9 );
M_.params( 10 ) = 0.95;
rhoa = M_.params( 10 );
M_.params( 12 ) = 0.02;
sigmaa = M_.params( 12 );
M_.params( 11 ) = 0.66;
rhopsi = M_.params( 11 );
M_.params( 13 ) = (-0.05);
sigmapsi = M_.params( 13 );
M_.params( 14 ) = 0.025;
delta = M_.params( 14 );
M_.params( 15 ) = 0.972;
sigma = M_.params( 15 );
M_.params( 16 ) = 1;
kappa = M_.params( 16 );
M_.params( 17 ) = 3.2706;
C = M_.params( 17 );
M_.params( 18 ) = 4.1817;
Y = M_.params( 18 );
M_.params( 19 ) = 0.0379;
Z = M_.params( 19 );
M_.params( 20 ) = 0.9110;
I = M_.params( 20 );
M_.params( 21 ) = 2.2781;
Ni = M_.params( 21 );
M_.params( 22 ) = 6.8323;
Nn = M_.params( 22 );
M_.params( 23 ) = 36.4417;
K = M_.params( 23 );
M_.params( 24 ) = 1.0094;
Qn = M_.params( 24 );
M_.params( 25 ) = 1.0220;
RKin = M_.params( 25 );
M_.params( 26 ) = 1.0129;
RKii = M_.params( 26 );
M_.params( 27 ) = 9.1104;
N = M_.params( 27 );
M_.params( 28 ) = 27.3313;
D = M_.params( 28 );
M_.params( 29 ) = 1.0070;
Qbar = M_.params( 29 );
M_.params( 30 ) = 0.0050;
MUi = M_.params( 30 );
M_.params( 31 ) = 4.2991;
PHIi = M_.params( 31 );
M_.params( 32 ) = 3.9369;
PHIn = M_.params( 32 );
M_.params( 33 ) = 0.5661;
OMEGAi = M_.params( 33 );
M_.params( 34 ) = 0.5207;
OMEGAn = M_.params( 34 );
M_.params( 35 ) = 0.5321;
NU = M_.params( 35 );
M_.params( 36 ) = 1.0034;
RKni = M_.params( 36 );
M_.params( 37 ) = 1.0125;
RKnn = M_.params( 37 );
M_.params( 38 ) = 0;
etac = M_.params( 38 );
M_.params( 39 ) = 0;
etauc = M_.params( 39 );
M_.params( 40 ) = 0;
etai = M_.params( 40 );
M_.params( 41 ) = 0;
etaomegai = M_.params( 41 );
M_.params( 42 ) = 0;
etaomegan = M_.params( 42 );
M_.params( 43 ) = 0;
etarkii = M_.params( 43 );
M_.params( 44 ) = 0;
etarkin = M_.params( 44 );
M_.params( 45 ) = 0;
etarkni = M_.params( 45 );
M_.params( 46 ) = 0;
etarknn = M_.params( 46 );
M_.params( 47 ) = 0;
nug = M_.params( 47 );
%
% INITVAL instructions
%
options_.initval_file = 0;
oo_.steady_state( 1 ) = 0;
oo_.steady_state( 2 ) = 0;
oo_.steady_state( 3 ) = 0;
oo_.steady_state( 4 ) = 0;
oo_.steady_state( 5 ) = 0;
oo_.steady_state( 6 ) = 0;
oo_.steady_state( 7 ) = 0;
oo_.steady_state( 8 ) = 0;
oo_.steady_state( 9 ) = 0;
oo_.steady_state( 10 ) = 0;
oo_.steady_state( 11 ) = 0;
oo_.steady_state( 12 ) = 0;
oo_.steady_state( 13 ) = 0;
oo_.steady_state( 14 ) = 0;
oo_.steady_state( 15 ) = 0;
oo_.steady_state( 16 ) = 0;
oo_.steady_state( 17 ) = 0;
oo_.steady_state( 18 ) = 0;
oo_.steady_state( 19 ) = 0;
oo_.steady_state( 20 ) = 0;
oo_.steady_state( 21 ) = 0;
oo_.steady_state( 22 ) = 0;
oo_.steady_state( 23 ) = 0;
oo_.steady_state( 24 ) = 0;
oo_.steady_state( 25 ) = 0;
oo_.steady_state( 26 ) = 0;
oo_.steady_state( 27 ) = 0;
oo_.steady_state( 28 ) = 0;
oo_.steady_state( 29 ) = 0;
oo_.steady_state( 30 ) = 0;
oo_.steady_state( 31 ) = 0;
oo_.steady_state( 32 ) = 0;
oo_.steady_state( 33 ) = 0;
oo_.steady_state( 34 ) = 0;
oo_.steady_state( 35 ) = 0;
oo_.steady_state( 36 ) = 0;
oo_.steady_state( 37 ) = 0;
oo_.steady_state( 38 ) = 0;
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
steady;
oo_.dr.eigval = check(M_,options_,oo_);
%
% SHOCKS instructions
%
M_.exo_det_length = 0;
M_.Sigma_e(1, 1) = M_.params(12)^2;
M_.Sigma_e(2, 2) = M_.params(13)^2;
var_list_ = char();
info = stoch_simul(var_list_);
save('LGK_results.mat', 'oo_', 'M_', 'options_');
if exist('estim_params_', 'var') == 1
  save('LGK_results.mat', 'estim_params_', '-append');
end
if exist('bayestopt_', 'var') == 1
  save('LGK_results.mat', 'bayestopt_', '-append');
end
if exist('dataset_', 'var') == 1
  save('LGK_results.mat', 'dataset_', '-append');
end
if exist('estimation_info', 'var') == 1
  save('LGK_results.mat', 'estimation_info', '-append');
end
if exist('dataset_info', 'var') == 1
  save('LGK_results.mat', 'dataset_info', '-append');
end
if exist('oo_recursive_', 'var') == 1
  save('LGK_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
