Next Article in Journal
Some Fixed Point Theorems for (ap)-Quasicontractions
Next Article in Special Issue
A New Application of Gauss Quadrature Method for Solving Systems of Nonlinear Equations
Previous Article in Journal
Permutation Symmetry in Coherent Electrons Scattering by Disordered Media
Previous Article in Special Issue
Generating Optimal Eighth Order Methods for Computing Multiple Roots
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Family of Derivative Free Optimal Fourth Order Methods for Computing Multiple Roots

1
Department of Mathematics, Sant Longowal Institute of Engineering and Technology, Longowal, Sangrur 148106, India
2
Department of Mathematics, Chandigarh University, 140413, Mohali, NH-95 Chandigarh-Ludhiana Highway, Punjab 140413, India
3
Department of Physics and Chemistry, Technical University of Cluj-Napoca, 400114 Cluj-Napoca, Romania
4
Institute of Doctoral Studies, Babeş-Bolyai University, 400084 Cluj-Napoca, Romania
*
Authors to whom correspondence should be addressed.
Symmetry 2020, 12(12), 1969; https://doi.org/10.3390/sym12121969
Submission received: 9 November 2020 / Revised: 25 November 2020 / Accepted: 26 November 2020 / Published: 28 November 2020
(This article belongs to the Special Issue Symmetry in Numerical Analysis and Numerical Methods)

Abstract

:
Many optimal order multiple root techniques, which use derivatives in the algorithm, have been proposed in literature. But contrarily, derivative free optimal order techniques for multiple root are almost nonexistent. By this as an inspirational factor, here we present a family of optimal fourth order derivative-free techniques for computing multiple roots of nonlinear equations. At the beginning the convergence analysis is executed for particular values of multiplicity afterwards it concludes in general form. Behl et. al derivative-free method is seen as special case of the family. Moreover, the applicability and comparison is demonstrated on different nonlinear problems that certifies the efficient convergent nature of the new methods. Finally, we conclude that our new methods consume the lowest CPU time as compared to the existing ones. This illuminates the theoretical outcomes to a great extent of this study.

1. Introduction

Construction of optimal higher-order methods, in the sense of Kung-Traub conjecture [1], free from the derivatives, is always required for the multiple roots of nonlinear function of the form χ ( x ) = 0 with multiplicity θ , i.e., χ ( j ) ( α ) = 0 , j = 0 , 1 , 2 , , θ 1 and χ ( θ ) ( α ) 0 . The well-known Newton’s method [2] is one of the simplest method for obtaining multiple roots of the nonlinear function, which is given by
x t + 1 = x t θ χ ( x t ) χ ( x t ) , t = 0 , 1 , 2 , .
Numerous higher order methods, have been developed in literature by Dong [3], Geum et al. [4], Hansen [5], Li et al. [6,7], Neta [8], Osada [9], Sharifi et al. [10], Sharma and Sharma [11], Zhou et al. [12], Victory and Neta [13], Agarwal et al. [14] and Soleymani et al. [15]. Such methods require the evaluations of derivatives. The without derivative methods are important in case where derivative χ of χ is very small or is costly to evaluate. One such without derivative method is the Traub-Steffensen method [16] which used
χ ( x t ) χ ( x t + b χ ( x t ) ) χ ( x t ) b χ ( x t ) , b R { 0 } ,
or
χ ( x t ) χ [ w t , x t ] ,
for the derivative χ in Newton method (1). Here w t = x t + b χ ( x t ) and χ [ w , x ] = χ ( w ) χ ( x ) w x is divided difference. Then method (1) takes the form of
x t + 1 = x t θ χ ( x t ) χ [ w t , x t ] .
Very recently, researchers have proposed some higher order derivative free methods. For example; Kumar et al. [17] have developed quadratically convergent method, Sharma et al. [18,19], Kumar et al. [20] and Behl et al. [21] developed fourth methods, and Sharma et al. [22] developed eighth order methods for computing the multiple solutions. The methods of [17,18,19,20,21,22] require two, three and four function evaluations per step and, therefore, according to Kung-Traub conjecture these possess optimal convergence [1]. Our main objective of this work is to develop derivative-free multiple root methods of high computational efficiency, which may attain a high convergence order using as small number of function evaluations as possible. Consequently, we develop a class of two-step derivative-free methods with fourth order of convergence. The presented scheme requires three function evaluations per step and, hence, it satisfy optimal criteria [1]. The methodology is based on the classical Traub-Steffensen method (2) and further modified by employing Traub-Steffensen-like iteration in the second step.

2. Construction of Method

