On High Order Barycentric Root-Finding Methods

To approximate a simple root of a real function f we construct a family of iterative maps, which we call Newton-barycentric functions, and analyse their convergence order. The performance of the resulting methods is illustrated by means of numerical examples.


INTRODUCTION
The classical Newton's iterative scheme for approximating the roots of an equation has been generalized by many authors in order to define iterative maps going from cubical to arbitrary orders of convergence (see for instance 1), (2), (3), (5), (6), (8), (9), (10), (11), (12), (13 , and references therein). In particular we mention the third-order method introduced in 13 and a family of high-order methods based on quadrature rules, recently described in 3 .
In that work the authors have introduced a recursive algorithm for constructing iterative maps. The first iterative function of this family is the classical Newton's method, which can be seen as the result of applying the rectangles rule to approximate a certain integral. As it is well-known, Newton's method has, in general, second order of convergence (when applied to the computation of simple roots). In the referred work 3 , the authors have shown that iterative maps of arbitrarily high order can be obtained if quadrature rules of higher degree are used (instead of the rectangles rule) to compute the mentioned integral. The properties of these methods were analysed and their convergence has been illustrated by several numerical examples.
In the present work we also construct a family of high order iterative maps, starting with the Newton's method, but we follow a different approach. All the iterative maps t to be considered have the common structure t(x) = x -[ϕ(x)] -1 f(x), where f is a real function with a simple zero at z. The function ϕ will be called a model function. This model function depends on f and on another function h which we name step function.
The family of iterative maps to be discussed is constructed by choosing a certain model function ϕ. We remark that our approach does not follow the traditional path for generating iterative maps by direct or inverse hyperosculatory interpolation (see, for instance 12 ) or Taylor expansions around a zero of f (see for instance 10 ).
Given a real -valued function f defined on an open set D ⊂ ℝ, we assume that f is sufficiently smooth in a neighbourhood of a simple zero z of f. In what follows we show how to construct a family of iterative maps t generating a sequence x i + 1 = t(x i ), i = 0, 1, ..., converging locally to z.

On High Order Barycentric Root-Finding Methods
Maps of higher order of convergence can produce accurate approximations with a smaller number of iterations. However, in general, they have expressions of increasing complexity, and so augmenting its computational cost. In this paper we show that at least for some Newton-barycentric methods the gain in accuracy compensates the computational cost. Like other root-finding methods these iterative maps can be easily extended to the case of multivariate functions. Numerical experiments have been carried out and the results are promising. The corresponding theory is still under construction.
The paper is organized as follows. In Section 2 we introduce the notion of a model function ϕ and prove that the iterative map t(x) = x -[ϕ(x)] -1 f(x) has a certain order of convergence (see Proposition 2.1). This is the main idea for the construction of recursive families of iterative functions.
In Section 3 we describe in detail one of these families, which we call the Newton-barycentric iterative maps. We begin by specifying the choice of the model function in the case of this family (Subsection 3.1). Then we describe the algorithm for obtaining iterative maps, based on this model function and we give the formulae of some of these maps (Subsection 3.2). In Section 4 we analyse the application of some Newton-barycentric formulas to some numerical examples. We finish in Section 5 with some conclusions.

RECURSIVE FAMILIES OF ITERATIVE MAPS
In this section we describe a procedure for constructing recursive families of higher order iterative methods. Proposition 2.1. Let z be a simple zero of a function f : D ⊂ ℝ ⟼ ℝ and ϕ a sufficiently smooth function in a neighbourhood of z, such that its derivatives ϕ ( i ) satisfy the j + 1 equalities (2.1) for j ≥ 0 a fixed integer. Then, for any initial value x 0 sufficiently close to z, the iterative process x k + 1 = t(x k ), k = 0, 1, ... , with

(2.2)
converges to z and its order of convergence is at least j + 2.
Proof. From (2.2) it is obvious that the zero z of f is a fixed point of the map t (that is, t(z) = z). Let us consider the function Δ t defined by Note that its derivatives are We now use induction on j to prove that the hypotheses in (2.1) imply that Δ t(z) = 0, Δ (1) t(z) = -1, Δ ( i ) t(z) = 0, for j ≥ 2, and consequently t has the referred order of convergence.

