var y c i n l k z;
varexo eps_z;
parameters alpha beta phi psi delta rho sigma kappa lambda ;
  alpha=0.33;
  beta=0.975;
  phi=2;
  psi=0.5;
  delta=0.06;
  rho=0.9;
  sigma=0.01;
  kappa=psi*(1-phi)-1;
  lambda=(1-psi)*(1-phi);
model(linear);
#yss=0.9;
#css=0.7;
#iss=0.2;
#nss=0.47;
#lss=0.53;
#kss=3.5;
#zss=1;
c-l=z+alpha*k(-1)-alpha*n;

kappa*c+lambda*l=kappa*c(+1)+lambda*l(+1)+((alpha*zss*((nss/kss)^(1-alpha)))/(alpha*zss*((nss/kss)^(1-alpha))+(1-delta)))*z(+1)+
((alpha*(1-alpha)*zss*((n/k)^(1-alpha)))/(alpha*zss*((nss/kss)^(1-alpha))+(1-delta)))*n(+1)-
((alpha*(1-alpha)*zss*((n/k)^(1-alpha)))/(alpha*zss*((nss/kss)^(1-alpha))+(1-delta)))*k;

y=z+alpha*k(-1)+(1-alpha)*n;

y=(css/(css+iss))*c+(iss/(css+iss))*i;

k=(iss/(iss+(1-delta)*kss))*i+((1-delta)*kss/(iss+(1-delta)*kss))*k(-1);

(nss/(nss+lss))*n+(lss/(nss+lss))*l=0;

z=rho*z(-1)+eps_z;

end;

steady;
check;

shocks;
var eps_z;
stderr 0.01;
end;

stoch_simul(order=1);



上述是我的dynare代码，我是按照结构方程计量经济学第二版第五章里面的第一个例子自己写的一个dynare程序，我首先把原来的非线性方程组进行对数线性化写成对数偏差的形式，但是最后会报错，错误形势如下：
错误使用 print_info (line 80)
The steady state contains NaN or Inf
出错 steady (line 92)
    print_info(info,options_.noprint, options_);
出错 newtest (line 136)
steady;
出错 dynare (line 180)
evalin('base',fname) ;


其中的yss表示的是水平值的稳态，是手动求出来的