Consider the following two-step iterative scheme θ 2 :
z t = x t θ χ ( x t ) χ [ w t , x t ] , x t + 1 = z t θ H ( s t , k t ) 1 2 s t χ ( x t ) χ [ w t , x t ] ,
where s t = χ ( z t ) χ ( x t ) θ , k t = χ ( z t ) χ ( w t ) θ and H : C 2 C is analytic in a neighborhood of ( 0 , 0 ) . The second step is weighted by the factor H ( s , k ) , so we can call it weight factor or more appropriately weight function.
In Theorems 1–3, we demonstrate that the presented iterative scheme (3) attains highest fourth-order of convergence, without adding any extra evaluation of function or its derivative.
Theorem 1.
Assume that χ : C C is an analytic function in a domain containing a multiple zero (say, α) with multiplicity θ = 2 . Suppose that the initial point x 0 is close enough to α, then the convergence order of the Formula (3) is at least 4, provided that H 00 = 0 , H 10 = 1 2 , H 01 = 1 2 , H 20 = H 02 2 H 11 , H 02 R and H 11 R , where H i j = i + j s i k j H ( s t , k t ) | ( s t = 0 , k t = 0 ) , for 0 i , j 2 .
Proof. 
Assume that the error at t-th stage is e t = x t α . Using the Taylor’s expansion of χ ( x t ) about α and keeping into mind that χ ( α ) = 0 , χ ( α ) = 0 and χ ( α ) 0 , we have
χ ( x t ) = χ ( α ) 2 ! e t 2 1 + B 1 e t + B 2 e t 2 + B 3 e t 3 + B 4 e t 4 + ,
where B n = 2 ! ( 2 + n ) ! χ ( 2 + n ) ( α ) χ ( α ) for n N .
Similarly χ ( w t ) about α , we have
χ ( w t ) = χ ( α ) 2 ! e w t 2 1 + B 1 e w t + B 2 e w t 2 + B 3 e w t 3 + B 4 e w t 4 + ,
where e w t = w t α = e t + b χ ( α ) 2 ! e t 2 1 + B 1 e t + B 2 e t 2 + B 3 e t 3 + B 4 e t 4 + .
Then the first step of (3) yields
e z t = z t α = 1 2 b χ ( α ) 2 + B 1 e t 2 1 16 ( b χ ( α ) ) 2 8 b χ ( α ) B 1 + 12 B 1 2 16 B 2 e t 3 + 1 64 ( ( b χ ( α ) ) 3 20 b χ ( α ) B 1 2 + 72 B 1 3 + 64 b χ ( α ) B 2 10 B 1 ( b χ ( α ) ) 2 + 16 B 2 + 96 B 3 ) e t 4 + O ( e t 5 ) .
Expanding χ ( z t ) about α , it follows that
χ ( z t ) = χ ( α ) 2 ! e z t 2 1 + B 1 e z t + B 2 e z t 2 + B 3 e z t 3 + B 4 e z t 4 + .
Using (4), (5) and (7) in s t and k t , after some simple calculations we have
s t = 1 2 b χ ( α ) 2 + B 1 e t 1 16 ( b χ ( α ) ) 2 6 b χ ( α ) B 1 + 16 ( B 1 2 B 2 ) e t 2 + 1 64 ( ( b χ ( α ) ) 3 22 b χ ( α ) B 1 2 + 4 29 B 1 3 + 14 b χ ( α ) B 2 2 B 1 3 ( b χ ( α ) ) 2 + 104 B 2 + 96 B 3 ) e t 3 + O ( e t 4 )
and
k t = 1 2 b χ ( α ) 2 + B 1 e t 1 16 3 ( b χ ( α ) ) 2 2 b χ ( α ) B 1 + 16 ( B 1 2 B 2 ) e t 2 + 1 64 ( 7 ( b χ ( α ) ) 3 + 24 b χ ( α ) B 2 14 b χ ( α ) B 1 2 + 116 B 1 3 2 B 1 11 ( b χ ( α ) ) 2 + 104 B 2 + 96 B 3 ) e t 3 + O ( e t 4 ) .
Taylor expansion of H ( s t , k t ) in the neighborhood of ( 0 , 0 ) is
H ( s t , k t ) H 00 + s t H 10 + k t H 01 + 1 2 s t 2 H 20 + s t k t H 11 + 1 2 k t 2 H 02 .
Using (4)–(10) in the second step of (3), then we have
e t + 1 = H 00 e t 1 2 H 00 + H 01 + H 10 1 b χ ( α ) 2 + B 1 e t 2 + n = 1 2 ψ n e t n + 2 + O ( e t 5 ) ,
where ψ n = ψ n ( b , B 1 , B 2 , B 3 , H 00 , H 10 , H 01 , H 20 , H 11 , H 02 ) .
We will get at least fourth order if we set coefficients of e t , e t 2 and e t 3 simultaneously equal to zero. Then, we have
H 00 = 0 , H 10 = 1 2 , H 01 = 1 2 , H 20 = H 02 2 H 11 .
Now using Equation (12) in (11), we have
e t + 1 = 1 16 b χ ( α ) 2 + B 1 b χ ( α ) ( 2 H 02 + 2 H 11 3 ) B 1 + 3 B 1 2 + ( b χ ( α ) ) 2 ( H 11 + H 02 1 ) 4 B 2 e t 4 + O ( e t 5 ) .
Thus, the theorem is proved. □
Theorem 2.
Using the hypotheses of Theorem 1, the order of convergence of scheme (3) for the case θ = 3 is at least 4, if H 00 = 0 , H 10 = 1 H 01 , H 20 = H 02 2 H 11 , H 01 R , H 02 R and H 11 R .
Proof. 
Keeping into mind that χ ( α ) = 0 , χ ( α ) = 0 , χ ( α ) = 0 , and χ ( α ) 0 , then we have
χ ( x t ) = χ ( α ) 3 ! e t 3 1 + B ¯ 1 e t + B ¯ 2 e t 2 + B ¯ 3 e t 3 + B ¯ 4 e t 4 + ,
where B ¯ n = 3 ! ( 3 + n ) ! χ ( 3 + n ) ( α ) χ ( α ) for n N .
Similarly, χ ( w t ) about α
χ ( w t ) = χ ( α ) 3 ! e w t 3 1 + B ¯ 1 e w t + B ¯ 2 e w t 2 + B ¯ 3 e w t 3 + B ¯ 4 e w t 4 + ,
where e w t = w t α = e t + b χ ( α ) 2 ! e t 3 1 + B ¯ 1 e t + B ¯ 2 e t 2 + B ¯ 3 e t 3 + B ¯ 4 e t 4 + .
Then the first step of (3) yields
e z t = z t α = B ¯ 1 3 e t 2 + b χ ( α ) 6 4 9 B ¯ 1 2 + 2 3 B ¯ 2 e t 3 + 16 27 B ¯ 1 3 + 1 9 B ¯ 1 ( 2 b χ ( α ) 13 B ¯ 2 ) + B ¯ 3 e t 4 + O ( e t 5 ) .
Expanding χ ( z t ) about α , it follows that
χ ( z t ) = χ ( α ) 3 ! e z t 3 1 + B ¯ 1 e z t + B ¯ 2 e z t 2 + B ¯ 3 e z t 3 + B ¯ 4 e z t 4 + .
Using (14), (15) and (17) in s t and k t , after some simple calculations we have
s t = B ¯ 1 3 e t + b χ ( α ) 6 5 9 B ¯ 1 2 + 2 3 B ¯ 2 e t 2 + 23 27 B ¯ 1 3 + B ¯ 1 b χ ( α ) 6 16 9 B ¯ 2 + B ¯ 3 e t 3 + O ( e t 4 )
and
k t = B ¯ 1 3 e t + b χ ( α ) 6 5 9 B ¯ 1 2 + 2 3 B ¯ 2 e t 2 + 23 27 B ¯ 1 3 + 2 9 B ¯ 1 b χ ( α ) 2 8 B ¯ 2 + B ¯ 3 e t 3 + O ( e t 4 ) .
Using (10) and (14)–(19) in the second step of (3), then we have
e k + 1 = H 00 e t 1 3 H 00 + H 01 + H 10 1 e t 2 + n = 1 2 ϕ n e t n + 2 + O ( e t 5 ) ,
where ϕ n = ϕ n ( b , B ¯ 1 , B ¯ 2 , B ¯ 3 , H 00 , H 10 , H 01 , H 20 , H 11 , H 02 ) .
If we set coefficients of e t , e t 2 and e t 3 simultaneously equal to zero. Then we have
H 00 = 0 , H 10 = 1 H 01 , H 20 = H 02 2 H 11 .
Now using Equation (21) in (20), we have
e t + 1 = B ¯ 1 27 3 2 b χ ( α ) ( H 01 1 ) + 2 B ¯ 1 2 3 B ¯ 2 e t 4 + O ( e t 5 ) .
Thus, the theorem is proved. □
Remark 1.
From above results we observe that the number of conditions on H i j is 4 , 3 corresponding to θ = 2 , 3 to obtain the fourth convergence order of the method (3). Their error equations also contain the term involving the parameter b. However, for the cases θ 4 , it has been seen that the error equation in each such case does not contain b term. We shall prove this fact in the next section.

