About this Journal Submit a Manuscript Table of Contents
The Scientific World Journal
Volume 2013 (2013), Article ID 837243, 9 pages
Research Article

An Iterative Solver in the Presence and Absence of Multiplicity for Nonlinear Equations

1Young Researchers and Elite Club, Zahedan Branch, Islamic Azad University, Zahedan, Iran
2Department of Mathematics, University of Venda, Private Bag X5050, Thohoyandou 0950, South Africa
3Department of Mathematics, Science and Letter Faculty, Kocaeli University, Umuttepe Campus, Kocaeli, Turkey

Received 17 June 2013; Accepted 1 October 2013

Academic Editors: S. Amat, M. Inc, Z. Mukandavire, S. Ostrovska, X. Qin, and S. A. Tersian

Copyright © 2013 Fazlollah Soleymani et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


We develop a high-order fixed point type method to approximate a multiple root. By using three functional evaluations per full cycle, a new class of fourth-order methods for this purpose is suggested and established. The methods from the class require the knowledge of the multiplicity. We also present a method in the absence of multiplicity for nonlinear equations. In order to attest the efficiency of the obtained methods, we employ numerical comparisons alongside obtaining basins of attraction to compare them in the complex plane according to their convergence speed and chaotic behavior.

1. Introduction

There are several methods for computing a zero with multiplicity of a nonlinear equation ; see, for example, [14]. Newton’s iterative method is only of first order unless it is modified to obtain the second order of convergence. This modification requires knowledge of the multiplicity. Traub in [5] has suggested the use of any method for or . Any such method will require higher derivatives than the corresponding one for simple zeros. In general, constructing iterative methods for finding multiple zeros is not an easy task and careful attention is required in their developments.

This work tries to deal with the solution of nonlinear equations when the multiplicity of the roots is available or can be approximated and also when is very high or not available. It then will be clear that for the choice , the corresponding method can easily be applied for simple zeros.

In this paper, we consider first the following class of methods for simple zeros (with , : where , , , and with the following error equation:

We then extend (1) for multiple zeros without lowering the convergence order or applying much more computational evaluations of the function or derivatives. Thus, the proposed scheme will be of order four for multiple zeros with two evaluations of the first-order derivative and one evaluation of the function, which clearly implies the consistency of the general solver with the optimality conjecture of Kung-Traub [6]. Furthermore, we will discuss some qualitative issues using the attraction basins as a criterion for comparison. To this end, we refer the readers to see [79] on obtaining a background on this matter regarding a rational map on the Riemann sphere, Julia set, Newton’s fractals, the definitions of dynamic and chaotic behaviors, and so forth.

The remaining sections of the paper are organized as follows. Section 2 will shortly report some of the existing methods for multiple roots using three functional evaluations per full step. It will be followed by Section 3 wherein the main contribution of the paper lies, as a general fourth-order solver in the presence and the absence of multiplicity for nonlinear equations. Section 4 will be devoted to the comparison of different methods based on basins of attraction to observe the chaotic behaviors of the corresponding fractals produced by each iterative method in the complex plane. Numerical reports will be furnished in Section 5 and the conclusion of the study has been presented in the last section.

2. A Brief Review

Here, we mention some of the famous schemes for multiple zeros using three functional evaluations per full cycle shortly. The old and well-known third-order method of Dong in [10] can be given as Dong has also proposed the following third-order method in [10]:

In 2010, J. R. Sharma and R. Sharma in [11] presented the following quartically convergent method: Recently in 2011, Zhou et al. in [12] presented the following optimal fourth-order scheme:

All the reported methods above require the multiplicity of the zeros and three functional evaluations per computing step to proceed. Clearly, only methods (5) and (6) are optimal in the sense of Kung-Traub for constructing multipoint methods without memory for solving nonlinear equations. Some other new developments for multiple zeros can be found in [13, 14].

3. A New General Solver for Multiple Zeros

Let us consider the following class for multiple zeros as a generalized form of (1) and by using weight function approach: We now show that this class is of fourth order for multiple roots under certain conditions on the weight function.

Theorem 1. Let be a multiple zero of a sufficiently differentiable function for an open interval with the multiplicity , which includes as an initial approximation of . Then, the class of methods without memory (7) is of optimal local order four, when and , where .

Proof. To find the asymptotic error constant of (7) where , we expand any terms of (7) around the multiple zero in the th iterate whence . Clearly, for such symbolic computations, a Computer Algebra System such as [15] needs to be used.
A Taylor expansion around yields and , so that, using algebraic software again, Now we easily have . Furthermore, we have by Taylor expansion Thus by using (10), we obtain And also, By using (9), (11), and (12), class (7) satisfies the error equation:
This shows that our class (7) reaches the local quartically convergence using three evaluations. The proof is complete.

The proposed class (7) uses only three functional evaluations per full cycle to achieve fourth order of convergence based on Theorem 1, which implied as its classical computational efficiency index that is the same as (5) and (6) and is higher than of methods (3) and (4).

A simple optimal and efficient method from class (7) could be presented in what follows: by choosing

Another optimal and efficient method from class (7) with a simplified error equation could be presented in what follows: by choosing This method has the simplified error equation:

An important challenge in the iterative methods based on the known multiplicity is to find the order of multiplicity correctly. However, most algorithms to determine the order of multiplicity may lead to mutually opposite requirements; for example,(i)Traub in [5] showed that when is very close to the multiple root of ;(ii)Lagouanelle in [16] proposed the following approximate formula: when once again is very close to the multiple root of .

The above-mentioned ways demand a very close approximation to calculate a multiplicity of high accuracy. On the other hand, to find a very close approximation to a multiple root, it is necessary to use precise multiplicity. Sometimes, both of the requirements cannot be attained at the same time.

Taking into account the opposite demands mentioned and additional calculations to find multiplicity, in those cases where we cannot provide an accurate multiplicity, it is sometimes better to apply a method which does not explicitly require the order of multiplicity, in spite of its lower computational efficiency arising from additional functional evaluations per iteration.

Hence, by applying the optimal class of methods (1) on the transformation , we can easily extend it for dealing with multiple roots in the absence of multiplicity, when high precision alongside high order is needed. For a simple weight function in (1) which is easily constructed by the following Mathematica command:


 DSolve[{G’’’ [t] == 0, G 1 == 1,

G 1 == −1/4, G’’ 1 == 7/4}, G[t], t];

weight[t_] = G[t]/. %[ 1 ]

We obtain for simple zeros and the following alternative for multiple zeros in the absence of multiplicity: wherein

Therefore, now we have an efficient method (22) of order four for finding the multiple roots too. Note that until now, we have distinguished two kinds of methods: those which deal with a known order of multiplicity and others, such as (22), with no information on multiplicity.

4. The Dynamical Behavior of the Methods

We here investigate the comparison of iterative schemes in the complex plane using basins of attraction. It is well known that a fixed point is a point of a function that does not change under some transformation. We further recall that if we regard the evolution of a dynamical system as a series of transformations, then there may or may not be a point which remains fixed under each transformation. The final state, that a dynamical system evolves towards, corresponds to an attracting fixed point of the evolution function, but the two concepts are not equivalent because not all fixed points attract the evolution of nearby points.

The aim herein is to use the basin of attraction as another way for comparing the iteration algorithms; see, for example, [17].

In this section, we take a rectangle and we assign a color to each point according to the multiple root at which the corresponding iterative method starting from converges, and we mark the point as black if the method does not converge. In this way, we distinguish the attraction basins by their colors for different methods.

The criteria, we have used in our Mathematica codes, are that the maximum number of iterations is 100. That is to say, if the method does not reach to the considered accuracy after 100 of its full iteration steps, we allocate the black color and also consider it as Not Convergent. The considered accuracy is the obtained residual of function to be less than .

We have used methods (3), (5), (6), and (14) for some complex functions having multiple roots with known multiplicity.

For the first test, we have taken the following function.

Test Problem 1. We have
Equation (24) has five roots , , , , and of the multiplicity five. Based on Figures 1 and 2 we can see that method (14), method (5), and method (3) are the good ones. Method (14) has larger basins in the areas corresponding to each root.

Figure 1: The methods (3) (a) and (5) (b) for Test Problem 1.
Figure 2: The methods (6) (a) and (14) (b) for Test Problem 1.

The second test problem is a polynomial as follows. Despite a few diverging points on the Julia set, we find that the basins of attraction for each root are largest for method (14).

Test Problem 2. We have
The zeros for this test are , , , and of the multiplicity four. The results are presented in Figures 3 and 4. The performance of the methods (5) and (14) are acceptable. Methods (3), (5), (6), and (14) are quite the same though (14) has larger basins again.

Figure 3: The methods (3) (a) and (5) (b) for Test Problem 2.
Figure 4: The methods (6) (a) and (14) (b) for Test Problem 2.

The third test problem is chosen as follows.

Test Problem 3. We have
The roots are , , , , and , , with the known multiplicity 2. The results are presented in Figures 5 and 6. The performance of the methods (5) and (14) are acceptable. The worse methods in this case are (3) and (6).

Figure 5: The methods (3) (a) and (5) (b) for Test Problem 3.
Figure 6: The methods (6) (a) and (14) (b) for Test Problem 3.

5. Numerical Reports

In this section, we exhibit numerical results showing the behavior of the methods in this paper with roots of known multiplicity . Iterative methods (3), (4), (5), (6), and (14) for multiple roots, which also require three functional evaluations, have been chosen. The list of test nonlinear functions including multiple zeros with their multiplicity is presented in Table 1.

Table 1: The test functions considered in this paper.

The results are summarized in Tables 2 and 3 after three full iterations for two different initial guesses. As they show, the novel scheme is comparable with all of the methods. All numerical instances were performed using 500 digits floating point arithmetic. We have computed the root of each test function for the initial guess .

Table 2: Results of comparisons for different methods after three full iterations.
Table 3: Results of comparisons for different methods after three full iterations.

The tests show that the considered optimal methods generate results of approximately the same accuracy. The formulas in this section have been employed on several other nonlinear equations. The experience shows that there is no clear winner among the optimal methods (5), (6), and (14) in the sense that in different situations different methods may be the winners.

Following the results here and the basins of attraction and also considering the optimal efficiency index of method (14), we could conclude that class (7) is efficient and can be considered as a good tool to find multiple zeros iteratively.

Note that the application of such iterative methods for finding generalized inverses could be considered for future researches in this trend of study; see, for example, [18, 19].

6. Concluding Remarks

We have demonstrated the performance of a new class of fourth-order methods. Convergence analysis has proved that the new methods obtained from the class preserve their order of convergence. We do not have to evaluate the second-order derivative of the functions in the presented methods when the multiplicity is known. We have also presented a quartically convergent method in the absence of multiplicity for nonlinear equations.

Due to the fact that the basin of attraction is a method to visually comprehend how an algorithm behaves as a function of the various starting points, thus the dynamical behavior of the proposed methods alongside the comparison with the existing methods has been given in Section 4 to clearly show the efficiency of the methods.

Conflict of Interests

The authors declare that they do not have any conflict of interests.


  1. E. Halley, “A new, exact and easy method of finding the roots of equations generally and that without any previous reduction,” Philosophical Transactions of the Royal Society, vol. 18, pp. 136–148, 1694. View at Publisher · View at Google Scholar
  2. E. Hansen and M. Patrick, “A family of root finding methods,” Numerische Mathematik, vol. 27, no. 3, pp. 257–269, 1976. View at Publisher · View at Google Scholar · View at Scopus
  3. B. Neta and A. N. Johnson, “High-order nonlinear solver for multiple roots,” Computers and Mathematics with Applications, vol. 55, no. 9, pp. 2012–2017, 2008. View at Publisher · View at Google Scholar · View at Scopus
  4. H. D. Victory Jr. and B. Neta, “A higher order method for multiple zeros of nonlinear functions,” International Journal of Computer Mathematics, vol. 12, no. 3-4, pp. 329–335, 1983. View at Scopus
  5. J. F. Traub, Iterative Methods for the Solution of Equations, Prentice Hall, New York, NY, USA, 1964.
  6. H. T. Kung and J. F. Traub, “Optimal order of one-point and multipoint iteration,” Journal of the ACM, vol. 21, no. 4, pp. 643–651, 1974. View at Scopus
  7. C. Getz and J. Helmstedt, Graphics with Mathematica: Fractals, Julia Sets, Patterns and Natural Forms, Elsevier, 2004.
  8. M. McClure, “Newton's method for complex polynomials,” A preprint version of a, “Mathematical graphics” column from Mathematica in Education and Research, pp. 1–15.
  9. E. W. Weisstein, “Newton's Method,” 2006, http://mathworld.wolfram.com/NewtonsMethod.html.
  10. C. Dong, “A family of multipoint iterative functions for finding multiple roots of equations,” International Journal of Computer Mathematics, vol. 21, no. 3-4, pp. 363–367, 1987. View at Publisher · View at Google Scholar
  11. J. R. Sharma and R. Sharma, “Modified Jarratt method for computing multiple roots,” Applied Mathematics and Computation, vol. 217, no. 2, pp. 878–881, 2010. View at Publisher · View at Google Scholar · View at Scopus
  12. X. Zhou, X. Chen, and Y. Song, “Constructing higher-order methods for obtaining the multiple roots of nonlinear equations,” Journal of Computational and Applied Mathematics, vol. 235, no. 14, pp. 4199–4206, 2011. View at Publisher · View at Google Scholar · View at Scopus
  13. M. Sharifi, D. K. R. Babajee, and F. Soleymani, “Finding the solution of nonlinear equations by a class of optimal methods,” Computers and Mathematics with Applications, vol. 63, no. 4, pp. 764–774, 2012. View at Publisher · View at Google Scholar · View at Scopus
  14. J. R. Sharma and R. Sharma, “Modified Chebyshev-Halley type method and its variants for computing multiple roots,” Numerical Algorithms, vol. 61, pp. 567–578, 2012. View at Publisher · View at Google Scholar · View at Scopus
  15. S. Wolfram, The Mathematica Book, Wolfram Media, 5th edition, 2003.
  16. J. L. Lagouanelle, “Sur une mtode de calcul de l'ordre de multiplicit des zros d'un polynme,” Comptes Rendus de l'Académie des Sciences A, vol. 262, pp. 626–627, 1966.
  17. F. Soleimani, F. Soleymani, and S. Shateyi, “Some iterative methods free from derivatives and their basins of attraction,” Discrete Dynamics in Nature and Society, vol. 2013, Article ID 301718, 10 pages, 2013. View at Publisher · View at Google Scholar
  18. F. Soleymani and P. S. Stanimirović, “A higher order iterative method for computing the Drazin inverse,” The Scientific World Journal, vol. 2013, Article ID 708647, 11 pages, 2013. View at Publisher · View at Google Scholar
  19. F. Toutounian and F. Soleymani, “An iterative method for computing the approximate inverse of a square matrix and the Moore-Penrose inverse of a non-square matrix,” Applied Mathematics and Computation, vol. 224, pp. 671–680, 2013. View at Publisher · View at Google Scholar