Next Article in Journal
A 6-Letter ‘DNA’ for Baskets with Handles
Previous Article in Journal
Certain Hermite–Hadamard Inequalities for Logarithmically Convex Functions with Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Derivative-Free Root Finding Methods Based on Inverse Interpolation

Centre for Advanced Studies in Pure and Applied Mathematics, Bahauddin Zakariya University, Multan 60800, Pakistan
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(2), 164; https://doi.org/10.3390/math7020164
Submission received: 12 December 2018 / Revised: 20 January 2019 / Accepted: 22 January 2019 / Published: 12 February 2019
(This article belongs to the Section Mathematics and Computer Science)

Abstract

:
Finding a simple root for a nonlinear equation f ( x ) = 0 , f : I R R has always been of much interest due to its wide applications in many fields of science and engineering. Newton’s method is usually applied to solve this kind of problems. In this paper, for such problems, we present a family of optimal derivative-free root finding methods of arbitrary high order based on inverse interpolation and modify it by using a transformation of first order derivative. Convergence analysis of the modified methods confirms that the optimal order of convergence is preserved according to the Kung-Traub conjecture. To examine the effectiveness and significance of the newly developed methods numerically, several nonlinear equations including the van der Waals equation are tested.

1. Introduction

In this paper, we present optimal derivative-free methods to solve a nonlinear equation of the form f ( x ) = 0 [1,2,3]. Multipoint iterative methods for this problem have been extensively studied in the last decade as they are computationally efficient than the one-point methods such as the methods of Newton, Halley and Laguerre. According to the conjecture of Kung and Traub [4], the order of convergence of any multipoint method requiring n + 1 evaluations cannot exceed the bound 2 n . The methods that satisfy this bound are called optimal methods. There is a vast literature on optimal multipoint methods, which are developed by using the famous one-step Newton method or the Steffensen method at the first step. The following is the iteration of Newton’s scheme to find a simple root α of a nonlinear equation f ( x ) = 0 , where f : I R R is a scalar function on an open interval I [1]:
x n + 1 = x n f ( x n ) f ( x n ) , n 0 .
For a background study of multi-point optimal methods, one may consult [2,3,5,6,7,8,9]. Steffensen’s iterative scheme is a well-known modification of Newton’s method obtained by using the approximation
f ( x n ) f ( x n ) f ( z n ) x n z n = f [ x n , z n ] ,
in the Newton’s scheme and is given as follows [10]:
x n + 1 = x n f ( x n ) f [ z n , x n ] , n 0 ,
where z n = x n + f ( x n ) . Both methods are quadratic in some neighborhood of α but Steffensen’s method has an advantage that it does not need the evaluation of the function’s derivative, which may be problematic and expensive to calculate for certain functions. To determine the computational efficiency of an iterative method, Ostrowski [11] defined the efficiency index as q 1 / n , where q is the convergence order and n is the number of functional evaluations per iterative step. In a recent paper [6], Cordero and Torregrosa conjectured an approximation of the first derivative:
f ( x n ) f [ x n , z n ] , z n = x n + γ f ( x n ) m , m q , γ R { 0 }
to transform a multipoint with derivative iterative method of order 2 q to a derivative-free method possessing the same order. Here, we use this conjecture by omitting the parameter γ to develop a family of n-point optimal derivative-free methods based on inverse interpolation. In Section 2, we present optimal derivative-free iterative methods based on inverse interpolation. The derivative-free forms of the methods are obtained using the conjecture of Cordero and Torregrosa [6] such that the convergence order is preserved. In Section 3, we consider some nonlinear equations and van der Waals equation for the numerical comparisons of presented methods with the existing methods of same kind.

2. Optimal Iterative Methods Based on Inverse Interpolation

In this section, we present optimal iterative methods based on inverse interpolation involving first derivative. The first derivative is further approximated by using the conjecture of Cordero and Torregrosa [6] to transform derivative methods into derivative-free methods such that the convergence order is preserved.

2.1. Optimal Two-Point Fourth Order Method

