Задача Коши в Mathcad для дифференциальных уравнений n-го порядка с одной неизвестной (обыкновенное дифференциальное уравнение - ОДУ) формулируется следующим образом. Найти решение дифференциального уравнения
в виде функции y=y(x), которая удовлетворяет заданным начальным условиям
где - заданные значение. Решение задачи Коши для обыкновенных дифференциальных уравнений второго и более высоких порядков можно свести к системе уравнений. Решение задачи Коши для ОДУ первого порядка в Mathcad с использованием различных функций приведено на листинге.
Наибольшее распространение для решения задачи Коши в Mathcad получил метод Рунге – Кутта. Суть метода состоит в последовательном отыскании искомого значения функции yi+1 по формуле
За h принимается достаточно малый шаг, с помощью которого весь интервал задачи Коши разбивается на дискретные точки, в которых и ищется решение. Погрешность результатов пропорциональна пятой степени шага (h5).
Геометрический смысл метода Рунге – Кутта состоит в следующем. Из очередной точки (xi,yi) выбирается направление (угол) α1, для которого tg(α1)=f(xi,yi). На этом направлении вычисляется точка с координатами Затем из точки (xi,yi) выбирается направление (угол) α2, для которого tg(α2)=f
На этом направлении вычисляется точка с координатами Далее из точки (xi,yi) выбирается направление (угол)α3, для которого
На этом направлении в Mathcad вычисляется точка с координатами После чего из точки (xi,yi) выбирается направление (угол) α4, для которого . Все четыре полученных направления усредняются в соответствии с формулой для расчета Δyi. На этом результирующем направлении и строится расчетная точка с координатами
Метод Рунге – Кутта благодаря высокой точности широко используется при численном решении дифференциальных уравнений и в частности в Mathcad. Существует несколько разновидностей данного метода, которые нашли свое отражение в рассмотренных выше функциях. На листинге можно не только сравнить результаты, полученные на основе различных функций, но и оценить эти результаты с позиций точности расчетов.
Путем сравнения результатов решения задачи, можно сделать вывод о точности решения задачи. Наиболее точный результат позволяет получить функция Rkadapt.