Formule van Newton-Cotes

Formule van Newton-Cotes voor $n=2$

Een formule van Newton-Cotes, genoemd naar de bedenkers Isaac Newton en Roger Cotes, is een formule voor de numerieke benadering van een integraal als gewogen som van functiewaarden in equidistante punten. De basisgedachte daarbij is de integrand te benaderen door een polynoom en de benaderende polynoom exact te integreren.

Definitie

Gegeven is de op het interval $[a,b]$ integreerbare functie $f.$ Het interval wordt opgedeeld in $n$ deelintervallen van gelijke lengte, de zogeheten stapgrootte $h=(b-a)/n$ door:

$a\le x_0<x_1<\dots <x_n\le b$

Het getal $n$ heet de orde van de benadering. De functie $f$ wordt benaderd met behulp van de bij de opdeling behorende Lagrange-polynomen ${\ell }_i$door:

$f(x)\approx P(x)=\sum _{i=0}^n{\ell }_i(x)f(x_i)$

De benadering voor de integraal van $f$ over $[a,b]$ is:

$\underset{a}{\overset{b}{\int }}f(x)dx\approx \underset{a}{\overset{b}{\int }}P(x)dx=\underset{a}{\overset{b}{\int }}\sum _{i=0}^n{\ell }_i(x)f(x_i)dx=\sum _{i=0}^n\underset{a}{\overset{b}{\int }}{\ell }_i(x)dxf(x_i)=\sum _{i=0}^n{w}_{i}f(x_i),$

dus als gewogen som van de functiewaarden in de deelpunten, met als gewichtsfactoren:

$w_i=\underset{a}{\overset{b}{\int }}{\ell }_i(x)\mathrm{d}x$

In deze definitie worden de eindpunten $a$ en $b$ in de benadering meegenomen: we spreken van een gesloten formule. Op analoge wijze kunnen ook zogenaamde open formules afgeleid worden, waarin de eindpunten $a$ en $b$ geen deelpunten zijn. Er zijn daarvoor verschillende methoden in gebruik. Sommige gebruiken de punten $x_1,\dots ,x_n$ zoals boven, andere kiezen

$x_1=a+\begin{array}{c}\frac{1}{2}\end{array}h,x_2=x_1+h,\dots ,x_n=x_{n-1}+h=b-\begin{array}{c}\frac{1}{2}\end{array}h$

In de tabel staan voor de gesloten formules de genormeerde gewichtsfactoren:

${\omega }_i=\frac{w_i}{b-a}=\frac{w_i}{nh}$

orde n	methode	gewichtsfactoren ωi
1	trapeziumregel	$\frac{1}{2}\frac{1}{2}$
2	regel van Simpson	$\frac{1}{6}\frac{4}{6}\frac{1}{6}$
3	3/8 - regel	$\frac{1}{8}\frac{3}{8}\frac{3}{8}\frac{1}{8}$
4	regel van Boole	$\frac{7}{90}\frac{32}{90}\frac{12}{90}\frac{32}{90}\frac{7}{90}$
5		$\frac{19}{288}\frac{75}{288}\frac{50}{288}\frac{50}{288}\frac{75}{288}\frac{19}{288}$
6		$\frac{41}{840}\frac{216}{840}\frac{27}{840}\frac{272}{840}\frac{27}{840}\frac{216}{840}\frac{41}{840}$

Opmerking

De regel van Boole wordt ook wel regel van Bode genoemd. Dit berust op een typefout uit het verleden waarbij Boole gelezen is als Bode.

De gewichtsfactoren voor hogere ordes zijn:

41 216 27 272 27 216 41 / 840

751 3577 1323 2989 2989 1323 3577 751 / 17280

989 5888 −928 10496 −4540 10496 −928 5888 989 /28350

2857 15741 1080 19344 5778 5778 19344 1080 15741 2857 / 89600

De oneven rijen hebben een lagere graad van algebraïsche nauwkeurigheid dan de even rijen. De hogere even rijen hebben negatieve coëfficiënten, wat de afrondingsfouten verhoogt. De regel van Simpson wordt verreweg het meest toegepast, vooral omdat hij eenvoudig is en toch redelijk nauwkeurig.

Voor een nauwkeurige benadering is het enerzijds van belang de orde van de benadering laag te kiezen en anderzijds dat de integrand niet te veel verandert over de deelintervallen. Deze eisen zijn tegenstrijdig, zodat men het integratiegebied wel opdeelt in subintervallen en in elk subinterval weer een lage orde benadering bepaalt. Voorbeelden hiervan zijn de regel van Milne en de regel van Weddle.

Berekening van de gewichten

Aan een voorbeeld is goed te zien hoe de gewichtsfactoren berekend worden. Zonder de algemeenheid te schaden kiezen we $a=0$ en $b=1.$ We bekijken het geval $n=3,$ dus $h=1/3.$

${\ell }_0(x)=-\frac{9}{2}(x-h)(x-2h)(x-1)$

${\ell }_1(x)=\frac{27}{2}x(x-2h)(x-1)$

${\ell }_2(x)=-\frac{27}{2}x(x-h)(x-1)$

${\ell }_3(x)=\frac{9}{2}x(x-h)(x-2h)$

$w_0=\underset{0}{\overset{1}{\int }}{\ell }_0(x)\mathrm{d}x=-\frac{9}{2}\underset{0}{\overset{1}{\int }}(x-h)(x-2h)(x-1)\mathrm{d}x=\frac{1}{8}$

$w_1=\underset{0}{\overset{1}{\int }}{\ell }_1(x)\mathrm{d}x=\frac{27}{2}\underset{0}{\overset{1}{\int }}x(x-2h)(x-1)\mathrm{d}x=\frac{3}{8}$

$w_2=\underset{0}{\overset{1}{\int }}{\ell }_2(x)\mathrm{d}x=-\frac{27}{2}\underset{0}{\overset{1}{\int }}x(x-h)(x-1)\mathrm{d}x=\frac{3}{8}$

$w_3=\underset{0}{\overset{1}{\int }}{\ell }_3(x)\mathrm{d}x=\frac{9}{2}\underset{0}{\overset{1}{\int }}x(x-h)(x-2h)\mathrm{d}x=\frac{1}{8}$

Voorbeeld

We benaderen de oppervlakte van een kwart cirkel boven de $x$-as, om 0 en met straal 1. We kiezen $n=4$ en berekenen:

$\underset{0}{\overset{1}{\int }}f(x)\mathrm{d}x=\underset{0}{\overset{1}{\int }}\sqrt{1-x^2}\mathrm{d}x\approx w_{0}f(0)+w_{1}f(0,25)+w_{2}f(0,5)+w_{3}f(0,75)+w_{4}f(1)=$

$\frac{1}{90}(7\cdot 1+32\cdot \frac{\sqrt{15}}{4}+12\cdot \frac{\sqrt{3}}{2}+32\cdot \frac{\sqrt{7}}{4}+7\cdot 0)=\frac{7+8\sqrt{15}+6\sqrt{3}+8\sqrt{7}}{90}=0,773$

Vergelijk het resultaat met de exacte waarde:

$\underset{0}{\overset{1}{\int }}\sqrt{1-x^2}\mathrm{d}x=\frac{\pi }{4}\approx 0,785$