To construct an optimal two-point method, we use the following quadratic polynomial,
x = R ( f ( x ) ) = a 1 + b 1 ( f ( x ) f ( x n ) ) + g 1 ( f ( x ) f ( x n ) ) 2 .
By substituting x = x n into Equation (5), we get:
a 1 = x n = R ( f ( x n ) ) .
Now by differentiating Equation (5) with respect to x, we get:
1 = R ( f ( x ) ) f ( x ) = ( b 1 + 2 g 1 ( f ( x ) f ( x n ) ) 2 ) f ( x ) .
Therefore,
b 1 = 1 f ( x n ) .
Now, by substituting x = y n and using Equations (6) and (7) in Equation (5), we get:
g 1 = 1 [ f ( y n ) f ( x n ) ] f [ y n , x n ] 1 [ f ( y n ) f ( x n ) ] f ( x n ) ,
where y n is the Newton’s iterate and f [ y , x ] = f ( y ) f ( x ) ( y x ) .
Hence, by using Equations (6)–(8) in Equation (5), we obtain the following two-point optimal fourth order method:
y n = x n f ( x n ) f ( x n ) , n 0 , x n + 1 = R ( 0 ) = y n + g 1 f ( x n ) 2 ,
where g 1 is given by Equation (8).
Theorem 1.
Let α be a simple root of f, where f : I R R is a sufficiently differentiable function in an open interval I. If x 0 is sufficiently close to α , then the iterative method in Equation (9) is fourth order convergent and possesses the following error relation:
e n + 1 = c 2 ( c 3 + 2 c 2 2 ) e n 4 + O ( e n 5 ) ,
where c j = f ( j ) ( α ) j ! f ( α ) , j 2 and e n = x n α .
Proof. 
By using Taylor’s expansions about α , we have
f ( x n ) = f ( α ) [ e n + c 2 e n 2 + c 3 e n 3 + c 4 e n 4 ] + O ( e n 5 )
and
f ( x n ) = f ( α ) [ 1 + 2 c 2 e n + 3 c 3 e n 2 + 4 c 4 e n 3 ] + O ( e n 4 )
By using above expressions in the first step of Equation (9), we get
y n α = c 2 e n 2 + ( 2 c 3 2 c 2 2 ) e n 3 + ( 3 c 4 7 c 3 c 2 + 4 c 2 3 ) e n 4 + O ( e n 5 ) .
Again, by using Taylor’s expansion, we have
f ( y n ) = c 2 f ( α ) e n 2 + ( 2 c 3 2 c 2 2 ) f ( α ) e n 3 + ( 3 c 4 7 c 3 c 2 + 4 c 2 3 ) f ( α ) e n 4 + O ( e n 5 ) .
Hence, we get the following error equation of the method in Equation (9) by using the above expressions in its second step:
e n + 1 = c 2 ( c 3 + 2 c 2 2 ) e n 4 + O ( e n 5 ) .
Thus, the proof is complete. □
Now, we modify the new two-point optimal scheme in Equation (9) to obtain a derivative-free method by using the conjecture in Equation (4) such that the optimal order is preserved.
Therefore, with the help of the approximation given in Equation (4), the iterative method in Equation (9) is modified as follows:
y n = x n f ( x n ) f [ z n , x n ] , z n = x n + f ( x n ) 2 , n 0 , x n + 1 = y n + h 1 f ( x n ) 2 ,
where
h 1 = 1 [ f ( y n ) f ( x n ) ] f [ y n , x n ] 1 [ f ( y n ) f ( x n ) ] f [ z n , x n ] .
Similar to Theorem 1, we can prove that the iterative method in Equations (11) and (12) has convergence order four with the following error equation:
e n + 1 = c 2 ( c 3 + f ( α ) 2 c 2 2 c 2 2 ) e n 4 + O ( e n 5 ) .

2.2. Optimal Three-Point Eighth Order Method

