我觉得方程写的很完美，代码也没有什么错误啊，为什么找不到稳态啊、是不是初始值的问题啊，求大神啊
var y,c,ch,cf,k,a,h;
varexo e;
parameters beta,sigma,alpha,delta,gama,s,rho,psi,eta,theta;
alpha=0.36;
sigma=0.95;
beta=0.99;
theta=0.4;
delta=0.025;
gama=0.2;
rho=0.95;
psi=0.3;
eta=0.3;
s=2;
model;
c=ch+s*cf;
a=rho*a(-1)+e;
k=(y-c)+(1-delta)*k(-1);
y=exp(a)*(k(-1)^alpha)*(h^(1-alpha));
psi*s=(1-psi)*(ch/cf)^(1/gama);
(1-delta+(h(+1)^(1-alpha))*exp(a(+1))*alpha*(k^(alpha-1)))*beta=(ch(+1)/ch)^(1/gama)*(c^(-sigma))/(c(+1)^(-sigma))*(c/c(+1))^(1/gama);
(ch^(1/gama))*theta*(h^(eta+alpha))=(1-alpha)*exp(a)*(k(-1)^alpha)*(c^((1/gama)-sigma))*psi;
end;
initval;
y=0.5;
c=1.1;
ch=0.1;
cf=0.5;
h=0.29;
k=1;
a=0;
e=0;

end;

shocks;
var e; stderr 0.009;
end;


stoch_simul;