%
% 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;
% Define global variables.
global M_ options_ oo_ estim_params_ bayestopt_ dataset_ dataset_info estimation_info ys0_ ex0_
options_ = [];
M_.fname = 'IRF_as_no_policy_eO_robust';
M_.dynare_version = '5.2';
oo_.dynare_version = '5.2';
options_.dynare_version = '5.2';
%
% Some global variables initialization
%
global_initialization;
M_.exo_names = cell(3,1);
M_.exo_names_tex = cell(3,1);
M_.exo_names_long = cell(3,1);
M_.exo_names(1) = {'eA'};
M_.exo_names_tex(1) = {'eA'};
M_.exo_names_long(1) = {'eA'};
M_.exo_names(2) = {'eP'};
M_.exo_names_tex(2) = {'eP'};
M_.exo_names_long(2) = {'eP'};
M_.exo_names(3) = {'eO'};
M_.exo_names_tex(3) = {'eO'};
M_.exo_names_long(3) = {'eO'};
M_.endo_names = cell(31,1);
M_.endo_names_tex = cell(31,1);
M_.endo_names_long = cell(31,1);
M_.endo_names(1) = {'qk'};
M_.endo_names_tex(1) = {'qk'};
M_.endo_names_long(1) = {'qk'};
M_.endo_names(2) = {'Phi_k'};
M_.endo_names_tex(2) = {'Phi\_k'};
M_.endo_names_long(2) = {'Phi_k'};
M_.endo_names(3) = {'eps1'};
M_.endo_names_tex(3) = {'eps1'};
M_.endo_names_long(3) = {'eps1'};
M_.endo_names(4) = {'eps2'};
M_.endo_names_tex(4) = {'eps2'};
M_.endo_names_long(4) = {'eps2'};
M_.endo_names(5) = {'theta'};
M_.endo_names_tex(5) = {'theta'};
M_.endo_names_long(5) = {'theta'};
M_.endo_names(6) = {'qb'};
M_.endo_names_tex(6) = {'qb'};
M_.endo_names_long(6) = {'qb'};
M_.endo_names(7) = {'Y'};
M_.endo_names_tex(7) = {'Y'};
M_.endo_names_long(7) = {'Y'};
M_.endo_names(8) = {'K'};
M_.endo_names_tex(8) = {'K'};
M_.endo_names_long(8) = {'K'};
M_.endo_names(9) = {'N'};
M_.endo_names_tex(9) = {'N'};
M_.endo_names_long(9) = {'N'};
M_.endo_names(10) = {'rk'};
M_.endo_names_tex(10) = {'rk'};
M_.endo_names_long(10) = {'rk'};
M_.endo_names(11) = {'w'};
M_.endo_names_tex(11) = {'w'};
M_.endo_names_long(11) = {'w'};
M_.endo_names(12) = {'Rb'};
M_.endo_names_tex(12) = {'Rb'};
M_.endo_names_long(12) = {'Rb'};
M_.endo_names(13) = {'I'};
M_.endo_names_tex(13) = {'I'};
M_.endo_names_long(13) = {'I'};
M_.endo_names(14) = {'C'};
M_.endo_names_tex(14) = {'C'};
M_.endo_names_long(14) = {'C'};
M_.endo_names(15) = {'BY'};
M_.endo_names_tex(15) = {'BY'};
M_.endo_names_long(15) = {'BY'};
M_.endo_names(16) = {'Lambda'};
M_.endo_names_tex(16) = {'Lambda'};
M_.endo_names_long(16) = {'Lambda'};
M_.endo_names(17) = {'Rk'};
M_.endo_names_tex(17) = {'Rk'};
M_.endo_names_long(17) = {'Rk'};
M_.endo_names(18) = {'spread'};
M_.endo_names_tex(18) = {'spread'};
M_.endo_names_long(18) = {'spread'};
M_.endo_names(19) = {'omegat'};
M_.endo_names_tex(19) = {'omegat'};
M_.endo_names_long(19) = {'omegat'};
M_.endo_names(20) = {'At'};
M_.endo_names_tex(20) = {'At'};
M_.endo_names_long(20) = {'At'};
M_.endo_names(21) = {'L'};
M_.endo_names_tex(21) = {'L'};
M_.endo_names_long(21) = {'L'};
M_.endo_names(22) = {'S'};
M_.endo_names_tex(22) = {'S'};
M_.endo_names_long(22) = {'S'};
M_.endo_names(23) = {'ppit'};
M_.endo_names_tex(23) = {'ppit'};
M_.endo_names_long(23) = {'ppit'};
M_.endo_names(24) = {'omgt'};
M_.endo_names_tex(24) = {'omgt'};
M_.endo_names_long(24) = {'omgt'};
M_.endo_names(25) = {'haircut'};
M_.endo_names_tex(25) = {'haircut'};
M_.endo_names_long(25) = {'haircut'};
M_.endo_names(26) = {'theta_original'};
M_.endo_names_tex(26) = {'theta\_original'};
M_.endo_names_long(26) = {'theta_original'};
M_.endo_names(27) = {'lambda'};
M_.endo_names_tex(27) = {'lambda'};
M_.endo_names_long(27) = {'lambda'};
M_.endo_names(28) = {'Phi_k_tilde'};
M_.endo_names_tex(28) = {'Phi\_k\_tilde'};
M_.endo_names_long(28) = {'Phi_k_tilde'};
M_.endo_names(29) = {'l_ppi'};
M_.endo_names_tex(29) = {'l\_ppi'};
M_.endo_names_long(29) = {'l_ppi'};
M_.endo_names(30) = {'l_omega'};
M_.endo_names_tex(30) = {'l\_omega'};
M_.endo_names_long(30) = {'l_omega'};
M_.endo_names(31) = {'V'};
M_.endo_names_tex(31) = {'V'};
M_.endo_names_long(31) = {'V'};
M_.endo_partitions = struct();
M_.param_names = cell(41,1);
M_.param_names_tex = cell(41,1);
M_.param_names_long = cell(41,1);
M_.param_names(1) = {'beta'};
M_.param_names_tex(1) = {'beta'};
M_.param_names_long(1) = {'beta'};
M_.param_names(2) = {'alpha'};
M_.param_names_tex(2) = {'alpha'};
M_.param_names_long(2) = {'alpha'};
M_.param_names(3) = {'A_ss'};
M_.param_names_tex(3) = {'A\_ss'};
M_.param_names_long(3) = {'A_ss'};
M_.param_names(4) = {'chi'};
M_.param_names_tex(4) = {'chi'};
M_.param_names_long(4) = {'chi'};
M_.param_names(5) = {'ksi'};
M_.param_names_tex(5) = {'ksi'};
M_.param_names_long(5) = {'ksi'};
M_.param_names(6) = {'ppi'};
M_.param_names_tex(6) = {'ppi'};
M_.param_names_long(6) = {'ppi'};
M_.param_names(7) = {'gamma'};
M_.param_names_tex(7) = {'gamma'};
M_.param_names_long(7) = {'gamma'};
M_.param_names(8) = {'delta'};
M_.param_names_tex(8) = {'delta'};
M_.param_names_long(8) = {'delta'};
M_.param_names(9) = {'eta'};
M_.param_names_tex(9) = {'eta'};
M_.param_names_long(9) = {'eta'};
M_.param_names(10) = {'eps_min'};
M_.param_names_tex(10) = {'eps\_min'};
M_.param_names_long(10) = {'eps_min'};
M_.param_names(11) = {'S_ss'};
M_.param_names_tex(11) = {'S\_ss'};
M_.param_names_long(11) = {'S_ss'};
M_.param_names(12) = {'pphi'};
M_.param_names_tex(12) = {'pphi'};
M_.param_names_long(12) = {'pphi'};
M_.param_names(13) = {'theta_ss'};
M_.param_names_tex(13) = {'theta\_ss'};
M_.param_names_long(13) = {'theta_ss'};
M_.param_names(14) = {'B_ss'};
M_.param_names_tex(14) = {'B\_ss'};
M_.param_names_long(14) = {'B_ss'};
M_.param_names(15) = {'eps1_ss'};
M_.param_names_tex(15) = {'eps1\_ss'};
M_.param_names_long(15) = {'eps1_ss'};
M_.param_names(16) = {'qk_ss'};
M_.param_names_tex(16) = {'qk\_ss'};
M_.param_names_long(16) = {'qk_ss'};
M_.param_names(17) = {'qb_ss'};
M_.param_names_tex(17) = {'qb\_ss'};
M_.param_names_long(17) = {'qb_ss'};
M_.param_names(18) = {'L_ss'};
M_.param_names_tex(18) = {'L\_ss'};
M_.param_names_long(18) = {'L_ss'};
M_.param_names(19) = {'Rb_ss'};
M_.param_names_tex(19) = {'Rb\_ss'};
M_.param_names_long(19) = {'Rb_ss'};
M_.param_names(20) = {'rk_ss'};
M_.param_names_tex(20) = {'rk\_ss'};
M_.param_names_long(20) = {'rk_ss'};
M_.param_names(21) = {'K_ss'};
M_.param_names_tex(21) = {'K\_ss'};
M_.param_names_long(21) = {'K_ss'};
M_.param_names(22) = {'N_ss'};
M_.param_names_tex(22) = {'N\_ss'};
M_.param_names_long(22) = {'N_ss'};
M_.param_names(23) = {'Y_ss'};
M_.param_names_tex(23) = {'Y\_ss'};
M_.param_names_long(23) = {'Y_ss'};
M_.param_names(24) = {'w_ss'};
M_.param_names_tex(24) = {'w\_ss'};
M_.param_names_long(24) = {'w_ss'};
M_.param_names(25) = {'I_ss'};
M_.param_names_tex(25) = {'I\_ss'};
M_.param_names_long(25) = {'I_ss'};
M_.param_names(26) = {'C_ss'};
M_.param_names_tex(26) = {'C\_ss'};
M_.param_names_long(26) = {'C_ss'};
M_.param_names(27) = {'Rk_ss'};
M_.param_names_tex(27) = {'Rk\_ss'};
M_.param_names_long(27) = {'Rk_ss'};
M_.param_names(28) = {'Lambda_ss'};
M_.param_names_tex(28) = {'Lambda\_ss'};
M_.param_names_long(28) = {'Lambda_ss'};
M_.param_names(29) = {'omega_bar'};
M_.param_names_tex(29) = {'omega\_bar'};
M_.param_names_long(29) = {'omega_bar'};
M_.param_names(30) = {'rho_a'};
M_.param_names_tex(30) = {'rho\_a'};
M_.param_names_long(30) = {'rho_a'};
M_.param_names(31) = {'rho_ppi'};
M_.param_names_tex(31) = {'rho\_ppi'};
M_.param_names_long(31) = {'rho_ppi'};
M_.param_names(32) = {'rho_omega'};
M_.param_names_tex(32) = {'rho\_omega'};
M_.param_names_long(32) = {'rho_omega'};
M_.param_names(33) = {'rho_b'};
M_.param_names_tex(33) = {'rho\_b'};
M_.param_names_long(33) = {'rho_b'};
M_.param_names(34) = {'phi'};
M_.param_names_tex(34) = {'phi'};
M_.param_names_long(34) = {'phi'};
M_.param_names(35) = {'stdA'};
M_.param_names_tex(35) = {'stdA'};
M_.param_names_long(35) = {'stdA'};
M_.param_names(36) = {'stdP'};
M_.param_names_tex(36) = {'stdP'};
M_.param_names_long(36) = {'stdP'};
M_.param_names(37) = {'stdO'};
M_.param_names_tex(37) = {'stdO'};
M_.param_names_long(37) = {'stdO'};
M_.param_names(38) = {'stdB'};
M_.param_names_tex(38) = {'stdB'};
M_.param_names_long(38) = {'stdB'};
M_.param_names(39) = {'Phi_k_tilde_ss'};
M_.param_names_tex(39) = {'Phi\_k\_tilde\_ss'};
M_.param_names_long(39) = {'Phi_k_tilde_ss'};
M_.param_names(40) = {'l_omega_ss'};
M_.param_names_tex(40) = {'l\_omega\_ss'};
M_.param_names_long(40) = {'l_omega_ss'};
M_.param_names(41) = {'l_ppi_ss'};
M_.param_names_tex(41) = {'l\_ppi\_ss'};
M_.param_names_long(41) = {'l_ppi_ss'};
M_.param_partitions = struct();
M_.exo_det_nbr = 0;
M_.exo_nbr = 3;
M_.endo_nbr = 31;
M_.param_nbr = 41;
M_.orig_endo_nbr = 31;
M_.aux_vars = [];
M_ = setup_solvers(M_);
M_.Sigma_e = zeros(3, 3);
M_.Correlation_matrix = eye(3, 3);
M_.H = 0;
M_.Correlation_matrix_ME = 1;
M_.sigma_e_is_diagonal = true;
M_.det_shocks = [];
M_.surprise_shocks = [];
M_.heteroskedastic_shocks.Qvalue_orig = [];
M_.heteroskedastic_shocks.Qscale_orig = [];
options_.linear = false;
options_.block = false;
options_.bytecode = false;
options_.use_dll = false;
M_.orig_eq_nbr = 31;
M_.eq_nbr = 31;
M_.ramsey_eq_nbr = 0;
M_.set_auxiliary_variables = exist(['./+' M_.fname '/set_auxiliary_variables.m'], 'file') == 2;
M_.epilogue_names = {};
M_.epilogue_var_list_ = {};
M_.orig_maximum_endo_lag = 1;
M_.orig_maximum_endo_lead = 1;
M_.orig_maximum_exo_lag = 0;
M_.orig_maximum_exo_lead = 0;
M_.orig_maximum_exo_det_lag = 0;
M_.orig_maximum_exo_det_lead = 0;
M_.orig_maximum_lag = 1;
M_.orig_maximum_lead = 1;
M_.orig_maximum_lag_with_diffs_expanded = 1;
M_.lead_lag_incidence = [
 0 5 36;
 0 6 37;
 0 7 38;
 0 8 39;
 0 9 40;
 0 10 0;
 0 11 0;
 1 12 0;
 0 13 0;
 0 14 41;
 0 15 0;
 0 16 0;
 0 17 0;
 0 18 0;
 0 19 0;
 0 20 42;
 0 21 0;
 0 22 0;
 0 23 0;
 2 24 0;
 0 25 0;
 0 26 0;
 3 27 43;
 4 28 44;
 0 29 0;
 0 30 0;
 0 31 0;
 0 32 0;
 0 33 0;
 0 34 0;
 0 35 45;]';