In [12], Neta and Petkovic proposed the following three-point optimal method based on inverse interpolation:
y n = x n f ( x n ) f ( x n ) , n 0 , w n = ψ 4 ( x n , y n ) , x n + 1 = y n + b 2 f ( x n ) 2 g 2 f ( x n ) 3 ,
where ψ 4 is a real function that should be chosen such that it provides the fourth order convergence of the sequence { x n } and it requires already computed values f ( x n ) , f ( x n ) and f ( y n ) . The iteration in Equation (9) is an example of such function. The values of g 2 and b 2 are given as:
g 2 = 1 [ f ( y n ) f ( x n ) ] [ f ( y n ) f ( w n ) ] f [ y n , x n ] 1 [ f ( w n ) f ( x n ) ] [ f ( y n ) f ( w n ) ] f [ w n , x n ] + 1 [ f ( w n ) f ( x n ) ] [ f ( y n ) f ( w n ) ] f ( x n ) 1 [ f ( y n ) f ( x n ) ] [ f ( y n ) f ( w n ) ] f ( x n ) ,
b 2 = 1 [ f ( y n ) f ( x n ) ] f [ y n , x n ] 1 f ( x n ) [ f ( y n ) f ( x n ) ] g 2 [ f ( y n ) f ( x n ) ] .
By using the approximation given by Equation (4) and the iterative scheme in Equations (11) and (12) at the second step of the three-point method in Equation (13), we obtain a new optimal derivative-free method as follows:
y n = x n f ( x n ) f [ z n , x n ] , z n = x n + f ( x n ) 3 , n 0 , w n = y n + h 2 f ( x n ) 2 , x n + 1 = y n + b 3 f ( x n ) 2 g 3 f ( x n ) 3 ,
where
h 2 = 1 [ f ( y n ) f ( x n ) ] f [ y n , x n ] 1 [ f ( y n ) f ( x n ) ] f [ z n , x n ] , g 3 = 1 [ f ( y n ) f ( x n ) ] [ f ( y n ) f ( w n ) ] f [ y n , x n ] 1 [ f ( w n ) f ( x n ) ] [ f ( y n ) f ( w n ) ] f [ w n , x n ] + 1 [ f ( w n ) f ( x n ) ] [ f ( y n ) f ( w n ) ] f [ z n , x n ] 1 [ f ( y n ) f ( x n ) ] [ f ( y n ) f ( w n ) ] f [ z n , x n ] , b 3 = 1 [ f ( y n ) f ( x n ) ] f [ y n , x n ] 1 f [ z n , x n ] [ f ( y n ) f ( x n ) ] g 3 [ f ( y n ) f ( x n ) ] .
Theorem 2.
Let α be a simple root of f, where f : I R R is a sufficiently differentiable function in an open interval I. For an initial approximation x 0 sufficiently close to α , the iterative methods given by Equations (16) and (17) is optimal eighth order convergent with the following error relation:
e n + 1 = c 2 2 ( 10 c 2 5 + 15 c 2 3 c 3 5 c 2 c 3 2 2 c 4 c 2 2 + 2 c 2 3 f ( α ) 3 + c 3 c 4 c 3 f ( α ) 3 c 2 ) e n 8 + O ( e n 9 ) ,
where c j = f ( j ) ( α ) j ! f ( α ) , j 2 and e n = x n α .
Proof. 
With the help of Taylor’s expansions, the proof is similar to the proof for Theorem 1 and those already taken in [6,8]. Hence, it is omitted. □

2.3. Optimal Four-Point Sixteenth Order Method

Neta and Petkovic [12] also presented an optimal four-point scheme with sixteenth order convergence using inverse interpolation requiring an evaluation of first derivative f ( x n ) at each step and four evaluations of the function.
We transform their four-point optimal scheme and use the scheme in Equations (11) and (12) at the second step to develop a new derivative-free optimal four-point method as follows:
y n = x n f ( x n ) f [ z n , x n ] , z n = x n + f ( x n ) 4 , n 0 , w n = y n + h 2 f ( x n ) 2 , t n = y n + b 3 f ( x n ) 2 g 3 f ( x n ) 3 , x n + 1 = y n + b 4 f ( x n ) 2 g 4 f ( x n ) 3 + g 5 f ( x n ) 4 ,
where h 2 , g 3 and b 3 are given in Equation (17) and the values of g 5 , g 4 and b 4 are given as:
g 5 = φ t φ w [ f ( t n ) f ( w n ) ] φ y φ w [ f ( y n ) f ( w n ) ] [ f ( t n ) f ( y n ) ] , g 4 = φ t φ w [ f ( t n ) f ( w n ) ] g 5 ( [ f ( t n ) f ( x n ) ] + [ f ( w n ) f ( x n ) ] ) , b 4 = φ t g 4 [ f ( t n ) f ( x n ) ] g 5 [ f ( t n ) f ( x n ) ] 2
where
φ t = 1 f [ t n , x n ] [ f ( t n ) f ( x n ) ] 1 f [ z n , x n ] [ f ( t n ) f ( x n ) ] , φ w = 1 f [ w n , x n ] [ f ( w n ) f ( x n ) ] 1 f [ z n , x n ] [ f ( w n ) f ( x n ) ] , φ y = 1 f [ y n , x n ] [ f ( y n ) f ( x n ) ] 1 f [ z n , x n ] [ f ( y n ) f ( x n ) ] .
Theorem 3.
Let α be a simple root of f, where f : I R R is a sufficiently differentiable function in an open interval I. For an initial approximation x 0 close enough to α , the iterative method given by Equations (18) and (19) is optimal sixteenth order convergent with the following error equation:
e n + 1 = c 2 4 ( c 4 c 2 f ( α ) 4 c 3 2 4 c 4 c 3 f ( α ) 4 c 2 3 4 c 5 c 2 2 c 4 c 3 3 c 4 c 3 4 + 20 c 2 7 c 5 24 c 2 5 c 4 2 + 15 c 3 5 c 2 180 c 3 4 c 2 3 + 715 c 3 3 c 2 5 + 20 f ( α ) 4 c 2 8 176 c 2 8 c 4 1250 c 2 7 c 3 2 + 980 c 3 c 2 9 280 c 2 11 5 c 2 2 f ( α ) 4 c 3 3 + 24 c 2 3 c 4 2 c 3 6 c 4 2 c 3 2 c 2 + 4 c 2 4 c 5 c 4 40 c 3 f ( α ) 4 c 2 6 5 c 5 c 3 3 c 2 + 25 c 3 2 c 2 3 c 5 + c 3 2 c 5 c 4 + 25 c 3 2 f ( α ) 4 c 2 4 260 c 4 c 3 2 c 2 4 40 c 3 c 2 5 c 5 + 63 c 4 c 3 3 c 2 2 + 380 c 2 6 c 4 c 3 + 4 c 4 f ( α ) 4 c 2 5 ) e n 16 + O ( e n 17 ) ,
where c j = f ( j ) ( α ) j ! f ( α ) , j 2 and e n = x n α .
Proof. 
The proof is similar to those already taken in [6,8] by using the Taylor’s expansions. □

