%-------------------------------------------------------------------------
%Programa para obtener la dinámica del modelo de Solow con crecimiento de
%la tecnología
%-------------------------------------------------------------------------
%Borro todo antes
clear
clc
%Parámetros del modelo
A0=1;
a=0.01;
rho=0.9;
alpha=0.2;
n=0.005;
delta=0.1;
s=0.15;
kbar=((s)/(a+n+a*n+delta))^(1/(1-alpha));
k0=kbar;
k(1)=k0;
theta(1)=1;
%Loop para la dinámica del capital en el modelo empezando con k0=kbar
for j=1:70
    k(j+1)=(((1-delta)*k(j))/((1+n)*(1+a)))+((s*(theta(j)*k(j)^(alpha)))/((1+n)*(1+a)));
    y(j)=theta(j)*k(j)^(alpha);
    Y(j)=(1+a)^(j)*(1+n)^(j)*y(j);
    theta(j+1)=theta(j)^(rho)*exp(0.007*randn(1));
    lY(j)=log(Y(j));
end
n=length(lY);
for i=1:n-1
    loggY(i)=(lY(i+1)-lY(i));
end
[T,C]=hpfilter(loggY);
cY=std(C)
plot(loggY,'g')
hold on
plot(T)
legend('Crecimiento de la economía','Tendencia')
title('La simulación de una economía')