var C L Y K R I W A G ;
varexo e ;
parameters  alpha beta tau delta gamma rho ;
alpha=0.35 ;
beta=0.97 ;
tau=0.225 ;
delta=0.06 ;
gamma=0.40 ;
rho=0.95 ;
model ;
I=Y-C-G;
K=I+(1-delta)*K(-1);
Y=A*(K(-1)^alpha)*(L^(1-alpha));
R=alpha*A*(K(-1)^(alpha-1))*(L^(1-alpha));
W=(1-alpha)*A*(K(-1)^alpha)*(L^(-alpha));
A=(1-rho)+rho*A(-1)+e;
C(+1)/C=beta*((1-tau)*R(+1)-delta+1);
C=(gamma/(1-gamma))*(1-tau)*W*(1-L);
G=tau*W*L+tau*R*K(-1);
end;
initval ;
C=0.918 ;
L=0.333 ;
K=12.669 ;
Y=1.235 ;
I=0.351;
R=0.002;
W=2;
A=1;
e=0;
end ;
steady ;
resid;

check;

shocks;

var e;
stderr 1;

end;

stoch_simul;