График решения в виде временныхзависимостей

Функция Rkadapt обычно дает более точный результат благодаря автоматическому изменению шага. Естественно, по скорости вычислений она проигрывает функции rkfixed, хотя и не всегда – если решение меняется медленно, это может привести к заметному уменьшению числа шагов. Таким образом, функция Rkadapt наиболее привлекательна для решения системдифференциальных уравнений, дающих медленно меняющееся решение. В последнем примере решение представлено в виде временных зависимостей.

Если решение системыдифференциальных уравнений имеет вид гладких функций, то вместо функции rkfixed, описанной ранее, целксообразно применять новую функцию:

Bulstoer (y, x1, n, F)

Она возвращает матрицу решения системы обыкновенных дифференциальных уравнений, правая часть которых (в виде первых производных неизвестных функций) записана в векторе F (х,у) при заданных в векторе у начальных условиях и при решении на интервале от х1 до х2 для nточек решения, не считая начальной точки.