De methode van Newton-Raphson
Prof.dr.ir. C. Vuik
Technische Universiteit Delft
e-mail: c.vuik@tudelft.nl
De methode van Newton-Raphson is een numerieke methode om het nulpunt
van een functie te bepalen. We zullen de methode uitleggen aan de hand
van het volgende probleem: gegeven de functie f bepaal het nulpunt
p zodanig dat f(p) = 0.
Stel p0 is een benadering van het nulpunt p. Maak nu de raaklijn
h(x) aan f door het punt (p0,f(p0)). De formule van h is:
h(x) = f(p0) + |
df dx
|
(p0) (x-p0). |
|
Om een nieuwe (betere) benadering te vinden van p bepalen we p1,
die het nulpunt is van h, dus h(p1) = 0. Invullen geeft:
0 = f(p0) + |
df dx
|
(p0) (p1-p0). |
|
Als we nu p1 naar links brengen en alle andere termen naar rechts,
dan krijgen we (onder de voorwaarde dat [(df)/(dx)](p0)
is ongelijk aan 0):
Zo kunnen we doorgaan om p2, p3,
... te vinden. De procedure
staat uitgelegd in het onderstaande plaatje.
Figuur: De Newton-Raphson methode
Voorbeeld
Stel we willen het positieve nulpunt bepalen van de functie f(x) = x2-2. Het
exacte antwoord in 6 decimalen is
1.41421. We nemen als startpunt
p0 = 1. De vergelijking van de raaklijn is nu
De nieuwe benadering van het nulpunt p1 is nu:
Zo doorgaande vinden we:
p1 = 1.50000
p2 = 1.41666
p3 = 1.41421
Hieruit blijkt dat de oplossing al na 3 stappen gevonden is.