Infinity Online Help - The problem of two bodies

The problem of two bodies

Data input for Infinity:

l1 := diff(x1(t),t)= x3(t)*x2(t)
l2 := diff(x2(t),t) = -x1(t)*x3(t)
l3 := diff(x3(t),t) = -0.51*x1(t)*x2(t)
SupT := 10
SupLocError := 1e-3
Complicated := 1
analysis({l1,l2,l3, x1(0) := 0, x2(0) := 1, x3(0) := 1}, {x1(t),x2(t),x3(t)});

Solutions received by Infinity:

Data input for Maple:

l1 := diff(x1(t), t) = x2(t)*x3(t);
l2 := diff(x2(t), t) = -x1(t)*x3(t);
l3 := diff(x3(t), t) = -0.51*x1(t)*x2(t);
sol := dsolve({l1,l2, l3, x1(0) = 0, x2(0) = 1, x3(0) = 1},{x1(t),x2(t),x3(t)}, type = numeric,method=taylorseries);
with(plots):
range0 := 0..10;
g1 := odeplot(sol,[t,x1(t)],range0,color = red):
g2 := odeplot(sol,[t,x2(t)],range0,color = green):
g3 := odeplot(sol,[t,x3(t)],range0,color = blue):
display(g1,g2,g3);

Solutions received by Maple:

Note: If you solved your particular task with the help of Infinity and you want to share your know-how, send us the solution. We will publish it at our site and put it into the help file.