%Programa para ver la regla de oro gráficamente desde el punto de vista 
%del equilibrio del estado estacionario
%Borro todo antes
clear
clc
%Parámetros del modelo
a=0.02;
theta=0.36;
n=0.05;
delta=0.1;
s1=0.25;
s2=0.45;
sgr=theta;
%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:6];
%Genero un vector y que va a reportarme el lado izquierdo de la ecuación
y=((1+a)*(1+n)-(1-delta))*x;
%Genero un vector z que me va a reportar el lado derecho
z1=(s1*(x.^(theta)));
z2=(s2*(x.^(theta)));
zgr=(sgr*(x.^(theta)));
%También para comparar el consumo en cada caso necesitamos tener un gráfico
%del ingreso
g=(x.^(theta));
plot(x,y)
hold on
plot(x,z1,'g')
hold on
plot(x,z2,'r')
hold on
plot(x,zgr,'k')
hold on
plot(x,g,'y')
hold off
legend('(\delta+n)k','s1Ak^\theta','s2Ak^\theta','sgrAk^\theta','Ingreso')
title('Regla de oro (2)')