%--------------------------------------------------------------------------
%Programa que encuentra el máximo del consumo en estado estacionario y me
%lo linkea al nivel de ahorro que reporta ese máximo
%--------------------------------------------------------------------------
%Borro todo antes
clear
clc
%Parámetros del modelo
a=0.02;
theta=0.36;
n=0.05;
delta=0.1;
%El cambio importante con respecto a los otros programas viene en el ahorro
%Acá el ahorro debe ser un vector que va de 0 a 1 y que nos va a mostrar
%donde se de un consumo máximo
s=[0:0.01:1];
kbar=((s)./(a+n+a*n+delta)).^(1/(1-theta));
ybar=kbar.^(theta);
cbar=(1-s).*ybar;
%Me creo una matriz que tenga dos columnas en una la del consumo y en otra
%la del ahorro. El signo .' hace referencia a la transpuesta de una matriz
%o un vector
A=[cbar.' s.'];
%Con n agarro el tamaño del vector cbar que me va a servir para generar un
%loop
n=length(cbar);
%Acá comienza el loop que me va a comparar los valores de cbar.' y me va a
%elegir el máximo
mayor=[A(1,1) A(1,2)];
for i=1:n-1
    if mayor(1,1)>A(i+1,1)
    else mayor=[A(i+1,1) A(i+1,2)];
    end
end
%Reporto el valor de ahorro que me indica el máximo de cbar
sgr=mayor(1,2)