Используйте функцию
Jacob, чтобы вычислить расширенный якобиан для обыкновенного дифференциального уравнения (ОДУ), а затем предоставить ее как входные данные для решателей
Stiffb и
Stiffr.
1. Определите систему четырех неизвестных:
Переменные yi являются функциями x.
2. Определите константы системы.
3. Определите начальные значения.
4. Определите векторную функцию D(x,y), соответствующую правой части системы.
5. Используйте функцию
augment, чтобы создать расширенный якобиан:
6. Вызовите функции Stiffb и Stiffr:
Возвращенные матрицы содержат 5 столбцов, соответствующих числу точек и решениям для четырех неизвестных.
7. Извлеките решения для четырех неизвестных из возвращенных матриц Stiffb и Stiffr:
8. Постройте графики и сравните возвращенные решения для каждого неизвестного из двух функций:
Графики показывают, что два решателя ОДУ возвращают идентичные решения.
