Главная
Infinity Online Help
Введение
О программе
Лицензионное соглашение
Как купить
Описание программы
Общее представление
Операции
Функции
Объкты
Переменные
Интерфейс
Графики
От простого к сложному
Примеры
Решения-утки
Задача Алейникова
Брюсселятор
Задача двух тел
Нелинейная автономная стационарная система с локально неустойчивым участком
Задача с несколькими разрывами второго рода
Задача с разрывом второго рода
Регистрация
Поддержка
  Задача с разрывом второго рода



Ввод данных для 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)})



Решения, полученные Infinity:



Ввод данных для 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);

Решения, полученные Maple:


Из вышеприведенных графиков видно, что Maple, в отличие от Infinity, не имеет средств для корректного преодоления разрывов второго рода в решениях.

Copyright © 1998-2004, MathRevolt. All rights reserved