Рассмотрим тот же самый пример: дифференциальное уравнение , частное решение, удовлетворяющее условию , промежуток и его разбиение на 10 частей
( – длина каждой части).
Цель усовершенствования состоит в том, чтобы приблизить «красные квадратики» ломаной к соответствующим «зелёным точкам» точного решения .
И идея модификации такова: отрезки должны быть параллельны касательным, которые проведены к графику функции не на левых краях, а «посерединке» интервалов разбиения. Что, естественно, улучшит качество приближения.
Алгоритм решения работает в том же русле, но формула, как нетрудно догадаться, усложняется:
, где
Плясать вновь начинаем от частного решения и сразу же находим 1-й аргумент «внешней» функции:
Далее следуют уже знакомые по предыдущему параграфу вычисления , после чего можно рассчитать 2-й аргумент «внешней» функции: .
Теперь находим нашего «монстра», который на поверку оказался не таким уж и страшным – обратите внимание, что это ТА ЖЕ функция , вычисленная в другой точке:
Умножаем результат на шаг разбиения:
|
|
Таким образом:
Алгоритм заходит на второй круг, не поленюсь, распишу его подробно:
рассматриваем пару и находим 1-й аргумент «внешней» функции:
Рассчитываем и находим её 2-й аргумент:
Вычислим значение:
и его произведение на шаг:
Таким образом:
Далее рассматриваем пару и т.д.
Вычисления разумно провести в Экселе (растиражировав формулы по той же схеме – см. видеоролик выше), а результаты свести в таблицу:
Числа целесообразно округлять до 4-5-6 знаков после запятой. Нередко в условии той или иной задачи есть прямое указание, с какой точностью следует проводить округление. Я подровнял сильно «хвостатые» значения до 6 знаков.
По результатам 2-го и 3-го столбцов (слева) построим ломаную , и для сравнения я снова приведу график точного решения :
Результат существенно улучшился! – красные квадратики практически «спрятались» за зелёными точками точного решения.