Решение обыкновенных дифференциальных уравнений (задача Коши), страница 2

Рис.  7.1. Применение функции odesolve для решения задачи с начальными                                                                      условиями для дифуравнения 31- го порядка (7.1) в Mathcad.                                                                                                                  tK = 25, m не задано, расчетная таблица выводится дополнительным вводом – рис.  Б

Правила применения функции odesolve(t, tК,число шагов m)

1.  Функции odesolve используется в решающем блоке с оператором Given (см. рис.  7.1).

2.  Переменная и число tК должна быть действительными.

3.  Число шагов m – необязательный параметр.

4.  Наивысшая производная дифуравнения dnx(t)/dtn должна быть линейной функцией.

5.  По умолчанию используется фиксированный шаг, но может устанавливаться и переменный (команда adaptive в меню от правой кнопки мышки на имени odesolve).

6.  Старшая производная должна быть линейной, низшие производные могут возводиться в степень и умножаться друг на друга.

7.  Ограничения – неравенства не допускаются; количество ограничений - равенств (записываемых с использованием логических операторов) должно быть не более n (n – порядок дифуравнения); они могут помещаться в любом месте решающего блока и должны быть независимы.

8.  Начальные условия должны быть простыми: у(0) = а1, у¢(0) = а2 и т.д.; сложные выражения типа у(0) +у¢ (0) = а3 не допускаются.

9.  Для записи дифуравнения можно использовать оператор дифференцирования d/dx , d2/dx2… и т.д. или применить запись с апострофами  y¢(x) , y¢¢(x)… (каждый апостроф вводится клавишами <Ctrl>+<F7> – см. рис. 7.1). В случае уравнения высокого порядка запись через оператор дифференцирования – предпочтительна (подсчет большого числа апострофов затруднителен - см. рис. 7.1).

10.   Для задачи с начальными условиями в начальной точке следует задавать у(0) и все ее n - 1 производные, причем, для производных следует применить запись с апострофами (запись с оператором дифференцирования для начальных условий Mathcad не воспринимает).