Initializing variables for steady state

Steady state
1

Hi Everyone! I am currently trying to run a model based on the.

The authors calibrated some values and used Bayesian estimation for the others, I am trying to initialize them just to arrive to a random steady state. I am getting NaN errors in all the residuals when running the steady function.
When I run resid 1, the Nan appear in a selection of variables but I could not find mistaked with my values. How do you initialize values in this case?

rbc.mod (5.6 KB)

2 
a_t=a_t(-1)+ e_at;
a_nt=a_nt(-1)+e_ant;
P_if=P_if(-1)+e_pif;
epsilon_c=epsilon_c(-1)+e_epsilonc;

implies that there are four unit roots and consequently infinitely many steady states. You cannot hope to simply find one numerically. You need to provide one analytical steady state.