var ct, cn, c, yn, p, pn, q; % En orden respectivo, consumo transable, cons no transable, cons total, prod no transable, ipc total, precio no transables y TCR (endógenas).

varexo b, r, pt, yt; % En orden respectivo, bonos, tasa int real, precio transables, dotación de transables (exógenas).

parameters beta, sigma, theta, kappa, gamma;
beta = 0.9787; % factor subjetivo de descuento en el consumo (Soto y Bergoeing)
sigma = 1.5; % elasticidad intertemporal entre consumo presente y futuro (Guajardo)
theta = 1.4; % elasticidad intratemporal de sustitución entre transables y no transables (arbitrario, pero la idea es q sea mayor q 1)
kappa = 1.7829; % factor de desutilidad del trabajo en la función de los hogares (Soto y Bergoeing)
gamma = 0.58; % ponderación del consumo transable dentro del consumo total (tesis L. Cancino)

model;
ct(+1)/ct = ((beta*(1+r(+1)))^sigma)*((pt(+1)/pt)^-theta)*((p/p(+1))^(sigma-theta)); % Euler (solución intertemporal)
cn/ct = (gamma/(1-gamma))*((pn/pt)^-theta); % Solución intratemporal
yn^((1+theta)/theta) = ((theta-1)/(theta*kappa))*(c^(-1/sigma))*(cn^(1/theta))*((pn/p)^-theta); % Trade off consumo total vs producción no transable
b(+1) = (1+r)*b+pn*yn+pt*yt-p*c; % Restricción presupuestaria
c = (((gamma^(1/theta))*(ct^((theta-1)/theta)))+((1-gamma)^(1/theta))*cn^((theta-1)/theta))^(theta/(theta-1)); % Índice de consumo agregado
p = (gamma*(pt^(1-theta))+(1-gamma)*(pn^(1-theta)))^(1/(1-theta)); % Índice total de precios
q = pt/pn; % Tipo de cambio real
end;

initval;
ct = 0.3254;
cn = 0.2356;
c = 0.5610;
yn = 0.2356;
p = 1;
pn = 1;
q = 1;
b = 0;
r = 0.0218;
pt = 1;
yt = 0.3254;
end;
steady;

check;

endval;
ct = 0.3254;
cn = 0.2356;
c = 0.5610;
yn = 0.2356;
p = 1;
pn = 1;
q = 1;
b = 2;
r = 0.04;
pt = 1;
yt = 0.3254;
end;
steady;

simul (periods=50);