Квадратурная формула Симпсона

При шаге , отрезок разбивается на N=40 равностоящих узлов.

Как и следовало ожидать, наименьшее количество равностоящих узлов N=40 получается при вычислении интеграла по квадратурной формуле Симпсона.

Апостериорная оценка погрешности методом Рунге. Автоматический выбор шага интегрирования.

Пусть - квадратурная формула, примененная на частичном отрезке и имеющая порядок m, то есть . Для формул прямоугольников и трапеций m = 3, а для формулы Симпсона m = 5. Проведем на каждом частичном отрезке все вычисления дважды, один раз с шагом , другой раз с шагом . Тогда справедлива оценка:

Если для заданного правая часть не превосходит , то получим:

то есть будет достигнута заданная точность .

Если же на каком-то из частичных отрезков эта оценка не будет выполняться, то шаг на этом отрезке надо измельчить еще в два раза и снова оценить погрешность. Заметим, что для некоторых функций такое измельчение может продолжаться слишком долго. Поэтому в соответствующей программе следует предусмотреть ограничение сверху на число измельчений, а также возможность увеличения .