Data input for 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)})

Solutions received by Infinity:

In the chosen scale the area containing the precise solution merges with the known precise solution.

Data input for 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);

Solutions received by Maple:


Note: If you solved your particular task with the help of Infinity and you want to share your know-how, send us the solution. We will publish it at our site and put it into the help file.