%%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
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
