var Yn Yt Kn Kt In It C Cn Ct f P Pn Pt Zn Zt;

varexo EZn EZt;

parameters a r B s c h dt v dn ROn ROt re; 

a=0.63;
r=0.38;
B = 0.9612;
s = 0.057;
c = 0.05; 
h =0.05;
ROt = 0.34; 
ROn=0.33;
dt=0.013;
dn=0.0075;
re=0.04;
v=1.7;

model;

%1 Produccion transables

Yt=(Zt)*(Kt(-1))^c;
%2 Produccion de no transables
Yn=(Zn)*(Kn(-1))^h;
%3 Inversion transables 
It=Kt-(1-dt)*Kt(-1);
%4Inversion  no transables
In=Kn-(1-dn)*Kn(-1);
%5 Consumo agregado
C=(a*((Ct)^(1-(1/r)))+(1-a)*((Cn)^(1-(1/r))))^(1/(1-(1/r)));
%6 Cambio en la riqueza neta de la economia domestica en en resto del mundo(re=rentabilidad extranjera) que es f y como lo definimos  
f=((Pt)/P)*(Yt-(Ct)-(It))+((Pn)/P)*(Yn-(Cn)-(In))+(1+re)*f(-1);
%7 Relacion de consumo de los bienes  tranasables sobre el consumo de los bienes no transablles  (numerario)
((Ct)/Cn)=(((1-a)/a)*((Pt)/Pn))^r;
%8 Relacion de consumo de los bienes tranasables sobre el total de consumo 
((Ct)/C) =(a^r)*((Pt/P)^(-r));

%9Relacion de consumo de los bienes no tranasables sobre el total de consumo 
((Cn)/C) =((1-a)^r)*(Pn/P)^(-r);

%10Nivel de precios de la economia 
P=((a^r)*((Pt)^(1-r))+((1-a)^r)*((Pn)^(1-r)))^(1/(1-r));

%11 Consumo futuro
C(+1)=((B*(1+re(+1)))^(1/s))*(C);

%12
c*(Zt(+1))*((Kt)^(c-1))+(1- dt)=(1+re(+1))*((P(+1)/P)/(Pt(+1)/Pt));

%13
h*(Zn(+1))*((Kn)^(h-1))+(1- dn)=(1+re(+1))*((P(+1)/P)/(Pn(+1)/Pn));

%14 Shock Como decirle al programa como "e" distrubuye normal.
  
Ln(Zt)=ROt*ln(Zt(-1))+EZt;  

Ln(Zn)=ROn*ln(Zn(-1))+EZn; 
 
end;

initval;

Kt=2.93;
Kn=1.26;
Yt=1.93;
Yn=1.36;
It=6.8;
In=4.66;
Ct=0.68;
Cn=0.77;
C=0.8;
Pt=1.67;
Pn=1;
P=1.5;
Zt=0.4;
Zn=0.4;
f=0.5;

end; 
resid;
steady;