=====================================
ODE: Orbit Equations
ID:  370ce4f580b4903af42480d2167c358f
=====================================
xd0 = x2
xd1 = x3
xd2 = x0+2.0*x3-par1*(x0+par0)/(sqrt((x0+par0)**2+x1**2)*((x0+par0)**2+x1**2))-par0*(x0-par1)/(sqrt((x0-par1)**2+x1**2)*((x0-par1)**2+x1**2))
xd3 = x1-2.0*x2-par1*x1/(sqrt((x0+par0)**2+x1**2)*((x0+par0)**2+x1**2))-par0*x1/(sqrt((x0-par1)**2+x1**2)*((x0-par1)**2+x1**2))
-------------------------------------
IVP: 124
-------------------------------------
x0 = 0.994
x1 = 0.0
x2 = 0.0
x3 = -2.00158510637908252240537862224
par0 = 0.012277471
par1 = 0.987722529