Root Finding Methods
There are not many general tools available to determine in advance whether a root-finding problem can be solved. For our purposes, the main goal is to obtain information about the existence of a root, and if a root does exist, it is important to estimate an interval in which it lies.
To locate such a root, two main methods can be used:
1. Graphical Method
One of the first attempts is to plot the graph of the function. Indeed, if the goal is to solve \( f(x) = 0 \) and the graph of \( f \) can be drawn, then the intersection of this graph with the x-axis (axis of \( x \)) gives an approximate location of a root.
There is nothing wrong with using this method, but it can be difficult to apply when the graph is complex. Sometimes, it is possible to write \( f = f_1 - f_2 \), where \( f_1 \) and \( f_2 \) are simple functions. Then, we look for the intersection points of their graphs, the x-coordinates of which are exactly the roots of \( f \).
Example 3
Consider the function defined by:
\[ f(x) = \frac{\log x - 1}{x}, \quad x \in \left] 0, +\infty \right[ \]
This function has a unique root in the interval \( \left] \frac{3}{2}, 2 \right[ \). Let \( f_1(x) = \log x \) and \( f_2(x) = 1 \), then solving \( f(x) = 0 \) is equivalent to \( f_1(x) = f_2(x) \).
The intersection point of the graphs of \( f_1 \) and \( f_2 \) corresponds to the desired root.
2. Intermediate Value Theorem Method
Another method is to verify that the function \( f \) is continuous on an interval \( [a, b] \). If two points \( a \) and \( b \) can be found such that:
\[ f(a) > 0 \quad \text{and} \quad f(b) < 0 \]
Then, by the Intermediate Value Theorem, there exists a point \( c \in (a, b) \) such that \( f(c) = 0 \).
Finding such points \( a \) and \( b \) may require intuition, trial and error, reasoning, or graphical tools.
Example
Consider the following polynomial:
\[ p(x) = x^4 - x^3 - x^2 + 7x - 10 \]
Here are the values of the polynomial for various \( x \):
| \( x \) | -4 | -3 | -2 | -1 | 0 | 1 | 2 | 3 | 4 |
|---|---|---|---|---|---|---|---|---|---|
| \( p(x) \) | 266 | 68 | 4 | -16 | -10 | -4 | 8 | 56 | 194 |
The table shows that the polynomial \( p \) changes sign between:
- \( x = -2 \) and \( x = -1 \) → a root in \( ]-2, -1[ \)
- \( x = 1 \) and \( x = 2 \) → a root in \( ]1, 2[ \)
Numerical Methods
In the following sections, we will present several numerical solution methods. All of these methods are iterative, meaning they start from an initial estimate and produce a sequence of approximations that is supposed to converge to the solution.
In some cases, this sequence converges to a limit. It must then be verified whether this limit is indeed a solution of the equation \( f(x) = 0 \). Another important question is the speed of convergence of the method.