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

Digits := 100
f(t) := 2*heaviside(t);
c := 1;
a := -3;
Complicated := 1
SupLocError := 1e-12
analysis({diff(x(t),t)+a*x(t) = f(t) + c*x(t)^2, x(0) := 0}, {x(t)}, {f(t)},{0},{2*(1-exp(t))/(exp(t) - 2)})

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

В выбранном масштабе область, содержащая точное решение, сливается с известным точным решением.

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

restart;
Digits := 50:
f(t) := 2;
c := 1;
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;
g1 := odeplot(sol1,[t,x(t)],range0,color = red):
g2 := plot(2*(1-exp(t))/(exp(t)-2),t=0..0.691317,color = green):
display(g1,g2);

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

Copyright © 1998-2004, MathRevolt. All rights reserved