Hi Everyone!

I would like to study **permanent shock** in my model.

For the steady state I use /matlab/ steady state file with the same name as model file (model_steadystate.m). And for temporary shocks there is no problem. Everything works properly.

But when it comes to permanent shock I use the /steady/ command twice. As I understand the first steady state comes from the /matlab/ file. But how to deliver the second steady state, in the /endval/ block. How to do this externaly like in the case of /initval/ block?

Unfortunalety I have many endogenous variables in my model and their steady state values are obtained numerically. So writing everything in the /endval/ block is very time consuming and difficult. Is it possible to implement in some way t**wo different steady states** for the /initval/ and /endval/ block?

Look at the attached example. Once you have defined a steady_state_model block, you need only to set the value of the exogenous variables in initval and endval. I have not tested it with the old format of a separate <mod_filename>_steadystate.m file, but I think that it should work as well.

Please, let me know.

Best

Michel

rbc_det4.mod (2.24 KB)

**Thank You very much.** From Your answer I finally have understood what is the idea of getting permanent shocks in /dynare/. I have used this in my case by writing external /matlab/ steady state file where I exploit not only /global M_/ but also /global oo_/ calling for the values form the main /dynare/ file with the whole model specification. It is very nice that I didn’t have to write any additional file or change anything in the source files.

Best regards

Karolina

PS. For anyone interested. Everythings goes similarly like in fs2000ns steadystate.m (in the models folder under UserGuide). So there is:

k = M_.params(1); and so on

but I use also

epsFE = oo_.exo_steady_state(1); and so on.

I am getting that error and I don’t know how to solve it. I have attached my code. Someone can help me?

??? Error using ==> print_info at 57

Impossible to find the steady state. Either the model doesn’t have a steady state, there are an infinity of steady

states, or the guess values are too far from the solution

Error in ==> steady at 92

print_info(info,options_.noprint);

Error in ==> paper4 at 322

steady;

Error in ==> dynare at 120

evalin(‘base’,fname) ;

Paper4.mod (2.98 KB)