3. Main Result

We shall prove the order of convergence of scheme (3) for θ 4 by the following theorem:
Theorem 3.
Using the hypotheses of Theorem 1, the order of convergence of scheme (3) for the cases θ 4 is at least 4, if H 00 = 0 , H 10 = 1 H 01 , H 20 = H 02 2 H 11 , H 01 R , H 02 R and H 11 R . Moreover, error in the scheme is given by
e t + 1 = 1 2 θ 3 ( 1 + θ ) B ¯ ¯ 1 3 2 θ B ¯ ¯ 1 B ¯ ¯ 2 e t 4 + O ( e t 5 ) .
Proof. 
Keeping into mind that χ ( j ) ( α ) = 0 , j = 0 , , θ 1 and χ ( θ ) ( α ) 0 , then developing χ ( x t ) about α in the Taylor’s series
χ ( x t ) = χ θ ( α ) θ ! e t θ 1 + B ¯ ¯ 1 e t + B ¯ ¯ 2 e t 2 + B ¯ ¯ 3 e t 3 + B ¯ ¯ 4 e t 4 + ,
where B ¯ ¯ n = θ ! ( θ + n ) ! χ ( θ + n ) ( α ) χ ( θ ) ( α ) for n N .
Also from the expansion of χ ( w t ) about α , it follows that
χ ( w t ) = χ θ ( α ) θ ! e w t θ 1 + B ¯ ¯ 1 e w t + B ¯ ¯ 2 e w t 2 + B ¯ ¯ 3 e w t 3 + B ¯ ¯ 4 e w t 4 + ,
where e w t = w t α = e t + β f θ ( α ) θ ! e t θ 1 + B ¯ ¯ 1 e t + B ¯ ¯ 2 e t 2 + B ¯ ¯ 3 e t 3 + B ¯ ¯ 4 e t 4 + .
From the first step of (3)
e z t = B ¯ ¯ 1 4 e t 2 + 1 16 4 B ¯ ¯ 2 3 B ¯ ¯ 1 2 e t 3 + 25 64 B ¯ ¯ 1 3 B ¯ ¯ 1 B ¯ ¯ 2 + 1 16 ( b χ ( 4 ) ( α ) + 12 B ¯ ¯ 3 ) e t 4 + O ( e t 5 ) , if θ = 4 . B ¯ ¯ 1 θ e t 2 + 1 θ 2 2 θ B ¯ ¯ 2 ( 1 + θ ) B ¯ ¯ 1 2 e t 3 + 1 θ 3 ( 1 + θ ) 2 B ¯ ¯ 1 3 θ ( 4 + 3 θ ) B ¯ ¯ 1 B ¯ ¯ 2 + 3 θ 2 B ¯ ¯ 3 e t 4 + O ( e t 5 ) , if θ 5 .
Expansion of χ ( z t ) around α yields
χ ( z t ) = χ ( θ ) ( α ) θ ! e z t θ 1 + B ¯ ¯ 1 e z t + B ¯ ¯ 2 e z t 2 + B ¯ ¯ 3 e z t 3 + B ¯ ¯ 4 e z t 4 + .
Using (23), (24) and (26) in the expressions of s t and k t , we have that
s k = B ¯ ¯ 1 4 e t + 1 8 4 B ¯ ¯ 2 3 B ¯ ¯ 1 2 e t 2 + 1 128 67 B ¯ ¯ 1 3 152 B ¯ ¯ 1 B ¯ ¯ 2 + 8 ( b χ ( 4 ) ( α ) + 12 B ¯ ¯ 3 ) e t 3 + O ( e t 4 ) , if θ = 4 . B ¯ ¯ 1 θ e t + 1 θ 2 2 θ B ¯ ¯ 2 ( 2 + θ ) B ¯ ¯ 1 2 e t 2 + 1 2 θ 3 ( 2 θ 2 + 7 θ + 7 ) B ¯ ¯ 1 3 2 θ ( 7 + 3 θ ) B ¯ ¯ 1 B ¯ ¯ 2 + 6 θ 2 B ¯ ¯ 3 e t 3 + O ( e t 4 ) , if θ 5
and
k t = B ¯ ¯ 1 4 e t + 1 8 4 B ¯ ¯ 2 3 B ¯ ¯ 1 2 e t 2 + 1 128 67 B ¯ ¯ 1 3 152 B ¯ ¯ 1 B ¯ ¯ 2 + 8 ( b χ ( 4 ) ( α ) + 12 B ¯ ¯ 3 ) e t 3 + O ( e t 4 ) , if θ = 4 . B ¯ ¯ 1 θ e t + 1 θ 2 2 θ B ¯ ¯ 2 ( 2 + θ ) B ¯ ¯ 1 2 e t 2 + 1 2 θ 3 ( 2 θ 2 + 7 θ + 7 ) B ¯ ¯ 1 3 2 θ ( 7 + 3 θ ) B ¯ ¯ 1 B ¯ ¯ 2 + 6 θ 2 B ¯ ¯ 3 e t 3 + O ( e t 4 ) , if θ 5 .
Inserting (10) and (23)–(28) in the second step of (3), it follows that
e t + 1 = H 00 e t + 1 θ ( H 00 + H 01 + H 10 1 ) B ¯ ¯ 1 e t 2 + n = 1 2 φ n e t n + 2 + O ( e t 5 ) .
where φ n = φ n ( b , B ¯ ¯ 1 , B ¯ ¯ 2 , B ¯ ¯ 3 , H 00 , H 10 , H 01 , H 20 , H 11 , H 02 ) , for θ = 4 and φ n = φ n ( B ¯ ¯ 1 , B ¯ ¯ 2 , B ¯ ¯ 3 , H 00 , H 10 , H 01 , H 20 , H 11 , H 02 ) , for θ 5 .
If the coefficients of e t , e t 2 and e t 3 vanish then we have
H 00 = 0 , H 10 = 1 H 01 , H 20 = H 02 2 H 11 .
Then, error Equation (29) is given by
e t + 1 = 1 2 θ 3 ( 1 + θ ) B ¯ ¯ 1 3 2 θ B ¯ ¯ 1 B ¯ ¯ 2 e t 4 + O ( e t 5 ) .
Thus, the theorem is proved. □
Remark2. 
This fourth order convergence rate is achieved by using only χ ( x t ) , χ ( w t ) and χ ( z t ) per iteration. Therefore, the scheme (3) is optimal by the Kung-Traub conjecture [1].
Remark3. 
Note that parameter b, which is utilized in w t , shows up just in the error equations of the cases θ = 2 , 3 yet not for θ 4 . We have seen that this parameter appears in the coefficients of e t 5 and higher order. However, we do not need such terms in order to show the required fourth order convergence.