M_.nstatic = 19;
M_.nfwrd   = 8;
M_.npred   = 2;
M_.nboth   = 2;
M_.nsfwrd   = 10;
M_.nspred   = 4;
M_.ndynamic   = 12;
M_.dynamic_tmp_nbr = [29; 13; 0; 0; ];
M_.model_local_variables_dynamic_tt_idxs = {
};
M_.equations_tags = {
  1 , 'name' , 'V' ;
  2 , 'name' , '2' ;
  3 , 'name' , '3' ;
  4 , 'name' , 'theta_original' ;
  5 , 'name' , 'eps1' ;
  6 , 'name' , 'eps2' ;
  7 , 'name' , '7' ;
  8 , 'name' , '8' ;
  9 , 'name' , 'rk' ;
  10 , 'name' , 'w' ;
  11 , 'name' , 'Rb' ;
  12 , 'name' , '12' ;
  13 , 'name' , '13' ;
  14 , 'name' , '14' ;
  15 , 'name' , '15' ;
  16 , 'name' , 'Lambda' ;
  17 , 'name' , 'BY' ;
  18 , 'name' , 'Rk' ;
  19 , 'name' , 'spread' ;
  20 , 'name' , 'omegat' ;
  21 , 'name' , '21' ;
  22 , 'name' , '22' ;
  23 , 'name' , 'lambda' ;
  24 , 'name' , 'haircut' ;
  25 , 'name' , 'theta' ;
  26 , 'name' , 'At' ;
  27 , 'name' , 'ppit' ;
  28 , 'name' , 'omgt' ;
  29 , 'name' , '29' ;
  30 , 'name' , 'l_ppi' ;
  31 , 'name' , 'l_omega' ;
};
M_.mapping.qk.eqidx = [2 3 6 12 13 22 30 31 ];
M_.mapping.Phi_k.eqidx = [2 3 5 18 29 30 31 ];
M_.mapping.eps1.eqidx = [2 5 7 12 13 29 30 31 ];
M_.mapping.eps2.eqidx = [2 4 6 12 13 21 ];
M_.mapping.theta.eqidx = [3 15 20 24 25 30 31 ];
M_.mapping.qb.eqidx = [7 11 15 ];
M_.mapping.Y.eqidx = [8 9 10 15 17 ];
M_.mapping.K.eqidx = [8 9 12 13 21 ];
M_.mapping.N.eqidx = [1 8 10 14 ];
M_.mapping.rk.eqidx = [2 9 12 13 18 29 ];
M_.mapping.w.eqidx = [10 14 ];
M_.mapping.Rb.eqidx = [11 19 30 31 ];
M_.mapping.I.eqidx = [12 15 ];
M_.mapping.C.eqidx = [1 14 15 16 ];
M_.mapping.BY.eqidx = [17 ];
M_.mapping.Lambda.eqidx = [2 7 16 29 30 31 ];
M_.mapping.Rk.eqidx = [18 19 ];
M_.mapping.spread.eqidx = [19 ];
M_.mapping.omegat.eqidx = [20 ];
M_.mapping.At.eqidx = [8 26 ];
M_.mapping.L.eqidx = [15 22 23 ];
M_.mapping.S.eqidx = [21 22 ];
M_.mapping.ppit.eqidx = [2 4 12 13 20 21 27 30 31 ];
M_.mapping.omgt.eqidx = [2 12 13 20 28 31 ];
M_.mapping.haircut.eqidx = [24 ];
M_.mapping.theta_original.eqidx = [4 25 ];
M_.mapping.lambda.eqidx = [15 23 25 ];
M_.mapping.Phi_k_tilde.eqidx = [29 ];
M_.mapping.l_ppi.eqidx = [30 ];
M_.mapping.l_omega.eqidx = [31 ];
M_.mapping.V.eqidx = [1 ];
M_.mapping.eA.eqidx = [26 ];
M_.mapping.eP.eqidx = [27 ];
M_.mapping.eO.eqidx = [28 ];
M_.static_and_dynamic_models_differ = false;
M_.has_external_function = false;
M_.state_var = [8 20 23 24 ];
M_.exo_names_orig_ord = [1:3];
M_.maximum_lag = 1;
M_.maximum_lead = 1;
M_.maximum_endo_lag = 1;
M_.maximum_endo_lead = 1;
oo_.steady_state = zeros(31, 1);
M_.maximum_exo_lag = 0;
M_.maximum_exo_lead = 0;
oo_.exo_steady_state = zeros(3, 1);
M_.params = NaN(41, 1);
M_.endo_trends = struct('deflator', cell(31, 1), 'log_deflator', cell(31, 1), 'growth_factor', cell(31, 1), 'log_growth_factor', cell(31, 1));
M_.NNZDerivatives = [131; -1; -1; ];
M_.static_tmp_nbr = [17; 5; 0; 0; ];
M_.model_local_variables_static_tt_idxs = {
};
load data_pam_haoning_robust
M_.params(1) = pam(1);
beta = M_.params(1);
M_.params(8) = pam(2);
delta = M_.params(8);
M_.params(2) = pam(3);
alpha = M_.params(2);
M_.params(3) = log(pam(4));
A_ss = M_.params(3);
M_.params(29) = pam(5);
omega_bar = M_.params(29);
M_.params(4) = pam(6);
chi = M_.params(4);
M_.params(5) = pam(7);
ksi = M_.params(5);
M_.params(9) = pam(8);
eta = M_.params(9);
sigma = pam(9);
M_.params(10) = pam(10);
eps_min = M_.params(10);
M_.params(6) = pam(11);
ppi = M_.params(6);
M_.params(7) = pam(12);
gamma = M_.params(7);
M_.params(14) = log(pam(13));
B_ss = M_.params(14);
M_.params(30) = 0.9;
rho_a = M_.params(30);
M_.params(31) = pam(15);
rho_ppi = M_.params(31);
M_.params(32) = pam(16);
rho_omega = M_.params(32);
M_.params(33) = pam(17);
rho_b = M_.params(33);
M_.params(34) = pam(18);
phi = M_.params(34);
M_.params(35) = pam(19);
stdA = M_.params(35);
M_.params(36) = M_.params(35);
stdP = M_.params(36);
M_.params(37) = M_.params(35);
stdO = M_.params(37);
M_.params(38) = pam(22);
stdB = M_.params(38);
M_.params(12) = pam(23);
pphi = M_.params(12);
load data_steady_state_as_haoning_robust
M_.params(16) = log(steady_state_as(1));
qk_ss = M_.params(16);
Phi_k_ss = log(steady_state_as(2));
M_.params(15) = steady_state_as(3);
eps1_ss = M_.params(15);
eps2_ss = steady_state_as(4);
M_.params(13) = steady_state_as(5);
theta_ss = M_.params(13);
M_.params(17) = log(steady_state_as(6));
qb_ss = M_.params(17);
M_.params(19) = steady_state_as(7);
Rb_ss = M_.params(19);
M_.params(20) = steady_state_as(8);
rk_ss = M_.params(20);
M_.params(21) = log(steady_state_as(9));
K_ss = M_.params(21);
M_.params(22) = log(steady_state_as(10));
N_ss = M_.params(22);
M_.params(23) = log(steady_state_as(11));
Y_ss = M_.params(23);
M_.params(24) = steady_state_as(12);
w_ss = M_.params(24);
M_.params(25) = log(steady_state_as(13));
I_ss = M_.params(25);
M_.params(26) = log(steady_state_as(14));
C_ss = M_.params(26);
omegat_ss = steady_state_as(15);
BY_ss = steady_state_as(16);
M_.params(27) = steady_state_as(17);
Rk_ss = M_.params(27);
M_.params(39) = log(M_.params(1)*(1+M_.params(9)/(1-M_.params(9))*(M_.params(15)/M_.params(10))^((-1)/M_.params(9)))*(M_.params(20)+(1-M_.params(8))*exp(Phi_k_ss)));
Phi_k_tilde_ss = M_.params(39);
M_.params(41) = (1-M_.params(8))*(M_.params(15)/M_.params(10))^((-1)/M_.params(9))*M_.params(9)/(1-M_.params(9))*M_.params(1)*100*M_.params(19)*(1-M_.params(13)^(1/M_.params(9)))-(1+M_.params(9)/(1-M_.params(9))*(M_.params(15)/M_.params(10))^((-1)/M_.params(9)))*M_.params(1)*100*M_.params(19)*M_.params(6)*exp(M_.params(16)-Phi_k_ss);
l_ppi_ss = M_.params(41);
M_.params(40) = exp(M_.params(16)-Phi_k_ss)*(1+M_.params(9)/(1-M_.params(9))*(M_.params(15)/M_.params(10))^((-1)/M_.params(9)))*M_.params(1)*100*M_.params(19)*M_.params(6)*(1-M_.params(29))+(1-M_.params(29))*(1-M_.params(8))*M_.params(9)/(1-M_.params(9))*(M_.params(15)/M_.params(10))^((-1)/M_.params(9))*M_.params(1)*100*M_.params(19)-(1-M_.params(8))*(M_.params(15)/M_.params(10))^((-1)/M_.params(9))*M_.params(9)/(1-M_.params(9))*M_.params(1)*100*M_.params(19)*(1-M_.params(13)^(1/M_.params(9)))*(1-M_.params(29));
l_omega_ss = M_.params(40);
M_.params(11) = (M_.params(6)+(1-M_.params(8))*(eps2_ss/M_.params(10))^((-1)/M_.params(9)))*exp(M_.params(21));
S_ss = M_.params(11);
M_.params(18) = M_.params(16)+M_.params(11);
L_ss = M_.params(18);
M_.params(28) = 1/exp(M_.params(26));
Lambda_ss = M_.params(28);
spread_ss = Rk_ss-Rb_ss;
V_ss = (C_ss-chi*((exp(N_ss))^(1+ksi))/(1+ksi))/(1-beta);
haircut_ss = 1-theta_ss;
%
% INITVAL instructions
%
options_.initval_file = false;
oo_.steady_state(1) = M_.params(16);
oo_.steady_state(2) = Phi_k_ss;
oo_.steady_state(3) = M_.params(15);
oo_.steady_state(4) = eps2_ss;
oo_.steady_state(5) = M_.params(13);
oo_.steady_state(6) = M_.params(17);
oo_.steady_state(7) = M_.params(23);
oo_.steady_state(8) = M_.params(21);
oo_.steady_state(9) = M_.params(22);
oo_.steady_state(10) = M_.params(20);
oo_.steady_state(11) = M_.params(24);
oo_.steady_state(12) = M_.params(19);
oo_.steady_state(13) = M_.params(25);
oo_.steady_state(14) = M_.params(26);
oo_.steady_state(15) = BY_ss;
oo_.steady_state(16) = M_.params(28);
oo_.steady_state(17) = M_.params(27);
oo_.steady_state(18) = spread_ss;
oo_.steady_state(19) = omegat_ss;
oo_.steady_state(20) = M_.params(3);
oo_.steady_state(31) = V_ss;
oo_.steady_state(21) = M_.params(18);
oo_.steady_state(22) = M_.params(11);
oo_.steady_state(23) = log(M_.params(6));
oo_.steady_state(24) = log(M_.params(29));
oo_.steady_state(25) = haircut_ss;
oo_.steady_state(26) = M_.params(13);
oo_.steady_state(28) = M_.params(39);
oo_.steady_state(30) = M_.params(40);
oo_.steady_state(29) = M_.params(41);
oo_.steady_state(27) = 0;
oo_.exo_steady_state(1) = 0;
oo_.exo_steady_state(2) = 0;
oo_.exo_steady_state(3) = 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
%
% SHOCKS instructions
%
M_.exo_det_length = 0;
M_.Sigma_e(3, 3) = (M_.params(37))^2;
options_.solve_algo = 2;
options_.steady.maxit = 100;
steady;
oo_.dr.eigval = check(M_,options_,oo_);
options_.hp_filter = 1600;
options_.irf = 20;
options_.nograph = true;
options_.order = 1;
options_.periods = 0;
var_list_ = {};
[info, oo_, options_, M_] = stoch_simul(M_, options_, oo_, var_list_);
impluse_as_eO =[I_eO K_eO Y_eO C_eO qk_eO eps1_eO eps2_eO Rb_eO omegat_eO spread_eO S_eO haircut_eO N_eO Phi_k_tilde_eO l_ppi_eO l_omega_eO];
save data_impluse_as_eO_robust_haoning impluse_as_eO;
delete *variables.m
delete *static.m
delete *asv


