close all
yi=[GEO(2):10:10*GEO(1)];
C=[];
for i=1:length(yi)

C(i)=-1/2*(sqrt(mu/yi(i))-wt*yi(i))^2+1/2*wt^2*yi(i)^2+pot(0,yi(i),0,GEO,mu);
end

plot(yi,C,'LineWidth',1.8)
title('Evolución de C para órbitas inicialmente circulares')
ylabel('C')
xlabel('Distancia inicial (Km)')
grid on
