Метод Рунге-Кутта с адаптивным шагом

Решим ту же систему дифференциальных уравнений методом Рунге-Кутта с адаптивным шагом.

Вектор правых частей уравнения и вектор начальных условий уже сформированы выше в примере решения системы дифференциальных уравнений методом Рунге-Кутта

Вызов функции rkadapt():

Результат работы функции – матрица X_rkа(таблица) значений, первый столбец которой содержит значения времени, второй искомые значения x₁(t), третий – x₂(t).

Для построения графика функций предварительно зададим интервал изменения индекса:

Рис. 3.3. Сравнение численного решения методом Рунге-Кутта с адаптивным шагом (штрихи) и аналитического (сплошные линии) решения системы дифференциальных уравнений