Steady state files for deterministic models

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 two 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.


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

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

Error in ==> paper4 at 322

Error in ==> dynare at 120
evalin(‘base’,fname) ;
Paper4.mod (2.98 KB)