2.4. n-Point Method of Optimal Order 2 n

In [2], Petkovic et al. proposed an n-point n-step iterative scheme of order 2 n that requires one evaluation of first derivative f ( x n ) and n evaluations of function at each step. We develop the following n-point n-step derivative-free method by using the approximation given in Equation (4):
ϕ 1 ( x k ) = x k f ( x k ) f [ z k , x k ] , z k = x k + f ( x k ) m , m n 0 , ϕ 2 ( x k ) = R 2 ( 0 ) , x k + 1 = ϕ n ( x n ) = R n ( 0 ) .
where R n is an inverse interpolating polynomial of degree n given as:
x = R n ( f ( x ) ) = d 0 + d 1 [ f ( x ) f ( x k ) ] + d 2 [ f ( x ) f ( x k ) ] 2 + + d n [ f ( x ) f ( x k ) ] n ,
with the conditions:
R n ( f ( x k ) ) = x k , R n ( f ( x k ) ) = 1 f [ z n , x n ] , R n ( f ( ϕ 1 ) ) = f ( ϕ 1 ) , . . . , R n ( f ( ϕ n ) ) = f ( ϕ n ) .
By the use of the conditions in Equation (23), the coefficients d 0 , d 1 , . . . , d n can be determined easily and hence we obtain an n-point derivative-free family of following form:
ϕ 1 ( x k ) = x k f ( x k ) f [ z k , x k ] , z k = x k + f ( x k ) m , m n 0 , ϕ 2 ( x k ) = R 2 ( 0 ) , x k + 1 = ϕ n ( x k ) = R n ( 0 ) = ϕ 1 ( x k ) + d 2 f ( x k ) 2 + + ( 1 ) n d n f ( x k ) n .
Theorem 4.
Let α be a simple root of f, where f : I R R is a sufficiently differentiable function in an open interval I. For an initial approximation x 0 sufficiently close to α , the n-point iterative method defined by Equation (24) is of optimal order 2 n .

3. Review of Some Four-Point Optimal Methods and Numerical Results

