Risky steady state for SGU 2003

[jpfeifer]

Please do not cross-post. See my reply at

3rd order perturbation closing devices Uribe

Setting shock_mat(1+burnin,strmatch('eps_sigma_a',M_.exo_names,'exact'))= 1; does not work, because there is no shock named eps_sigma_a in your model.