Infinity Online Help - The problem with a discontinuity of the second kind

The problem with a discontinuity of the second kind

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

Solutions received by Infinity:

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

Solutions received by Maple:

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.