Запишем формулу Симпсона в общем виде:

(1)

(2)

Погрешность формулы Симпсона - метода четвертого порядка:

, (3)

Так как метод Симпсона позволяет получить высокую точность, если не слишком велика. В противном случае метод второго порядка может дать большую точность.

Например, для функции форма трапеции при для дает точный результат , тогда как по формуле Симпсона получаем

3. Геометрическая иллюстрация

На отрезке длиной 2h строится парабола, проходящая через три точки ,. Площадь под параболой, заключенная между осью OX и прямыми, принимают равной интегралу.

Особенностью применения формулы Симпсона является тот факт, что число разбиений отрезка интегрирования - четное.

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

(4)

Это формула Симпсона «трех восьмых».

Для произвольного отрезка интегрирования формула (4) может быть «продолжена»; при этом число частичных отрезков должно быть кратно трем ( точек).

, m=2,3, . (5)

- целая часть

Можно получить формулы Ньютона-Котеса старших порядков :

(6)

- количество отрезков разбиения;

- степень используемого полинома;

- производная -го порядка в точке ;

- шаг разбиения.

В таблице 1 выписаны коэффициенты . Каждая строка соответствует одному набору промежутков узлами для построения многочлена k-ой степени. Чтобы воспользоваться этой схемой для большего количества наборов (например, при k=2 и n=6), нужно «продолжить» коэффициенты, а затем сложить их.

Таблица 1:

Алгоритм оценки погрешности формул трапеции и Симпсона можно записать в виде: (7),

где - коэффициент, зависящий от метода интегрирования и свойств подынтегральной функции;

h - шаг интегрирования;

p - порядок метода.