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

В методе Симпсона в каждой части деления подынтегральная функция аппроксимируется квадратичной параболой a₀x²+a₁x+a₂. В результате вся кривая подынтегральной функции на участке [a,b] заменяется кусочно-непрерывной линией, состоящей из отрезков квадратичных парабол. Приближенное значение интеграла I равно сумме площадей под квадратичными параболами.

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

L_k=2h.

В результате количество частей деления N2=n/2. Тогда n в методе Симпсона всегда четное число.

Определим площадь S₁ на участке [x₀, x₂] (рис.12.2).

Исходя из геометрического смысла определенного интеграла, площадь S₁ равна определенному интегралу от квадратичной параболы на участке [x₀, x₂]:

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

На основании этого условия строим систему линейных уравнений:

Решая эту систему, найдем коэффициенты параболы.

В результате имеем: ..

Для участка [x₂, x₄]: ..

:::::::::::::::::::

Для участка [x_i-1, x_i+1]: .,

где .

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

где

N2 - количество частей деления.

Точность метода Симпсона имеет порядок (h³/h⁴).

Схема алгоритма метода Симпсона представлена на рис.12.7.

Рис. 12.7. Схема алгоритма Симпсона (с автоматическим выбором шага)

Формула Симпсона (также Ньютона-Симпсона ^[1]) относится к приёмам численного интегрирования. Получила название в честь британского математика Томаса Симпсона (1710—1761).

Суть приёма заключается в приближении подынтегральной функции на отрезке интерполяционным многочленом второй степени , то есть приближение графика функции на отрезке параболой. Метод Симпсона имеет порядок погрешности 4 и алгебраический порядок точности4.