%%Radio de Hill para un mismo C con diferentes inclinaciones
close all
Csuf=Cs
in=[0:1:80]*pi/180;
Yin=[]
for i=1:length(in)
    YhINC = fzero(@(yi) HillRadINC(yi,Csuf,in(i)),100e8);
    Yin(i)=YhINC;
end
semilogy(in*180/pi,Yin,'LineWidth',1.8)

str=sprintf('Evolución del radio de Hill para órbitas inclinadas, C=%g',C);
title(str);
ylabel('R_{Hill}')
xlabel('Inclinación (º)')
grid on