1) Wstęp
Ekstrapolacja Riehardsona polega na obliczaniu dokładniejszej wartości całki na podstawie dwóch wartości wyznaczonych numerycznie dla różnych wielkości podprzedziałów, na jakie podzielony jest przedział całkowania. Metodę można użyć wtedy, gdy dysponujemy wzorem na błąd metody całkowania.
Zastosowanie ekstrapolacji Riehardsona do metody trapezów przybliżonego całkowania Błąd metody trapezów wyraża wzór
(1)
12
Niech (/*!) będzie przybliżoną wartością całki wyznaczoną metodą trapezów, a błędem obliczenia całki przy wielkości podprzedziału . Wówczas dokładną wartość całki można wyznaczyć z równania
/s /,(/!,) +R,(h1)s/1(h1) + ch1: (b-a)f ' ,
gdzie c =
12
jest stałą ( przy założeniu, że f jest niezależne od wielkości
podprzedziału h). Analogicznie, gdy/2( ^2) Jest przybliżoną wartością całki wyznaczoną metodą trapezów, a ft2(h2) błędem obliczenia całki przy wielkości podprzedziału h2, można napisać
/ = /2(h2) + ą(h2) = /2(h2) + ch* (3)
Po przyrównaniu prawych stron równań (2) oraz (3) otrzymujemy /,(/!,) + ctf S /2(/lJ + c/l2 (4)
Z równania (4) wyznaczamy c
C =
hf-hi
Po podstawieniu prawej strony równania (5) do równania (3) dostajemy
(6)
-1
Błąd zależności (6) jest rzędu o(hĄ), podczas g(h/ błąd metody trapezów jest rzędu o(h2).
Do obliczenia pewnej wielkości stosuje się metodę numeryczną z parametrem h. Wynikiem jej działania jest F(h). wartością dokładną jest F(0). Trudności obliczeniowe rosną gdy h maleje.
2