var y c k l w r z b Te;

varexo e;

parameters beta teta mu phi r_1 rho delta sigma;

beta = 0.49;
teta = 0.23;
mu = 0.6;
phi = 0.75;
r_1 = 0.95;
rho = 0.6;
delta = 0.6; 
sigma = 0.7;

model;

c + b(+1) + k(+1) -(1-delta)*k = w*l+(1+r_1)*b + r*k;               % Eq 0  --- Restricción Presupuestaria 

c/(1-l) = (mu/(1-mu))/w;                                            % eq 1 

(c/(1-l))^(-mu) = beta*(c(+1)/(1-l(+1)))^(-mu)*(1-delta+r(+1));     % eq 2 

(c/(1-l))^(-mu)=beta*(c(+1)/(1-l(+1)))^(-mu)*(1+r_1);              % eq 3


% La firma 

y = z*(l^teta)*(k^(1-teta));                                        % eq 4   función de Utilidad de la firma

(1-Te)*teta*z*(l^(teta-1))*(k^(1-teta)) = w;                        % eq 5

(1-Te)*(1-teta)*z*(l^teta)*(k^(-teta))  = r;                        % eq 6


%  Shokc de productividad

log(z(+1)) = rho*log(z)+ e;                                         % eq 7


% Uso del consumo 

c+k(+1)-(1-delta)*k+ Te*y = y;                                      % eq 8 



end;

initval;

y = 1.7;
c = 0.71;
k = 0.5;
l = 0.3;
w = 2.07;
r = 0.03;
z = 1;
Te = 0.5;


end;

steady(solve_algo=4,maxit=1000);

check;

shocks;

var e = sigma^2;

end;

stoch_simul(periods=1000);