%--------------------------------------------------------------------------
%Ahora vamos a buscar otra forma de obtener el máximo de consumo y cuál es
%el 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.'];
B=[cbar.' s.'];
n=length(cbar);
%Ahora lo que hago es ordenar este vector de mayor a menor
C=sort(A,1,'descend');
%Ahora me genero un loop que me compare el primer valor del sort que se que
%es el maximo de consumo con la matriz B y que me reporte el valor de
%ahorro de este máximo
mayor=[B(1,1) B(1,2)];
for i=1:n
    if C(1)==B(i,1)
        mayor=[B(i,1) B(i,2)];
    else
    end
end
sgr=mayor(1,2)