Next Article in Journal
Similarity Measures of q-Rung Orthopair Fuzzy Sets Based on Cosine Function and Their Applications
Next Article in Special Issue
On the Semilocal Convergence of the Multi–Point Variant of Jarratt Method: Unbounded Third Derivative Case
Previous Article in Journal
Systems of Variational Inequalities with Nonlinear Operators
Previous Article in Special Issue
Optimal Fourth, Eighth and Sixteenth Order Methods by Using Divided Difference Techniques and Their Basins of Attraction and Its Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Higher Order Chebyshev-Halley-Type Family of Iterative Methods for Multiple Roots

1
Department of Mathematics, King Abdualziz University, Jeddah 21589, Saudi Arabia
2
Instituto Universitario de Matemática Multidisciplinar, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
3
Fac. de Ciencias Económicas, Universidad Laica “Eloy Alfaro de Manabí”, Manabí 130214, Ecuador
4
Departamento de Matemática Aplicada, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(4), 339; https://doi.org/10.3390/math7040339
Submission received: 15 January 2019 / Revised: 21 March 2019 / Accepted: 22 March 2019 / Published: 9 April 2019
(This article belongs to the Special Issue Iterative Methods for Solving Nonlinear Equations and Systems)

Abstract

:
The aim of this paper is to introduce new high order iterative methods for multiple roots of the nonlinear scalar equation; this is a demanding task in the area of computational mathematics and numerical analysis. Specifically, we present a new Chebyshev–Halley-type iteration function having at least sixth-order convergence and eighth-order convergence for a particular value in the case of multiple roots. With regard to computational cost, each member of our scheme needs four functional evaluations each step. Therefore, the maximum efficiency index of our scheme is 1.6818 for α = 2 , which corresponds to an optimal method in the sense of Kung and Traub’s conjecture. We obtain the theoretical convergence order by using Taylor developments. Finally, we consider some real-life situations for establishing some numerical experiments to corroborate the theoretical results.

1. Introduction

One important field in the area of computational methods and numerical analysis is to find approximations to the solutions of nonlinear equations of the form:
f ( x ) = 0 ,
where f : D C C is the analytic function in the enclosed region D , enclosing the required solution. It is almost impossible to obtain the exact solution in an analytic way for such problems. Therefore, we concentrate on obtaining approximations of the solution up to any specific degree of accuracy by means of an iterative procedure, of course doing it also with the maximum efficiency. In [1], Kung and Traub conjectured that a method without memory that uses n + 1 functional evaluations per iteration can have at most convergence order p = 2 n . If this bound is reached, the method is said to be optimal.
For solving nonlinear Equation (1) by means of iterations, we have the well-known cubically-convergent family of Chebyshev–Halley methods [2], which is given by:
x n + 1 = x n 1 + 1 2 L f ( x n ) 1 α L f ( x n ) f ( x n ) f ( x n ) , α R ,
where L f ( x n ) = f ( x n ) f ( x n ) { f ( x n ) } 2 . A great variety of iterative methods can be reported in particular cases. For example, the classical Chebyshev’s method [1,3], Halley’s method [1,3], and the super-Halley method [1,3] can be obtained if α = 0 , α = 1 2 , and α = 1 , respectively. Despite the third-order convergence, the scheme (2) is considered less practical from a computational point of view because of the computation of the second-order derivative.
For this reason, several variants of Chebyshev–Halley’s methods free from the second-order derivative have been presented in [4,5,6,7]. It has been shown that these methods are comparable to the classical third-order methods of the Chebyshev–Halley-type in their performance and can also compete with Newton’s method. One family of these methods is given as follows:
y n = x n f ( x n ) f ( x n ) , x n + 1 = x n 1 + f ( y n ) f ( x n ) α f ( y n ) f ( x n ) f ( x n ) , α R
We can easily obtain some well-known third-order methods proposed by Potra and Pták [4] and Sharma [5] (the Newton-secant method (NSM)) for α = 0 and α = 1 . In addition, we have Ostrowski’s method [8] having optimal fourth-order convergence, which is also a special case for α = 2 . This family is important and interesting not only because of not using a second- or higher order derivative. However, this scheme also converges at least cubically and has better results in comparison to the existing ones. Moreover, we have several higher order modifications of the Chebyshev–Halley methods available in the literature, and some of them can be seen in [9,10,11,12].
In this study, we focus on the case of the multiple roots of nonlinear equations. We have some fourth-order optimal and non-optimal modifications or improvements of Newton’s iteration function for multiple roots in the research articles [13,14,15,16,17]. Furthermore, we can find some higher order methods for this case, but some of them do not reach maximum efficiency [18,19,20,21,22,23]; so, this topic is of interest in the current literature.
We propose a new Chebyshev–Halley-type iteration function for multiple roots, which reaches a high order of convergence. Specifically, we get a family of iterative methods with a free parameter α , with sixth-order convergence. Therefore, the efficiency index is 6 1 / 4 , and for α = 2 , this index is 8 1 / 4 , which is the maximum value that one can get with four functional evaluations, reaching optimality in the sense of Kung and Traub’s conjecture. Additionally, an extensive analysis of the convergence order is presented in the main theorem.
We recall that ξ C is a multiple root of the equation f ( x ) = 0 , if it is verified that:
f ( ξ ) = 0 , f ( ξ ) = 0 , , f ( m 1 ) ( ξ ) = 0 a n d f ( m ) ( ξ ) 0 ,
the positive integer ( m 1 ) being the multiplicity of the root.
We deal with iterative methods in which the multiplicity must be known in advance, because this value, m, is used in the iterative expression. However, we point out that these methods also work when one uses an estimation of the multiplicity, as was proposed in the classical study carried out in [24].
Finally, we consider some real-life situations that start from some given conditions to investigate and some standard academic test problems for numerical experiments. Our iteration functions here are found to be more comparable and effective than the existing methods for multiple roots in terms of residual errors and errors among two consecutive iterations, and also, we obtain a more stable computational order of convergence. That is, the proposed methods are competitive.

