Общая формула Симпсона и ее остаточный член

Пусть n=2m есть четное число и - значения функции для равноотстоящих точек с шагом . Применяя формулу Симпсона к каждому удвоенному промежутку длины 2h, будем иметь .

Следовательно, .

Отсюда получаем общую формулу Симпсона:

Введя обозначения , формулу можно записать в более простом виде:

Если функция непрерывно дифференцируема до четвертого порядка, то ошибка формулы Симпсона на каждом удвоенном промежутке дается формулой:

, где .

Суммируя все эти ошибки, получим остаточный член общей формулы Симпсона в виде:

непрерывна на отрезке [ a,b ], поэтому найдется точка такая, что .

Следовательно

, (8.9)

где .

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

, отсюда , т.е. h имеет порядок . Говорят, что степень точности метода Симпсона равна четырем

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

Предполагая, что на отрезке [ a,b ] производная меняется медленно, в силу формулы (8.9), получаем приближенное выражение для искомой ошибки

, где коэффициент M будем считать постоянным на промежутке интегрирования. Пусть и - приближенные значения интеграла , полученные по формуле Симпсона соответственно с шагом h и H=2h. Имеем: и . Отсюда