Some Special Cases

Based on the forms of function H ( s , k ) that satisfy the conditions of Theorems 1–3. Then, we get a new optimal family of order fourth as follows:
z t = x t θ χ ( x t ) χ [ w t , x t ] , x t + 1 = z t θ s t ( 1 H 01 ) + k t H 01 1 2 s t 2 ( H 02 + 2 H 11 ) + s t k t H 11 + 1 2 k t 2 H 02 1 2 s t χ ( x t ) χ [ w t , x t ] .
( 1 )
For H 01 = 1 2 , H 02 = 0 and H 11 = 0 in expression (32), we have
x t + 1 = z t θ s t + k t 2 ( 1 2 s t ) χ ( x t ) χ [ w t , x t ] .
It is important to note that the above method (33) is Behl et al. method [21]. This shows that Behl et al. method [21] is the special case of our family (32).
( 2 )
If H 01 = 1 2 , H 02 = 0 and H 11 = 1 in expression (32), we have
x t + 1 = z t θ s t 2 s t 2 + k t + 2 s t k t 2 ( 1 2 s t ) χ ( x t ) χ [ w t , x t ] .
( 3 )
If H 01 = 1 2 , H 02 = 1 and H 11 = 0 in expression (32), we get
x t + 1 = z t θ s t + s t 2 + k t k t 2 2 ( 1 2 s t ) χ ( x t ) χ [ w t , x t ] .
( 4 )
Let H 01 = 1 2 , H 02 = 1 and H 11 = 1 in expression (32), we obtain
x t + 1 = z t θ s t + s t 2 + k t 2 s t k t + k t 2 2 ( 1 2 s t ) χ ( x t ) χ [ w t , x t ] .
( 5 )
Let H 01 = θ 1 2 , H 02 = 0 and H 11 = 0 in expression (32), we have
x t + 1 = z t θ ( 3 θ ) s t + ( θ 1 ) k t 2 ( 1 2 s t ) χ ( x t ) χ [ w t , x t ] .
In above each case z t = x t θ χ ( x t ) χ [ w t , x t ] . For future reference the proposed methods (33), (34),(35), (36) and (37) are denoted by BM, NM1, NM2, NM3 and NM4, respectively.

