%Este es un programa de Matlab que simula 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.1;
gamma=0.8;
%sigma es el desvio estándar del shock
sigma=0.1;
%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 que está determinada
%totalmente por la ecuación intertemporal del capital y la evolución del
%shock
for i=1:100
    k(i+1)=((beta/(1+beta))*(1-theta))*A(i)*k(i)^(theta);
    A(i+1)=1-gamma+gamma*A(i)+sigma*randn(1);
    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
%Aca vemos los comovimientos y desvios estándares de la economía
desvioY=std(Y)
desvioI=std(I)
A=[Y.' I.'];
corrYI=corr(A)
%Acá usamos la función de calculo de la utilidad
u=util(cyoung,cold,beta);
subplot(3,1,1),plot(A)
title('La evolución de la tecnología')
subplot(3,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('La simulación de la economía')
subplot(3,1,3),plot(u)
title('La utilidad de las diferentes generaciones')