The steady state has NaNs or Inf

[jpfeifer]

The problem is

Yss=Yhss*Phss/Pss+Ylss*Plss/Pss+Gss;
CGss=ksiss*Yss;
Gss=(beta-1+tauNss*(Whss*Nhss+Wlss*Nlss)+tauKss*(rhss*Khss+rlss*Klss)+tauCss*Pss*Css+tauHss*Phss*Yhss+tauLss*Plss*Ylss-Phss*CGss)/Pss;

You are using Gss before it is known.
See

Difference between different types of moments

Your parameter initializations are not recursive. Try using F9 to execute the initializations before the model-block to see that zeta uses psi_nu in its definition, but psi_nu is only defined later on.

as well as

Difference between different types of moments

Put resid(1); before the steady command. Your model is supposed to be linear, but given initial values of 0, resid(1) returns: [quote]Residuals of the static equations: Equation number 1 : 0 Equation number 2 : 0 Equation number 3 : 0 Equation number 4 : 0 Equation number 5 : 0 Equation number 6 : 0 Equation number 7 : 0 Equation number 8 : -0.03475 Equation number 9 : -0.03475 Equation number 10 : 0 Equation number 11 : 0 Equation number 12 : 0 Equation number 13 : 0 Equation …