4. Numerical Results

In order to validate of theoretical results that have been proven in previous sections, the new methods BM, NM1, NM2, NM3 and NM4 are checked numerically by imposing them on some nonlinear equations. Moreover, they are also compared with some existing derivative free optimal fourth order methods. We consider, for example, the methods by Sharma et al. [18,19] and Kumar et al. [20]. The methods are expressed as follows:
Method by Sharma et al. [18] (SK1):
z t = x t θ χ ( x t ) χ [ w t , x t ] , x t + 1 = z t s t + ( θ 1 ) k t + θ s t 2 + θ s t k t χ ( x t ) χ [ w t , x t ] .
Method by Sharma et al. [19] (SK2):
z t = x t θ χ ( x t ) χ [ w t , x t ] , x t + 1 = z t θ 2 h + 3 θ h 2 2 1 + 1 d t χ ( x t ) χ [ w t , x t ] ,
where
h = s t 1 + s t , and d t = χ ( w t ) χ ( x t ) θ .
Method by Kumar et al. [20] (KM):
z t = x t θ χ ( x t ) χ [ w t , x t ] , x t + 1 = z t ( θ + 2 ) s t 1 2 s t χ ( x t ) χ [ w t , x t ] + χ [ w t , z t ] .
Computational work is compiled in the programming software, e.g., Mathematica [23]. Performance of the new methods is tested by selecting value of the parameter b = 0.01 . The tabulated results obtained by the methods for each problem include (i) number of iterations ( t ) required to obtain the solution using the stopping criterion | x t + 1 x t | + | χ ( x t ) | < 10 100 , (ii) estimated error | x t + 1 x t | in the first three iterations, (iii) calculated convergence order (CCO) and (iv) elapsed time (CPU time in seconds), which is measured by the command “TimeUsed[ ]” (Table 1. The calculated convergence order (CCO) to confirm the theoretical convergence order is calculated by the formula (see [24])
CCO = log | ( x t + 2 α ) / ( x t + 1 α ) | log | ( x t + 1 α ) / ( x t α ) | , for each t = 1 , 2 ,
From the computed results in Table 2, we can observe the good convergence behavior of the proposed methods like that of existing methods. This also explains stable nature of the methods. It is also clear that the approximations to the solutions by the proposed methods have greater or equal accuracy than those computed by existing methods. We display the value 0 of | x t + 1 x t | at the stage when stopping criterion | x t + 1 x t | + | χ ( x t ) | < 10 100 has been satisfied. From the calculation of computational order of convergence shown in each table, we verify the fourth order of convergence. The efficient nature of presented methods can be verified by the fact that the amount of CPU time consumed by the methods is less than that of the time taken by existing methods. This conclusion is also confirmed by similar numerical experiments on many other different problems.

5. Conclusions

In the foregoing study, we have proposed a family of fourth order derivative-free numerical methods for solving nonlinear equations with multiple roots of known multiplicity. Analysis of the convergence has been carried out, which proves the order four under standard assumptions of the function whose zeros we are looking for. In addition, our proposed scheme also satisfies the Kung-Traub hypothesis of optimal order of convergence. Some special cases have been discussed. These are employed to solve nonlinear equations including those arising in practical problems. The new methods are compared with existing techniques of same order. We conclude the work with a remark that derivative-free methods are good alternatives to Newton-type schemes in the cases when derivatives are expensive to compute or difficult to obtain.

Author Contributions

Methodology, S.K. and D.K.; Formal analysis, J.R.S.; Investigation, S.K.; Data Curation, D.K.; Conceptualization, L.J.; and Writing—review and editing, S.K. and D.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We would like to express our gratitude to the anonymous reviewers for their help with the publication of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kung, H.T.; Traub, J.F. Optimal order of one-point and multipoint iteration. J. Assoc. Comput. Mach. 1974, 21, 643–651. [Google Scholar] [CrossRef]
  2. Schröder, E. Über unendlich viele Algorithmen zur Auflösung der Gleichungen. Math. Ann. 1870, 2, 317–365. [Google Scholar] [CrossRef] [Green Version]
  3. Dong, C. A family of multipoint iterative functions for finding multiple roots of equations. Int. J. Comput. Math. 1987, 21, 363–367. [Google Scholar] [CrossRef]
  4. 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]
  5. Hansen, E.; Patrick, M. A family of root finding methods. Numer. Math. 1977, 27, 257–269. [Google Scholar] [CrossRef]
  6. Li, S.; Liao, X.; Cheng, L. A new fourth-order iterative method for finding multiple roots of nonlinear equations. Appl. Math. Comput. 2009, 215, 1288–1292. [Google Scholar]
  7. Li, S.G.; Cheng, L.Z.; Neta, B. Some fourth-order nonlinear solvers with closed formulae for multiple roots. Comput Math. Appl. 2010, 59, 126–135. [Google Scholar] [CrossRef] [Green Version]
  8. Neta, B. New third order nonlinear solvers for multiple roots. App. Math. Comput. 2008, 202, 162–170. [Google Scholar] [CrossRef] [Green Version]
  9. Osada, N. An optimal multiple root-finding method of order three. J. Comput. Appl. Math. 1994, 51, 131–133. [Google Scholar] [CrossRef] [Green Version]
  10. Sharifi, M.; Babajee, D.K.R.; Soleymani, F. Finding the solution of nonlinear equations by a class of optimal methods. Comput. Math. Appl. 2012, 63, 764–774. [Google Scholar] [CrossRef] [Green Version]
  11. Sharma, J.R.; Sharma, R. Modified Jarratt method for computing multiple roots. Appl. Math. Comput. 2010, 217, 878–881. [Google Scholar] [CrossRef]
  12. Zhou, X.; Chen, X.; Song, Y. Constructing higher-order methods for obtaining the multiple roots of nonlinear equations. J. Comput. Appl. Math. 2011, 235, 4199–4206. [Google Scholar] [CrossRef] [Green Version]
  13. Victory, H.D.; Neta, B. A higher order method for multiple zeros of nonlinear functions. Int. J. Comput. Math. 1983, 12, 329–335. [Google Scholar] [CrossRef]
  14. Agarwal, R.P.; Gala, S.; Ragusa, M.A.; Salimi, M. A regularity criterion in weak spaces to Boussinesq equations. Mathematics 2020, 8, 920. [Google Scholar] [CrossRef]
  15. Soleymani, F.; Babajee, D.K.R.; Lotfi, T. On a numerical technique for finding multiple zeros and its dynamics. J. Egypt. Math. Soc. 2013, 21, 346–353. [Google Scholar] [CrossRef]
  16. Traub, J.F. Iterative Methods for the Solution of Equations; Chelsea Publishing Company: New York, NY, USA, 1982. [Google Scholar]
  17. Kumar, D.; Sharma, J.R.; Argyros, I.K. Optimal one-point iterative function free from derivatives for multiple roots. Mathematics 2020, 8, 709. [Google Scholar] [CrossRef]
  18. Sharma, J.R.; Kumar, S.; Jäntschi, L. On a class of optimal fourth order multiple root solvers without using derivatives. Symmetry 2019, 11, 1452. [Google Scholar] [CrossRef] [Green Version]
  19. Sharma, J.R.; Kumar, S.; Jäntschi, L. On Derivative Free Multiple-Root Finders with Optimal Fourth Order Convergence. Mathematics 2020, 8, 1096. [Google Scholar] [CrossRef]
  20. Kumar, S.; Kumar, D.; Sharma, J.R.; Cesarano, C.; Agarwal, P.; Chu, Y.M. An Optimal Fourth Order Derivative-Free Numerical Algorithm for Multiple Roots. Symmetry 2020, 12, 1038. [Google Scholar] [CrossRef]
  21. Behl, R.; Alharbi, S.K.; Mallawi, F.O.; Salimi, M. An optimal derivative-free Ostrowski’s scheme for multiple roots of nonlinear equations. Mathematics 2020, 8, 1809. [Google Scholar] [CrossRef]
  22. Sharma, J.R.; Kumar, S.; Argyros, I.K. Development of optimal eighth order derivative-free methods for multiple roots of nonlinear equations. Symmetry 2019, 11, 766. [Google Scholar] [CrossRef] [Green Version]
  23. Wolfram, S. The Mathematica Book, 4th ed.; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
  24. Weerakoon, S.; Fernando, T.G.I. A variant of Newton’s method with accelerated third-order convergence. Appl. Math. Lett. 2000, 13, 87–93. [Google Scholar] [CrossRef]
  25. Douglas, J.M. Process Dynamics and Control; Prentice Hall: Englewood Cliffs, NJ, USA, 1972. [Google Scholar]
  26. Behl, R.; Zafar, F.; Alshormani, A.S.; Junjua, M.U.D.; Yasmin, N. An optimal eighth-order scheme for multiple zeros of unvariate functions. Int. J. Comput. Meth. 2018. [Google Scholar] [CrossRef]
  27. Bradie, B. A Friendly Introduction to Numerical Analysis; Pearson Education Inc.: New Delhi, India, 2006. [Google Scholar]
  28. Hoffman, J.D. Numerical Methods for Engineers and Scientists; McGraw-Hill Book Company: New York, NY, USA, 1992. [Google Scholar]
  29. Zeng, Z. Computing multiple roots of inexact polynomials. Math. Comput. Lett. 2004, 74, 869–903. [Google Scholar] [CrossRef] [Green Version]