Here, we use the transformation in Equation (4) to modify optimal four-point iterative scheme of Sharifi et al. [13] denoted by SL16 and optimal four-point method by Geum and Kim [14] denoted by GK16. Table 1 presents the original and modified forms of four-points methods where
K 1 ( u n ) = 1 + ρ u n + ( 9 + 5 ρ / 2 ) u n 2 1 + ( ρ 2 ) u n + ( 4 + ρ / 2 ) u n 2 , K 2 ( u n , v n , w n ) = 1 + 2 u n + ( 2 + ξ ) w n 1 v n + ξ w n , K 3 ( u n , v n , w n , t n ) = 1 + 2 u n + ( 2 + ξ ) v n w n 1 v n 2 w n t n + 2 ( 1 + ξ ) v n w n 1 2 u n w n [ 6 + 12 u n + u n 2 ( 24 11 ρ ) + u n 3 ( 11 ρ 2 66 ρ + 136 ) + 4 ξ ] + ( 2 u n ( ξ 2 2 ξ 9 ) 4 ξ 6 ) w n 2
and
L 1 ( u n ) = 1 + 2 u n + 5 u n 2 6 u n 3 , L 2 ( u n , v n , w n ) = 1 + 2 u n + 4 w n + 6 u n 2 + v n L 3 ( u n , v n , w n , p n , q n , t n ) = 1 + 6 u n 2 + 2 u n v n 3 + v n + 4 w n 4 w n 2 + u n w n + 6 u n 2 w n + 2 u n 3 w n 10 u n w n 2 + t n + 2 q n + 8 p n + 2 u n t n + 2 v n w n + 6 u n 2 t n 4 v n 2 w n + 24 u n 4 w n ,
are weight functions,
u n = f ( y n ) f ( x n ) , v n = f ( r n ) f ( y n ) , w n = f ( r n ) f ( x n ) , t n = f ( s n ) f ( r n ) , p n = f ( s n ) f ( x n ) , q n = f ( s n ) f ( y n )
and ρ , ξ are free parameters. Here, we have chosen ρ = 2 , ξ = 2 .
A family of optimal derivative-free iterative methods of arbitrary high order by using polynomial interpolation presented by Cordero et al. [7] is given as follows:
y 0 = x k , k 0 , y 1 = y 0 + f ( y 0 ) , x k + 1 = y k + 1 = y j f ( y j ) p j ( y j ) , j = 1 , 2 , . . . , n ,
where p j is the polynomial that interpolates f in y 0 , y 1 , , y j .
Now, we test all the discussed optimal with- and without-derivative methods using different types of nonlinear equations. We employed multi-precision arithmetic with 4000 significant decimal digits in the programming package of Maple 16 (Waterloo Maple Inc., Waterloo, ON, Canada) to obtain a high accuracy and to avoid the loss of significant digits.
In chemistry, several nonlinear systems can be found, for example in the investigation of stability of chemical reactions. Here, our concern is to deal with the solution of Van der Waals equation for finding the volume of a gas. Van der Waals equation is given by
( P + a n 2 V 2 ) ( V n b ) = n R T ,
which relates the pressure (P), volume (V) and temperature (T) of a gas. Here, n represents the number of moles of gas present and R is the universal gas constant. The term involving the parameter a corrects the pressure for intermolecular attractive forces, while the term involving b is a correction for that portion of the volume of gas that is not compressible due to the intrinsic volume of the gas molecules. We have to find the solution of a nonlinear equation in V to determine the volume V of the gas in terms of the remaining parameters as follows:
P V 3 ( n b P + n R T ) V 2 + a n 2 V a n 2 b = 0 .
Suppose that one mole of chlorine gas has a pressure of 2 atmospheres and a temperature of 313 K. For chlorine gas, a = 6.29 atm·L 2 /mol 2 and b = 0.0562 L/mol. The universal gas constant has the value R = 0.08206 atm·L/mol·K. Hence, we obtain following equation, which is cubic in V:
f 1 ( V ) = 2 V 3 25.79718 V 2 + 6.29 V 0.353498 ,
with a root α = 0.087113811 .
We compare the convergence behavior of the modified methods MGK16, MSL16 and Equations (18) and (19) denoted by MNP16 with their with derivative versions, i.e., GK16, SL16 and (NP16) [12], respectively, and the family (Equation (26)) denoted by (CT16) for n = 4 , using f 1 ( V ) and the nonlinear functions given in Table 2. The error between approximated and exact root ( x n α ) and the computational order of convergence (coc) for first three iterations of various methods is displayed in Table 3, Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9, where E ( i ) denotes E × 10 i , which means that the absolute error for the corespondent method is zero up to i 1 decimal places. It is observed that the computational order of convergence (coc) supports the theoretical order of convergence. The formula to compute the computational order of convergence (coc) is given by [15]:
coc log ( x n + 1 α ) / ( x n α ) log ( x n α ) / ( x n 1 α ) .
It was observed from the numerical experiments that the modified four-point methods MGK16, MSL16 and MNP16 are comparable and competitive to the methods GK16, SL16, NP16 and CT16. Especially, for the case of f 7 , all the with-derivative four-point methods failed to converge when the initial guess was taken far from the required root, while the modified derivative-free four-points methods provided remarkably fast convergence. We see that schemes GK16 and CT16 failed to converge to the required root for the case of f 4 .

4. Conclusions

