Data input for Infinity:

Digits := 50
l1 := 9.533*y1(t) + diff(y1(t),t) = 2 + y1(t)^2*y2(t)
l2 := -8.533*y1(t) + diff(y2(t), t) = - y1(t)^2*y2(t)
SupT := 10
SupLocError := 1e-15
Complicated := 1
analysis({l1,l2, y1(0) := 1, y2(0) := 4.2665}, {y1(t), y2(t)})

Solutions received by Infinity:


Data input for Maple:

SupT := 10;
l1 := 9.533*y1(t) + diff(y1(t),t) = 2 + y1(t)^2*y2(t):
l2 := -8.533*y1(t) + diff(y2(t), t) = -y1(t)^2*y2(t):
sol := dsolve({l1,l2,y1(0) = 1, y2(0) = 4.2665},{y1(t),y2(t)}, type = numeric, abserr = Float(1, -15));
with(plots):
range0 := 0..SupT;
g1 := odeplot(sol,[t,y1(t)],range0,color = green, numpoints = 2000):
g2 := odeplot(sol,[t,y2(t)],range0,color = red, numpoints = 2000):
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.