New Fourth-order Schroder-type Methods for Finding Zeros of Nonlinear Equations Having Unknown Multiplicity

In this paper we define two new fourth -order Schroder -type methods for finding zeros of nonlinear equations having unknown multiplicity. In terms of computational cost the new iterative methods requires six evaluations of functions per iteration. It is proved tha t the new methods have a convergence of order four. Numerical comparisons are included to demonstrate exceptional convergence speed of the proposed methods.


Introduction
Solving nonlinear equations is one of the most important problems in numerical analysis. In this paper, we consider iterative methods to find a multiple root α of multiplicity m, i.e., α In recent years, some modifications of the Newton-type methods for multiple roots have been proposed and analysed [3,[6][7][8]. We have found that there are not many methods known to handle the case when multiplicity is unknown, hence we present a new method for finding multiple zeros of a nonlinear equation, based on four evaluations of the function per iteration. In this study we have developed a new iterative method for finding multiple roots of nonlinear equations of a higher order than the existing well-known iterative methods [4][5][6]. Furthermore, in this paper we have shown further development of the second-order methods [4][5][6]. In addition, the new iterative method has a better efficiency index than the classical Schroder method and Noor et al. methods [4]. Hence, the proposed third-order method is significantly better when compared with these established methods.

Basic Definitions
In order to establish the order of convergence of the new third-order method, we use the following definitions.
is the error equation. If the error equation exists, then p is the order of convergence of the iterative method [1,3,9]. Definition 3 Let r be the number of function evaluations of the method. The efficiency of the method is measured by the concept of efficiency index and defined as , r p (3) where p is the order of convergence of the method [1,3,9]. Definition 4 Suppose that 2 1 , and n x are three successive iterations closer to the root α . Then the computational order of convergence may be approximated by where ( ) ( ),

Development of the Method and Convergence Analysis
In this section we define a new third-order iterative method for finding multiple roots of a nonlinear equation. In order to construct the new iterative method we require one function and three derivatives evaluations.
The scheme is given as where ∈  n , 0 x is the initial guess and provided that the denominator of (5) is not equal to zero. We have found that the formula (5) can also be expressed as where ( ) ( ) α Then the iterative method defined by scheme (5) has third-order convergence.

Proof
Let α be a multiple root of multiplicity m of a sufficiently differentiable function ( ) Using the Taylor series expansion of ( ) ( ) Substituting appropriate expressions in (5), Simplifying (18), we obtain the asymptotic error constant The expression (19) establishes the asymptotic error constant for the third-order of convergence for the new Schroder-type method defined by (5).

Formulas for Approximating the Multiplicity
In this section, we state some of the recently presented formulas to approximate the multiplicity m of the method given in [7]. First we state the classical Lagouanelle formula [2] for approximating the multiplicity m given as and the new improved formula for approximating the multiplicity m expressed as In [7] Thukral presented new alternative formulas for approximating the multiplicity m which are expressed as , The approximations of m are obtained by these formulas are based on the estimates of (5) and the errors of the above approximants are given by and the performance of these formulas are displayed in the table 4 and 5.

The Established Methods
For the purpose of comparison, four iterative methods presented in [4] are considered. Since these methods are well established, the essential formulas are used to calculate the approximate solution of the given nonlinear equations and thus compare the effectiveness of the new third-order method.
The first of the methods is the classical Schroder second-order method presented in [5] and is expressed as where 1 2 , t t are given by (8).
The other four second-order methods are by Noor et al. and are presented in [4], and the general form of the method is expressed as Further details of construction of the above formulas are given [4].

Numerical Results
In this section, we shall present the numerical results obtained by employing the iterative methods (5), (29)-(33) to solve some nonlinear equations with unknown multiplicity m. To demonstrate the performance of the new higher order iterative methods, we use ten particular nonlinear equations. The consistency and stability of the results are determined by examining the convergence of the new iterative methods. The findings are generalised by illustrating the effectiveness of the higher order methods for determining the multiple root of a nonlinear equation. Consequently, the errors obtained by each of the methods are displayed in the following tables. In fact, the errors displayed are of absolute value and insignificant approximations by the various methods have been omitted in the following tables. The numerical computations listed in the tables were performed on an algebraic system called Maple.
The new third-order method requires four function evaluations and has the order of convergence three. To determine the efficiency index of the new method, definition 3 will be used. Hence, the efficiency index of the new method given by (5) is 4

1.316 ≈
whereas the efficiency index of the second-order methods (29)-(34) is given by 3 We can see that the efficiency index of the new third-order method has a better efficiency index than the second-order methods. The test functions, multiplicity m and their exact root α are displayed in table 1. The difference between the root α and the approximation n x for test functions with initial guess 0 x are displayed in table 2. Table 2 shows the absolute errors obtained by each of iterative methods described, we see that the new third order method is producing better results than the established methods. Furthermore, the computational order of convergence (COC) are displayed in table 3. From the table 3, we observe that the COC perfectly coincides with the theoretical result.
In addition, the difference between the multiplicity m and the approximation  m based on the present method (5) are displayed in table 4. In table 4 we see a remarkable precision of the multiplicity m and there is no significant difference between the Lagouanelle formula (20) and the recently introduced formulas (22)-(26). Furthermore, in table 5 we display the difference between the multiplicity m and the approximation  m based on the Schroder method (29). As before, we find that there is no significant difference between the formulas, they are all performing equally. As expected, we see a remarkable precision of the multiplicity m based on new formula (5) approximants are significant better than the approximants of the Schroder method (29). In fact, n x is calculated by using the same total number of function evaluations (TNFE) for all methods.

Conclusions
In this paper, a new third-order iterative method for solving nonlinear equations with multiple roots has been introduced. Convergence analysis proves that the new method preserve the order of convergence. The prime motive for presenting the new method was to establish a higher order of convergence method than the existing second order methods [4]. The effectiveness of the new third-order method is examined by showing the accuracy of the multiple roots of several nonlinear equations. After an extensive experimentation, it can be concluded that the convergence of the new third-order method is remarkably fast. The main purpose of demonstrating the new method for different types of nonlinear equations was purely to illustrate the accuracy of the approximate solution, the stability of the convergence, the consistency of the results and to determine the efficiency of the new iterative method. We have shown numerically, and verified, that the new iterative method converge to the order three. Empirically, we have found that the approximants obtained by the new third order method is producing better precision of the multiplicity m than the approximants of classical Schroder second order method. Finally, we have constructed a new higher order iterative method with better efficiency index than the existing second order methods, but unfortunately the new method is not of optimal order, hence further investigation is essential.