Regula falsi

Regula falsi of de methode van regula falsi is een algoritme uit de numerieke wiskunde om de nulpunten van een continue functie te bepalen. Het algoritme convergeert trager dan de methode van Newton-Raphson, maar is stabieler. De methode maakt gebruik van iteraties van het gezochte punt en combineert eigenschappen van halveringsmethode en de secant-methode.

Methode

Animatie

De regula falsi berekent opeenvolgende benaderingen in de vorm van een interval $[a_{n},\,b_{n}]$ waarin een nulpunt van de functie $f$ ligt door recursie en begint met twee punten in de buurt van het nulpunt waarvan de functiewaarden tegengestelde tekens hebben, die dus onder en boven de $x$-as liggen, dus gegeven een functie $f(x)$ met twee waarden $a_{0}$ en $b_{0}$ zodat $a_{0}<b_{0}$ en $f(a_{0})\ f(b_{0})<0$.

Het nulpunt $c_{n}$ bevindt zich daarna steeds ergens tussen de twee punten $a_{n}$ en $b_{n}$. Vervolgens wordt het snijpunt met de $x$-as van de lijn, die de punten op de grafiek van de functie bij de twee vorige punten $(a_{n},f(a_{n}))$ en $(b_{n},f(b_{n}))$ met elkaar verbindt, bepaald.

$c_{n}=b_{n}-f(b_{n}){\frac {b_{n}-a_{n}}{f(b_{n})-f(a_{n})}}$

Uit het teken van de functiewaarde in dit snijpunt wordt bepaald in welk interval, gevormd door dit punt en een van de vorige, het nulpunt ligt. Zo wordt het interval waarin zich het nulpunt bevindt steeds verkleind, dus zodat

$c_{n}=a_{n}-{\frac {a_{n}-b_{n}}{f(a_{n})-f(b_{n})}}f(a_{n})$

en

$a_{n+1}=c_{n},\ b_{n+1}=b_{n}\quad$ de ondergrens wordt opgeschoven als $f(a_{n})\ f(c_{n})>0$

of

$a_{n+1}=a_{n},\ b_{n+1}=c_{n}\quad$ de bovengrens wordt teruggeschoven als $f(b_{n})\ f(c_{n})>0$