
Data input for Infinity:
Digits := 200
f(t) := 2*heaviside(t); c := 2; a := -3;
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})
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 := 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);
Solutions received by Maple:
Precise solutions:
plot(tan(t+Pi/4)-1, t = 0..10, y = -10..10);
From the above-given plots it is clear that, unlike
Infinity, Maple doesn’t possess any means of overcoming the
discontinuity of the second kind correctly.
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.
Copyright © 1998-2004,
MathRevolt.
All rights reserved
|