A novel Newton-Raphson style root finding algorithm

Komi Agbalenyo,Vincent Cailliez,Jonathan Cailliez

2023-09-02

Abstract:Many problems in applied mathematics require root finding algorithms. Unfortunately, root finding methods have limitations. Firstly, regarding the convergence, there is a trade-off between the size of it's domain and it's rate. Secondly the numerous evaluations of the function and its derivatives penalize the efficiency of high order methods. In this article, we present a family of high order methods, that require few functional evaluations ( One for each step plus one for each considered derivative at the start of the method), thus increasing the efficiency of the methods.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the trade - off between convergence and efficiency in existing root - finding algorithms (such as the Newton - Raphson method). Specifically, although high - order convergence methods increase the convergence speed, they require more calculations of the function and its derivatives, thus reducing the practical efficiency. In addition, the convergence domains of traditional methods are limited for certain functions. For example, for the arctangent function (atan), when the initial guess value is greater than 1.39, the traditional method fails to converge.

To solve these problems, the paper proposes a new class of high - order root - finding algorithms. These algorithms maintain high - order convergence while requiring fewer function evaluations (only one function evaluation per iteration and a one - time calculation of the required derivatives at the start of the method). This not only improves the efficiency of the method but also expands its convergence domain.

### Main contributions
1. **Improve efficiency**: The new method significantly improves efficiency by reducing the number of function evaluations in each iteration.
2. **Expand convergence domain**: The new method can converge within a wider range of initial guess values.
3. **Theoretical derivation**: The paper elaborately derives second - order, third - order and fourth - order convergence methods and shows how to generalize them to any order.

### Key formulas
- **Second - order convergence method**:
\[
U_{n + 1}=U_n-\frac{g(U_n)}{g'(l)}
\]

- **Third - order convergence method**:
\[
U_{n + 1}=U_n-\frac{g(U_n)}{g'(l)}+\frac{g(U_n)^2g''(l)}{2g'(l)^3}
\]

- **Fourth - order convergence method**:
\[
U_{n + 1}=U_n-\frac{g(U_n)}{g'(l)}+\frac{g(U_n)^2g''(l)}{2g'(l)^3}-\frac{g(U_n)^3[3g''(l)^2 - g'(l)g'''(l)]}{6g'(l)^5}
\]

Through these improvements, the method proposed in the paper shows higher efficiency and a wider convergence domain than existing methods in numerical experiments.