Решение жесткого ОДУ методом 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) (для двух последних функций).
