function EOt=pplanetaLPCPfun(t,EO)
EOt=zeros(1:6);

T0=(t/86400)/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);
    EO(1)
    myfun = @(E,M) M-E+EO(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);
TA=EO(6)*180/pi
