Of all the methods to find the root of a function f(x) = 0, the Regula Falsi method is the oldest one. Being a closed bracket method, it is similar in many ways to the bisection method. Here, the algorithm of regula falsi method has been presented along with its flowchart and features.
Regula falsi method is also known by the name of false position method. Interpolation is the approach of this method to find the root of nonlinear equations by finding new values for successive iterations. In this method, unlike the secant method, one interval always remains constant.
The overall programming effort has been made easy for this method with the simple and easy-to-understand algorithm and flowchart given below. The convergence is always guaranteed in this method and is of first order. The false position method may be slow, but it is found superior to the bisection method in many ways. The iterative formula used here is:
[highlight color=”yellow”]x = [x0*f(x1) – x1*f(x0)] / (f(x1) – f(x0))[/highlight]
Features of Regula Falsi Method:
- No. of initial guesses – 2
- Type – closed bracket
- Convergence – linear
- Rate of convergence – slow
- Accuracy – good
- Approach – interpolation
- Programming effort – easy
Regula Falsi Method Algorithm:
- Read values of x0, x1 and e
*Here x0 and x1 are the two initial guesses
e is the degree of accuracy or the absolute error i.e. the stopping criteria*
- Computer function values f(x0) and f(x1)
- Check whether the product of f(x0) and f(x1) is negative or not.
If it is positive take another initial guesses.
If it is negative then goto step 5.
x = [x0*f(x1) – x1*f(x0)] / (f(x1) – f(x0))
- Check whether the product of f(x1) and f(x) is negative or not.
If it is negative, then assign x0 = x;
If it is positive, assign x1 = x;
- Check whether the value of f(x) is greater than 0.00001 or not.
If yes, goto step 5.
If no, goto step 8.
*Here the value 0.00001 is the desired degree of accuracy, and hence the stopping criteria.*
- Display the root as x.
Regula Falsi Method Flowchart:
The regula falsi (false position) method algorithm and flowchart given above are not exactly the same, only the approach to the method is same. A little modification to the iteration formula has been done in the flowchart. The two initial guesses should be of opposite nature for this method.