Ввод данных для Infinity:

l1 := 0.1*diff(x1(t),t)-x1(t)-x2(t) = -x1(t)^3/3
l2 := -2*a*x1(t) + diff(x2(t),t) = -2*x1(t)^3
a := 1.2
SupT := 10
Complicated := 0
SupLocError := 1e-6
analysis({l1,l2, x1(0):=2, x2(0):=eval(2^3/3-a*2)},{x1(t),x2(t)},{0},{0}, {exp(log(2)-2*t), ((exp(log(2)-2*t))^3)/3-1.2*exp(log(2)-2*t)});

Решения, полученные Infinity с флагом Complicated = 0:



Решения, полученные Infinity с флагом Complicated = 1:

Digits := 100
l1:= 0.1*diff(x1(t),t)-x1(t)-x2(t) = -1/3*x1(t)^3;
l2:= -2*a*x1(t) + diff(x2(t),t) = -2*x1(t)^3;
a:=1.2;
SupT := 20
SupLocError := 1e-40
Complicated := 1
analysis({l1,l2, x1(0):=2, x2(0):=eval(2^3/3-a*2)},{x1(t),x2(t)},{0},{0}, {exp(log(2)-2*t), ((exp(log(2)-2*t))^3)/3-1.2*exp(log(2)-2*t)});

Ввод данных для Maple:

x10 := 2: a := 1+2*mu: mu := 0.1:
ode1 := mu*diff(x1(t),t)-x1(t)-x2(t) = -x1(t)^3/3:
ode2 := -2*a*x1(t)+diff(x2(t),t) = -2*x1(t)^3:
sol := dsolve({ode1,ode2,x1(0) = x10, x2(0) = x10^3/3 - a*x10},{x1(t),x2(t)}, type =
numeric,method=taylorseries);
with(plots):
range0 := 0..10;
g1 := odeplot(sol,[t,x1(t)],range0,color= green, numpoints = 2000):
g2 := odeplot(sol,[t,x2(t)],range0,color = red, numpoints = 2000):
display(g1,g2);

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

Точные решения
plot([exp(log(2)-2*t), ((exp(log(2)-2*t))^3)/3-1.2*exp(log(2)-2*t)], t = 0..10);

Copyright © 1998-2004, MathRevolt. All rights reserved