%-------------------------------------------------------------------------
%Programa para obtener la dinámica del modelo de Solow sin crecimiento de
%la tecnología
%-------------------------------------------------------------------------
%Parámetros del modelo
A=1;
theta=0.36;
n=0.05;
delta=0.1;
s=0.25;
k0=0.5;
k(1)=k0;
T=100;
%Calculo el estado estacionario para luego verificar
kbar=((s*A)/(n+delta))^(1/(1-theta))
%Loop para la dinámica del capital en el modelo
for j=1:200
    k(j+1)=(((1-delta)*k(j))/(1+n))+((s*A*(k(j)^(theta)))/(1+n));
end
%Ahora quiero usar otra forma de ver el estado estacionario
%Recordemos que (delta+n)k=sAk^(theta)
%Genero un vector para k que voy a nombrar x
x=[0:0.01:4];
%Genero un vector y que va a reportarme el lado izquierdo de la ecuación
y=(delta+n)*x;
%Genero un vector z que me va a reportar el lado derecho
z=(s*A*(x.^(theta)));
subplot(2,1,1), plot(k,'b')
subplot(2,1,2), plot(x,y)
hold on
plot(x,z)
hold off