var a b  y wA wB c  k kG I IG l lA lB m n i w  r;
//
varexo epsm epsn epsi;
//-----------------------------------------------
// 2. Declaration of Parameters
//-----------------------------------------------
parameters
tauK
tauL
G
chi
beta
alphaA
betaA
alphaB
betaB
gamma
theta
delta
deltag
rhom
rhon

rhoi
sigmam
sigman
sigmai;
//-----------------------------------------------
// 3. Parameter Calibration
//-----------------------------------------------
tauK=0.8;
tauL=0.8;
G = 0.1;
beta = 0.9992;
alphaA = 0.21;
betaA = 0.79;
alphaB =0.31;
betaB =0.42;
gamma =0.27;
theta=0.5;
delta = 0.5;
deltag=0.05;
rhom = 0.9;
rhon =0.9;
//rhog = 0.9;
rhoi=0.9;
chi = 0.5;
sigmam = 0.001;
sigman = 0.001;
sigmai=0.001;
// parameters used for initial value calculation
//-----------------------------------------------
// 4. The MODEL
//-----------------------------------------------
model;
(1/(c)) = beta*(1/(c(+1)))*(1+(1-tauK)*r-delta);
//c = (1-tauL)*w*l +(1-tauK)*r*k-I;
chi*(c)/(1-l) = (1-tauL)*w;
a = exp(m)*(kG^alphaA)*(lA^betaA);
b = exp(n)*(k^alphaB)*(lB^betaB);
y = (theta*b^gamma+(1-theta)*a^gamma)^(1/gamma);
y=c+I+IG+G;
c+I+IG+G =(theta*b^gamma+(1-theta)*a^gamma)^(1/gamma);
l =lA+lB;
I = k(+1)-(1-delta)*k;
IG = k(+1)-(1-deltag)*k;
wA*lA = betaA*a;
wB*lB = betaB*b;
w*l= wA*lA+wB*lB;
r*k= alphaB*b;
G + exp(i)*IG = w*tauL*l + r*tauK*k;

exp(m)= exp((m(-1))^rhom)*exp(epsm);
exp(n)= exp((n(-1))^rhon)*exp(epsn);
exp(i)= exp((i(-1))^rhoi)*exp(epsi);
//m = rhom*m(-1)+epsm;
//i = rhoi*i(-1)+epsi;
//n = rhon*n(-1)+epsn;
end;
resid(1);
//-----------------------------------------------
// 5. Initial guesses for steady-state computation
//-----------------------------------------------
initval;
k=0.20;
kG=0.08;
I=0.3;
IG=0.1;
l=0.8;
lA=0.2;
lB=0.6;
m=0.05;
n=0.05;
i=0.05;
//end;
//steady_state_model;
a=exp(m)*(kG^alphaA)*(lA^betaA);
b=exp(n)*(k^alphaB)*(lB^betaB);
y=(theta*b^gamma+(1-theta)*a^gamma)^(1/gamma);
//c=0.0176676;
c=(1-tauL)*(betaA*exp(m)*(kG^alphaA)*(lA^betaA)+betaB*exp(n)*(k^alphaB)*(lB^betaB))*(1-l)/l;
w=betaA*exp(m)*(kG^alphaA)*(lA^betaA)+betaB*exp(n)*(k^alphaB)*(lB^betaB);
wA=betaA*exp(m)*(kG^alphaA)*(lA^betaA)/lA;
wB=betaB*exp(n)*(k^alphaB)*(lB^betaB)/lB;
//r= 0.798342;
r=alphaB*exp(n)*(k^alphaB)*(lB^betaB)/k;
//epsm=0;
//epsn=0;
//epsi=0;
end;
//-----------------------------------------------
// 6. Specification of shocks
//-----------------------------------------------
shocks;
var epsm = sigmam^2;
var epsn = sigman^2;
var epsi = sigmai^2;
//var epsm; stderr 0.014;
//var epsn; stderr 0.014;
//var epsi;stderr 0.015;
end;

steady;
check;
solve_algo=3;

stoch_simul ( order=1);

运行结果提示如下错误：

错误使用 print_info (line 74)
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

出错 steady (line 92)
    print_info(info,options_.noprint, options_);

出错 try4 (line 272)
steady;

出错 dynare (line 180)
evalin('base',fname) ;

想请教各位是什么原因出差呢？是方程错误还是其他情况呢？

非常感谢