Table 1. Following problems are considered in this paper.
Table 1. Following problems are considered in this paper.
ProblemsRootMultiplicityInitial Guess
Isothermal continuous stirred tank reactor problem [25]
χ 1 ( x ) = x 4 + 11.50 x 3 + 47.49 x 2 + 83.06325 x + 51.23266875 −2.852−2.7
Van der Waals problem [26]
χ 2 ( x ) = x 3 5.22 x 2 + 9.0825 x 5.2675 1.7522
Planck law radiation problem [27]
χ 3 ( x ) = e x 1 + x 5 3 4.9651142317…35.5
Manning problem for isentropic supersonic flow [28]
χ 4 ( x ) = [ tan 1 5 2 tan 1 ( x 2 1 ) + 6 ( tan 1 x 2 1 6
tan 1 1 2 5 6 ) 11 63 ] 4 1.8411294068…41.2
Standard test problem [20]
χ 5 ( x ) = x ( x 2 + 1 ) ( 2 e x 2 + 1 + x 2 1 ) cosh 3 π x 2 i51.2i
Clustering problem [29]
χ 6 ( x ) = ( x 2 ) 15 ( x 4 ) 5 ( x 3 ) 10 ( x 1 ) 20 1200.7
Table 2. Numerical results for problems.
Table 2. Numerical results for problems.
Methodst | x 2 x 1 | | x 3 x 2 | | x 4 x 3 | CCOCPU-Time
Problem χ 1 ( x )
SK14 5.02 × 10 3 4.91 × 10 12 4.84 × 10 48 4.0000.0812
SK24 5.02 × 10 3 5.00 × 10 12 5.34 × 10 48 4.0000.0853
KM4 5.02 × 10 3 4.83 × 10 12 4.83 × 10 48 4.0000.0798
BM4 5.02 × 10 3 4.84 × 10 12 4.48 × 10 48 4.0000.0788
NM14 5.02 × 10 3 4.85 × 10 12 4.54 × 10 48 4.0000.0778
NM24 5.02 × 10 3 4.82 × 10 12 4.41 × 10 48 4.0000.0779
NM34 5.02 × 10 3 4.84 × 10 12 4.48 × 10 48 4.0000.0784
NM44 5.02 × 10 3 4.84 × 10 12 4.48 × 10 48 4.0000.0783
Problem χ 2 ( x )
SK16 3.03 × 10 2 1.26 × 10 3 5.30 × 10 8 4.0000.0724
SK26 3.40 × 10 2 2.14 × 10 3 6.88 × 10 7 4.0000.0942
KM5 2.25 × 10 2 2.69 × 10 4 2.37 × 10 11 4.0000.0704
BM5 2.34 × 10 2 3.43 × 10 4 9.30 × 10 11 4.0000.0692
NM15 2.34 × 10 2 3.44 × 10 4 9.34 × 10 11 4.0000.0654
NM25 2.34 × 10 2 3.43 × 10 4 9.26 × 10 11 4.0000.0472
NM35 2.34 × 10 2 3.43 × 10 4 9.30 × 10 11 4.0000.0494
NM45 2.34 × 10 2 3.43 × 10 4 9.30 × 10 11 4.0000.0502
Problem χ 3 ( x )
SK13 5.56 × 10 6 1.32 × 10 25 04.0000.4962
SK23 6.34 × 10 6 2.70 × 10 25 04.0000.4726
KM3 4.93 × 10 6 6.76 × 10 26 04.0000.4137
BM3 4.91 × 10 6 6.62 × 10 26 04.0000.4232
NM13 4.91 × 10 6 6.62 × 10 26 04.0000.4062
NM23 4.91 × 10 6 6.61 × 10 26 04.0000.4204
NM33 4.91 × 10 6 6.62 × 10 26 04.0000.4247
NM43 4.94 × 10 6 6.86 × 10 26 04.0000.4251
Problem χ 4 ( x )
SK15 2.88 × 10 1 2.21 × 10 2 3.24 × 10 9 4.0003.3120
SK25 2.73 × 10 1 1.97 × 10 2 2.84 × 10 9 4.0003.2642
KM5 3.11 × 10 1 2.60 × 10 2 4.32 × 10 9 4.0003.3230
BM5 3.11 × 10 1 2.60 × 10 2 4.31 × 10 9 4.0003.2114
NM15 3.11 × 10 1 2.60 × 10 2 4.31 × 10 9 4.0003.1423
NM25 3.11 × 10 1 2.60 × 10 2 4.31 × 10 9 4.0003.1876
NM35 3.11 × 10 1 2.60 × 10 2 4.31 × 10 9 4.0003.2591
NM45 3.11 × 10 1 2.60 × 10 2 4.32 × 10 9 4.0002.9642
Problem χ 5 ( x )
SK14 7.14 × 10 5 5.13 × 10 18 1.36 × 10 70 4.0000.5691
SK24 7.93 × 10 5 1.16 × 10 17 5.21 × 10 69 4.0000.5724
KM4 6.43 × 10 5 2.07 × 10 18 2.22 × 10 72 4.0000.5772
BM4 6.66 × 10 5 2.38 × 10 18 3.91 × 10 72 4.0000.5547
NM14 6.65 × 10 5 2.37 × 10 18 3.84 × 10 72 4.0000.5462
NM24 6.67 × 10 5 2.39 × 10 18 3.98 × 10 72 4.0000.5531
NM34 6.66 × 10 5 2.38 × 10 18 3.91 × 10 72 4.0000.5684
NM44 6.12 × 10 5 1.70 × 10 18 1.00 × 10 72 4.0000.5642
Problem χ 6 ( x )
SK14 9.74 × 10 3 5.21 × 10 8 4.57 × 10 29 4.0000.1377
SK25 1.39 × 10 2 4.13 × 10 7 3.65 × 10 25 4.0000.1421
KM4 3.41 × 10 3 1.50 × 10 10 5.63 × 10 40 4.0000.1324
BM4 3.42 × 10 3 1.51 × 10 10 5.86 × 10 40 4.0000.1257
NM14 3.41 × 10 3 1.51 × 10 10 5.83 × 10 40 4.0000.1246
NM24 3.42 × 10 3 1.51 × 10 10 5.89 × 10 40 4.0000.1098
NM34 3.42 × 10 3 1.51 × 10 10 5.86 × 10 40 4.0000.1249
NM44 3.35 × 10 3 1.40 × 10 10 4.34 × 10 40 4.0000.0914
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kumar, S.; Kumar, D.; Sharma, J.R.; Jäntschi, L. A Family of Derivative Free Optimal Fourth Order Methods for Computing Multiple Roots. Symmetry 2020, 12, 1969. https://doi.org/10.3390/sym12121969

AMA Style

Kumar S, Kumar D, Sharma JR, Jäntschi L. A Family of Derivative Free Optimal Fourth Order Methods for Computing Multiple Roots. Symmetry. 2020; 12(12):1969. https://doi.org/10.3390/sym12121969

Chicago/Turabian Style

Kumar, Sunil, Deepak Kumar, Janak Raj Sharma, and Lorentz Jäntschi. 2020. "A Family of Derivative Free Optimal Fourth Order Methods for Computing Multiple Roots" Symmetry 12, no. 12: 1969. https://doi.org/10.3390/sym12121969

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