Метод Симпсона

Отрезок интегрирования разбивается на четное число 2n равных частей длины h=(b-a)/2n.

a=x₀ < x₁ < …< x_2n-1 < x_2n = b

Рассматривается пара соседних участков и через три точки кривой с координатами (x₀,y₀), (x₁,y₁), (x₂,y₂) проводится парабола с осью, параллельной оси Oy.

Ее уравнение y=Ax² +Bx +C. Площадь криволинейной трапеции на участке [x₀,x₂] заменяется площадью криволинейной трапеции, ограниченной параболой.

Если вынести за скобку общий множитель x₂-x₀ и привести к общему знаменателю, получится

(1)

Неизвестные коэффициенты A, B, C находятся из условия, что при значениях x равных x₀, x₁,x₂, функция f(x) принимает соответственно значения y₀, y₁, y₂. Если взять , то условия можно записать

(2)

Второе равенство умножается на четыре, все три равенства складываются, получается

(3)

Что совпадает с квадратной скобкой в правой части равенства (1). Если подставить (3) в (1), и заметить, что x₂ –x₁=2h, (h=(b-a)/(2n)), то получится

(4)

Для каждой следующей пары участков получится такая же формула

(5)

Суммируя равенства вида (4) и (5) по всем участкам можно получить

или, если обозначить N=2n

.

Это и есть формула Симпсона. Ее называют также формулой парабол.

Переменная с=1 для нечетных i и –1 для четных, т.о. реализуется коэффициент (-1) i -1.

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