Using dynare4.4.3 for slove two-country DSGE problems, I have already got the steady value, and I also doing the linear, but both two solutions(1.got the steady value then input the model；2.linear the model）have the same error! What does these mean? and how to solve? Thanks very much!!
Warning: Matrix is singular to working precision.
In stochastic_solvers at 188
In resol at 137
In stoch_simul at 88
In liukun1212 at 376
In dynare at 180
Warning: Matrix is singular to working precision.
In stochastic_solvers at 190
In resol at 137
In stoch_simul at 88
In liukun1212 at 376
In dynare at 180
Error using eig
Input to EIG must not contain NaN or Inf.
Error in stochastic_solvers (line 192)
dr.eigval = eig(transition_matrix(dr));
Error in resol (line 137)
[dr,info] = stochastic_solvers(dr,check_flag,M,options,oo);
Error in stoch_simul (line 88)
[oo_.dr,info,M_,options_,oo_] = resol(0,M_,options_,oo_);
Error in liukun1212 (line 376)
info = stoch_simul(var_list_);
Error in dynare (line 180)
evalin(‘base’,fname) ;
and the mod below:
var ch yh ih kh wh lh dh qh chh chf ihh ihf shh shf cah nfah vah zh cf yf if kf wf lf df qf cff cfh iff ifh sff sfh caf nfaf vaf zf ;
varexo epszh epszf;
parameters sigma kappa alpha delta phi beta omega rhoh rhof etaa etab etac etad muc muy mui muk muw mul mud muq much mucf muih muif mush musf ;
sigma = 2.00; //elasticity of substitution (consumption)
kappa = 0.34; //captial weight
alpha = 0.85; //home bias
delta = 0.06; //depreciation rate
phi = 2.00; //elasticity of substitution between home and foreign
beta = 0.96; //discount factor
omega = 0.50; //elasticity of labor supply w.r.t. real wage
rhoh = 0.50; //persistence of home
rhof = 0.30; //persistence of foreign
etaa = 1-beta+deltabeta-deltakappabeta;
etab = 1-beta+deltabeta;
etac = sigma+omega-kappasigma-kappaomega;
etad = (1-kappa)^(1/(sigma+omega)); //temporary variable
//steady
muc=(etaa^(omega/(sigma+omega)))((kappabeta)^((kappa+omegakappa)/etac))etad/(etab^((omega+kappa)/etac));
muy=(etaa^(-sigma/(sigma+omega)))((kappabeta)^((kappa+omegakappa)/etac))etad/(etab^((kappa-sigma+kappasigma+kappaomega)/etac));
mui=delta*(etaa^(-sigma/(sigma+omega)))((kappabeta)^((kappa+sigma+omega-kappasigma)/etac))etad/(etab^((kappa+omega)/etac));
muk=(etaa^(-sigma/(sigma+omega)))((kappabeta)^((kappa+sigma+omega-kappasigma)/etac))etad/(etab^((kappa+omega)/etac));
muw=(1-kappa)((kappabeta)/etab)^(kappa/(1-kappa));
mul=(etaa^(-sigma/(sigma+omega)))((kappabeta)^(kappa*(1-sigma)/etac))etad/(etab^((kappa-sigma)/etac));
mud=(1-beta)(etaa^(-sigma/(sigma+omega)))((kappabeta)^((kappa+sigma+omega-kappasigma)/etac))etad/(beta(etab^((kappa+omega)/etac)));
muq=(etaa^(-sigma/(sigma+omega)))((kappabeta)^((kappa+sigma+omega-kappasigma)/etac))etad/(etab^((kappa+omega)/etac));
much=alpha(etaa^(omega/(sigma+omega)))((kappabeta)^((kappa+omegakappa)/etac))etad/(etab^((omega+kappa)/etac));
mucf=(1-alpha)(etaa^(omega/(sigma+omega)))((kappabeta)^((kappa+omegakappa)/etac))etad/(etab^((omega+kappa)/etac));
muih=alphadelta*(etaa^(-sigma/(sigma+omega)))((kappabeta)^((kappa+sigma+omega-kappasigma)/etac))etad/(etab^((kappa+omega)/etac));
muif=(1-alpha)delta(etaa^(-sigma/(sigma+omega)))((kappabeta)^((kappa+sigma+omega-kappasigma)/etac))etad/(etab^((kappa+omega)/etac));
mush=0.5+0.5((1-kappa)(2alpha-1))/(1-kappa(2alpha-1));
musf=0.5-0.5((1-kappa)(2alpha-1))/(1-kappa*(2*alpha-1));
model(linear);
//consumption CES home
ch(-1) = alpha*chh(-1)+(1-alpha)*chf(-1); // 1
//consumption CES foreign
cf(-1) = alpha*cff(-1)+(1-alpha)*cfh(-1); // 2
//production function home
yh(-1) = kappa*kh(-1)+(1-kappa)*lh(-1)+zh(-1); // 3
//production function foreign
yf(-1) = kappa*kf(-1)+(1-kappa)*lf(-1)+zf(-1); // 4
//dynamic capital home
kh = (1-delta)kh(-1)+deltaih(-1); // 5
//dynamic capital foreign
kf = (1-delta)kf(-1)+deltaif(-1); // 6
//investment index home
ih(-1) = alpha*ihh(-1)+(1-alpha)*ihf(-1); // 7
//investment index foreign
if(-1) = alpha*iff(-1)+(1-alpha)*ifh(-1); // 8
//wage and output home
wh(-1)+lh(-1) = yh(-1); // 9
//wage and output foreign
wf(-1)+lf(-1) = yf(-1); // 10
//dividend appartment home
yh(-1) =(muddh(-1)+muiih(-1))/(mud+mui); // 11
//dividend appartment foreign
yf(-1) = (muddf(-1)+muiif(-1))/(mud+mui); // 12
//firm optimal home
kh-sigmach(-1)+sigmach = ((kappamuy)/(kappamuy+1-delta))*yh(-1); // 13
//firm optimal foreign
kf-sigmacf(-1)+sigmacf = ((kappamuy)/(kappamuy+1-delta))*yf(-1); // 14
//budget consumption home
mucch(-1)+shhmushmuq+shfmusf*muq = (wh(-1)+lh(-1))muwmul+shh(-1)mushmuq+(shh(-1)+dh(-1))mushmud+shf(-1)musfmuq+(shf(-1)+df(-1))musfmud; //15
//budget consumption foreign
muccf(-1)+sffmushmuq+sfhmusf*muq = (wf(-1)+lf(-1))muwmul+sff(-1)mushmuq+(sff(-1)+df(-1))mushmud+sfh(-1)musfmuq+(sfh(-1)+dh(-1))musfmud; //16
//euler of consumption home
sigmach-sigmach(-1) = (muqqh+muddh)/(muq+mud); //17
//euler of consumption foreign
sigmacf-sigmacf(-1) = (muqqf+muddf)/(muq+mud); //18
//euler of equity share home/foreign
sigma*(ch(-1)-cf(-1))=sigma*(ch-cf); //19
//euler labor home
wh(-1)-omegalh(-1) = sigmach(-1); //20
//euler labor foreign
wf(-1)-omegalf(-1) = sigmacf(-1); //21
//market clean production home
muchchh(-1)+mucfcfh(-1)+muihihh(-1)+muififh(-1) = muy*yh(-1); //22
//market clean production foreign
muchcff(-1)+mucfchf(-1)+muihiff(-1)+muifihf(-1) = muy*yf(-1); //23
//market clean equity
mushshh(-1)+musfsfh(-1) = mushsff(-1)+musfshf(-1); //24
//NFA expression home
nfah(-1)=(musfmuq/muy)(shf(-1)+qf(-1))-(musfmuq/muy)(sfh(-1)+qh(-1)); //25
//NFA expression foreign
nfaf(-1) = -nfah(-1); //26
//CA expression home
cah(-1) = nfah-nfah(-1)-vah(-1); //27
//CA expression foreign
caf(-1) = cah(-1); //28
//VA expression home
vah(-1) = (musfmuq/muy)(qf(-1)-qf)-(musfmuq/muy)(qh-qh(-1)); //29
//VA expression foreign
vaf(-1) = -vah(-1); //30
//consumption chh and chf home
chh(-1) = ch(-1); //31
//consumption cff and cfh foreign
cff(-1) = cf(-1); //32
//investment ihh and ihf home
ihh(-1) = ih(-1); //33
//investment ihh and ihf foreign
iff(-1) = if(-1); //34
//technology shock home
zh(-1) = rhoh*zh(-2)+epszh; //35
//technology shock foreign
zf(-1) = rhof*zf(-2)+epszf; //36
end;
shocks;
var epszh; stderr 0.0210;
var epszf; stderr 0.0150;
corr epszh, epszf = 0.25;
end;
steady;
stoch_simul(order=1,irf=20);