Fifth-Order Iterative Method for Solving Multiple Roots of the Highest Multiplicity of Nonlinear Equation

A three-step iterative method with fifth-order convergence as a new modification of Newton’s method was presented. This method is for finding multiple roots of nonlinear equation with unknown multiplicity m whose multiplicity m is the highest multiplicity. Its order of convergence is analyzed and proved. Results for some numerical examples show the efficiency of the new method.


Introduction
This paper addresses the problem of multiple roots x * of nonlinear equation f (x) = 0 with unknown multiplicity m whose multiplicity m is the highest multiplicity, where f : [a, b] ⊂ R → R is a nonlinear differential function on [a, b].In case the multiplicity m is given explicitly, there are many iterative methods established via various techniques (see [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15] for more details).If the multiplicity m is not known explicitly, Traub [16] utilized a simple transformation F (x) = f (x)/f (x) instead of f (x) for computing a multiple root of f (x) = 0.In this case, the aim of solving a multiple root is reduced to that of solving a simple root of the transformed equation f (x) = 0, and thus any iterative method can be used to preserve the original convergence order.However, Newton's method for this transformed equation requires evaluations of the derivatives f (x) and f (x).In order to avoid the evaluations of these derivatives with the multiplicity m unknown, for multiple roots, King [17] proposed the secant method which does not use the function F = f /f , but rather use . Wu and Fu [18] further used and transformed the problem of solving multiple roots of f (x) = 0 into that of solving simple root of f (x) = 0. Actually, they established the following iteration formulae: where p ∈ R, |p| < ∞.So, the sequence {x n } produced by the iteration Formulae (1) is at least quadratically convergent for multiple roots.Moreover, Wu et al. [19] defined a function where m is the multiplicity, and employed the modified Steffensen's method (see [20,21]) to compute the approximate solution of the equation f (x) = 0, where h n (> 0) is the step size of iteration and |t| < ∞.Parida and Gupta [22] suggested another transformation where δ = sign(f , and transform the task of solving multiple zeros of f into that of solving simple zero of F .In this case, they utilized a quadratically convergent derivative free Newton-like iterative method: where the parameter p should be chosen such that the denominator is the largest in magnitude.Yun [23] suggested a new transformation of f (x) as took ε such that max and proposed an iterative method as follows: Recently, for the transformed equation K(x) = 0 with a simple root, Yun [24] proposed a Steffensen-type iterative formula where In this paper we construct a new modified Newton's method.We will present the proof that the method is three-step iterative method with fifth-order convergence for nonlinear equations of multiple roots with unknown multiplicity m, whose multiplicity m is the highest multiplicity and without requiring the use of the second derivative.

Iterative Method with Fifth-Order Convergence for Solving Multiple Roots
We consider the simple transformation (see [16,25]): and use a Newton-like iterative method: In order to avoid computing the first derivatives of function F (x n ), F (y n ) and F (z n ), we approximate them as follows: (See [25][26][27][28] for the detail discussions of Equations ( 11)-( 13) respectively.)Substituting the approximations of F (x n ), F (y n ) and F (z n ) given by Equations ( 11)- (13) in Equation (10), we establish the following new iterative method: We give the following convergence theorem for the proposed method Equation ( 14) as follows.
Theorem 1. Suppose that f ∈ C 1 (D)(D ⊆ R → R) has multiple roots x * ∈ D, whose multiplicity m is the highest multiplicity of nonlinear equation, where D is an open interval.If the initial point x 0 is sufficiently close to x * , the iterative method defined by ( 14) has fifth-order convergence.
Proof.Without loss of generality, we assume that f (x) has two multiple roots where x * is a multiple root of Equation ( 15) with multiplicity m and x 1 is a multiple root of Equation ( 15) with multiplicity n (m > n), h(x) is a continuous function with h(x * ) = 0 and h(x 1 ) = 0.According to Equation ( 15), we have Dividing Equation ( 15) by Equation ( 16), we get From Equation ( 17), we can see that the problem of computing multiple roots of f (x) = 0 can be reduced to the equivalent problem of computing simple root x * of F (x) = 0.
The proof is completed.We further consider how to find the highest multiplicity of the root x * in the iterative method.If x n is the n-th iteration computed by an iterative method applied to f , then from Equation (9), we have where Consequently, when the iteration becomes closer to the root x * , we can estimate its multiplicity by computing In the practical computing root x * process, some iteration number is no more than two by using Equation (14).According to this case, we compute the root x * by using Equation ( 14), then we select the initial value near the root x * , and we use Newton iterative method to compute the highest multiplicity.Therefore, m is approximately the reciprocal of the divided difference of f for successive iteration x n and x n+1 .

Numerical Results
We employ the proposed modification of Newton's method with three-step Equation ( 14) (MNM) to solve some nonlinear equations.All the computations were done by using Visual C++ 6.0 and were satisfied the condition such that |f In order to show the effectiveness of our iterative method, we provide at least three different initial iterative values.From different initial iterative values, they can be convergent to the same iterative solution whose multiplicity is the highest multiplicity of nonlinear equation.We used the following test functions and obtained the approximate zeros x * round up to the 17-th decimal place: , x 0 = −12.0,23.0, 1.9, m = 7, n = 4, x * = 2.236067977499790 ), x 0 = 11, −11, 1.8, 0.5, m = 5, n = 2, x * = 2.6457513110645906 The computational results indicate that our proposed iterative method can converge to multiple roots whose multiplicity is the highest multiplicity of nonlinear equation.In the next section, for the special case which the multiplicity of the roots of nonlinear equation is a single multiplicity, we present the analysis result for comparison with previous methods.

Comparison with Previous Methods
In this section, we use the proposed modification of Newton's method with three-step (14) (MNM) (the Formula (14) in our paper) to solve some nonlinear equations which the multiplicity of the roots is a single multiplicity, and compare them with King's method [17] (KM, (4), ( 13) )), the high-order convergence iteration methods without employing derivatives given in [18] (WFM, ( 6) )), the improved method for finding multiple roots and it ąŕs multiplicity of nonlinear equations given in [22] (PGM, (( 6), (11)) , where δ = sign(f (11)), the derivative free iterative method for finding multiple roots of nonlinear equations given in [23] (YM, ((7), (10)) Hε(2x k −x k−1 )−Hε(x k−1 ) (10)), transformation methods for finding multiple roots of nonlinear equations [24] (YM, ((10), (12)) , where f ε (x) ≈ εf (x)f (x) (10), ), as well as Newton's first order method(NM) All the computations were done by using Visual C++ 6.0.Since King's method and Yun's method require two starting values, we have used x 1 = x 0 − 0.1.We used the following test functions and obtained the approximate zeros x * round up to the 17th decimal place: Note that we used NC in Table 1 to mean that the method does not converge to the root.And these methods can converge to root by using closer initial values.The computational results in Table 1 demonstrate that our proposed iterative method (MNM) requires less number of iterations than those of KM, WFM, PGM, YM, and NM.Therefore, it is significant and applicable and can compete with other existing methods.
Table 1.Comparison of various aiterative methods for the case of single multiplicity.

Efficiency of Iterative Methods
In the following we compare the efficiency of methods mentioned in Section 1 whose the multiplicity of the roots of nonlinear equation is a single multiplicity.We consider the definition of efficiency index [29][30][31] as EFF = r 1 θ , where r is the order of the method and θ is number of function (and derivatives) evaluations per iteration required by the method.These results are presented in Table 2.In Table 2, we listed the methods according to decreasing order of efficiency index, and multiplicity m of roots with all methods under the row of King's method (including King's method) are unknown.For unknown multiplicity m, our new method comes second, next to King's method.

Table 2 .
Comparison the efficiency of various iterative methods for the case of single multiplicity.