Матричные вычисления в Mathcad
d7c8102a

Решение жесткого ОДУ методом RADAUS (продолжение листингов 9 8 и 9 9)



Рисунок 9.14. Решение жесткого ОДУ методом RADAUS (продолжение листингов 9.8 и 9.9)



В заключение приведем соответствующие встроенные функции, которые применяются для решения жестких систем ОДУ не на всем интервале, а только в одной заданной точке t1.

  •  radau(y0, t0, t1, асc, F, k, s) — метод RADAUS.
  •  stiffb(y0,t0,t1,acc,F, J, k, s) — метод Булирша—Штера.
  •  stiffr (y0,t0, t1, acc, F, J, k, s) — метод Розенброка.


Имена этих функций пишутся со строчной буквы, а их действие и набор параметров полностью аналогичны рассмотренным нами ранее для функций, относящихся к решению в заданной точке нежестких систем (см. разд. 9.3.2). Отличие заключается в специфике применяемого алгоритма и необходимости задания матричной функции якобиана J(t,y) (для двух последних функций).

 


Содержание раздела