oo_.time = toc(tic0);
disp(['Total computing time : ' dynsec2hms(oo_.time) ]);
if ~exist([M_.dname filesep 'Output'],'dir')
    mkdir(M_.dname,'Output');
end
save([M_.dname filesep 'Output' filesep 'IRF_as_no_policy_eO_robust_results.mat'], 'oo_', 'M_', 'options_');
if exist('estim_params_', 'var') == 1
  save([M_.dname filesep 'Output' filesep 'IRF_as_no_policy_eO_robust_results.mat'], 'estim_params_', '-append');
end
if exist('bayestopt_', 'var') == 1
  save([M_.dname filesep 'Output' filesep 'IRF_as_no_policy_eO_robust_results.mat'], 'bayestopt_', '-append');
end
if exist('dataset_', 'var') == 1
  save([M_.dname filesep 'Output' filesep 'IRF_as_no_policy_eO_robust_results.mat'], 'dataset_', '-append');
end
if exist('estimation_info', 'var') == 1
  save([M_.dname filesep 'Output' filesep 'IRF_as_no_policy_eO_robust_results.mat'], 'estimation_info', '-append');
end
if exist('dataset_info', 'var') == 1
  save([M_.dname filesep 'Output' filesep 'IRF_as_no_policy_eO_robust_results.mat'], 'dataset_info', '-append');
end
if exist('oo_recursive_', 'var') == 1
  save([M_.dname filesep 'Output' filesep 'IRF_as_no_policy_eO_robust_results.mat'], 'oo_recursive_', '-append');
end
if ~isempty(lastwarn)
  disp('Note: warning(s) encountered in MATLAB/Octave code')
end
