本人在求解DSGE模型的时候遇到了困难，现已将模型设置好，但是Dynare代码出现错误如下，恳请前辈们慷慨相助。。。或者联系本人QQ350641328或791935570，定有重谢！！！本人的程序如下所示：
var a b y c  k kG I IG l lA lB m n i w wA wB 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
rhog
rhoi
sigmam
sigman
sigmai;
//-----------------------------------------------
// 3. Parameter Calibration
//-----------------------------------------------
tauK=0.002;
tauL=0.002;
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.02;
deltag=0.05;
rhom = 0.9;
rhon =0.9;
rhog = 0.9;
rhoi=0.9;
chi = 1;
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;
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(i) = exp((i(-1))^rhoi)*exp(epsi);
exp(n) = exp((n(-1))^rhon)*exp(epsn);
end;
resid(1);
//-----------------------------------------------
// 5. Initial guesses for steady-state computation
//-----------------------------------------------
initval;
k=20;
kG=8;
I=0.3;
IG=0.1;
l=0.8;
lA=0.2;
lB=0.6;
m=0.001;
n=0.001;
i=0.001;
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=(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=alphaB*exp(n)*(k^alphaB)*(lB^betaB)/k;
end;
//-----------------------------------------------
// 6. Specification of shocks
//-----------------------------------------------
shocks;
var epsm = sigmam^2;
var epsn = sigman^2;
var epsi = sigmai^2;
end;
steady;
solve_algo=1;

stoch_simul ( order=1);

本人的错误如下所示：??? Error using ==> print_info at 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

Error in ==> steady at 92
    print_info(info,options_.noprint, options_);

Error in ==> try3 at 277
steady;

Error in ==> dynare at 162
evalin('base',fname) ;