Метод Гауса

Для отримання підвищеної точності для чисельного інтегрування користуються формулою Гауса

, (8.13)

в якої не фіксуються не тільки вузли інтерполяції x1, x2,…,xn, а й квадратурні коефіцієнти С1,…,Сn. При цьому Zn невідомих величин x1,…,xn; С1,…,Сn визначається із умови, що формула є точною у випадку будь-якого багаточлена 2n-1.

Таким чином, для будь-якого багаточлена (2n-1)-й степеню

(8.14)

повинна виконуватися рівність:

. (8.15)

Багаточлен f(x) степеню 2n-1 можна показати у вигляді

f(x)=F(x)Q(x)+R(x) , (8.16)

де F(x)-шуканий багаточлен n-ї степені, а Q(x) та R(X)- відповідно частинне ділення f(x) на F(x) та залишок від цього ділення, степінь багаточленів Q(x) та R(x) не перевищують (2n-1).

Вираз для F(x) можна записати так:

(8.17)

тут величини x1,…,xn - шукані абсциси формули Гаусса, а А12,…,Аn - постійні.

Оскільки шукана функція F(x) у вузлах x1,…,xn перетворюється на нуль, то

(8.18)

(8.19)

Але для багаточлена R(x) степеню не вище n-1 також повинна виконуватися рівність:

(8.20)

Bіднімаючи (8.20) від (8.19),отримаємо

Із останнього відношення можна визначити шукану функцію F(x). Оскільки ця рівність справедлива для якого-небудь багаточлена Q(x) степеню n-1

, (8.21)

то при будь-яких коефіцієнтах маємо таку систему рівнянь:

(8.22)

Підставляючи в (8.22) вирази для F(x) із формули (7.2 ) та інтегруючи, отримаємо для визначення коефіцієнтів А1, … Аn систему n рівнянь

(8.23)

з яких видно, що А1357=…=0 та, отже, шуканий многочлен має вигляд:

(8.24)

Відмітимо, що при парному n корені рівняння F(x)=0 попарно рівні за абсолютним значенням, але протилежні за знаком, а при непарному n коренем є також і х=0.

Визначивши із системи (8.23) коефіцієнти Аі (і=1,2,…,n), складемо рівняння F(x)=0 та знайдемо його корені х1,…,хn, тобто шукані абсциси формули Гаусса, а потім обчислимо коефіцієнти Сi (і=1,2,…,n) за формулою

(8.25)

Приклад 1. Побудувати квадратурну формулу Гаусса для випадку n=2 на відрізку інтегрування [-1, 1].

Розв‘язок. Загальний вигляд квадратурної формули Гаусса при n=2 та заданих межах інтегрування:

,

де підлягають визначенню квадратурні коефіцієнти с1 та с2, а також абсцис х1 та х2.

Для визначення абсцис складемо многочлен , коефіцієнти А1 та А2 якого знайдемо із системи вигляду (8.22)

підстановкою багаточлена F(x) у систему. Маємо

; ,

тобто А1=0, А2 =-1/3. Тоді , звідки

та

Коефіцієнти с1 та с2 обчислимо за формулою (8.25)

;

Отже,

Для обчислення інтеграла загального вигляду слід замінити змінні

, (i=1,2,….,n) (8.26)

Тоді формула Гаусса прийме вигляд

Значення квадратурних коефіцієнтів Гаусса сі (і=1,2...n) та абсцис хі (і=1,2...n) наведені в таблиці 8.3.

Таблиця 8.3 – Значення квадратурних коефіцієнтів Гаусса

n=1 x1=0,5 с1=2
n=2 -x1=x2=0,577350 с12=1
n=3 -x1=x3=0,774597, x2=0 с1=с3=0,555555, с2=0,888889
n=4 -x1=x4=0,861136, -x2=x3=0,339981 с14=0,347855, с23=0,652145
n=5 -x1=x5=0,906180, -x2=x4=0,538470 с1=с5=0,236927, с2=с4=0,478629
n=6 -x1=x6=0,932470, -x2=x5=0,661210, -x3=x4=0,238620 с1=с6=0,171324, с2=с5=0,360761, с3=с4=0,467914
n=7 -x1=x7=0,949108, -x2=x6=0,741531, -x3=x5=0,405845, x4=0 с1=с7=0,129485, с2=с6=0,279705, с3=с5=0,381830, с4=0,417960
n=8 -x1=x8=0,960290, -x2=x7=0,796666, -x3=x6=0,525532, x4=0,183434 с1=с8=0,101228, с2=с7=0,222381, с3=с6=0,313707, с4=с5=0,362684,