2. Construction of the Higher Order Scheme

In this section, we present the new Chebyshev–Halley-type methods for multiple roots of nonlinear equations, for the first time. In order to construct the new scheme, we consider the following scheme:
y n = x n m f ( x n ) f ( x n ) , z n = x n m 1 + η 1 α η f ( x n ) f ( x n ) , x n + 1 = z n H ( η , τ ) f ( x n ) f ( x n ) ,
where the function:
H ( η , τ ) = η τ β ( α 2 ) 2 η 2 ( η + 1 ) + τ 3 + τ 2 ( η + 1 ) ( τ + 1 )
with:
η = f ( y n ) f ( x n ) 1 m , τ = f ( z n ) f ( y n ) 1 m , β = m ( α ( α + 2 ) + 9 ) η 3 + η 2 ( α ( α + 3 ) 6 τ 3 ) + η ( α + 8 τ + 1 ) + 2 τ + 1 ,
where α R is a free disposable variable. For m = 1 , we can easily obtain the scheme (3) from the first two steps of the scheme (4).
In Theorem 1, we illustrate that the constructed scheme attains at least sixth-order convergence and for α = 2 , it goes to eighth-order without using any extra functional evaluation. It is interesting to observe that H ( η , τ ) plays a significant role in the construction of the presented scheme (for details, please see Theorem 1).
Theorem 1.
Let us consider x = ξ to be a multiple zero with multiplicity m 1 of an analytic function f : C C in the region containing the multiple zero ξ of f ( x ) . Then, the present scheme (4) attains at least sixth-order convergence for each α, but for a particular value of α = 2 , it reaches the optimal eighth-order convergence.
Proof. 
We expand the functions f ( x n ) and f ( x n ) about x = ξ with the help of a Taylor’s series expansion, which leads us to:
f ( x n ) = f ( m ) ( ξ ) m ! e n m 1 + c 1 e n + c 2 e n 2 + c 3 e n 3 + c 4 e n 4 + c 5 e n 5 + c 6 e n 6 + c 7 e n 7 + c 8 e n 8 + O ( e n 9 ) ,
and:
f ( x n ) = f m ( ξ ) m ! e n m 1 ( m + ( m + 1 ) c 1 e n + ( m + 2 ) c 2 e n 2 + ( m + 3 ) c 3 e n 3 + ( m + 4 ) c 4 e n 4 + ( m + 5 ) c 5 e n 5 + ( m + 6 ) c 6 e n 6 + ( m + 7 ) c 7 e n 7 + ( m + 8 ) c 8 e n 8 + O ( e n 9 ) ) ,
respectively, where c k = m ! ( m 1 + k ) ! f m 1 + k ( ξ ) f m ( ξ ) , k = 2 , 3 , 4 , 8 and e n = x n ξ is the error in the nth iteration.
Inserting the above expressions (5) and (6) into the first substep of scheme (4) yields:
y n ξ = c 1 m e n 2 + i = 0 5 ϕ i e n i + 3 + O ( e n 9 ) ,
where ϕ i = ϕ i ( m , c 1 , c 2 , , c 8 ) are given in terms of m , c 2 , c 3 , , c 8 , for example ϕ 0 = 1 m 2 2 m c 2 ( m + 1 ) c 1 2 and ϕ 1 = 1 m 3 3 m 2 c 3 + ( m + 1 ) 2 c 1 3 m ( 3 m + 4 ) c 1 c 2 , etc.
Using the Taylor series expansion and the expression (7), we have:
f ( y n ) = f ( m ) ( ξ ) e n 2 m [ c 1 m m m ! + ( 2 m c 2 ( m + 1 ) c 1 2 ) c 1 m m e n m ! c 1 + c 1 m 1 + m 1 2 m ! c 1 3 { ( 3 + 3 m + 3 m 2 + m 3 ) c 1 4 2 m ( 2 + 3 m + 2 m 2 ) c 1 2 c 2 + 4 ( m 1 ) m 2 c 2 2 + 6 m 2 c 1 c 3 } e n 2 + i = 0 5 ϕ ¯ i e n i + 3 + O ( e n 9 ) ] .
We obtain the following expression by using (5) and (8):
η = c 1 e n m + 2 m c 2 ( m + 2 ) c 1 2 m 2 e n 2 + θ 0 e n 3 + θ 1 e n 4 + θ 2 e n 5 + O ( e n 6 ) ,
where θ 0 = ( 2 m 2 + 7 m + 7 ) c 1 3 + 6 m 2 c 3 2 m ( 3 m + 7 ) c 1 c 2 2 m 3 , θ 1 = 1 6 m 4 12 m 2 ( 2 m + 5 ) c 1 c 3 + 12 m 2 ( ( m + 3 ) c 2 2 2 m c 4 ) 6 m ( 4 m 2 + 16 m + 17 ) c 1 2 c 2 + ( 6 m 3 + 29 m 2 + 51 m + 34 ) c 1 4 and θ 2 = 1 24 m 5 12 m 2 ( 10 m 2 + 43 m + 49 ) c 1 2 c 3 24 m 3 ( ( 5 m + 17 ) c 2 c 3 5 m c 5 ) + 12 m 2 ( 10 m 2 + 47 m + 53 ) c 2 2 2 m ( 5 m + 13 ) c 4 c 1 4 m ( 30 m 3 + 163 m 2 + 306 m + 209 ) c 1 3 c 2 + ( 24 m 4 + 146 m 3 + 355 m 2 + 418 m + 209 ) c 1 5 .
With the help of Expressions (5)–(9), we obtain:
z n ξ = ( α 2 ) c 1 2 m 2 e n 3 + i = 0 4 ψ i e n i + 4 + O ( e n 9 ) ,
where ψ i = ψ i ( α , m , c 1 , c 2 , , c 8 ) are given in terms of α , m , c 2 , c 3 , , c 8 with the first two coefficients explicitly written as ψ 0 = 1 2 m 3 2 α 2 10 α + ( 7 4 α ) m + 11 c 1 3 + 2 m ( 4 α 7 ) c 1 c 2 and ψ 1 = 1 6 m 4 6 α 3 + 42 α 2 96 α + ( 29 18 α ) m 2 + 6 ( 3 α 2 14 α + 14 ) m + 67 c 1 4 + 12 m 2 ( 5 3 α ) c 1 c 3 + 12 m 2 ( 3 2 α ) c 2 2 + 12 m 3 α 2 + 14 α + ( 5 α 8 ) m 14 c 1 2 c 2 .
By using the Taylor series expansion and (10), we have:
f ( z n ) = f ( m ) ( ξ ) e n 3 m ( α 2 ) c 1 2 m 2 m m ! + i = 1 5 ψ ¯ i e n i + O ( e n 6 ) .
From Expressions (8) and (11), we further have:
τ = ( α 2 ) c 1 m e n + 2 α 2 + 8 α + ( 2 α 3 ) m 7 c 1 2 + 2 m ( 3 2 α ) c 2 2 m 2 e n 2 + γ 1 e n 3 + γ 2 e n 4 + O ( e n 5 ) ,
where γ 1 = 1 3 m 3 3 α 3 + 18 α 2 30 α + ( 4 3 α ) m 2 + 3 ( 2 α 2 7 α + 5 ) m + 11 c 1 3 + 3 m 2 ( 4 3 α ) c 3 + 3 m ( 4 α 2 + 14 α + 3 α m 4 m 10 ) c 1 c 2 and γ 2 = 1 24 m 4 24 m 2 6 α 2 + 20 α + ( 4 α 5 ) m 14 c 1 c 3 + 12 m 2 ( 8 α 2 + 24 α + 4 α m 5 m 13 ) c 2 2 + 2 m ( 5 4 α ) c 4 12 m 12 α 3 66 α 2 + 100 α + 2 ( 4 α 5 ) m 2 + ( 20 α 2 + 64 α 41 ) m 33 c 1 2 c 2 + 24 α 4 + 192 α 3 492 α 2 + 392 α + 6 ( 4 α 5 ) m 3 + ( 72 α 2 + 232 α 151 ) m 2 + 6 ( 12 α 3 66 α 2 + 100 α 33 ) m + 19 c 1 4 .
By using Expressions (9) and (12), we obtain:
H ( η , τ ) = ( α 2 ) c 1 2 m 2 e n 2 + λ 1 e n 3 + λ 2 e n 4 + O ( e n 5 )
where λ 1 = c 1 2 m 3 c 1 2 2 α 2 + 8 α + ( 4 α 7 ) m 7 + 2 ( 7 4 α ) c 2 m and λ 2 = 1 6 m 3 c 1 4 6 α 3 + 36 α 2 66 α + ( 29 18 α ) m 2 + 3 ( 6 α 2 22 α + 17 ) m + 34 + 12 ( 5 3 α ) c 3 c 1 m 2 + 12 ( 3 2 α ) c 2 2 m 2 + 6 c 2 c 1 2 m 6 α 2 + 22 α + 2 ( 5 α 8 ) m 17 .
Now, we use the expressions (5)–(13) in the last substep of Scheme (4), and we get:
e n + 1 = i = 1 3 L i e n i + 5 + O ( e n 9 ) ,
where L 1 = ( α 2 ) 2 c 1 3 m 6 c 1 2 α 2 α + m 2 α 2 + 4 α 17 m 3 2 c 2 ( m 1 ) m , L 2 = ( α 2 ) c 1 2 12 c 2 c 1 2 m 10 α 3 24 α 2 39 α + ( 16 α 27 ) m 2 ( 10 α 3 + 27 α 2 262 α + 301 ) m + 91 + 12 c 3 c 1 m 2 ( 4 α + ( 4 α 7 ) m + 8 ) + 12 c 2 2 m 2 ( 12 α + 4 ( 3 α 5 ) m + 21 ) + c 1 4 24 α 4 + 168 α 3 156 α 2 662 α + ( 52 α 88 ) m 3 ( 60 α 3 + 162 α 2 1616 α + 1885 ) m 2 + 2 ( 18 α 4 12 α 3 711 α 2 + 2539 α 2089 ) m + 979 and L 3 = c 1 24 m 8 24 c 2 c 3 c 1 m 3 ( 42 α 2 146 α + 125 ) m 6 ( 7 α 2 26 α + 24 ) 24 c 2 3 m 3 24 α 2 + 84 α + ( 24 α 2 80 α + 66 ) m 73 + 12 c 3 c 1 3 m 2 2 ( 15 α 4 63 α 3 5 α 2 + 290 α 296 ) + ( 54 α 2 190 α + 165 ) m 2 + ( 30 α 4 28 α 3 + 968 α 2 2432 α + 1697 ) m + 12 c 1 2 m 2 c 2 2 80 α 4 304 α 3 226 α 2 + 1920 α + 2 ( 81 α 2 277 α + 234 ) m 2 + ( 80 α 4 112 α 3 + 2712 α 2 6410 α + 4209 ) m 1787 4 ( α 2 ) c 4 m ( 3 α + ( 3 α 5 ) m + 6 ) 2 c 2 c 1 4 m 3 ( 96 α 5 804 α 4 + 1504 α 3 + 2676 α 2 10612 α + 8283 ) + 4 ( 177 α 2 611 α + 521 ) m 3 3 ( 220 α 4 + 280 α 3 7556 α 2 + 18400 α 12463 ) m 2 + 4 ( 108 α 5 234 α 4 4302 α 3 + 22902 α 2 38593 α + 20488 ) m + c 1 6 48 α 6 480 α 5 + 996 α 4 + 5472 α 3 29810 α 2 + 50792 α + ( 276 α 2 956 α + 818 ) m 4 + ( 360 α 4 448 α 3 + 12434 α 2 30518 α + 20837 ) m 3 + ( 432 α 5 1236 α 4 16044 α 3 + 92306 α 2 161292 α + 88497 ) m 2 + ( 168 α 6 + 888 α 5 + 5352 α 4 55580 α 3 + 173290 α 2 224554 α + 97939 ) m 29771 .
It is noteworthy that we reached at least sixth-order convergence for all α . In addition, we can easily obtain L 1 = L 2 = 0 by using α = 2 .
Now, by adopting α = 2 in Expression (14), we obtain:
e n + 1 = A 0 12 c 3 c 1 m 3 12 c 2 c 1 2 m ( 3 m 2 + 30 m 1 ) + 12 c 2 2 m 2 ( 2 m 1 ) + c 1 4 ( 10 m 3 + 183 m 2 + 650 m 3 ) 24 m 8 e n 8 + O ( e n 9 ) ,
where A 0 = ( c 1 3 ( m + 1 ) 2 c 1 c 2 m ) . The above Expression (15) demonstrates that our proposed Scheme (4) reaches eighth-order convergence for α = 2 by using only four functional evaluations per full iteration. Hence, it is an optimal scheme for a particular value of α = 2 according to the Kung–Traub conjecture, completing the proof. □

