u1'=v1 
v1'=u1+2*v2-muhat*(u1+mu)/d1-mu*(u1-muhat)/d2
u2'=v2 
v2'=u2-2*v1-muhat*u2/d1-mu*u2/d2
d1=((u1+mu)^2+u2^2)^(3/2)
d2=((u1-muhat)^2+u2^2)^(3/2)
par mu=.012277471
!muhat =1-mu
init u1=.994,u2=0,v1=0,v2=-2.0015851063790825
@ total=20,meth=5,tol=1e-6,atol=1e-6
@ xp=u1,yp=u2,xlo=-1.5,ylo=-1.5,xhi=1.0001,yhi=1.5
done