Ввод данных для 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)});

Решения, полученные Infinity:



Ввод данных для 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);

Решения, полученные Maple:

Copyright © 1998-2004, MathRevolt. All rights reserved