4.8.3  Die Trapezformel: Stückweise Approximation mit P1.

Man kann nun versuchen, die Trapezformel zu verbessern, indem man f(x) auf Δk = [xk,xk+1] durch das Lagrange-Polynom ersten Grades

P1(x) = x xk+1 xk xk+1f(xk) + x xk xk+1 xkf(xk+1)

ersetzt, k = 0,,n 1. Man sieht leicht, dass

xkxk+1 P1(x)dx = h f(xk+1) + f(xk) 2 .

Damit ergibt sich abf(x)dx = k=0n1xkxk+1 f(x)dx k=0n1xkxk+1 P1(x)dx k=0n1h f(xk) + f(xk+1) 2

und folglich

abf(x)dx b a n f(a) + f(b) 2 + f(x1) + + f(xk1) =: ST a,b(f).

Für die Fehlerabschätzung setzen wir voraus, dass f C2([a,b], ). Dann gilt nach (4.46)

f(x) P1(x) = f(2)(θ x) 2! (x xk)(x xk+1),x Δk,

für ein geeignetes θx Δk. Als Differenz zweier integrierbarer Funktionen ist die rechte Seite auf Δk integrierbar und R(Δk) := xkxk+1 f(x)dx xkxk+1 P1(x)dx = 1 2!xkxk+1 f(θx)(x xk)(x xk+1)dx.

Dabei gilt |R(Δk)| 1 2! max xΔk|f(x)|xkxk+1 (xk+1 x)(x xk)dx 1 12h3 max xΔk|f(x)|.

Daraus folgt abf(x)dx S T a,b(f) k=0n1|R(Δ k)| 1 12 b a n 3 n max x[a,b]|f(x)| 1 12 (b a)3 n2 max x[a,b]|f(x)|.

Wir merken an, dass diese Abschätzung nicht besser ist als die Fehlerabschätzung für die Rechteckformel! Die Trapezformel liefert also keinen wesentlichen Vorteil gegenüber der Rechteckformel.