
Ввод данных для Infinity:
restart
Digits := 300
f(t) := heaviside(t)
l1 := diff(x1(t),t) = 1/x2(t)+f(t);
l2 := diff(x2(t),t) = x1(t) + x2(t)^2;
SupT := 10
SupLocError := 1e-20
Complicated := 1
analysis({l1,l2, x1(0) := 1.1, x2(0) := 0.1}, {x1(t), x2(t)}, {f(t)})
Решения, полученные Infinity:
Ввод данных для Maple:
l1 := diff(x1(t),t) = 1/x2(t) + 1:
l2 := diff(x2(t),t) = x1(t) + x2(t)^2:
sol := dsolve({l1,l2,x1(0) = 1.1, x2(0) = 0.1},{x1(t),x2(t)}, type =
numeric,method=taylorseries);
with(plots):
range0 := 0..5;
g1 := odeplot(sol,[t,x1(t)],range0,color = red, numpoints = 2000):
g2 := odeplot(sol,[t,x2(t)],range0,color = green, numpoints = 2000):
display(g1,g2);
Решения, полученные Maple:
Из вышеприведенных графиков видно, что Maple, в отличие от Infinity,
не имеет средств для корректного преодоления разрывов второго рода в
решениях.
Copyright © 1998-2004,
MathRevolt.
All rights reserved
|