Efficient iterative methods for finding simultaneously all the multiple roots of polynomial equation

Two new iterative methods for the simultaneous determination of all multiple as well as distinct roots of nonlinear polynomial equation are established, using two suitable corrections to achieve a very high computational efficiency as compared to the existing methods in the literature. Convergence analysis shows that the orders of convergence of the newly constructed simultaneous methods are 10 and 12. At the end, numerical test examples are given to check the efficiency and numerical performance of these simultaneous methods.


Introduction
A wide range of theoretical and practical problems arise in various fields of mathematical, economical, physical, and engineering sciences which can be formulated as a polynomial equation of degree n with arbitrary real or complex coefficient: f (x) = x n + a n-1 x n-1 + · · · + a 0 = n j=1 (xζ j ) = (xζ i ) where ζ 1 · · · ζ n denote all the simple or complex roots of (1). Approximating all roots of the nonlinear polynomial equation using simultaneous methods has a lot of applications in sciences and engineering because simultaneous iterative methods are less time consuming since they can be implemented for parallel processing as well. Further details about their convergence properties, computational efficiency, and parallel processing may be found in  and the references cited there in. The main objective of this paper is to develop simultaneous methods which have a higher convergence order and are more efficient as compared to the existing methods. A very high computational efficiency is achieved by using two suitable corrections [26,27] with convergence orders equal to ten and twelve with a minimal number of function evaluations in each step.

Construction of simultaneous methods for multiple roots
Consider two-step fourth-order Newton's method [26] for finding multiple roots of nonlinear equation (1) where σ is the multiplicity of exact root, say ζ , of (1). We would like to convert (2) into a simultaneous method for extracting all the distinct as well as multiple roots of (1). We use the third-order Dong et al. method [26] as a correction to increase the efficiency and convergence order requiring no additional evaluation of the function: Suppose that the nonlinear polynomial equation (1) has n roots. Then This implies This gives The multiple root equation (5) can be written as Replacing x j by x * j in (6), we have where x * j = u j using (3) .
Using (7) in the first step of (2), we have , k = 0, 1, . . . , Thus we have constructed a new simultaneous method (8) abbreviated as MNS10M for extracting all distinct as well as multiple roots of polynomial equation (1).

Convergence analysis
In this section, the convergence analysis of a family of two-step simultaneous methods (8) given in a form of the following theorem is presented.
Theorem 1 Let ζ 1 , . . . , ζ n be simple roots of (1). If x (0) 1 , . . . , x (0) n are the initial approximations of the roots respectively and sufficiently close to the actual roots, then the order of convergence of method (8) equals ten.
Proof Let i = x iζ i , i = y iζ i , and i = z iζ i be the errors in x i , y i , and z i approximations respectively. Consider the first step of (8), which is Then, obviously, for distinct roots, we have Thus, for multiple roots, we have from (8) If it is assumed that absolute values of all errors j (j = 1, 2, 3, . . .) are of the same order as, say, | j | = O| |, then from (9) we have From the second equation of (8), we get which shows that the convergence order of method (8) is ten. Hence we have proved the theorem.

Improvement of efficiency and convergence order
To improve the convergence order of method (8) from 10 to 12, using same function evaluation, we use (7), i.e., where Z * j is a fourth-order method [27]. Using (11) in the first step of (2), we have Thus we have constructed a new simultaneous method (12), abbreviated as MNS12M for extracting all multiple roots of polynomial equation (1). For multiplicity unity, we used method (12) for determing all the distinct roots of (1), abbreviated as MNS12D.

Convergence analysis
In this section, the convergence analysis of a family of two-step simultaneous methods (12) is given in a form of the following theorem.
n are the initial approximations of the roots respectively and sufficiently close to the actual roots, then the order of convergence of method (12) equals twelve.
Proof Let i = x iζ i , i = y iζ i , and i = z iζ i be the errors in x i , y i , and z i approximations respectively. Consider the first step of (12), which is Then, obviously, for distinct roots, we have Thus, for multiple roots, we have from (6) , If it is assumed that absolute values of all errors j (j = 1, 2, 3, . . .) are of the same order as, say, | j | = O| |, then from (13) we have From the second equation of (12), we have If it is assumed that absolute values of all errors j (j = 1, 2, 3, . . .) are of the same order as, say, | j | = O| |, then we have which shows that the convergence order of method (12) is twelve. Hence we have proved the theorem.

Computational analysis
Here we compare the computational efficiency and convergence behavior of the Petkovic et al. [28] method (abbreviated as PJM10D) and the new simultaneous iterative methods (8) and (12). As presented in [28], the efficiency of an iterative method can be estimated using the efficiency index given by where D is the computational cost and r is the order of convergence of the iterative method.
The number of addition and subtraction, multiplications, and divisions per iteration for all n roots of a given polynomial of degree m is denoted by AS m , M m , and D m . The computational cost can be approximated as and thus (15) becomes Applying (17) and by data given in Table 1, we calculate the percentage ratio ρ( (8), (X)) and ρ( (12), (X)) [28] given by ρ (12), (X) = EF (12) EF(X) where X is the Petkovic method PJM10D. These ratios are graphically displayed in Fig. 1(a), (b), (c). It is evident from Fig. 1(a), (b), (c) that the new methods (8) and (12) are more efficient as compared to the Petkovic method PJM10D. We also calculate the CPU execution time, as all the calculations are done using Maple 18 on (Processor Intel(R) Core(TM) i3-3110m CPU@2.4 GHz with 64-bit operating system. We observe that CPU times of the methods MMS10M and MNS12M are less than those of PJM10D, showing the dominant efficiency of our methods (8) and (12) as compared to them.

Numerical results
Here some numerical examples are considered in order to demonstrate the performance of our family of two-step tenth-order simultaneous methods, namely MNS10M (8) and MNS12M (12). We compare our family of methods with the Petkovic et al. [28] method of convergence of order ten for finding all distinct roots of (1) (abbreviated as PJM10D). All the computations are performed using Maple 15 with 64 digits floating point arithmetic. We take ∈= 10 -30 as a tolerance and use the following stopping criteria for estimating the roots: where e i represents the absolute error of function values in (i) Numerical test examples from [10,28,29] are provided in Tables 2, 3, and 4. In all tables, CO represents the convergence order, n represents the number of iterations, and CPU represents execution time in seconds. All calculations are done using Maple 15 on (Processor Intel(R) Core(TM) i3-3110m CPU@2.4 GHz with 4 GB (3.89 GB USABLE)) with
The initial approximations have been taken as The initial approximations have been taken as (0) x 2 = 0.9, x 3 = 1.8, x 4 = 2.9,

Results and discussion
From Tables 2-4 and from Fig. 1(a)-(c), we conclude that • Our methods MNS10D and MNS12D are more efficient as compared to PJM10D in terms of the number of iterations and CPU time. • Our methods MNS10M and MNS12M are applicable for multiple as well as distinct roots, whereas PJM10D is applicable for distinct roots only.

Conclusion
We have developed here two simultaneous two-step methods of order ten and twelve, namely MNS10D, MNS10M, MNS12D, and MNS12M for determination of all the distinct as well as multiple roots of nonlinear polynomial equation (1). From Tables 1-4, we observed that our methods are very effective and more efficient as compared to the existing method PJM10D [28].