3. Numerical Experiments

In this section, we illustrate the efficiency and convergence behavior of our iteration functions for particular values α = 0 , α = 1 , α = 1.9 , and α = 2 in Expression (4), called O M 1 , O M 2 , O M 3 , and O M 4 , respectively. In this regards, we choose five real problems having multiple and simple zeros. The details are outlined in the examples (1)–(3).
For better comparison of our iterative methods, we consider several existing methods of order six and the optimal order eight. Firstly, we compare our methods with the two-point family of sixth-order methods proposed by Geum et al. in [18], and out of them, we pick Case 4c, which is mentioned as follows:
y n = x n m f ( x n ) f ( x n ) , m > 1 , x n + 1 = y n m + a 1 u n 1 + b 1 u n + b 2 u n 2 × 1 1 + c 1 s n f ( y n ) f ( y n ) ,
where:
a 1 = 2 m 4 m 4 16 m 3 + 31 m 2 30 m + 13 ( m 1 ) 4 m 2 8 m + 7 , b 1 = 4 2 m 2 4 m + 3 ( m 1 ) 4 m 2 8 m + 7 , b 2 = 4 m 2 8 m + 3 4 m 2 8 m + 7 , c 1 = 2 ( m 1 ) , u n = f ( y n ) f ( x n ) 1 m , s n = f ( y n ) f ( x n ) 1 m 1 ,
called G M 1 .
In addition, we also compare them with one more non-optimal family of sixth-order iteration functions given by the same authors of [19], and out of them, we choose Case 5YD, which is given by:
y n = x n m f ( x n ) f ( x n ) , m 1 , w n = x n m u n 2 2 u n 1 u n 1 5 u n 2 f ( x n ) f ( x n ) , x n + 1 = x n m u n 2 2 u n 1 5 u n 2 u n + v n 1 f ( x n ) f ( x n ) ,
where u n = f ( y n ) f ( x n ) 1 m and v n = f ( w n ) f ( x n ) 1 m , and this method is denoted as G M 2 .
Moreover, we compare our methods with the optimal eighth-order iterative methods proposed by Zafar et al. [21]. We choose the following two schemes out of them:
y n = x n m f ( x n ) f ( x n ) , w n = y n m u n 6 u n 3 u n 2 + 2 u n + 1 f ( x n ) f ( x n ) , x n + 1 = w n m u n v n ( 1 + 2 u n ) ( 1 + v n ) 2 w n + 1 A 2 P 0 f ( x n ) f ( x n )
and:
y n = x n m f ( x n ) f ( x n ) , w n = y n m u n 1 5 u n 2 + 8 u n 3 1 2 u n f ( x n ) f ( x n ) , x n + 1 = w n m u n v n ( 1 + 2 u n ) ( 1 + v n ) 3 w n + 1 A 2 P 0 ( 1 + w n ) f ( x n ) f ( x n ) ,
where u n = f ( y n ) f ( x n ) 1 m , v n = f ( w n ) f ( y n ) 1 m , w n = f ( w n ) f ( x n ) 1 m , and these iterative methods are denoted in our tables as Z M 1 and Z M 2 , respectively.
Finally, we demonstrate their comparison with another optimal eighth-order iteration function given by Behl et al. [22]. However, we consider the following the best schemes (which was claimed by them):
y n = x n m f ( x n ) f ( x n ) , z n = y n m f ( x n ) f ( x n ) h n ( 1 + 2 h n ) , x n + 1 = z n + m f ( x n ) f ( x n ) t n h n 1 t n 1 2 h n h n 2 + 4 h n 3 2 k n
and:
y n = x n m f ( x n ) f ( x n ) , z n = y n m f ( x n ) f ( x n ) h n ( 1 + 2 h n ) , x n + 1 = z n m f ( x n ) f ( x n ) t n h n 1 t n 1 + 9 h n 2 + 2 k n + h n ( 6 + 8 k n ) 1 + 4 h n ,
with h n = f ( y n ) f ( x n ) 1 m , k n = f ( z n ) f ( x n ) 1 m t n = f ( z n ) f ( y n ) 1 m , which are denoted B M 1 and B M 2 , respectively.
In order to compare these schemes, we perform a numerical experience, and in Table 1 and Table 2, we display the difference between two consecutive iterations | x n + 1 x n | , the residual error in the corresponding function | f ( x n ) | , and the computational order of convergence ( ρ ) (we used the formula given by Cordero and Torregrosa [25]:
ρ ln ( x k + 1 x k / x k x k 1 ) ln ( x k x k 1 / x k 1 x k 2 )
We make our calculations with several significant digits (a minimum of 3000 significant digits) to minimize the round-off error. Moreover, the computational order of convergence is provided up to five significant digits. Finally, we display the initial guess and approximated zeros up to 25 significant digits in the corresponding example where an exact solution is not available.
All computations have been performed using the programming package Mathematica 11 with multiple precision arithmetic. Further, the meaning of a ( ± b ) is shorthand for a × 10 ( ± b ) in the numerical results.
Example 1.
Population growth problem:
The law of population growth is defined as follows:
d N ( t ) d t = γ N ( t ) + η ,
where N(t) = the population at time t, η = the fixed/constant immigration rate, and γ = the fixed/constant birth rate of the population. We can easily obtain the following solution of the above differential equation:
N ( t ) = N 0 e γ t + η γ ( e γ t 1 ) ,
where N 0 is the initial population.
For a particular case study, the problem is given as follows: Suppose a certain population contains 1,000,000 individuals initially, that 300,000 individuals immigrate into the community in the first year, and that 1,365,000 individuals are present at the end of one year. Find the birth rate ( γ ) of this population.
To determine the birth rate, we must solve the equation:
f 1 ( x ) = 1365 1000 e x 300 x ( e x 1 ) .
wherein x = γ and our desired zero of the above function f 1 is 0.05504622451335177827483421. The reason for considering the simple zero problem is to confirm that our methods also work for simple zeros. We choose the starting point as x 0 = 0.5 .
Example 2.
The van der Waals equation of state:
P + a 1 n 2 V 2 ( V n a 2 ) = n R T ,
explains the behavior of a real gas by introducing in the ideal gas equations two parameters, a 1 and a 2 , specific for each gas. The determination of the volume V of the gas in terms of the remaining parameters requires the solution of a nonlinear equation in V,
P V 3 ( n a 2 P + n R T ) V 2 + a 1 n 2 V a 1 a 2 n 2 = 0 .
Given the constants a 1 and a 2 of a particular gas, one can find values for n , P , and T, such that this equation has three simple roots. By using the particular values, we obtain the following nonlinear function:
f 2 ( x ) = x 3 5.22 x 2 + 9.0825 x 5.2675 .
which has three zeros; out of them, one is the multiple zero α = 1.75 of multiplicity two, and the other is the simple zero α = 1.72 . Our desired root is α = 1.75 , and we chose x 0 = 1.8 as the initial guess.
Example 3.
Eigenvalue problem:
For this, we choose the following 8 × 8 matrix:
A = 12 12 36 12 0 0 12 8 148 129 397 147 12 6 109 74 72 62 186 66 8 4 54 36 32 24 88 36 0 0 24 16 20 13 45 19 8 6 13 10 120 98 330 134 8 24 90 60 132 109 333 115 12 6 105 66 0 0 0 0 0 0 0 4 .
The corresponding characteristic polynomial of this matrix is as follows:
f 3 ( x ) = ( x 4 ) 3 ( x + 4 ) ( x 8 ) ( x 20 ) ( x 12 ) ( x + 12 ) .
The above function has one multiple zero at α = 4 of multiplicity three. In addition, we consider x 0 = 2.7 as the starting point.
Example 4.
Let us consider the following polynomial equation:
f 4 ( z ) = ( x 1 ) 3 1 50 .
The desired zero of the above function f 4 is α = 2 with multiplicity of order 50, and we choose initial guess x 0 = 2.1 for this problem.

4. Conclusions

We presented an eighth-order modification of the Chebyshev–Halley-type iteration scheme having optimal convergence to obtain the multiple solutions of the scalar equation. The proposed scheme is optimal in the sense of the classical Kung–Traub conjecture. Thus, the efficiency index of the present methods is E = 8 4 1.682 , which is better than the classical Newton’s method E = 2 2 1.414 . Finally, the numerical experience corroborates the theoretical results about the convergence order, and moreover, it can be concluded that our proposed methods are highly efficient and competitive.

Author Contributions

R.B. and E.M. have contribute to the theoretical results and F.C. and D.A. have carried out the numerical experience.

Funding

This research was partially supported by Ministerio de Economía y Competitividad under Grant MTM2014-52016-C2-1-2-P and by the project of Generalitat Valenciana Prometeo/2016/089.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Traub, J.F. Iterative Methods for the Solution of Equations; Prentice-Hall: Englewood Cliffs, NJ, USA, 1964. [Google Scholar]
  2. Gutiérrez, J.M.; Hernández, M.A. A family of Chebyshev–Halley type methods in Banach spaces. Bull. Aust. Math. Soc. 1997, 55, 113–130. [Google Scholar] [CrossRef]
  3. Kanwar, V.; Singh, S.; Bakshi, S. Simple geometric constructions of quadratically and cubically convergent iterative functions to solve nonlinear equations. Numer. Algorithms 2008, 47, 95–107. [Google Scholar] [CrossRef]
  4. Potra, F.A.; Pták, V. Nondiscrete Induction and Iterative Processes; Research Notes in Mathematics; Pitman: Boston, MA, USA, 1984; Volume 103. [Google Scholar]
  5. Sharma, J.R. A composite third order Newton-Steffensen method for solving nonlinear equations. Appl. Math. Comput. 2005, 169, 242–246. [Google Scholar]
  6. Argyros, I.K.; Ezquerro, J.A.; Gutiérrez, J.M.; Hernández, M.A.; Hilout, S. On the semilocal convergence of efficient Chebyshev-Secant-type methods. J. Comput. Appl. Math. 2011, 235, 3195–3206. [Google Scholar] [CrossRef]
  7. Xiaojian, Z. Modified Chebyshev–Halley methods free from second derivative. Appl. Math. Comput. 2008, 203, 824–827. [Google Scholar] [CrossRef]
  8. Ostrowski, A.M. Solutions of Equations and System of Equations; Academic Press: New York, NY, USA, 1960. [Google Scholar]
  9. Amat, S.; Hernández, M.A.; Romero, N. A modified Chebyshev’s iterative method with at least sixth order of convergence. Appl. Math. Comput. 2008, 206, 164–174. [Google Scholar] [CrossRef]
  10. Kou, J.; Li, Y. Modified Chebyshev–Halley method with sixth-order convergence. Appl. Math. Comput. 2007, 188, 681–685. [Google Scholar] [CrossRef]
  11. Li, D.; Liu, P.; Kou, J. An improvement of Chebyshev–Halley methods free from second derivative. Appl. Math. Comput. 2014, 235, 221–225. [Google Scholar] [CrossRef]
  12. Sharma, J.R. Improved Chebyshev–Halley methods with sixth and eighth order convergence. Appl. Math. Comput. 2015, 256, 119–124. [Google Scholar] [CrossRef]
  13. Neta, B. Extension of Murakami’s high-order non-linear solver to multiple roots. Int. J. Comput. Math. 2010, 87, 1023–1031. [Google Scholar] [CrossRef]
  14. Zhou, X.; Chen, X.; Song, Y. Constructing higher-order methods for obtaining the multiple roots of nonlinear equations. J. Comput. Math. Appl. 2011, 235, 4199–4206. [Google Scholar] [CrossRef] [Green Version]
  15. Hueso, J.L.; Martínez, E.; Teruel, C. Determination of multiple roots of nonlinear equations and applications. J. Math. Chem. 2015, 53, 880–892. [Google Scholar] [CrossRef]
  16. Behl, R.; Cordero, A.; Motsa, S.S.; Torregrosa, J.R. On developing fourth-order optimal families of methods for multiple roots and their dynamics. Appl. Math. Comput. 2015, 265, 520–532. [Google Scholar] [CrossRef]
  17. Behl, R.; Cordero, A.; Motsa, S.S.; Torregrosa, J.R.; Kanwar, V. An optimal fourth-order family of methods for multiple roots and its dynamics. Numer. Algorithms 2016, 71, 775–796. [Google Scholar] [CrossRef]
  18. Geum, Y.H.; Kim, Y.I.; Neta, B. A class of two-point sixth-order multiple-zero finders of modified double-Newton type and their dynamics. Appl. Math. Comput. 2015, 270, 387–400. [Google Scholar] [CrossRef] [Green Version]
  19. Geum, Y.H.; Kim, Y.I.; Neta, B. A sixth-order family of three-point modified Newton-like multiple-root finders and the dynamics behind their extraneous fixed points. Appl. Math. Comput. 2016, 283, 120–140. [Google Scholar] [CrossRef] [Green Version]
  20. Behl, R.; Cordero, A.; Motsa, S.S.; Torregrosa, J.R. An eighth-order family of optimal multiple root finders and its dynamics. Numer. Algorithms 2018, 77, 1249–1272. [Google Scholar] [CrossRef]
  21. Zafar, F.; Cordero, A.; Motsa, S.S.; Torregrosa, J.R. Optimal iterative methods for finding multiple roots of nonlinear equations using free parameters. J. Math. Chem. 2018, 56, 1884–1901. [Google Scholar] [CrossRef]
  22. Behl, R.; Alshomrani, A.S.; Motsa, S.S. An optimal scheme for multiple roots of nonlinear equations with eighth-order convergence. J. Math. Chem. 2018, 56, 2069–2084. [Google Scholar] [CrossRef]
  23. Behl, R.; Zafar, F.; Alshomrani, A.S.; Junjuaz, M.; Yasmin, N. An optimal eighth-order scheme for multiple zeros of univariate functions. Int. J. Comput. Methods 2018, 1843002. [Google Scholar] [CrossRef]
  24. McNamee, J.M. A comparison of methods for accelerating convergence of Newton’s method for multiple polynomial roots. ACM Signum Newsl. 1998, 33, 17–22. [Google Scholar] [CrossRef]
  25. Cordero, A.; Torregrosa, J.R. Variants of Newton’s method using fifth-order quadrature formulas. Appl. Math. Comput. 2007, 190, 686–698. [Google Scholar] [CrossRef]
Table 1. Comparison on the basis of the difference between two consecutive iterations | x n + 1 x n | for the functions f 1 f 4 .
Table 1. Comparison on the basis of the difference between two consecutive iterations | x n + 1 x n | for the functions f 1 f 4 .
fn OM 1 OM 2 OM 3 OM 4 GM 1 GM 2 ZM 1 ZM 2 BM 1 BM 2
f 1 12.3 (−3)8.4 (−4)9.3 (−5)3.5 (−5)*3.6 (−5)1.6 (−4)2.3 (−4)7.6 (−5)3.7 (−5)
22.0 (−16)9.0 (−20)8.8 (−28)2.0 (−37)*1.4 (−29)4.2 (−31)8.9 (−30)2.6 (−34)5.0 (−37)
39.7 (−95)1.3 (−115)6.4 (−166)2.5 (−295)*5.4 (−173)1.0 (−243)5.5 (−233)5.4 (−270)5.7 (−292)
ρ 5.99976.00006.00018.0000*6.00008.00008.00008.00008.0000
f 2 11.3 (−3)8.2 (−4)4.0 (−3)3.5 (−4)9.5 (−4)3.9 (−4)3.9 (−4)4.1 (−3)2.7 (−4)2.6 (−4)
22.5 (−10)4.2 (−12)6.4 (−16)8.7 (−18)2.7 (−11)1.0 (−14)5.2 (−17)9.8 (−17)1.1 (−18)1.4 (−19)
32.0 (−50)8.7 (−62)6.5 (−87)1.5 (−126)2.0 (−56)3.9 (−78)5.9 (−120)1.2 (−117)6.3 (−134)1.0 (−141)
ρ 5.97575.99286.02147.99635.98365.99757.99457.99417.99718.0026
f 3 19.1 (−5)3.6 (−5)8.0 (−6)6.0 (−6)8.5 (−5)4.8 (−5)4.9 (−6)5.2 (−6)2.0 (−6)1.8 (−6)
21.8 (−28)1.4 (−31)9.8 (−38)2.0 (−47)1.0 (−28)5.0 (−31)6.0 (−48)1.0 (−47)1.5 (−51)2.8 (−52)
31.2 (−170)4.4 (−190)3.3 (−229)2.5 (−379)3.1 (−172)5.8 (−187)2.7 (−383)2.3 (−381)1.4 (−412)1.3 (−418)
ρ 6.00006.00006.00008.00006.00006.00008.00008.00008.00008.0000
f 4 12.4 (−5)7.1 (−6)4.2 (−7)1.4 (−7)1.8 (−5)2.0 (−7)4.8 (−7)6.5 (−7)1.9 (−7)6.3 (−8)
21.5 (−26)1.7 (−30)3.9 (−40)6.7 (−54)1.1 (−26)1.8 (−41)5.7 (−49)8.4 (−48)8.0 (−53)4.2 (−57)
37.5 (−154)3.2 (−178)2.6 (−438)1.7 (−424)6.6 (−154)1.0 (−245)2.2 (−384)6.6 (−375)9.6 (−416)5.9 (−169)
ρ 6.00006.00006.00008.00006.00006.00008.00008.00008.00002.2745
* means that the corresponding method does not work.
Table 2. Comparison on the basis of residual errors | f ( x n ) | for the functions f 1 f 4 .
Table 2. Comparison on the basis of residual errors | f ( x n ) | for the functions f 1 f 4 .
fn OM 1 OM 2 OM 3 OM 4 GM 1 GM 2 ZM 1 ZM 2 BM 1 BM 2
f 1 12.71.01.1 (−1)4.2 (−2)*4.4 (−2)1.9 (−1)2.7 (−1)9.2 (−2)4.4 (−2)
22.4 (−13)1.1 (−16)1.1 (−24)2.4 (−34)*1.7 (−26)5.1 (−28)1.1 (−26)3.2 (−31)6.0 (−34)
31.2 (−91)1.6 (−112)7.8 (−163)3.0 (−292)*5.4 (−173)1.2 (−240)6.7 (−230)6.5 (−267)7.0 (−289)
f 2 15.0 (−8)2.1 (−8)4.8 (−9)3.6 (−9)2.8 (−8)4.6 (−9)4.6 (−9)5.1 (−9)2.3 (−9)2.0 (−9)
21.8 (−21)5.3 (−25)1.2 (−32)2.3 (−36)2.2 (−23)3.2 (−30)8.0 (−35)2.9 (−34)3.4 (−38)5.9 (−40)
31.2 (−101)2.2 (−124)1.3 (−174)6.9 (−254)1.2 (−113)4.6 (−157)1.1 (−240)4.3 (−236)1.2 (−268)3.1 (−284)
f 3 14.9 (−8)3.1 (−9)3.1 (−11)1.4 (−11)4.1 (−8)7.4 (−9)7.8 (−12)9.1 (−12)5.2 (−13)3.6 (−13)
23.9 (−79)1.8 (−88)6.1 (−107)4.9 (−136)7.1 (−80)8.0 (−87)1.4 (−137)6.9 (−137)2.1 (−148)1.5 (−150)
31.0 (−505)5.6 (−564)2.4 (−681)1.1 (−1131)1.9 (−510)1.2 (−554)1.3 (−1143)7.5 (−1138)1.9 (−1231)1.3 (−1249)
f 4 11.2 (−207)2.7 (−234)1.1 (−295)3.3 (−319)3.5 (−214)1.0 (−311)6.6 (−293)2.3 (−286)1.8 (−313)6.2 (−337)
21.9 (−1268)2.6 (−1465)3.8 (−1947)1.6 (−2635)1.9 (−1274)9.8 (−2014)3.4 (−2389)9.4 (−2331)9.8 (−2582)1.1 (−2795)
34.2 (−7633)2.3 (−8851)7.5 (−11856)6.1 (−21166)6.0 (−7636)7.3 (−12226)1.6 (−19159)7.1 (−18686)8.8 (−20728)3.4 (−8388)

Share and Cite

MDPI and ACS Style

Behl, R.; Martínez, E.; Cevallos, F.; Alarcón, D. A Higher Order Chebyshev-Halley-Type Family of Iterative Methods for Multiple Roots. Mathematics 2019, 7, 339. https://doi.org/10.3390/math7040339

AMA Style

Behl R, Martínez E, Cevallos F, Alarcón D. A Higher Order Chebyshev-Halley-Type Family of Iterative Methods for Multiple Roots. Mathematics. 2019; 7(4):339. https://doi.org/10.3390/math7040339

Chicago/Turabian Style

Behl, Ramandeep, Eulalia Martínez, Fabricio Cevallos, and Diego Alarcón. 2019. "A Higher Order Chebyshev-Halley-Type Family of Iterative Methods for Multiple Roots" Mathematics 7, no. 4: 339. https://doi.org/10.3390/math7040339

APA Style

Behl, R., Martínez, E., Cevallos, F., & Alarcón, D. (2019). A Higher Order Chebyshev-Halley-Type Family of Iterative Methods for Multiple Roots. Mathematics, 7(4), 339. https://doi.org/10.3390/math7040339

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop