%Este es un programa de Matlab que grafica el impulso respuesta
%la economía de ciclos reales del handout de la clase 12 de George
clc
clear
%Parámetros 
beta=0.98;
theta=0.36;
delta=0.10;
gamma=0.9;
%sigma es el desvio estándar del shock
sigma=0.2;
%Usamos el estado estacionario del capital
kbar=((beta/(1+beta))*(1-theta))^(1/(1-theta));
%Indicamos que la economía parte del estado estacionario tanto de capital
%como de el shock de tecnología
A(1)=1;
k(1)=kbar;
%Acá comienza el loop que simula la economía separando entre el período 
%que se da el shockque está determinada totalmente por la ecuación 
%intertemporal del capital y la evolución del shock
for i=1:100
    if i==2
    k(i+1)=((beta/(1+beta))*(1-theta))*A(i)*k(i)^(theta);
    A(i+1)=1-gamma+gamma*A(i)+sigma;
    I(i)=k(i+1)-(1-delta)*k(i);
    Y(i)=A(i)*k(i)^(theta);
    C(i)=Y(i)+(1-delta)*k(i)-k(i+1);
    w(i)=(1-theta)*A(i)*k(i)^(theta);
    r(i+1)=theta*A(i+1)*k(i+1)^(theta-1);
    cyoung(i)=w(i)-k(i+1);
    cold(i)=r(i+1)*k(i+1)+(1-delta)*k(i+1);
    else
    k(i+1)=((beta/(1+beta))*(1-theta))*A(i)*k(i)^(theta);
    A(i+1)=1-gamma+gamma*A(i);
    I(i)=k(i+1)-(1-delta)*k(i);
    Y(i)=A(i)*k(i)^(theta);
    C(i)=Y(i)+(1-delta)*k(i)-k(i+1);
    w(i)=(1-theta)*A(i)*k(i)^(theta);
    r(i+1)=theta*A(i+1)*k(i+1)^(theta-1);
    cyoung(i)=w(i)-k(i+1);
    cold(i)=r(i+1)*k(i+1)+(1-delta)*k(i+1);
    end
end
subplot(2,1,1),plot(A)
title('La evolución de la tecnología')
subplot(2,1,2),plot(k,'g')
hold on
plot(I)
hold on
plot(Y,'r')
hold on
plot(C,'y')
hold off
legend('Capital Fisico','Inversion','Ingreso','Consumo')
title('Impulso Respuesta')