(2.3)
and applying the derivative operator to this equation, we have

(2.4)
Since Δt(z) = 0, and ϕ satisfies (2.1) with i = 0, it follows that f (1) which means that the iterative process generated by t has local order of convergence p at least 2. That is, p ≥ j + 2.
For an integer m ≥ 2, assume that and (2.5) 3) and the Leibniz's rule for the derivatives of the product, we have Thus, by the induction hypotheses, we obtain Hence the iterative map t m + 1 has local order of convergence p ≥ m + 2 and the proof is complete.
Remark 2.1. The well-known result on the local order of convergence of the Newton's map t( immediately from Proposition 2.1. It is enough to see that (2.1) is verified for j = 0, i.e. ϕ (0) (z) = f (1) (z), and so t has local order of convergence at least 2.
With the aim of analysing in detail the process of creating iterative maps, proposed in Proposition 2.1, we will now introduce the definitions of model function and step function. A function ϕ is called a model function of degree j, at x = z, if it satisfies the j + 1 conditions (2.1), but does not satisfy a similar condition for i = j + 1.
A function h is called a step function of degree j at x = z (or simply a step function) if it satisfies the following j + 1 equalities: but does not satisfy the condition corresponding to i = j + 1.
In the same way, let

THE NEWTON-BARYCENTRIC MAPS
We now consider a family of iterative maps, based on the model function ϕ k , defined as follows: ϕ k is the inner product of a constant vectorial function U k and a function V k depending only on the first derivative f (1) evaluated at x + i h(x), for i = 0, ..., k.
Namely, we take where the choice of a will be discussed below, and

(3.2)
where h is a step function of degree k. If one proves that ϕ k = 〈a, V k (x)〉 is a model function of degree k then, by Proposition 2.1, the respective process t k (x) = x -[ϕ k (x)] -1 f(x) has order of convergence at least k + 2.

Choice of the model function
The next proposition shows that ϕ k is a model function of degree k if and only if U k (x) = a is the unique solution of a non homogeneous linear system. Moreover, this solution represents the barycentric coordinates of ϕ k in a basis defined by the components of V k . The name Newton-barycentric maps reflects this property, as well as the fact that the recursive process starts with the Newton's method. Proposition 3.1. Let f be a function satisfying the hypotheses of Proposition 2.1, h a step function of degree k, k ≥ 0 a fixed integer and ϕ k = 〈U k , V k 〉, with U k and V k defined by (3.1) and (3.2). That is,

(3.3)
Then, the derivative of order k of V k , evaluated at x = z, is V k ( k ) (z) = D k R k , where D k and R k are the following (k + 1)×(k + 1)  = (a 0 , a 1 , a 2 , ..., a k ) is the (unique) solution of the linear system

(3.5)
In particular, this solution satisfies the equality

(3.6)
(ii) If f ( i ) (z) ≠ 0 for i = 1, ..., k, then the function U k = a represents the (normalized) barycentric coordinates of the model function ϕ k relative to the basis, Moreover, the iterative process generated by t k (x) = x -[ϕ k (x)] -1 f(x) has order of convergence at least k + 2.
Proof. For i = 0, 1, ..., k the derivatives of order i of V k , evaluated at x = z, are: So, the equalities (3.5) hold.
For (i), since ϕ k = 〈a, V k (x)〉, it is straightforward to verify that the conditions (2.1) for ϕ k to be a model function are equivalent to the system R k a = b. So, U k = a must be a solution of this system. As R k is nonsingular, this is the unique solution. Furthermore, since z is a simple zero of f, the equality (3.6) holds because it is just the first equation of the system R k a = b.

(3.7)
the only solution is α = 0. Differentiating (3.7) and evaluating at x = z, we obtain the homogeneous linear system which admits only the solution α = 0 since both the diagonal matrix and R k are nonsingular.
The last assertion follows from Proposition 2.1 since by item (i) ϕ k is a model function of degree k. □ The expressions for the first five barycentric maps are shown in Table 1.

