
Ввод данных для Infinity:
Digits := 200
f(t) := 2*heaviside(t);
c := 1;
a := -2;
Complicated := 1
SupLocError := 1e-15
analysis({diff(x(t),t)+a*x(t) = f(t) + c*x(t)^2, x(0) := 0}, {x(t)},
{f(t)},{0},{tan(t+Pi/4)-1})
Решение, полученное Infinity:
В выбранном масштабе область, содержащая точное
решение, сливается с известным точным решением.
Ввод данных для Maple:
restart;
Digits := 50:
f(t) := 2;
c := 2;
a := -3;
sol1 := dsolve({diff(x(t),t)+a*x(t) = f(t) + c*x(t)^2, x(0) = 0}, {x(t)},
type =
numeric,method=lsode, abserr = Float(1, -20));
with(plots):
range0 := 0..10;
odeplot(sol1,[t,x(t)],range0,color = green, numpoints = 2000);
Решения, полученные с помощью Maple:

Точные решения:
plot(tan(t+Pi/4)-1, t = 0..10, y = -10..10);
Из вышеприведенных графиков видно, что Maple, в
отличие от Infinity, не имеет средств для корректного преодоления разрывов
второго рода в решениях
Copyright © 1998-2004,
MathRevolt.
All rights reserved
|