AU=149.598E6;
EOTiE=[1.675733036789938E-02,1.000010743023400E+00*AU,2.632995338813666E-03*pi/180,1.768849465159507E+02*pi/180,2.861361233053117E+02*pi/180,3.567168760884612E+02*pi/180];
VarTi=[-0.00004392,0.00000562*AU,-0.01294668*pi/180,0*pi/180,0.32327364*pi/180,35999.37244981*pi/180];
%
wm=EOTiE(5)+EOTiE(4);
wm=2*pi*(wm/2/pi-floor(wm/2/pi));
e=EOTiE(1);
E=2*atan(tan(EOTiE(6)/2)*sqrt((1-e)/(1+e)));
M=E-e*sin(E);
EOTiE(5)=wm;
EOTiE(6)=M+wm;
EOTi=[EOTiE VarTi];
format long
%%Manera 1
EOt=pplanetaLPCPfun(0,EOTi);
EOt(2)=EOt(2)/AU;
EOt(3:6)=EOt(3:6)*180/pi;
TiLPIni=EOt(6)
%Manera 2
EO=EOTi;
M=EO(6)-EO(5);
M=2*pi*(M/2/pi-floor(M/2/pi))
EO(5)=EO(5)-EO(4);
EOt(1)
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));
EO(2)=EO(2)/AU;
EO(3:6)=EO(3:6)*180/pi;
TiDIRIni=EO(6)
%Lo que deberia salir
TiInicial=[3.567168760884612E+02]
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
