%Función que propaga la posición de un planeta(TIERRA o JUPITER) a un instante t desde la
%época dados los elementos orbitales en la época
function EOt=pplanetaLPCPmod(t,EO)
EOt=zeros(1,6);
DeltaT=t/86400;
T0=DeltaT/36525;
EO(1:6)=EO(1:6)+EO(7:12)*T0;
M=EO(6)-EO(5);
M=2*pi*(M/2/pi-floor(M/2/pi));
EO(5)=EO(5)-EO(4);
myfun = @(E,M) M-E+EOt(1)*sin(E);
OPTIONS=optimset('TolFun',1e-10,'TolX',1e-10,'MaxIter',5000,'MaxFunEvals',5000);
X = fzero(@(E) myfun(E,M),pi,OPTIONS);
TAepoch=2*atan(sqrt((1+EOt(1))/(1-EOt(1)))*tan(X/2));
EO(6)=2*pi*(TAepoch/2/pi-floor(TAepoch/2/pi));
EOt=EO(1:6);
end
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
