var K I C q N Y W R;
varexo tau_n tau_k G;
//Hay tres variables exógenas que se determinan en el modelo
// tau_n  impuesto distorsionador del ingreso laboral
// tau_k  impuesto distorsionador de la renta de capital

parameters beta alpha delta eta sigma phi zeta;    //parámetros de base
beta = 0.98;  // factor de descuento, tasa anual
alpha = 0.35; // participación del capital en la producción
delta = 0.05; // tasa de depreciación, anual
eta = 1.0;    // elasticidad de la oferta de trabajo Frisch
sigma = 0.2;  // elasticidad de sustitución intertemporal
phi = 0.01;   // curvatura de la función de ajuste del costo
zeta = 0.5;               

model;                          
                                                // función de acumulación del capital
Y=(K(-1)^alpha)*(N^(1-alpha));                                      // función de producción
W=(1-alpha)*Y/N; 					                                // condición de max. del salario 
R=alpha*(Y/K(-1));                                                  // condición de max. del interés
Y=C+I+G+K(-1)*(phi/2)*((I/K(-1)-delta)^2);                          // condición de equilibrio
G=tau_n*N*W+tau_k*K(-1)*(R-delta-(phi/2)*((I/K(-1)-delta)^2))-T;    // blance presupuestal del gobierno: gasto=recaudación de impuesto por ingreso labora+renta del capital-transferencia
zeta*N^(1/eta)=(C^(-1/phi))*W*(1-tau_n);                            // "q" es el valor sombra de una unidad adicional de capital
q=beta*(((C(+1))^(-1/phi))*((1-tau_k(+1))*R(+1)+tau_k(+1)*delta+phi*((1-tau_k(+1))/2)*(((I(+1)/K))^2)-delta^2))+(1-delta)*q(+1));
q=(C^(-1/phi))*(1+phi*((I/K(-1))-1)*(1-tau_k));
end;

initval;
R=(1/0.98-0.183*0.1-(1-0.1))/(1-0.183);
N=1/3;
K=300;
Y=K^0.35*N^(1-0.35);
G=0.2*Y;
I=0.1*K;
C=0.8*Y-I;
tau_n = 0.362;
tau_k = 0.183;
end;

steady(aolve_algo=1);
check;

//shocks;
//var tau_n; stderr tau_n;
//var tau_k; stderr tau_k;
//end;
//steady;
//check;