Исправленный метод Эйлера решения задачи Коши для ОДУ 1 порядка

В усовершенствованном методе Эйлера усреднялись наклоны касательных, т.е. производные от искомой функции. В модифицированном методе происходит усреднение точек (рис.7.8).

Прямая L1 есть касательная к истинной кривой y=y(x) в точке (xm, ym). Ее наклон к оси OX равен углу Исправленный метод Эйлера решения задачи Коши для ОДУ 1 порядка - student2.ru , для которого

Исправленный метод Эйлера решения задачи Коши для ОДУ 1 порядка - student2.ru ,

или в силу (7.2):

Исправленный метод Эйлера решения задачи Коши для ОДУ 1 порядка - student2.ru .

Исправленный метод Эйлера решения задачи Коши для ОДУ 1 порядка - student2.ru

Рис.7.8. Модифицированный метод Эйлера

Прямая L2 есть касательная к решению уравнения (7.2) в точке Исправленный метод Эйлера решения задачи Коши для ОДУ 1 порядка - student2.ru , являющей­ся пересечением L1 c прямой x = xm+h/2. Наклон L2 равен углу Исправленный метод Эйлера решения задачи Коши для ОДУ 1 порядка - student2.ru , для которого

Исправленный метод Эйлера решения задачи Коши для ОДУ 1 порядка - student2.ru .

Прямая L параллельна прямой L2 и проходит через точку (xm, ym), а ее пересечение c пря­мой x = xm+1 и определяет окончательное значение ym+1 решения уравнения в тоске xm+1. Уравнение прямой L можно записать в виде

Исправленный метод Эйлера решения задачи Коши для ОДУ 1 порядка - student2.ru ,

где Исправленный метод Эйлера решения задачи Коши для ОДУ 1 порядка - student2.ru .

Поэтому

Исправленный метод Эйлера решения задачи Коши для ОДУ 1 порядка - student2.ru   (7.12)

Выражение (7.12) есть вычислительная формула модифицированного метода Эйлера.

Он также согласуется с разложением в ряд Тейлора с точностью до h2. Блок-схема этого алгоритма аналогична предыдущей и отличается лишь формулой в блоке «ордината следующей точки».

Наши рекомендации