\require{AMSmath}
WisFaq - de digitale vraagbaak voor wiskunde en wiskunde onderwijs


Printen

Regel van Simpson ; praktische uitwerking

Geg. : I = 0ò1 ex dx (= e-1)

Ik loop vast bij het uitschrijven van beanderingen mbv de regel van Simpson:

I*2 = 1/6(f(0)+4f(1/2)+f(1))

I*4 = 1/12(f(0)+4f(?)+2f(?)+f(?)) 1,718318

Ik kan van alle soorten benaderingen m.b.v. de bijbehorende formules niet bepalen of ik moet starten met f(a) en ook met f(b) eindigen of dat ik net een stap daarvoor moet eindigen.

Dit zelfde geldt voor de trapezium / rechthoeksregel.

In de formules staat dan wel bijv. f(x0) f(x1) f(xn) bijvoorbeeld, maar welke functiewaarden horen daar nu bij?

Maikel
Student hbo - donderdag 27 oktober 2005

Antwoord

Bij I*2 benader je de oppervlakte onder de grafiek van f op het interval [0;1] met 1/6(f(0)+4f(1/2)+f(1))*(1-0). (Die 1-0 is de breedte van het deelinterval.)
Wanneer je het interval [0,1] in twee gelijke delen opsplitst [0,1/2] en [1/2,1] kun je de regel van Simpson op beide deelintervallen toepassen.
Je krijgt dan:
[0,1/2]: 1/6(f(0)+4*f(1/4)+f(1/2)*(1/2-0)=1/12(f(0)+4*f(1/4)+f(1/2)
[1/2,1]:1/6(f(1/2)+4*f(3/4)+f(1)*(1-1/2)=1/12(f(1/2)+4*f(3/4)+f(1)

Neem je deze twee samen dan krijg je:
1/12(f(0)+4f(1/4)+f(1/2)+f(1/2)+4f(3/4)+f(1))=
1/12(f(0)+4f(1/4)+2f(1/2)+4f(3/4)+f(1))

Zou je [0;1] in vier deelintervallen verdelen en op elk deel de regel van Simpson toepassen dan krijg je dus
[0;1/4]: 1/24(f(0)+4f(1/8)+f(1/4))
[1/4;1/2]: 1/24(f(1/4)+4f(3/8)+f(1/2))
[1/2;3/4]: 1/24(f(1/2)+4f(5/8)+f(3/4))
[3/4;1]: 1/24(f(3/4)+4f(7/8)+f(1))

Tezamen:
1/24(f(0)+4f(1/8)+2f(1/4)+4f(3/8)+2f(1/2)+4f(5/8)+2f(3/4)+4f(7/8)+f(1)).

Op een soortgelijke manier kun je uitzoeken hoe het bij de andere benaderingsalgoritmes werkt.


donderdag 27 oktober 2005

©2001-2024 WisFaq