Construction of the recursive family of iterative maps
We recall that a model function ϕ depends on a certain step function h. Now, for each model function ϕ j included in the definition of we use a step function which is defined recursively by h j (x) = t j -1 (x) -x. The starter t 0 will be taken to be the Newton's map t 0 ( . The next proposition shows that the iterative map t m , defined recursively in (3.8), has local order of convergence at least m + 2.
Proposition 3.2. Let z be a simple zero of a given function f and t 0 the Newton's map For a given natural number m ≥ 1, define recursively the step function h m and the iterative map t m by

(3.8)
where ϕ j is constructed using h j as step function and ϕ j is a barycentric-type map, given by Proof. It is only necessary to prove that each function h j is a step function of degree j and the statement follows from Proposition 3.1.
Let m ≥ 1 be an integer. As h m (0) (z) = 0, h m (1) (z) = -1 and for any integer i such that 2 ≤ i ≤ m , we have h m ( i ) (z) = 0, and so h m is a step function of degree m. □ The Newton-barycentric maps t k of arbitrary degree k are thus completely defined by Proposition 3.2. Let us remark that the same procedure to construct recursive families of iterative methods can be applied using other types of model and step functions. Moreover the idea of this method is easily extendable to the case of multivariate functions. Multidimensional analogs of the Newtonbarycentric maps have been implemented and applied to the solution of systems of nonlinear equations and the numerical results obtained so far are promising.

NUMERICAL EXAMPLES
Besides the theoretical interest of this new family of iterative functions, some of the above described methods are of practical interest. Though they are not optimal in the sense of the Kung and Traub's conjecture 7 , they are more efficient than the Newton's method, in the sense that with the same number of function evaluations a more accurate result can be produced. Concerning for example the Newton's method (iterative function t 0 ) and the method with the iterative function t 1 , we know that the first one requires 2 function evaluations at each iteration, while the second one requires 3. This means that two iterations of the second method require as many function evalutions (6) as three iterations of the first one. However the result produced by two iterations of the second method is in general much more accurate than the one produced by 3 iterations of the Newton's method. This happens, because the second method has at least convergence order three, and it follows that t 1 ○ t 1 (the composition of t 1 with itself) has at least convergence order 9; on the other hand, since the Newton method has in general convergence order 2, the method t 0 ○ t 0 ○ t 0 has just convergence order 8. Concerning the method t 2 (fourth order of convergence) each iteration requires 5 function evaluations; however as we shall see in the examples below, one single iteration of this method is often sufficient to obtain a result with an error less than 10 -10 (which can only be obtained with three iterates of the Newton's method).
Let us present some numerical results that illustrate the above properties. In all the cases we have applied the methods with the iterative functions t 0 , t 1 and t 2 , starting with a certain initial approximation x 0 and performing enough iterations to obtain an approximation with absolute error less than10 -10 . Then we compare the error of the last approximation, the number of iterations and the number of function evaluations by each method. Moreover we have obtained computational estimates of the convergence order, which confirm the theoretical predictions. All the computations were carried out in a personal computer using Mathematica 14 . The results are displayed in Tables 2, 3, 4 and they illustrate the advantage of using the methods with the iterative functions t 1 and t 2 for approximating the roots of real functions. f(x) = x 3 + 4x 2 -10, the initial approximation is x 0 = 1, the solution is z = 1.3652300134141, with 12 significant digits (see this example in 13 ).  Table 3: f(x) = cos( x) -x, the initial approximation is x 0 = 0.1, the solution is z = 0.739085133215, with 12 significant digits.

CONCLUSIONS
In this work we have introduced a family of high order iterative methods for the numerical solution of nonlinear equations. This family starts with the Newton's method as the basis of the recurrence process. Then each member of the family is build by a well defined procedure, forming a sequence of increasing convergence order.
The numerical examples presented in Section 4 demonstrate that the first iterative functions of these family (in particular t 1 and t 2 ) offer a good alternative to the Newton's method, taking account their accuracy and number of function evaluations. We remark that these methods can be easily combined with algorithms for the separation of real roots, such as recently described in 4 , providing an effective tool for detection and computation of real roots.