function tf=ameliore(t0)
d1=1236.853086;
d0=0.827361;
p2=2*pi;
L=p2*frac(0.374897+1325.552410*t0);
LS=p2*frac(0.993133+99.997361*t0);
D=p2*(frac(0.5+d0+d1*t0)-0.5);
F=p2*frac(0.259086+1342.227825*t0);
dlm=22640*sin(L)-4586*sin(L-2*D)+2370*sin(2*D)+769*sin(2*L)-...
    668*sin(LS)-412*sin(2*F)-212*sin(2*L-2*D)-...
	206*sin(L+LS-2*D)+192*sin(L+2*D)-165*sin(LS-2*D)-...
	125*sin(D)-110*sin(L+LS)+148*sin(L-LS)-55*sin(2*F-2*D);
dls=6893*sin(LS)+72*sin(2*LS);
DLambda=D/p2+(dlm-dls)/1296000;
tf=t0-DLambda/d1;