In this paper, we have developed a family of optimal derivative-free root finding methods of arbitrary high order based on inverse interpolation by applying the conjecture of Cordero and Torregrosa. Some existing derivative based methods are modified using this conjecture. Convergence analysis is studied for the proposed optimal methods. Finally, numerical tests are provided that support the theoretical results. It was observed that the modified derivative free methods can compete and work better than their with-derivative versions. Especially, the proposed derivative-free methods provided remarkably fast convergence for the case of f 7 , even when the initial guess was taken far from the required root, while the derivative based methods failed to converge.

Author Contributions

Methodology, M.-u.-D.J.; writing original draft preparation, M.-u.-D.J.; investigation, M.-u.-D.J.; writing review and editing, F.Z., N.Y.; supervision, F.Z., N.Y.

Funding

This research received no external funding.

Acknowledgments

The authors are thankful to the respected reviewers for their useful comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ortega, J.M.; Rheiboldt, W.C. Iterative Solution of Nonlinear Equations in Several Variables; Academic Press: New York, NY, USA, 1970. [Google Scholar]
  2. Petković, M.S.; Neta, B.; Petković, L.D.; Džunić, J. Multipoint Methods for Solving Nonlinear Equations; Elsevier: Amsterdam, The Netherlands, 2013. [Google Scholar]
  3. Traub, J.F. Iterative Methods for the Solution of Equations; Prentice-Hall: Englewood Cliffs, NJ, USA, 1964. [Google Scholar]
  4. Kung, H.T.; Traub, J.F. Optimal order of one-point and multipoint iteration. J. Assoc. Comput. Math. 1974, 235, 643–651. [Google Scholar] [CrossRef]
  5. Behl, R.; Argyros, I.K.; Motsa, S.S. A new highly efficient and optimal family of eighth order methods for solving nonlinear equations. Appl. Math. Comput. 2016, 282, 175–186. [Google Scholar] [CrossRef]
  6. Cordero, A.; Torregrosa, J.R. Low-complexity root-finding iteration functions with no derivatives of any order of convergence. J. Comput. Appl. Math. 2015, 275, 502–515. [Google Scholar] [CrossRef] [Green Version]
  7. Cordero, A.; Hueso, J.L.; Martinez, E.; Torregrosa, J.R. Generating optimal derivative-free iterative methods for nonlinear equations by using polynomial interpolation. Math. Comput. Model. 2013, 57, 1950–1956. [Google Scholar] [CrossRef]
  8. Zafar, F.; Yasmin, N.; Akram, S.; Junjua, M.D. A general class of derivative-free optimal root finding methods based on rational interpolation. Sci. World J. 2015, 2015, 934260. [Google Scholar] [CrossRef] [PubMed]
  9. Zheng, Q.; Li, J.; Huang, F. An optimal Steffensen-type family for solving nonlinear equations. Appl. Math. Comput. 2011, 217, 9592–9597. [Google Scholar] [CrossRef]
  10. Steffensen, I.F. Remarks on iteration. Skand. Actuar. J. 1933, 16, 64–72. [Google Scholar] [CrossRef]
  11. Ostrowski, A.M. Solution of Equations and Systems of Equations; Academic Press: New York, NY, USA, 1960. [Google Scholar]
  12. Neta, B.; Petković, M.S. Construction of optimal order nonlinear solvers using inverse interpolation. Appl. Math. Comput. 2010, 217, 2448–2455. [Google Scholar] [CrossRef] [Green Version]
  13. Sharifi, S.; Salimi, M.; Siegmund, S.; Lotfi, T. A new class of optimal four-point methods with convergence order 16 for solving nonlinear equations. Math. Comput. Simul. 2016, 119, 69–90. [Google Scholar] [CrossRef] [Green Version]
  14. Geum, Y.H.; Kim, Y.I. A biparametric family of optimally convergent sixteenth order multipoint methods with their fourth-step weighting function as a sum of a rational and a generic two-variable function. J. Comput. Appl. Math. 2011, 235, 3178–3188. [Google Scholar] [CrossRef]
  15. 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]
