Ввод данных для 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)})
Решения, полученные Infinity:
Ввод данных для 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);
Решения, полученные Maple:
