On the Local Convergence of a Third Order Family of Iterative Processes

Efficiency is generally the most important aspect to take into account when choosing an iterative method to approximate a solution of an equation, but is not the only aspect to consider in the iterative process. Another important aspect to consider is the accessibility of the iterative process, which shows the domain of starting points from which the iterative process converges to a solution of the equation. So, we consider a family of iterative processes with a higher efficiency index than Newton’s method. However, this family of proecsses presents problems of accessibility to the solution x∗. From a local study of the convergence of this family, we perform an optimization study of the accessibility and obtain iterative processes with better accessibility than Newton’s method.


Introduction
In general, the roots of a nonlinear equation cannot be expressed in a closed form and this problem is commonly carried out by applying iterative methods.So, starting from one or several initial approximations of a solution x * of Equation ( 1), a sequence {x n } of approximations is constructed so that it converges to x * .We can get the sequence {x n } in different ways, depending on the iterative method that is applied.The best-known iterative scheme is Newton's method, Observe that we need the operator F to be differentiable Fréchet in order to apply Newton's method.
Efficiency is generally the most important aspect to take into account when choosing an iterative method to approximate a solution of Equation (1) (see [1][2][3][4]).The classic efficiency index, defined by Ostrowski in [5], provides a balance between the order of convergence ρ and the number of functional evaluations d: I = ρ 1 d .However, it is interesting to note that if Equation (1) represents a system of n nonlinear equations, the computational cost of evaluating the operators F and F is not similar, as it happens in the case of scalar equations.So, for one evaluation of F , n functional evaluations are required.Meanwhile, the evaluation of the associated jacobian matrix F requires n 2 functional evaluations, so that the evaluations of F and F cannot be considered in the same way.Therefore, this efficiency index must be modified to consider the number of functional evaluations in order for it to be a good efficiency measurement for an iterative process in the multivariate case.
Taking into account this extension of the efficiency index from the Chebyshev method with cubical convergence, Ezquerro and Hernández define in [6] a family of multipoint iterative processes with cubical convergence under the usual conditions required of the operator F .This family of iterative processes is given by and has better index of efficiency than Newton's method Equation (2).In this paper, we focus our attention on the analysis of the local convergence of sequence Equation (3) in Banach spaces under weak convergence conditions.In fact, we only consider conditions for the first derivative of the operator F , as in the case of the Newton method [7].
In any case, the efficiency index is not the only aspect to compare iterative processes.Another important aspect to take into account is the domain of accessibility associated to the iterative process, which is defined as the domain of starting points from which the iterative process converges to a solution of the equation.The location of starting approximations, from which the iterative methods converge to a solution of the equation, is a difficult problem to solve.This location is from the study of the convergence that is made of the iterative process.The analysis of the accessibility of an iterative process from the local study of the convergence is based on demanding conditions to the solution x * , from certain conditions on the operator F , and provide the so-called ball of convergence of the iterative process, that shows the accessibility to x * from the initial approximation x 0 belonging to the ball.This study requires conditions on the operator F .As we indicate previously, family Equation (3) has higher efficiency index than the Newton method, but it presents problems of accessibility to the solution x * of Equation ( 1) if we consider a semilocal study of the convergence (see [6]).
The paper is organized as follows.In Section 2, a local convergence result is provided.In Section 3, from the previous local convergence result, we observe that the size of the ball of convergence of family Equation (3) depends on the parameter p ∈ (0, 1] and, then, we do an optimization study.Moreover, we realize a comparative study from the point of view of the accessibility or, equivalently, from the ball of convergence with the Newton method.Finally, in the last section, a numerical test confirms the theoretical results obtained.