Table 1. Four-point methods and their modifications.
Table 1. Four-point methods and their modifications.
Original Iterative MethodModified Iterative Method
GK16:MGK16:
y n = x n f ( x n ) f ( x n ) , n 0 , y n = x n f ( x n ) f [ z n , x n ] , z n = x n + f ( x n ) 4 , n 0 ,
r n = y n K 1 ( u n ) f ( y n ) f ( x n ) , r n = y n K 1 ( u n ) f ( y n ) f [ z n , x n ] ,
s n = r n K 2 ( u n , v n , w n ) f ( r n ) f ( x n ) , s n = r n K 2 ( u n , v n , w n ) f ( r n ) f [ z n , x n ] ,
x n + 1 = s n K 3 ( u n , v n , w n , t n ) f ( s n ) f ( x n ) . x n + 1 = s n K 3 ( u n , v n , w n , t n ) f ( s n ) f [ z n , x n ] .
SL16:MSL16:
y n = x n f ( x n ) f ( x n ) , n 0 , y n = x n f ( x n ) f [ z n , x n ] , z n = x n + f ( x n ) 4 , n 0 ,
r n = y n L 1 ( u n ) f ( y n ) f ( x n ) , r n = y n L 1 ( u n ) f ( y n ) f [ z n , x n ] ,
s n = r n L 2 ( u n , v n , w n ) f ( r n ) f ( x n ) , s n = r n L 2 ( u n , v n , w n ) f ( r n ) f [ z n , x n ] ,
x n + 1 = s n L 3 ( u n , v n , w n , p n , q n , t n ) f ( s n ) f ( x n ) . x n + 1 = s n L 3 ( u n , v n , w n , p n , q n , t n ) f ( s n ) f [ z n , x n ] .
Table 2. Test functions.
Table 2. Test functions.
Test FunctionsExact Root α x 0
f 2 ( x ) = ( 2 + x 3 ) cos ( π x 2 ) + log ( x 2 + 2 x + 2 ) 1 0.93
f 3 ( x ) = x 2 e x + x cos 1 x 3 + 1 1.5650602 1.25
f 4 ( x ) = x e x + log ( 1 + x + x 4 ) 0 0.5
f 5 ( x ) = ( x 1 ) ( x + 1 + log ( 2 + x + x 2 ) ) 1 1.05
f 6 ( x ) = 20 x 5 x 2 + 1 2 0.42767729 0.38
f 7 ( x ) = e sin ( 8 x ) 4 x 0.34985721 7
Table 3. Numerical results for f 1 , x 0 = 0 .
Table 3. Numerical results for f 1 , x 0 = 0 .
ErrorGK16SL16NP16CT16MGK16MSL16MNP16
x 1 α 2.81 ( 5 ) 3.02 ( 4 ) 5.28 ( 5 ) 1.73 ( 4 ) 9.73 ( 6 ) 2.87 ( 4 ) 4.36 ( 5 )
x 2 α 6.12 ( 54 ) 4.70 ( 38 ) 1.17 ( 49 ) 2.96 ( 40 ) 2.44 ( 61 ) 2.05 ( 38 ) 5.52 ( 51 )
x 3 α 1.42 ( 832 ) 6.75 ( 580 ) 4.03 ( 764 ) 1.68 ( 612 ) 5.89 ( 951 ) 1.19 ( 585 ) 2.36 ( 785 )
coc16.0016.0016.0016.0016.0016.0016.00
Table 4. Numerical results for f 2 , x 0 = 0.93 .
Table 4. Numerical results for f 2 , x 0 = 0.93 .
ErrorGK16SL16NP16CT16MGK16MSL16MNP16
x 1 α 1.44 ( 10 ) 2.63 ( 9 ) 1.83 ( 10 ) 6.42 ( 10 ) 1.65 ( 10 ) 2.59 ( 9 ) 1.83 ( 10 )
x 2 α 8.08 ( 147 ) 1.61 ( 127 ) 2.77 ( 145 ) 9.99 ( 136 ) 6.82 ( 147 ) 1.23 ( 127 ) 2.58 ( 145 )
x 3 α 7.29 ( 2327 ) 6.48 ( 2019 ) 1.98 ( 2302 ) 1.18 ( 2148 ) 4.69 ( 2327 ) 8.22 ( 2021 ) 6.18 ( 2303 )
coc16.0016.0016.0016.0016.0016.0016.00
Table 5. Numerical results for f 3 , x 0 = 1.25 .
Table 5. Numerical results for f 3 , x 0 = 1.25 .
ErrorGK16SL16NP16CT16MGK16MSL16MNP16
x 1 α 5.84 ( 13 ) 7.09 ( 9 ) 2.49 ( 11 ) 4.10 ( 6 ) 1.37 ( 11 ) 1.57 ( 8 ) 1.46 ( 11 )
x 2 α 2.58 ( 201 ) 7.05 ( 135 ) 5.79 ( 176 ) 2.16 ( 89 ) 2.07 ( 179 ) 2.53 ( 129 ) 4.07 ( 180 )
x 3 α 5.49 ( 3215 ) 6.51 ( 2151 ) 4.04 ( 2810 ) 7.88 ( 1422 ) 1.46 ( 2864 ) 5.25 ( 2062 ) 5.01 ( 2877 )
coc16.0016.0016.0016.0016.0016.0016.00
Table 6. Numerical results for f 4 , x 0 = 0.5 .
Table 6. Numerical results for f 4 , x 0 = 0.5 .
ErrorGK16SL16NP16CT16MGK16MSL16MNP16
x 1 α Diverges 2.9732 1.42 ( 4 ) Diverges 1.20 ( 9 ) 1.51 ( 7 ) 6.69 ( 10 )
x 2 α Diverges 7.74 ( 1 ) 1.65 ( 67 ) Diverges 3.50 ( 149 ) 1.24 ( 112 ) 2.43 ( 152 )
x 3 α Diverges 4.59 ( 6 ) 1.85 ( 1074 ) Diverges 9.05 ( 2382 ) 4.90 ( 1794 ) 2.35 ( 2431 )
cocDiverges8.9416.00Diverges16.0016.0016.00
Table 7. Numerical results for f 5 , x 0 = 1.05 .
Table 7. Numerical results for f 5 , x 0 = 1.05 .
ErrorGK16SL16NP16CT16MGK16MSL16MNP16
x 1 α 1.02 ( 23 ) 6.99 ( 25 ) 1.29 ( 23 ) 3.40 ( 21 ) 4.19 ( 20 ) 8.96 ( 23 ) 2.69 ( 21 )
x 2 α 2.44 ( 370 ) 8.60 ( 390 ) 8.99 ( 369 ) 1.83 ( 327 ) 4.78 ( 310 ) 3.43 ( 356 ) 7.83 ( 330 )
x 3 α 00 1 ( 3999 ) 0000
coc16.0016.0016.0016.0016.0016.0016.00
Table 8. Numerical results for f 6 , x 0 = 0.38 .
Table 8. Numerical results for f 6 , x 0 = 0.38 .
ErrorGK16SL16NP16CT16MGK16MSL16MNP16
x 1 α 1.65 ( 5 ) 1.28 ( 3 ) 8.41 ( 11 ) 1.70 ( 11 ) 1.10 ( 5 ) 1.04 ( 3 ) 2.78 ( 11 )
x 2 α 2.04 ( 66 ) 2.68 ( 34 ) 1.59 ( 151 ) 1.86 ( 161 ) 8.28 ( 69 ) 1.04 ( 35 ) 5.53 ( 160 )
x 3 α 5.69 ( 1041 ) 5.39 ( 525 ) 4.62 ( 2403 ) 7.53 ( 2561 ) 8.66 ( 1079 ) 1.60 ( 547 ) 3.22 ( 2539 )
coc16.0016.0016.0016.0016.0016.0016.00
Table 9. Numerical results for f 7 , x 0 = 7 .
Table 9. Numerical results for f 7 , x 0 = 7 .
ErrorGK16SL16NP16CT16MGK16MSL16MNP16
x 1 α DivergesDivergesDiverges 1.98 ( 2 ) 4.92 ( 3 ) 3.00 ( 2 ) 1.50 ( 2 )
x 2 α DivergesDivergesDiverges 3.89 ( 12 ) 3.72 ( 26 ) 5.71 ( 11 ) 3.31 ( 17 )
x 3 α DivergesDivergesDiverges 1.20 ( 168 ) 5.72 ( 398 ) 8.21 ( 155 ) 9.46 ( 225 )
cocDivergesDivergesDiverges16.1216.0816.4916.20

Share and Cite

MDPI and ACS Style

Junjua, M.-u.-D.; Zafar, F.; Yasmin, N. Optimal Derivative-Free Root Finding Methods Based on Inverse Interpolation. Mathematics 2019, 7, 164. https://doi.org/10.3390/math7020164

AMA Style

Junjua M-u-D, Zafar F, Yasmin N. Optimal Derivative-Free Root Finding Methods Based on Inverse Interpolation. Mathematics. 2019; 7(2):164. https://doi.org/10.3390/math7020164

Chicago/Turabian Style

Junjua, Moin-ud-Din, Fiza Zafar, and Nusrat Yasmin. 2019. "Optimal Derivative-Free Root Finding Methods Based on Inverse Interpolation" Mathematics 7, no. 2: 164. https://doi.org/10.3390/math7020164

APA Style

Junjua, M. -u. -D., Zafar, F., & Yasmin, N. (2019). Optimal Derivative-Free Root Finding Methods Based on Inverse Interpolation. Mathematics, 7(2), 164. https://doi.org/10.3390/math7020164

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