Regel van Simpson

Uit Wikipedia, de vrije encyclopedie
Ga naar: navigatie, zoeken
De regel van Simpson kan worden afgeleid door benadering van de integrand (blauw) met de kwadratische interpolant (rood).

De regel van Simpson is een benaderingsformule om de numerieke waarde van een integraal te berekenen. De regel, die ontwikkeld is door Thomas Simpson, is een speciaal geval van een formule van Newton-Cotes.

De regel van Simpson benadert de integraal van een functie over het interval als de integraal van de tweedegraadspolynoom die de functie benadert en in de eindpunten en het midden van het interval met de functie overeenkomt, d.w.z. en . De polynoom is door deze eis bepaald. Er gelden de volgende drie lineaire vergelijkingen:

waaruit

De integraal van is dan

dus een gewogen som van de functiewaarden in de eindpunten en het midden van het interval.

Met de regel van Simpson wordt de integraal dus benaderd door:

Voor een goede benadering is het van belang dat de variatie van de integrand over het interval niet te groot is. In praktijk wordt daarom het interval verdeeld in een aantal subintervallen, en wordt op elk subinterval de regel van Simpson toegepast. Als de intervallen tussen de functiewaarden van gelijke lengte zijn, wordt de benadering:

Als het interval gehalveerd wordt, verkleint de fout met een factor 32. De regel van Simpson is daarmee efficiënter dan de meer eenvoudige trapeziumregel. Deze beide methoden zijn ingebouwd in de methode van Romberg.