Local Convergence
As we indicated above, the local convergence results for iterative methods require conditions on the operator F and the solution x * of Equation ( 1).Note that a local result provides what we call ball of convergence and denote by B(x * , R).From the value R, the ball of convergence gives information about the accessibility of the solution x * .
An interesting local result for Newton's method Equation ( 2) is given in [7] by Dennis and Schnabel, where the following conditions are required: (C 1 ) Let x * be a solution of Equation ( 1) and exist r > 0, so that B(x * , r) ⊂ Ω and the operator Now, our aim is to establish the local convergence study for family Equation (3), with cubical convergence (see [6]), from conditions (C 1 ) and (C 2 ).So, we obtain a local convergence result in the same initial conditions that for the Newton method, which is a second order method.
Firstly, we provide two technical lemmas, where we obtain some results about the operator F and the sequences {x n }, {y n } and {z n }, which define family Equation (3).Lemma 1.Let us suppose that conditions (C 1 )-(C 2 ) are satisfied.Let us assume that exists S, a positive real number, with S ≤ r.If βKS < 1, then, for all x ∈ B(x * , S), the operator [F (x)] −1 exists and Taking into account the Banach lemma on invertible operators and conditions (C 1 )-(C 2 ), the proof follows.
Lemma 2. Fixed n ∈ N, we assume that x n , z n ∈ B(x * , r) and there exists Γ n .Then, we obtain the following decompositions Proof.From the Taylor series of F (x * ), (i n ) and (ii n ) are easily obtained.Now, by considering iterative scheme Equation (3), it follows (iii n+1 ): Now, we are able to obtain a result of local convergence to family Equation (3).
Theorem 1.Let F : Ω ⊆ X → Y be a nonlinear continuously differentiable operator on a non-empty open convex domain Ω of a Banach space X with values in a Banach space Y .Suppose that conditions (C 1 ) and (C 2 ) are satisfied.Then, fixed p ∈ (0, 1], there exists R > 0, such that the sequence {x n }, given by the iterative process of family Equation (3), corresponding to the value of p prefixed, is well-defined and converges to a solution x * of the equation F (x) = 0, from every point x 0 ∈ B(x * , R).
Proof.Firstly, fixed p ∈ (0, 1], we define two auxiliary scalar functions: with t ∈ R and p ∈ (0, 1], which are strictly increasing in R + . Secondly, let R < min{r, R}, where R = ζp (1+ζp)βK and ζ p is the positive real root of the equation g(t) − 1 = 0.As x 0 ∈ B(x * , R), with R < r and βK R < βKR < 1, then by Lemma 1 the operator Γ 0 = [F (x 0 )] −1 exists and is such that So, y 0 is well-defined and, from Lemma 2 (i 0 ), we obtain On the other hand, taking into account that g(2) Moreover, z 0 is well-defined and z 0 ∈ B(x * , R), since, from item (ii 0 ) of Lemma 2, it follows as a consequence of f (δ) < 1.Furthermore, x 1 is well-defined and, from item (iii 1 ) of Lemma 2, we obtain that x 1 ∈ B(x * , R), since Now, as g(ζ p ) = 1 and δ < ζ p , it follows that g(δ) < 1 and then Following now an inductive procedure on n ∈ N, we have: Therefore, x n − x * < g(δ) n x 0 − x * , for all n ∈ N, and consequently lim n→+∞ As we have just seen in the previous result, the equation g(t)−1 = 0 determines the value of ζ p which allows us to obtain the radius of the ball of convergence.
Concerning the uniqueness of the solution x * , we have the following result.
. By the Banach lemma, we have that J is invertible, since The proof is complete.

On the Accessibility
A procedure to determine the accessibility of an iterative process to x * is to estimate the ball of convergence of the method from a local convergence result.In this way, from the local convergence result obtained in the previous section, we study the accessibility of family Equation (3).To start the comparative study of accessibility of the different iterative processes of family Equation (3), we can see in Figure1 that the function g(t) − 1 is strictly increasing on R + , tends to infinite and g(0) < 1.Therefore, for each fixed p, there is a positive real root ζ p of the equation g(t) − 1 = 0.In addition, as we can see in Table 1, the value of this root can be increased by increasing the value of p. So, we can conclude that the accessibility of the iterative processes of family Equation (3) increase by increasing the value of p.On the other hand, we want to compare the accessibility of the iterative processes of family Equation (3) with that of the Newton method.For this, we consider the result of convergence obtained by Dennis and Schnabel in [7].They prove that under conditions (C 1 ) and (C 2 ), for any starting point belonging , the Newton method is convergent.
As the radii of the balls of convergence of family Equation (3) and the Newton method are given under conditions (C 1 ) and (C 2 ), it is clear that the radius of the ball of convergence of the Newton method is larger than that of family Equation (3), for p ∈ (0, 0.896037), while for p ∈ (0.896037..., 1], the iterative processes of family Equation (3) have improved accessibility to the Newton method, coinciding both for p = 0.896037... (see Table 1).Therefore, we have obtained, for p ∈ (0.896037..., 1], iterative processes with better accessibility and efficiency than the Newton method.

Application
Next, we illustrate the previous result with the following example given in [7].We choose the max-norm.Let F : R 3 → R 3 be defined as F (x, y, z) = (x, y 2 + y, e z − 1).It is obvious that x * = (0, 0, 0) is a solution of the system.

1 Figure 1 .
Figure 1.Graph of the function g(t)-1 for different values of p.

Table 1 .
Positive real roots of the equation g(t) − 1 = 0 and the values of R for different values of p in Equation (3).