* Данный текст распознан в автоматическом режиме, поэтому может содержать ошибки
Численное решение ДУ в системе Mathcad
761
пустимый интервал между точками, в которых ищется решение (шаг интегриро вания). Обратите внимание, что функции, начинающиеся с малой буквы, дают ре шения только для конечной точки. Различаются функции также и методом реше ния. Если решение расходится, что может привести к переполнению разрядной сетки чисел, то прежде всего надо попытаться уменьшить шаг интегрирования. Матрица функция Якоби J, фигурирующая в этих функциях, имеет размер n ? (n+1) и представляется в виде:
Чаще всего эту матрицу несложно вычислить вручную или с помощью сим вольного дифференцирования. Фрагмент документа Mathcad с простым примером решения жесткой системы из двух дифференциальных уравнений показан на рис. 9.37. Для решения системы характерно наличие двух резко отличающихся стадий – крутой и пологой (или быстрой и медленной, если под независимой переменной подразумевать время).
Рис. 9.37. Пример решения жесткой системы дифференциальных уравнений