Next Article in Journal
Combining Optimization and Simulation for Designing a Robust Short-Sea Feeder Network
Next Article in Special Issue
Improved Sliding Mode Finite-Time Synchronization of Chaotic Systems with Unknown Parameters
Previous Article in Journal
Graphs from Features: Tree-Based Graph Layout for Feature Analysis
Previous Article in Special Issue
On the Solutions of Second-Order Differential Equations with Polynomial Coefficients: Theory, Algorithm, Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a Family of Jarratt-Like Sixth-Order Iterative Methods for Solving Nonlinear Systems with Their Basins of Attraction

Department of Mathematics, Dankook University, Cheonan 330-714, Korea
*
Author to whom correspondence should be addressed.
Algorithms 2020, 13(11), 303; https://doi.org/10.3390/a13110303
Submission received: 31 October 2020 / Revised: 16 November 2020 / Accepted: 17 November 2020 / Published: 20 November 2020

Abstract

:
We develop a family of three-step sixth order methods with generic weight functions employed in the second and third sub-steps for solving nonlinear systems. Theoretical and computational studies are of major concern for the convergence behavior with applications to special cases of rational weight functions. A number of numerical examples are illustrated to confirm the convergence behavior of local as well as global character of the proposed and existing methods viewed through the basins of attraction.
MSC:
65H05; 65H99; 41A25; 65B99

1. Introduction

Since exact solutions for nonlinear equations are rarely available, we usually resort to their numerical solutions. To locate the desired numerical roots, many authors [1,2,3,4,5,6,7,8,9] have developed high-order iterative methods including optimal eighth-order ones [10,11,12,13,14,15].
This paper is devoted to devise a class of sixth-order iterative root-finders for nonlinear systems by employing a three-step weighted Jarratt-like method below:
y n = x n γ · f ( x n ) 1 f ( x n ) , γ R , z n = x n T f ( s ) · f ( x n ) 1 f ( x n ) , x n + 1 = z n L f ( s ) · f ( x n ) 1 f ( z n ) ,
where s = f ( x n ) 1 f ( y n ) , γ is a parameter to be determined later and T f , L f : C C are weight functions being analytic [16,17,18] in a neighborhood of 1. Note that Scheme (1) uses two functional values as well as two derivatives. We are certainly able to introduce generic weight functions using one derivative and three functional values to develop general optimal eighth-order methods that covers the existing ones for the zero of a given scalar function. However, expanding such approach to a nonlinear system requires different weight functions. For unified analysis to be performed in both scalar and vector functions, we aim to develop a family of Jarratt-like sixth-order iterative methods by maintaining the same form of weight functions with two derivatives as well as two functional values. This extension to nonlinear systems is the main strength of this paper.
The robustness of the current analysis presented here covers most existing studies on higher-order root-finders using two derivatives and two function values for both scalar and vector equations. The results of Theorem 1 give us not only fairly generic scalar function solvers, but also some advantage of extending to a nonlinear system with any finite dimension. Such an extension is evidently characterized by Theorem 2 to be studied in this analysis.
Our major aim is not only to design a class of sixth-order methods by fully specifying the algebraic structure of generic weight functions T f ( s ) and L f ( s ) , but also to investigate their basins of attraction behind the extraneous fixed points [19] when applied to polynomials. The last sub-step of (1) in the form of weighted Newton’s method is clearly more convenient in dealing with extraneous fixed points which are the roots of the weight function T f ( s ) + L f ( s ) · f ( z ) f ( x ) .
The extraneous fixed points may lead us to attractive, indifferent, repulsive and chaotic orbits via the related basins of attraction.
Section 2 investigates the main theorem regarding the convergence behavior with the desired forms of weight functions, while Section 3 deals with special cases of weight functions that can cover many of the existing studies using two derivatives and two functional evaluations. Section 4 discusses the computational and long-term orbit behavior of the proposed iterative methods regarding scalar functions. Section 5 presents numerical experiments in a d-dimensional Euclidean space by solving a system of nonlinear vector equations f : R d R d encountered in a real life with d { 3 , 4 , 9 , 10 } . In addition, computational efficiency is addressed with issues related to the accuracy and applicability of the proposed methods. Concluding remarks are stated in Section 6.

2. Main Theorem

The main theorem for a nonlinear scalar equation will be pursued and extended later in Section 5 to a system of nonlinear vector equations:
Theorem 1.
Assume that f : D C C has a simple root α and is analytic in a region D containing α. Let θ j = f ( j ) ( α ) j ! f ( α ) for j = 2 , 3 , . Let x 0 be an initial guess chosen in a sufficiently small neighborhood of α. Let T f , L f : C C be analytic in a neighborhood of 1. Let T i = 1 i ! d i d s i T f ( s ) | s = 1 and L i = 1 i ! d i d s i L f ( s ) | s = 1 for 0 i 5 . If T 0 = 1 , T 1 = 1 2 γ , L 0 = 1 , L 1 = 1 γ , then iterative scheme (1) defines a family of fifth-order methods for γ R { 0 } . If we add further constraints either with { γ = 2 3 , T 2 = 9 8 , | T 3 | < , | L 2 | < } or with { γ = 1 , L 2 = 3 2 , | L 3 | < , | T 2 | < } , then iterative scheme (1) reduces to a family of sixth-order methods satisfying the error equation below. For n = 0 , 1 , 2 , ,
e n + 1 = 1 243 2 ( 27 + 8 L 2 ) θ 2 2 + 9 θ 3 · ( 135 + 64 T 3 ) θ 2 3 27 θ 2 θ 3 + 3 θ 4 e n 6 + O ( e n 7 ) , i f γ = 2 3 , 4 ( 1 2 T 2 ) ( 9 + 4 L 3 + 2 T 2 ) θ 2 5 + 4 ( 1 + L 3 + 3 T 2 ) θ 2 3 θ 3 5 4 θ 2 θ 3 2 e n 6 + O ( e n 7 ) , i f γ = 1 .
Proof. 
Taylor series expansion of f ( x n ) about α up to sixth-order terms with f ( α ) = 0 leads us to the following:
f ( x n ) = f ( α ) [ e n + θ 2 e n 2 + θ 3 e n 3 + θ 4 e n 4 + θ 5 e n 5 + θ 6 e n 6 + O ( e n 7 ) ] .
It follows that
f ( x n ) = f ( α ) [ 1 + 2 θ 2 e n + 3 θ 3 e n 2 + 4 θ 4 e n 3 + 5 θ 5 e n 4 + 6 θ 6 e n 5 + O ( e n 6 ) ] .
For brevity of notation, we denote e n by e, unless otherwise specified from now on. Symbolic computation of Mathematica [20] yields:
y n = x n γ f ( x n ) f ( x n ) = α + ( 1 γ ) e + γ θ 2 e 2 2 γ ( θ 2 2 θ 3 ) e 3 + Y 4 e 4 + Y 5 e 5 + Y 6 e 6 + O ( e 7 ) ,
where Y 4 = γ ( 4 θ 2 3 7 θ 2 θ 3 + 3 θ 4 ) , Y 5 = 2 γ ( 4 θ 2 4 10 θ 2 2 θ 3 + 3 θ 3 2 + 5 θ 2 θ 4 2 θ 5 ) , Y 6 = γ ( 16 θ 2 5 52 θ 2 3 θ 3 + 33 θ 2 θ 3 2 + 28 θ 2 2 θ 4 17 θ 3 θ 4 13 θ 2 θ 5 + 5 θ 6 ) .
In view of the fact that f ( y n ) = f ( x n ) | e n ( y n α ) , we get:
f ( y n ) = f ( α ) [ 1 2 ( γ 1 ) θ 2 e + [ 2 γ θ 2 2 + 3 ( γ 1 ) 2 θ 3 ] e 2 + Σ i = 3 5 D i e i + O ( e 6 ) ] ,
where D i = D i ( γ , θ 2 , θ 3 , , θ 6 ) for 3 i 5 . Hence, we have:
s = f ( y n ) f ( x n ) = 1 2 γ θ 2 e + 3 γ [ 2 θ 2 2 + ( 2 + γ ) θ 3 ] ) e 2 + Σ i = 3 5 E i e i + O ( e 6 ) ,
where E i = E i ( γ , θ 2 , θ 3 , , θ 6 ) for 3 i 5 .
Noting that s = 1 + O ( e ) and f ( x n ) f ( x n ) = O ( e ) , we need a Taylor expansion of T f ( s ) about s = 1 up to fifth-order terms:
T f ( s ) = T 0 + T 1 ( s 1 ) + T 2 ( s 1 ) 2 + T 3 ( s 1 ) 3 + T 4 ( s 1 ) 4 + T 5 ( s 1 ) 5 + O ( e 6 ) .
Thus, we find
z n = x n T f ( s ) · f ( x n ) f ( x n ) = α + ( 1 T 0 ) e + ( T 0 + 2 T 1 γ ) θ 2 e 2 +
[ 2 T 0 ( θ 2 2 θ 3 ) γ ( 8 T 1 θ 2 2 + 4 T 2 γ θ 2 2 + 3 T 1 ( 2 + γ ) θ 3 ) ] e 3 + Σ i = 4 6 W i e i + O ( e 7 ) ,
where W i = W i ( γ , θ 2 , θ 3 , , θ 6 , T 0 , , T 5 ) for 4 i 6 .
In view of the fact that f ( z n ) = f ( x n ) | e n ( z n α ) , we get:
f ( z n ) = f ( α ) [ ( 1 T 0 ) e + [ ( 1 T 0 ) 2 θ 2 + ( T 0 + 2 T 1 γ ) θ 2 ] e 2 + Σ i = 3 6 F i e i + O ( e 7 ) ] ,
where F i = F i ( γ , θ 2 , θ 3 , , θ 6 , T 0 , , T 5 ) for 3 i 6 . Noting that s = O ( 1 ) and f ( z n ) f ( x n ) = O ( e ) , we need a Taylor expansion of L f ( s ) about s = 1 up to fifth-order terms:
L f ( s ) = L 0 + L 1 ( s 1 ) + L 2 ( s 1 ) 2 + L 3 ( s 1 ) 3 + L 4 ( s 1 ) 4 + L 5 ( s 1 ) 5 + O ( e 6 ) .
Hence, we have:
x n + 1 = z n L f ( s ) · f ( z n ) f ( x n ) = α + ( L 0 1 ) ( T 0 1 ) e
+ [ ( T 0 2 L 1 T 0 γ + 2 ( L 1 + T 1 ) γ L 0 ( 1 + T 0 + T 0 2 + 2 T 1 γ ) ] θ 2 e 2 + Σ i = 3 6 G i e i + O ( e 7 ) ,
where G i = G i ( γ , θ 2 , θ 3 , , θ 6 , T 0 , , T 5 , L 0 , , L 5 ) for 3 i 6 .
By taking
T 0 = 1 , T 1 = 1 2 γ , L 0 = 1 , L 1 = 1 γ ,
we further obtain
x n + 1 = α + [ 4 ( 3 + 2 L 2 γ 2 ) ( 1 + 2 T 2 γ 2 ) θ 2 4 + ( 12 + 15 γ + 4 ( L 2 + 3 T 2 ) γ 2 6 ( L 2 + 2 T 2 ) γ 3 ) θ 2 2 θ 3
+ 3 2 ( γ 1 ) ( 3 γ 2 ) θ 3 2 ] e 5 + G 6 e 6 + O ( e 7 ) .
From ( γ 1 ) ( 3 γ 2 ) = 0 , we find two sets of relations:
{ γ = 2 3 , T 2 = 9 8 } , { γ = 1 , L 2 = 3 2 }
for vanishing the fifth-order term in (14). Hence, we eventually have two sets of relations:
{ γ = 2 3 , T 0 = 1 , T 1 = 3 4 , T 2 = 9 8 , L 0 = 1 , L 1 = 3 2 } ,
{ γ = 1 , T 0 = 1 , T 1 = 1 2 , L 0 = 1 , L 1 = 1 , L 2 = 3 2 } ,
which are substituted into G 6 in (14) and lead us to the desired relation (2) with { | T 3 | < , | L 2 | < } when γ = 2 3 , or with { | L 3 | < , | T 2 | < } when γ = 1 . □
Remark 1.
The fifth-order expansion of the weight functions is considered due to the fact that s = 1 + O ( e ) and f ( z n ) f ( x n ) = O ( e ) or f ( x n ) f ( x n ) = O ( e ) . However, the result of Theorem 1 shows that z = O ( e k ) , f ( z n ) f ( x n ) = O ( e k ) , f o r k { 3 , 4 } with γ { 1 , 2 3 } . Hence, we require only at most the third-order expansion for both weight functions to achieve the desired sixth-order convergence. This favorable fact is used below to establish the corresponding theorem for a family of sixth-order systems of nonlinear equations.

3. Special Cases of Weight Functions

Theorem 1 clearly covers the existing case study in R shown in [21]
γ = 2 3 , T f ( s ) = 1 + 1 s 2 ( β λ + λ s ) H ( 1 s ) , λ R , β R { 0 } L f ( s ) = 1 2 ( 5 3 s ) ,
where H ( t ) = 1 + 3 β 2 4 β t + 9 β 2 3 β 4 λ + 2 8 β 2 t 2 . Note that T f in (15) is a cubic-order rational weight function. It also covers the existing case study in R shown in [2]
γ = 2 3 , T f ( s ) = 1 + 3 ( 1 s ) 4 1 + 3 β ( 1 s ) 4 3 α ( 1 s ) , α , β R , L f ( s ) = 1 2 ( 5 3 s ) ,
where T f in (15) is a second-order rational weight function.
Although the result of Theorem 1 allows an infinite number of other forms of weight functions T f ( s ) and L f ( s ) , we are specifically interested only in second-order rational forms for the both weight functions T f ( s ) and L f ( s ) .
T f ( s ) = a i + b i s + g i s 2 c i + d i s + h i s 2 , L f ( s ) = p i + q i s + k i s 2 r i + σ i s + τ i s 2 , for   1 i 2 ,
where the desired coefficients are to be determined based on the results of Theorem 1.
Two cases can be considered. The first case can be given as follows:
  • Case A: { γ = 2 3 , T 0 = 1 , T 1 = 3 4 , L 0 = 1 , L 1 = 3 2 , T 2 = 9 8 , | T 3 | < , | L 2 | < }
    b 1 = 1 3 ( 9 a 1 31 g 1 + 16 h 1 ) , c 1 = 2 a 1 + 6 g 1 3 h 1 , d 1 = 2 3 ( 9 a 1 23 g 1 + 11 h 1 ) , T 3 = 9 ( 9 a 1 23 g 1 + 13 h 1 ) 16 ( 3 a 1 7 g 1 + 4 h 1 ) , p 1 = 1 3 ( 7 k 1 5 q 1 + 2 σ 1 + 4 τ 1 ) , r 1 = 1 3 ( 4 k 1 2 q 1 σ 1 + τ 1 ) , L 2 = 3 ( 2 k 1 + 3 σ 1 + 4 τ 1 ) 4 ( 2 k 1 + q 1 σ 1 2 τ 1 ) . Hence, the desired form of (17) becomes:
    T f ( s ) = a 1 + 1 3 ( 9 a 1 31 g 1 + 16 h 1 ) s + g 1 s 2 2 a 1 + 6 g 1 3 h 1 + 2 3 ( 9 a 1 23 g 1 + 11 h 1 ) s + h 1 s 2 , L f ( s ) = 1 3 ( 7 k 1 5 q 1 + 2 σ 1 + 4 τ 1 ) + q 1 s + k 1 s 2 1 3 ( 4 k 1 2 q 1 σ 1 + τ 1 ) + σ 1 s + τ 1 s 2 ,
    where a 1 , g 1 , h 1 , k 1 , q 1 , σ 1 , and τ 1 are free parameters. Similarly, we can find the second case:
  • Case B: { γ = 1 , T 0 = 1 , T 1 = 1 2 , L 0 = 1 , L 1 = 1 , T 2 = 3 2 , | T 2 | < , | L 3 | < }
    b 2 = 1 3 ( a 2 + 2 d 2 5 g 2 + 4 h 2 ) , c 2 = 1 3 ( 2 a 2 d 2 2 g 2 + h 2 ) , T 2 = 9 d 2 2 + 6 h 2 ( a 2 7 g 2 + 2 h 2 ) 6 d 2 ( 3 g 2 + 2 h 2 ) 8 ( a 2 + d 2 g 2 + 2 h 2 ) 2 , p 2 = 1 3 ( 7 k 2 + σ 2 + 2 τ 2 ) , r 2 = 1 3 ( 2 k 2 σ 2 τ 2 ) , q 2 = 1 3 ( 8 k 2 + σ 2 + 4 τ 2 ) , L 3 = 3 ( 3 σ 2 + 4 τ 2 ) 4 ( k 2 + σ 2 + τ 2 ) . Hence, the desired form of (17) becomes:
    T f ( s ) = a 2 + 1 3 ( a 2 + 2 d 2 5 g 2 + 4 h 2 ) s + g 2 s 2 1 3 2 a 2 d 2 2 g 2 + h 2 + d 2 s + h 2 s 2 , L f ( s ) = 1 3 ( 7 k 2 + σ 2 2 τ 2 ) + 1 3 ( 8 k 2 + σ 2 + 4 τ 2 ) s + k 2 s 2 1 3 ( 2 k 2 σ 2 τ 2 ) + σ 2 s + τ 2 s 2 ,
    where a 2 , d 2 , g 2 , h 2 , k 2 , σ 2 , and τ 2 are free parameters.
Although numerous forms of weight functions T f ( s ) and L f ( s ) satisfying (18) or (19) are applicable, we are specifically interested to the following forms:
  • Case A: γ = 2 3
    T f ( s ) { 1 2 ( 3 s + 1 3 s 1 ) , 1 8 ( 5 + 3 s 2 ) , 1 32 6 s + 1 + 125 6 s 1 , 23 8 3 s + 9 8 s 2 } , L f ( s ) { 1 4 ( 3 s + 1 3 s 1 ) 2 , 1 2 ( 3 s 1 ) , 2 3 s 1 , 2 s 5 s 3 , 1 4 ( 7 3 s 2 ) , 4 s 2 7 s 2 3 , 1 4 ( 1 + 3 s 2 ) , 4 1 + 3 s 2 , 5 3 s 2 , 2 ( 2 s ) ( 1 + s ) } .
Remark 2.
Existing studies with combinations of weight functions
( T f , L f ) { 1 2 ( 3 s + 1 3 s 1 ) , 1 4 ( 3 s + 1 3 s 1 ) 2 , 1 2 ( 3 s + 1 3 s 1 ) , 2 3 s 1 , 1 8 ( 5 + 3 s 2 ) , 1 2 ( 3 s 1 ) , 1 2 ( 3 s + 1 3 s 1 ) , 1 2 ( 3 s 1 ) } can be found, respectively, in [3,8,9,22]. These existing methods are denoted by E M 1 , E M 2 , E M 3 , and E M 4 , respectively, for use below. Note that selecting α = β = 2 in (16) readily yields T f ( s ) = 2 s 3 s 1 .
  • Case B: γ = 1
    T f ( s ) { 2 s 3 s 1 , 2 s 2 5 s 3 , 1 2 ( 3 s ) , 1 2 s ( 3 s 5 ) , 1 + s 2 s , 3 s 1 2 s 2 , 2 1 + s , 5 s 3 + s , 2 s ( s 3 ) , 1 4 ( 3 + 1 s 2 ) , 4 s 2 5 s 2 1 , 4 s 2 + 3 , 1 4 ( 5 s 2 ) , ( 3 s + 1 4 s ) 2 , ( s 5 ) 2 16 , 4 ( s 3 s 5 ) 2 , ( 4 s 5 s 1 ) 2 , 16 ( s + 3 ) 2 , 1 4 ( s + 3 s + 1 ) 2 } , L f ( s ) { 7 8 s + 3 s 2 2 , 1 2 ( 1 + 1 s 2 ) , 1 2 s 2 + 3 s , s + 1 3 s 1 , 1 32 4 s + 1 + 81 4 s 1 , 1 49 32 7 s + 96 7 s 3 } .
Remark 3.
Existing studies with combinations of weight functions
( T f , L f ) { 1 + s 2 s , 7 8 s + 3 s 2 2 , 2 1 + s , s + 1 3 s 1 , 1 + s 2 s , 1 2 ( 1 + 1 s 2 ) } can be found, respectively, in [7,22,23]. These existing methods are denoted by E M 5 , E M 6 and E M 7 , respectively, for use below.
In view of (20) and (21), we can select a total of 154 special pairs of second-order rational weight functions ( T f ( s ) , L f ( s ) ) . Excluding known studies, the following pairs of weight functions ( T f ( s ) , L f ( s ) ) may be of great interest to us. The corresponding methods to such pairs of weight functions ( T f ( s ) , L f ( s ) ) are denoted by L K i for 1 i 10 , respectively, and indicated on the right of (22) and (23).
  • Case A: γ = 2 3
    ( T f ( s ) , L f ( s ) ) = ( 1 2 ( 3 s + 1 3 s 1 ) , 2 s 5 s 3 ) , L K 1 ( 1 2 ( 3 s + 1 3 s 1 ) , 5 3 s 2 ) , L K 2 ( 1 8 ( 5 + 3 s 2 ) , 2 3 s 1 ) , L K 3 ( 1 8 ( 5 + 3 s 2 ) , 5 3 s 2 ) , L K 4 ( 23 8 3 s + 9 8 s 2 , 5 3 s 2 ) , L K 5 .
    One should be aware that Method L K 2 can be found by taking α = β = 2 from (16).
  • Case B: γ = 1
    ( T f ( s ) , L f ( s ) ) = ( 2 s 3 s 1 , s + 1 3 s 1 ) , L K 6 ( 3 s 2 , s + 1 3 s 1 ) , L K 7 ( 1 + s 2 s , s + 1 3 s 1 ) , L K 8 ( 2 1 + s , 1 2 ( 1 + 1 s 2 ) ) , L K 9 ( 5 s 3 + s , s + 1 3 s 1 ) , L K 10 .

4. Computational Experiments on Local and Global Convergence

For computational experiments, we first deal with local convergence of methods (1) for a variety of test functions along with the existing methods EM1–EM7; then we discuss global convergence underlying extraneous fixed points via basins of attraction. Numerical experiments for 17 methods EM1– EM7 and LK1–LK10 were implemented with Mathematica with 300 and 140 digits of minimum number of precision for scalar and vector equations, respectively.
Definition 1.
(Computational Convergence Order (COC) and Approximated Computational Convergence Order (ACOC)) Assume that theoretical asymptotic error constant η = lim n | e n | | e n 1 | p and convergence order p 1 are known. Define p n = | log | e n / η | log | e n 1 | | as the computational convergence order. Note that lim n p n = p . Approximated computational convergence order p ¯ n is defined as p ¯ n = | log ( | x n x n 1 | / | x n 1 x n 2 | ) log ( | x n 1 x n 2 | / | x n 2 x n 3 | ) | requiring knowledge at four points x n , x n 1 , x n 2 , x n 3 .

4.1. Local Convergence

Table 1 lists test scalar functions to measure the convergence behavior of proposed scheme (1). Computed values of x n is listed with up to 15 significant digits for proper readability. The error bound ϵ = 1 2 × 10 180 is assigned for scalar equations.
According to Table 2, sixth-order convergence is clearly seen. The values of computational asymptotic error constant agree up to 10 significant digits with η . It appears that the computational convergence order well approaches 6. In Table 3, we compare numerical errors | x n α | of proposed Methods LK1LK10 with those of existing Methods EM1EM7. The least errors within the prescribed error bound are highlighted in bold face. According to the comparison, Methods LK1 and LK8 display slightly better convergence for most test functions, while other remaining methods exhibit similar convergence.
A close inspection of the asymptotic error constant η ( α , θ i , T f , L f ) | x n + 1 α | | x n α | 6 addresses one’s attention to the local convergence dependent on f ( x ) , x 0 , α , T f and L f . Accordingly, the convergence of one method is hardly expected to be always better than the others.

4.2. Global Convergence

We usually locate a zero α of f ( x ) by means of a fixed point ξ of iterative maps:
x n + 1 = W f ( x n ) , n = 0 , 1 , .
In general, W f might possess other extraneous fixed points ξ α . It is well known that extraneous fixed points may result in attractive, indifferent or repulsive cycles as well as other periodic orbits influencing global convergence. Combining proposed methods (1) with maps (24), we find:
x n + 1 = W f ( x n ) = x n f ( x n ) f ( x n ) H f ( x n ) ,
where H f ( x n ) = T f ( s ) + L f ( s ) · f ( z n ) f ( x n ) can be regarded as a weight function of the classical Newton’s method. We are interested in the dynamics [24,25,26,27,28,29] of maps (1) underlying their extraneous fixed points of associated with their basins of attraction.
Good initial guesses for the numerical solutions of methods (1) can be determined from the basins of attraction which exhibits convergence of global character. Table 4 features statistical data including the average number of iterations per point, CPU time (in seconds), and number of points requiring 40 iterations. In the following examples, we take a 6 × 6 square centered at the origin containing all the zeros of the given test functions. We then take 601 × 601 equally spaced points in the square as initial points for methods (1). We color the point based on the root it converged to. This way we can figure out if the method converged within the maximum number of iteration and if it converged to the root closer to the initial point.
Figure 1, Figure 2, Figure 3 and Figure 4 present basins of attraction of 17 iterative maps when applied to various polynomials.
Example 1.
As a first example, we have taken a quadratic polynomial with all real roots:
p 1 ( z ) = ( z 2 1 ) .
Clearly the roots are ± 1 . Basins of attraction for all the listed methods are given in Figure 1. Consulting Table 4, we find that Methods LK1 and LK8 use the fewest iterations per point on average (AvgCon), and they also have the fewest black points. Other remaining methods have AvgCon ranging from 3.43 to 4.83. The fastest methods are EM2 with 49.406 s and EM6 with 48.032 s. Observe that Methods LK5 and LK6 exhibit more chaotic nature along the imaginary axis than others.
Example 2.
In our second example, we have taken a cubic polynomial:
p 2 ( z ) = ( z 3 + 4 z 2 10 ) .
Basins of attraction are given in Figure 2. We now consult Table 4 to find that the methods with the fewest AvgCon are EM1 with 4.2698 and EM6 with 3.9447 iterations. All the others require between 3.97 and 8.44. In terms of CPU time in seconds, the fastest are EM2 (129.172 s) and EM6 (128.797 s) and the slowest are LK5 (306.672 s) and EM5 (172.548 s). The methods having the most black points are LK5 with 75,147 and LK7 with 16,293, while the methods having the fewest are LK1 with 418 points and EM6 with no points. Method LK5 displays most chaotic nature near the basin boundaries axis, followed by Method LK4. Methods LK1, LK8, LK9 and EM6 present better stability than others.
Example 3.
As a third example, we have taken another cubic polynomial:
p 3 ( z ) = ( z 3 z ) .
Now, all the roots are real. The basins for this example are plotted in Figure 3. Based on Table 4 we find that the methods displaying the lowest AvgCon are EM3 with 4.3713 and LK8 with 4.2768. The fastest methods are EM2 with 131.750 s and EM6 with 129.766 s, while the slowest are LK1 with 229.297 s and LK3 with 164.906 s. The methods having the fewest black points are LK1 with 3000 and EM2 with 192. Methods LK4 and EM3 reveal the most chaotic nature, followed by Methods LK3 and LK5. Methods LK8 and EM3 are more stable than others.
Example 4.
As a fourth example, we have taken a quartic polynomial:
p 4 ( z ) = ( z 4 1 ) .
The basins are given in Figure 4. In terms of AvgCon, EM1 with 5.1189 and EM4 with 4.4454 are the best, while LK5 with 6.6932 and EM1 with 5.4899 are the worst. The fastest are EM2 with 110.578 s and EM6 with 100.078 s, while the slowest are LK4 with 541.970 s and EM1 with 285.126 s. The methods having the fewest black points are LK5 with 4 and EM2 with 0, while the methods having the most black points are LK4 with 8444 and LK1 with 5762. Methods EM1, EM3, EM5, LK3, LK4 and LK7 are more chaotic than the others, while Methods EM6 and LK6 are more stable. Method LK1 is of somewhat different stability character.

5. Extension to a Family of the Sixth-Order Methods for Nonlinear Systems of Equations

Let f : D C d C d with d N have a zero α D and be holomorphic in a neighborhood of α . Taylor expansion of f ( x n ) about α easily gives:
f ( x ( n ) ) = f ( α ) ( e ( n ) + c 2 e ( n ) 2 + + c m e ( n ) m ) + O ( | | e ( n ) m + 1 | | ) , n = 0 , 1 , ,
where e ( n ) = x ( n ) α and c j = 1 j ! f ( α ) 1 f ( j ) ( α ) for j 2 . For notational convenience, we drop the subscript n of e ( n ) and x ( n ) for the time being. We observe that f ( α ) and f ( α ) 1 are d × d matrices, with c j e ( n ) j C d . From (30), we find that the truncated f ( x ) defines a polynomial in e with matrix coefficients (independent of x). Hence, it is easily seen that
f ( x ) = f ( α ) ( I + 2 c 2 e + 3 c 3 e 2 + + m c m e m 1 ) + O ( | | e | | m ) , n = 0 , 1 , ,
where I is the d × d identity matrix. The inverse of f ( x ( n ) ) can be found by identifying B = ( 2 c 2 e + 3 c 3 e 2 + + m c m e m 1 from the relation
( I B ) 1 = I + B + B 2 + B 3 + , with   | | B | | < 1 .
Consequently, we find that:
f ( x ) 1 = ( I + X 1 e + X 2 e 2 + X 3 e 3 + X 4 e 4 + X 5 e 5 + X 6 e 6 ) f ( α ) 1 + O ( | | e 7 | | ) ,
where X 1 = 2 c 2 , X 2 = ( 4 c 2 2 3 c 3 ) , X 3 = 8 c 2 3 4 c 4 + 6 ( c 2 c 3 + c 3 c 2 ) , X 4 = 16 c 2 4 + 9 c 3 2 5 c 5 + 8 c 2 c 4 + 8 c 4 c 2 12 ( c 2 2 c 3 + c 3 c 2 2 + c 2 c 3 c 2 ) , X 5 = 32 c 2 5 6 c 6 + 2 ( 5 c 2 c 5 + 6 c 3 c 4 + 6 c 4 c 3 + 5 c 5 c 2 ) + 24 ( c 2 3 c 3 + c 3 c 2 3 + c 2 c 3 c 2 2 + c 2 2 c 3 c 2 ) 2 ( 9 c 2 c 3 2 + 8 c 2 2 c 4 + 9 c 3 2 c 2 + 8 c 4 c 2 2 + 8 c 2 c 4 c 2 + 9 c 3 c 2 c 3 ) , X 6 = 64 c 2 6 + 16 c 4 2 + 12 c 2 c 6 20 c 2 2 c 5 + 15 c 3 c 5 + 15 c 5 c 3 + 12 c 6 c 2 48 ( c 2 4 c 3 + c 3 c 2 4 + c 2 c 3 c 2 3 + c 2 2 c 3 c 2 2 + c 2 3 c 3 c 2 ) 24 ( c 2 c 3 c 4 + c 3 c 2 c 4 ) 3 ( 9 c 3 3 + 8 c 2 c 4 c 3 + 8 c 4 c 2 c 3 ) 4 ( 5 c 5 c 2 2 + 5 c 2 c 5 c 2 + 6 ( c 3 c 4 c 2 + c 4 c 3 c 2 ) ) + 4 ( 9 c 2 2 c 3 2 + 8 c 2 3 c 4 + 9 c 3 2 c 2 2 + 8 c 4 c 2 3 + 9 c 2 c 3 2 c 2 + 8 c 2 c 4 c 2 2 + 8 c 2 2 c 4 c 2 + 9 c 3 c 2 2 c 3 + 9 c 2 c 3 c 2 c 3 + 9 c 3 c 2 c 3 c 2 ) .
Additional computations show that:
f ( x ) 1 f ( x ) = e c 2 e 2 + ( 2 c 2 2 2 c 3 ) e 3 + ( 4 c 2 3 3 c 4 + 4 c 2 c 3 + 3 c 3 c 2 ) e 4 + j = 5 6 A j e j + O ( | | e | | 7 ) with A j = A j ( c 2 , c 3 , c 4 , c 5 , c 6 ) .
e y = y α = x α γ f ( x ) 1 f ( x ) = e ( 1 γ ) + γ c 2 e 2 2 γ ( c 2 2 c 3 ) e 3 + γ ( 4 c 2 3 + 3 c 4 4 c 2 c 3 3 c 3 c 2 ) e 4 + j = 5 6 B j e j + O ( | | e | | 7 ) with B j = B j ( γ , c 2 , c 3 , c 4 , c 5 , c 6 ) .
We find f ( y ( n ) ) = f ( x ) | e > e y and s = f ( x ) 1 f ( y ) = I 2 γ c 2 e + 3 γ ( 2 c 2 2 + c 3 ( γ 2 ) ) e 2 + j = 3 5 C j e j + O ( | | e | | 6 ) with C j = C j ( γ , c 2 , c 3 , c 4 , c 5 , c 6 ) .
( s I ) 2 = γ 2 4 c 2 2 e 2 + e 3 ( 24 c 2 3 6 ( 2 + γ ) c 2 c 3 6 ( 2 + γ ) c 3 c 2 ) ) e 3 + j = 4 5 D j e j + O ( | | e | | 6 ) with D j = D j ( γ , c 2 , c 3 , c 4 , c 5 , c 6 ) .
( s I ) 3 = 2 γ 2 [ 4 c 2 3 e 3 γ + ( 12 ( 1 + 2 γ ) c 2 4 + 6 ( 2 + γ ) γ c 2 2 c 3 + 6 ( 2 + γ ) γ c 3 c 2 2 + 6 ( 2 + γ ) γ c 2 c 3 c 2 ) e 4 + E 5 e 5 ] ) + O ( | | e | | 6 ) with E 5 = E 5 ( γ , c 2 , c 3 , c 4 , c 5 , c 6 ) .
Theorem 1 suggests us to use T f and L f as at most third-degree matrix polynomials in ( s I ) .
e z = z α = x α T f ( s ) f ( x ) 1 f ( x ) = e ( T 0 I + T 1 ( s I ) + T 2 ( s I ) 2 + T 3 ( s I ) 3 ) f ( x ) 1 f ( x ) = e ( 1 T 0 ) + c 2 e 2 ( T 0 + 2 T 1 γ ) + ( c 3 ( 2 T 0 3 T 1 ( 2 + γ ) γ ) 2 c 2 2 ( T 0 + 2 γ ( 2 T 1 + T 2 γ ) ) ) e 3 + j = 4 6 F j e j ] + O ( | | e | | 7 ) with E j = F j ( γ , c 2 , c 3 , c 4 , c 5 , c 6 , T 0 , T 1 , T 2 ) .
Equating the coefficients of the first and second-order terms in e z yields:
T 0 = 1 , T 1 = 1 / ( 2 γ ) .
We obtain: f ( z ) = f ( x ) | e > e z and f ( x ) 1 f ( z ) = ( 1 T 0 ) e + ( 2 + 3 T 0 + 2 T 1 γ ) c 2 e 2 + ( ( 2 + 2 T 0 + 3 T 0 2 T 0 3 3 T 1 ( 2 + γ ) γ ) c 3 4 ( 1 + 2 T 0 + 3 T 1 γ + T 2 γ 2 ) ) c 2 2 e 3 + j = 4 6 G j e j + O ( | | e | | 7 ) , with G j = G j ( γ , c 2 , c 3 , c 4 , c 5 , c 6 , T 0 , T 1 , T 2 ) .
e ( n + 1 ) = x ( n + 1 ) α = z ( n ) α L f ( s ) f ( x ( n ) ) 1 f ( z ( n ) ) = e z ( L 0 I + L 1 ( s I ) + L 2 ( s I ) 2 + L 3 ( s I ) 3 ) f ( x ( n ) ) 1 f ( z ( n ) ) = ( L 0 1 ) ( T 0 1 ) e + ( T 0 2 L 1 T 0 γ + 2 ( L 1 + T 1 ) γ + L 0 ( 2 3 T 0 2 T 1 γ ) ) c 2 e 2 + j = 3 6 H j e j ] + O ( | | e 7 | | ) with H j = H j ( γ , c 2 , c 3 , c 4 , c 5 , c 6 , T 0 , T 1 , T 2 , L 0 , L 1 , L 2 ) .
Now, we annihilate the first five coefficients of the terms up to the fifth-order terms of e ( n + 1 ) with the use of (34) by taking the set of coefficients below:
{ γ = 2 / 3 , T 0 = 1 , T 1 = 3 / 4 , T 2 = 9 / 8 , L 0 = 1 , L 1 = 3 / 2 } ,
{ γ = 1 , T 0 = 1 , T 1 = 1 / 2 , L 0 = 1 , L 1 = 1 , L 2 = 3 / 2 } .
The discussions thus far lead us to the following theorem for nonlinear systems of equations.
Theorem 2.
Let f : Ω C d C d with d N have a simple root α and be sufficiently Frétchet differentiable in Ω containing α. Let x ( 0 ) be an initial guess chosen close to α. Let T f , L f : C d × d C d × d be matrix functions sufficiently Frétchet differentiable in a neighborhood of I, being defined by:
T f ( s ) = T 0 I + T 1 ( s I ) + T 2 ( s I ) 2 + T 3 ( s I ) 3 + O ( | | ( s I ) 4 | | ) ,
L f ( s ) = L 0 I + L 1 ( s I ) + L 2 ( s I ) 2 + L 3 ( s I ) 3 + O ( | | ( s I ) 4 | | ) with T i , L i C for i = 0 , 1 , 2 , 3 and s = f ( x n ) 1 f ( x n ) . If { γ = 2 3 , T 0 = 1 , T 1 = 3 4 , T 2 = 9 8 , L 0 = 1 , L 1 = 3 2 , | T 3 | < , | L 2 | < } or { γ = 1 , T 0 = 1 , T 1 = 1 2 , L 0 = 1 , L 1 = 1 , L 2 = 3 2 , | L 3 | < , | T 2 | < } are given, then iterative scheme (1) reduces to a family of sixth-order methods satisfying the error equation below. For n = 0 , 1 , 2 , ,
e ( n + 1 ) = [ 2 ( 27 8 L 2 ) ( 135 + 64 T 3 ) 243 c 2 5 + 2 ( 27 8 L 2 ) 81 c 2 2 c 4 ( 135 + 64 T 3 ) 27 c 3 c 2 3 c 3 c 4 9 + c 3 2 c 2 2 ( 27 8 L 2 ) 9 c 2 2 c 3 c 2 ] e ( n ) 6 + O ( | | e ( n ) 7 | | ) , i f γ = 2 3 . [ 4 c 2 5 ( 1 + 2 T 2 ) ( 9 + 4 L 3 + 2 T 2 ) + ( 6 + 12 T 2 ) c 3 c 2 3 + ( 1 2 T 2 ) c 2 c 3 c 2 2 + 1 4 ( c 2 c 3 2 + 4 ( 9 + 4 L 3 + 2 T 2 ) c 2 3 c 3 6 c 3 c 2 c 3 ) ] e ( n ) 6 + O ( | | e ( n ) 7 | | ) , i f γ = 1 .
where c j = 1 j ! f ( α ) 1 f ( j ) ( α ) for j = 2 , 3 , .
Equation (37) clearly reduces to (2) for a scalar function by identifying c i with θ i . In what follows, we employ several test examples for the zeros of vector-valued functions to verify the convergence behavior claimed here. In terms of Euclidean norm | | | | , we display the error sizes for e k = | | x ( k + 1 ) x ( k ) | | , residual error | | f ( x ( k + 1 ) ) | | and ACOC using the error criterion of | | x ( k + 1 ) x ( k ) | | < 10 140 within 20 iterations.
Test Example 1
We consider a nonlinear algebraic vector equation f : R 3 R 3 defined by f ( x ) = 0 with x = ( x 1 , x 2 , x 3 ) T as follows:
π ( x 1 2 + x 2 2 / 2 ) 3 x 3 = 0 , x 1 2 + x 2 / 2 + 2 cos x 3 = 0 , x 1 x 2 cos x 2 · sin ( 2 x 3 ) 2 = 0 .
The exact solution is given by x = ( 1 , 2 , π ) T . We try to solve (38) with an initial guess vector x ( 0 ) = ( 0.8 , 1.8 , 3.0 ) T by method (1), and find the results in Table 5. We observe that ACOC approaches up to 6, which is the theoretical order of convergence.
Test Example 2
We consider a nonlinear ODE boundary-value problem given below:
2 y ( x ) y ( x ) + y ( x ) 2 + 4 y ( x ) 2 = 0 , y ( π 6 ) = 1 / 4 , y ( π 2 ) = 1 .
The exact solution is found to be y = ( sin x ) 2 . With the use of the central finite difference method, the first and second derivatives are approximated by:
y ( x ) y n + 1 y n 1 2 h , y ( x ) y n + 1 2 y n + y n 1 h 2 ,
where y n = y ( x n ) , h = 1 N ( π 2 π 6 ) = π 3 N , N is the number of divisions of the interval [ π 6 , π 2 ] . It can be shown that y ( x + h ) = y ( x ) + O ( h 3 ) , y ( x ) = O ( h 2 ) and y ( x ) = O ( h 2 ) in view of Taylor expansion of y ( x + h ) about x. This discretization yields the algebraic equations with 6 unknowns y 0 , y 1 , y 2 , y 3 , y 4 , y 5 :
y j 1 2 16 ( h 2 1 ) y j 2 + y j 1 ( 8 y j 2 y j + 1 ) 8 y j y j + 1 + y j + 1 2 = 0 ,
for j = 0 , 1 , , N 1 , with boundary conditions y 0 = y ( x 0 ) = y ( π 6 ) = 1 4 , and y N = y ( x N ) = y ( π 2 ) = 1 . Further computation after selecting N = 5 gives us a nonlinear algebraic vector equation f : R 4 R 4 defined by f ( y ) = 0 with y = ( y 1 , y 2 , y 3 , y 4 ) T of the form:
1 / 16 16 ( 1 + h 2 ) y 1 2 + 1 / 4 ( 8 y 1 2 y 2 ) 8 y 1 y 2 + y 2 2 = 0 , y 1 2 16 ( 1 + h 2 ) y 2 2 + y 1 ( 8 y 2 2 y 3 ) 8 y 2 y 3 + y 3 2 = 0 , y 2 2 16 ( 1 + h 2 ) y 3 2 + y 2 ( 8 y 3 2 y 4 ) 8 y 3 y 4 + y 4 2 = 0 , 1 + y 3 2 2 y 3 ( 1 + 4 y 4 ) 8 y 4 16 ( 1 + h 2 ) y 4 2 = 0 .
After solving (42) with an initial guess vector y ( 0 ) = ( 0.6 , 0.7 , 0.8 , 0.9 ) T by a typical method LK1, we find the results in Table 6 and Figure 5. It is seen that ACOC approaches up to 6, which is the theoretical order of convergence.
The errors | ( sin x ) 2 y i | for 1 i 4 ) at the internal nodes are, respectively, given by:
0.0045808024173838738376216319522 , 0.00731167760933327109201071704851 , 0.0073683157995758596940525729079 , 0.00474700015993753712029421111822 .
As a remark, we should note that the numerical solution by the central finite-difference methods is accurate within the range of Δ y ( x ) = O ( h 3 ) with h = 5 π 3 N = 0.00918704 .
Test Example 3
A two-dimensional nonlinear reaction-diffusion equation for a concentration u ( x , t ) of the substance under consideration in a bounded domain Ω R 2 with continuous boundary Ω is represented by an initial boundary value problem:
u t d Δ u = u ( a u ) in Ω × ( 0 , ) , u = g on Ω × ( 0 , ) ,
where d > 0 is a diffusion coefficient, a is a positive constant, g is continuous on Ω , and Δ is the Laplacian operator. For brevity of analysis, let d = 1 , a = 1 , and Ω = [ 0 , 1 ] × [ 0 , 1 ] (i.e., unit square region). We are interested in steady state solutions to (43), which lead us to elliptic partial differential equations with Dirichlet boundary conditions as follows:
u x x + u y y = u ( u 1 ) in [ 0 , 1 ] × [ 0 , 1 ] , u ( x , 0 ) = u ( x , 1 ) = x ( x 1 ) 2 + 1 , u ( 0 , y ) = u ( 1 , y ) = y ( y 1 ) 2 + 1 .
By using central divided differences with step h = 1 / 4 in each component of the space vector, we discretize (44) into a nonlinear system of equations with 25 nodes, 9 of which constitute interior nodal variables x 1 , x 2 , , x 9 in Ω , while the remaining 16 nodes are boundary nodes. As a result, we obtain a nonlinear algebraic vector equation f : R 9 R 9 defined by:
f ( x ) = A x + h 2 Ψ ( x ) b = 0 , with x = ( x 1 , x 2 , , x 9 ) T ,
where A = B I 0 I B I 0 I B , B = 4 h 2 1 0 1 4 h 2 1 0 1 4 h 2 , I is the identity matrix of size 3 × 3 , Ψ ( x ) = ( x 1 2 , x 2 2 , , x 9 2 ) T and b = ( 29 16 , 7 8 , 29 16 , 7 8 , 0 , 7 8 , 29 16 , 7 8 , 29 16 ) T .
We solve (38) with an initial guess vector x ( 0 ) = ( 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 ) T by a typical method LK1, and find the results in Table 7. It is evident that ACOC reaches up to 6, being the theoretical order of convergence. As can be seen in Table 7, the methods with γ = 2 / 3 appear to converge more quickly and better than those with γ = 1 .
Interior 16 nodal values of the steady-state solution of u ( x , t ) are illustrated with adjacent nodal points connected by straight lines in Figure 6.
Test Example 4
A d-dimensional nonlinear equation f : R d R d with d = m a x ( 4 , n N ) is given by:
0 = x i cos ( x i j = 1 4 x j ) , i = 1 , 2 , , d N .
The above nonlinear system is described in [4]. Selecting n = 10 , we find d = 10 and solve (46) in R 10 with an initial guess vector x ( 0 ) = ( 0.75 , 0.75 , , 0.75 ) T for the desired root x = ( 0.5149332 , 0.5149332 , , 0.5149332 ) T R 10 in Table 8. It is evident that ACOC reaches up to 6, being the theoretical order of convergence. As can be seen in Table 8, the methods with γ = 2 / 3 appear to converge more quickly and better than those with γ = 1 .

Computational Efficiency

The computational efficiency of an iterative method is defined by an efficiency index E = ρ 1 / d [30], with ρ as the order of convergence and d as the number of functional evaluations per iteration. We require n scalar functions for each f and n 2 for each f . The concept of the efficiency index E applied to a nonlinear system of vector equations has been extended to treat the concept of computational efficiency by using C E = ρ 1 / ( d + o p ) [4], where o p is the number of operations associated with products and quotients. Suppose that n is the size of the matrix needed in the nonlinear system of vector equations. Matrix inversion requires n 3 n 3 product-quotient operations and L U -decomposition technique for solving linear systems requires n 2 product-quotient operations, including the n 2 product-quotient operations related to matrix multiplication by a vector. Note that each method treated here follows the three set of linear systems (1) and has one matrix inverse f ( x n ) 1 . Consequently, the number of functional evaluations plus product-quotient operations d + o p becomes 2 n + 2 n 2 + n 3 n 3 + 6 n 2 = n 3 + 18 n 2 + 5 n 3 , which gives us the computational efficiency C E = 6 3 n 3 + 18 n 2 + 5 n for each listed method.
Many real-life application problems include ones related to: interval arithmetic benchmark, neurophysiology, chemical equilibrium, kinematic application, combustion application, and economics modeling, whose studies are described in [31]. The methods used therein are based on second-order Newton-like approach which may be more efficient in real-life problems in terms of speed and computational cost. On the other hand, our proposed family of sixth-order methods (1) is much more accurate than Newton-like methods, but has more complexities owing to the high-order formulation and require more CPU time to get the desired solution.
One certainly has to acknowledge that determining a better method than the other one should be avoided through solving a function with a randomly chosen initial guess vector and comparing the number of convergent iterations.

6. Conclusions

A family of Jarratt-like iterative methods for scalar and vector equations is developed and its convergence properties are theoretically established through Theorems 1 and 2. Computational aspects applied to various test equations agree well with the convergence behavior claimed in the theory developed. Global convergence behavior of the listed methods is illustrated for typical polynomials based upon their basins of attraction. The basins of attraction suggest selecting members of the iterative methods (1) give better convergence.
We will focus our future study on extending the current approach with different weight functions to the development of higher-order iterative root-finders.

Author Contributions

Formal analysis, Y.I.K.; Investigation, M.-Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Amat, S.; Busquier, S.; Plaza, S. Dynamics of the King and Jarratt iterations. Aeq. Math. 2005, 69, 212–223. [Google Scholar] [CrossRef]
  2. Bahl, A.; Cordero, A.; Sharma, R.; Torregrosa, J.R. A novel bi-parametric sixth order iterative scheme for solving nonlinear systems and its dynamics. Appl. Math. Comput. 2019, 357, 147–166. [Google Scholar] [CrossRef]
  3. Cordero, A.; Hueso, J.L.; Martinez, E.; Torregrosa, J.R. A modified Newton-Jarratt’s composition. Numer. Algor. 2010, 55, 87–99. [Google Scholar] [CrossRef]
  4. Capdevila, R.R.; Cordero, A.; Torregrosa, J.R. A New Three-Step Class of Iterative Methods for Solving Nonlinear Systems. Mathematics 2019, 7, 1221. [Google Scholar] [CrossRef] [Green Version]
  5. Cordero, A.; Torregrosa, J.R.; Vassileva, M.P. Three-step iterative methods with optimal eighth-order convergence. J. Comput. Appl. Math. 2011, 235, 3189–3194. [Google Scholar] [CrossRef] [Green Version]
  6. Narang, M.; Bhatia, S.; Kanwar, V. New two-parameter Chebyshev-Halley-like family of fourth and sixth-order methods for systems of nonlinear equations. Appl. Math. Comput. 2014, 248, 215–224. [Google Scholar] [CrossRef]
  7. Parhi, S.K.; Gupta, D.K. A sixth order method for nonlinear equations. Appl. Math. Comput. 2008, 203, 50–55. [Google Scholar] [CrossRef]
  8. Soleymani, F.; Lotfi, T.; Bakhtiari, P. A multi-step class of iterative methods for nonlinear systems. Optim. Lett. 2014, 8, 1001–1015. [Google Scholar] [CrossRef]
  9. Wang, X.; Kou, J.; Gu, C. Semilocal convergence of a sixth-order Jarratt method in Banach spaces. Numer. Algor. 2011, 57, 441–456. [Google Scholar] [CrossRef]
  10. Bi, W.; Wu, Q.; Ren, H. A new family of eighth-order iterative methods for solving nonlinear equations. Appl. Math. Comput. 2009, 214, 236–245. [Google Scholar] [CrossRef]
  11. Džunić, J.; Petković, M.S.; Petković, L.D. A family of optimal three-point methods for solving nonlinear equations using two parametric functions. Appl. Math. Comput. 2011, 217, 7612–7619. [Google Scholar] [CrossRef]
  12. Li, Z.; Jiang, J.; Hong, L.; Sun, J.Q. A subspace expanding technique for global zero finding of multi-degree-of-freedom nonlinear systems. Appl. Math. Mech. 2020, 41, 769–784. [Google Scholar] [CrossRef]
  13. Liu, L.; Wang, X. Eighth-order methods with high efficiency index for solving nonlinear equations. Appl. Math. Comput. 2010, 215, 3449–3454. [Google Scholar] [CrossRef]
  14. Sharifi, S.; Ferrara, M.; Salimi, M.; Siegmund, S. New modification of Maheshwari’s method with optimal eighth order convergence for solving nonlinear equations. Open Math. 2016, 14, 443–451. [Google Scholar] [CrossRef] [Green Version]
  15. Xiong, F.R.; Schütze, O.; Ding, Q.; SunXiao, J.Q. Finding zeros of nonlinear functions using the hybrid parallel cell mapping method. Commun. Nonlinear Sci. Numer. Simul. 2016, 334, 23–37. [Google Scholar] [CrossRef]
  16. Ahlfors, L.V. Complex Analysis; McGraw-Hill Book, Inc.: New York, NY, USA, 1979. [Google Scholar]
  17. Hörmander, L. An Introduction to Complex Analysis in Several Variables; North-Holland Publishing Company: Amsterdam, The Netherlands, 1973. [Google Scholar]
  18. Shabat, B.V. Introduction to Complex Analysis PART II, Functions of Several Variables; American Mathematical Society: Providence, RI, USA, 1992. [Google Scholar]
  19. Vrscay, E.R.; Gilbert, W.J. Extraneous Fixed Points, Basin Boundaries and Chaotic Dynamics for shröder and König rational iteration Functions. Numer. Math. 1988, 52, 1–16. [Google Scholar] [CrossRef]
  20. Wolfram, S. The Mathematica Book, 5th ed.; Wolfram Media: Champaign, IL, USA, 2003. [Google Scholar]
  21. Magre nan, Á.A. A new tool to study real dynamics: The convergence plane. Appl. Math. Comput. 2016, 275, 394–403. [Google Scholar]
  22. Xiao, X.-Y.; Yin, H.-W. Accelerating the convergence speed of iterative methods for solving nonlinear systems. Appl. Math. Comput. 2018, 333, 8–19. [Google Scholar] [CrossRef]
  23. Lotfi, T.; Bakhtiari, P.; Cordero, A.; Mahdiani, K.; Torregrosa, J.R. Some new efficient multipoint iterative methods for solving nonlinear systems of equations. Int. J. Comput. Math. 2015, 92, 1921–1934. [Google Scholar] [CrossRef] [Green Version]
  24. Andreu, C.; Cambil, N.; Cordero, A.; Torregrosa, J.R. A class of optimal eighth-order derivative-free methods for solving the Danchick-Gauss problem. Appl. Math. Comput. 2014, 232, 237–246. [Google Scholar] [CrossRef]
  25. Argyros, I.K.; Magre nán, A.Á. On the convergence of an optimal fourth-order family of methods and its dynamics. Appl. Math. Comput. 2015, 252, 336–346. [Google Scholar] [CrossRef]
  26. Chun, C.; Neta, B. Comparison of several families of optimal eighth order methods. Appl. Math. Comput. 2016, 274, 762–773. [Google Scholar] [CrossRef] [Green Version]
  27. Cordero, A.; García-Maimó, J.; Torregrosa, J.R.; Vassileva, M.P.; Vindel, P. Chaos in King’s iterative family. Appl. Math. Lett. 2013, 26, 842–848. [Google Scholar] [CrossRef] [Green Version]
  28. Geum, Y.H.; Kim, Y.I.; Magre nán, Á.A. A biparametric extension of King’s fourth-order methods and their dynamics. Appl. Math. Comput. 2016, 282, 254–275. [Google Scholar] [CrossRef]
  29. 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]
  30. Traub, J.F. Iterative Methods for the Solution of Equations; Chelsea Publishing Company: Chelsea, VT, USA, 1982. [Google Scholar]
  31. Grosan, C.; Abraham, A. A New Approach for Solving Nonlinear Equations Systems. IEEE Trans. Syst. Man Cybern. Part A Syst. Hum. 2008, 38, 698–714. [Google Scholar] [CrossRef]
Figure 1. Basins of attraction of the listed methods, for the roots of the polynomial z 2 1 .
Figure 1. Basins of attraction of the listed methods, for the roots of the polynomial z 2 1 .
Algorithms 13 00303 g001
Figure 2. Basins of attraction of the listed methods, for the roots of the polynomial z 3 + z 2 + z + 1 .
Figure 2. Basins of attraction of the listed methods, for the roots of the polynomial z 3 + z 2 + z + 1 .
Algorithms 13 00303 g002
Figure 3. Basins of attraction of the listed methods, for the roots of the polynomial z ( z 3 + 1 ) .
Figure 3. Basins of attraction of the listed methods, for the roots of the polynomial z ( z 3 + 1 ) .
Algorithms 13 00303 g003
Figure 4. Basins of attraction of the listed methods, for the roots of the polynomial z ( z 2 + i 8 ) .
Figure 4. Basins of attraction of the listed methods, for the roots of the polynomial z ( z 2 + i 8 ) .
Algorithms 13 00303 g004
Figure 5. An ODE boundary value problem for Example 2.
Figure 5. An ODE boundary value problem for Example 2.
Algorithms 13 00303 g005
Figure 6. Steady state solution of the reaction-diffusion equation for Example 3.
Figure 6. Steady state solution of the reaction-diffusion equation for Example 3.
Algorithms 13 00303 g006
Table 1. Test functions f i ( x ) with zeros α and initial guesses x 0 .
Table 1. Test functions f i ( x ) with zeros α and initial guesses x 0 .
i f i ( x ) α x 0
1 sin x log ( 1 + x 2 ) 0 0.01
2 3 + sin x x 2 1.97932014655621 2.0
3 2 x π + cos x · log ( x 2 + 1 ) π 2 1.53
4 2 x 3 + e x 2 + sin x 2 0.719549366870672 0.73
5 x 3 x 3 cos ( π x 6 ) + 1 ( x 2 + 1 ) 11 5 + 4 3 2 1.87
6 e x 3 + 1 x 5 + 7 cos ( x 3 + 1 ) 1 1 + i 3 2 0.52 + 0.85 i
7 x log x x + x 2 1 1.05
Here log z(z C ) represents a principal analytic branch with − π ≤ Im(logz) < π.
Table 2. Convergence for test functions f 1 ( x ) f 4 ( x ) with typically selected Methods EM1, LK1, EM5, LK6.
Table 2. Convergence for test functions f 1 ( x ) f 4 ( x ) with typically selected Methods EM1, LK1, EM5, LK6.
Method fn x n | f ( x n ) | | x n α | | e n / e n 1 6 | η p n
00.010.009899840.0100000
1 1.33986049407934 × 10 12 1.339 × 10 12 1.339 × 10 12 1.3398604941.2962962965.99282
EM1 f 1 2−7.50000879616187 × 10 72 7.500 × 10 72 7.500 × 10 72 1.296296296 6.00000
3−2.30714514140106 × 10 427 2.307 × 10 427 2.307 × 10 427
02.00.09070260.0206799
11.97932014655603 7.783 × 10 13 1.786 × 10 13 0.0022845037840.002483362146.02152
LK1 f 2 21.97932014655621 3.520 × 10 79 8.081 × 10 80 0.002483362140 6.00000
31.97932014655621 0.0 × 10 299 3.015 × 10 328
01.530.03239610.0407963
11.57079629958335 2.058 × 10 8 2.721 × 10 8 5.9023757917.1905181066.06171
EM5 f 3 21.57079632679490 2.208 × 10 45 2.919 × 10 45 7.190518106 6.00000
31.57079632679490 3.367 × 10 267 4.450 × 10 267
00.730.03180410.0104506
10.719549366862969 2.311 × 10 11 7.703 × 10 12 5.9130124096.1206425656.00757
LK6 f 4 20.719549366870672 3.837 × 10 66 1.278 × 10 66 6.120642565 6.00000
30.719549366870672 1.20621 × 10 327 1.507 × 10 327
Table 3. Local convergence of selected methods for various test functions.
Table 3. Local convergence of selected methods for various test functions.
γ Method | x n α | f ( x ) ; x 0
f 1 ; 0.01 f 2 ; 2.0 f 3 ; 1.53 f 4 ; 0.73 f 5 ; 1.87 f 6 ; 0.52 + 0.85 i f 7 ; 1.05
EM1 | x 1 α | 1.33e-12 * 4.03e-135.07e-91.64e-123.13e-59.14e-102.26e-9
| x 2 α | 7.50e-722.30e-771.99e-502.49e-712.59e-262.72e-542.34e-53
EM2 | x 1 α | 2.54e-127.48e-131.11e-84.50e-123.92e-58.81e-103.89e-9
| x 2 α | 6.61e-701.75e-755.43e-482.97e-681.63e-252.16e-541.11e-51
EM3 | x 1 α | 5.88e-121.68e-123.05e-81.49e-115.62e-58.08e-108.09e-9
| x 2 α | 2.26e-675.13e-736.77e-451.34e-642.73e-241.28e-541.94e-49
EM4 | x 1 α | 4.17e-121.20e-121.89e-88.28e-124.89e-58.49e-106.03e-9
| x 2 α | 2.05e-684.97e-742.37e-462.14e-669.32e-251.76e-542.45e-50
2 / 3 LK1 | x 1 α | 6.33e-131.78e-136.13e-93.26e-121.37e-59.89e-106.46e-10
| x 2 α | 3.58e-748.08e-808.66e-503.13e-699.64e-304.55e-544.72e-57
LK2 | x 1 α | 7.48e-122.10e-123.32e-81.561e-116.43e-58.34e-101.00e-8
| x 2 α | 1.20e-662.51e-721.29e-441.86e-648.09e-241.78e-549.18e-49
LK3 | x 1 α | 3.59e-121.04e-121.79e-88.13e-124.50e-58.39e-105.22e-9
| x 2 α | 7.27e-691.80e-741.55e-461.87e-664.76e-251.57e-548.82e-51
LK4 | x 1 α | 1.05e-112.93e-125.35e-82.82e-117.39e-57.94e-101.34e-8
| x 2 α | 1.32e-652.59e-713.71e-431.17e-622.37e-231.29e-547.28e-48
LK5 | x 1 α | 3.58e-119.46e-121.94e-71.24e-101.27e-46.41e-103.85e-8
| x 2 α | 6.72e-629.48e-683.57e-394.05e-581.74e-213.79e-551.24e-44
EM5 | x 1 α | 2.02e-123.88e-132.72e-82.23e-112.60e-52.17e-91.88e-9
| x 2 α | 1.16e-701.99e-772.91e-452.25e-632.11e-261.19e-511.16e-53
EM6 | x 1 α | 1.38e-123.93e-132.88e-98.25e-131.33e-52.32e-91.96e-9
| x 2 α | 9.18e-721.94e-773.98e-522.26e-731.08e-281.49e-519.28e-54
EM7 | x 1 α | 4.19e-138.51e-145.45e-93.56e-121.17e-52.03e-94.68e-10
| x 2 α | 2.00e-754.73e-823.20e-505.72e-694.62e-296.12e-526.03e-58
LK6 | x 1 α | 3.93e-121.12e-121.81e-87.70e-125.75e-51.91e-95.60e-9
| x 2 α | 1.36e-683.03e-741.65e-461.27e-662.27e-244.21e-521.41e-50
1LK7 | x 1 α | 7.75e-132.18e-131.10e-81.25e-112.17e-52.59e-91.60e-9
| x 2 α | 1.73e-733.02e-797.41e-484.21e-652.85e-273.33e-512.03e-54
LK8 | x 1 α | 2.27e-134.60e-142.11e-91.07e-128.59e-62.08e-92.65e-10
| x 2 α | 2.82e-776.39e-844.14e-531.29e-724.42e-307.24e-521.11e-59
LK9 | x 1 α | 3.38e-129.73e-132.33e-81.20e-111.99e-52.25e-94.43e-9
| x 2 α | 4.73e-691.11e-741.02e-452.98e-652.16e-271.22e-512.77e-51
LK10 | x 1 α | 1.36e-123.81e-132.49e-95.51e-121.87e-52.45e-92.08e-9
| x 2 α | 8.46e-721.55e-772.54e-521.31e-671.12e-272.21e-511.37e-53
* 1.33e-12 ≡ 1.33 × 10−12.
Table 4. Global convergence behavior of selected methods for various examples.
Table 4. Global convergence behavior of selected methods for various examples.
γ ExConvMethod
EM1EM2EM3EM4LK1LK2LK3LK4LK5
Ex1CPU55.98549.40659.70458.01651.75052.78159.29759.49958.985
NConpts000000000
Conpts360,000360,000359,964359,964360,000359,996360,000359,856351,604
AvgCon3.59563.64213.88833.88833.33673.95413.73594.11724.8282
Bkpts0036360401448396
Ex2CPU147.140129.172181.671169.453262.750148.141186.343192.281306.672
NConpts00003113000032030
2 / 3 Conpts348,170357,746343,698357,078328,452354,802348,498337,638252,823
AvgCon4.26984.49814.65984.66014.36914.81864.50584.91458.4395
Bkpts11,830225416,3022922418519811,50222,36275,147
Ex3CPU154.452131.750149.125161.625229.297149.266181.828147.829159.766
NConpts000060000
Conpts349,039354,026226,092349,416356,994344,930308,029226322287215
AvgCon4.42754.54294.37134.64757.73634.79424.48934.47664.9749
Bkpts10,9615974133,90810,584300015,07051,971133,67872,785
Ex4CPU211.297110.578391.172142.797504.171150.000303.328541.97146.984
NConpts0000182,324000118
Conpts357,420359,660354,742359,618177,562358,856356,208351,556359,878
AvgCon5.11895.16825.37575.27905.87425.37895.26015.55506.6932
Bkpts258034052583821141144379284444
EM5EM6EM7LK6LK7LK8LK9LK10-
Ex1CPU57.67148.03254.53254.29750.64154.10958.28156.328
NConpts0014800000
Conpts359,924360,000359,852354,908359,984360,000360,000359,856-
AvgCon3.99653.43663.59323.90173.50503.39353.74463.4631
Bkpts7600509216000
Ex2CPU172.548128.797159.328147.407144.687148.328172.500162.219
NConpts00000000
1 Conpts343,983360,000354,417356,718343,707359,474359522358554-
AvgCon4.88803.94474.44694.47694.28343.97764.572634.1854
Bkpts16,01705583328216,2935264781446
Ex3CPU152.344129.766151.719138.062138.062152.187164.906155.969
NConpts00006000
Conpts334,256359,808350,497337,990344,846357,584354,886353,352-
AvgCon4.81714.28234.51544.33294.38034.27684.70964.3668
Bkpts25744192950322,01015,154241651146648
Ex4CPU285.126100.078197.062102.141270.438128.39140.844131.531
NConpts00000000
Conpts355,406360,000357,822359,320354,238359,774359,820359,518-
AvgCon5.48994.78715.15234.44544.87864.83665.32924.8444
Bkpts4594021786805762226180482
Conv, convergent behavior; CPU, processing CPU time in seconds; NConpts, number of points whose each orbit is non-convergent but bounded; Conpts, number of points whose each orbit is convergent; AvgCon, average number of iterations for convergence per point; Bkpts, number of points whose each orbit tends to infinity within 40 iterations.
Table 5. Convergence results of Test Example 1.
Table 5. Convergence results of Test Example 1.
γ MTConvk
1234
| | x ( k + 1 ) x ( k ) | | 9.950514 × 10 5 1.605920 × 10 24 6.683987 × 10 144 -
EM1 | | f ( x ( k + 1 ) ) | | 2.411043 × 10 4 2.777255 × 10 24 1.003785 × 10 143 -
ACOC--6.031727-
| | x ( k + 1 ) x ( k ) | | 1.874364 × 10 4 1.263701 × 10 23 1.441060 × 10 138 -
EM2 | | f ( x ( k + 1 ) ) | | 1.728447 × 10 4 1.309191 × 10 23 1.564793 × 10 138 -
ACOC--5.995603-
| | x ( k + 1 ) x ( k ) | | 6.184026 × 10 4 1.132727 × 10 20 4.735571 × 10 120 -
EM3 | | f ( x ( k + 1 ) ) | | 5.710069 × 10 4 2.808419 × 10 20 6.200273 × 10 120 -
ACOC--5.937616-
| | x ( k + 1 ) x ( k ) | | 4.268556 × 10 4 1.600843 × 10 21 3.291193 × 10 125 -
EM4 | | f ( x ( k + 1 ) ) | | 3.530786 × 10 4 2.797766 × 10 21 5.120856 × 10 125 -
ACOC--5.950154-
| | x ( k + 1 ) x ( k ) | | 2.188288 × 10 4 2.036417 × 10 22 4.127188 × 10 132 -
2 / 3 LK1 | | f ( x ( k + 1 ) ) | | 3.544655 × 10 4 1.523420 × 10 22 1.380288 × 10 131 -
ACOC--6.083509-
| | x ( k + 1 ) x ( k ) | | 1.099663 × 10 3 1.535614 × 10 18 2.400517 × 10 108 -
LK2 | | f ( x ( k + 1 ) ) | | 1.373880 × 10 3 1.832089 × 10 18 4.662356 × 10 108 -
ACOC--6.045514-
| | x ( k + 1 ) x ( k ) | | 2.705259 × 10 4 1.464141 × 10 22 5.102270 × 10 132 -
LK3 | | f ( x ( k + 1 ) ) | | 2.328357 × 10 4 1.521995 × 10 22 5.744550 × 10 132 -
ACOC--5.992229-
| | x ( k + 1 ) x ( k ) | | 1.558340 × 10 3 2.184297 × 10 17 3.479449 × 10 101 -
LK4 | | f ( x ( k + 1 ) ) | | 2.171191 × 10 3 2.361842 × 10 17 4.473893 × 10 101 -
ACOC--6.048919-
| | x ( k + 1 ) x ( k ) | | 6.447430 × 10 3 3.754373 × 10 13 6.253864 × 10 75 -
LK5 | | f ( x ( k + 1 ) ) | | 1.431032 × 10 2 3.975474 × 10 13 9.442193 × 10 75 -
ACOC--6.036081-
| | x ( k + 1 ) x ( k ) | | 1.266127 × 10 1 1.334229 × 10 5 1.581618 × 10 27 8.488895 × 10 136
EM5 | | f ( x ( k + 1 ) ) | | 2.066671 × 10 1 1.227262 × 10 5 3.346691 × 10 27 1.247900 × 10 135
ACOC--5.5128914.937955
| | x ( k + 1 ) x ( k ) | | 4.042375 × 10 3 1.704335 × 10 12 6.648600 × 10 61 -
EM6 | | f ( x ( k + 1 ) ) | | 9.149067 × 10 3 1.610892 × 10 12 1.329968 × 10 60 -
ACOC--5.163562-
| | x ( k + 1 ) x ( k ) | | 1.664555 × 10 3 4.117280 × 10 15 4.812524 × 10 74 -
EM7 | | f ( x ( k + 1 ) ) | | 3.401916 × 10 3 5.017450 × 10 15 8.231433 × 10 74 -
ACOC--5.077437-
| | x ( k + 1 ) x ( k ) | | 1.159241 × 10 2 1.598320 × 10 10 3.471780 × 10 51 -
1LK6 | | f ( x ( k + 1 ) ) | | 2.345525 × 10 2 1.432171 × 10 10 5.191559 × 10 51 -
ACOC--5.173088-
| | x ( k + 1 ) x ( k ) | | 3.541294 × 10 1 9.084211 × 10 3 1.371646 × 10 11 1.161833 × 10 54
LK7 | | f ( x ( k + 1 ) ) | | 3.847192 × 10 1 9.088413 × 10 3 3.418907 × 10 11 1.280641 × 10 54
ACOC--5.5447764.882879
| | x ( k + 1 ) x ( k ) | | 3.163740 × 10 1 1.064594 × 10 4 4.230722 × 10 21 6.046883 × 10 104
LK8 | | f ( x ( k + 1 ) ) | | 1.573963 × 10 4 5.304441 × 10 21 9.960412 × 10 104 2.258926 × 10 154
ACOC--4.7223415.051280
| | x ( k + 1 ) x ( k ) | | 3.157359 × 10 1 7.508417 × 10 4 3.371253 × 10 16 2.325497 × 10 79
LK9 | | f ( x ( k + 1 ) ) | | 1.530191 × 10 3 3.643990 × 10 16 4.154766 × 10 79 3.237105 × 10 154
ACOC--4.7061025.115202
| | x ( k + 1 ) x ( k ) | | 3.156074 × 10 1 9.353115 × 10 4 1.945715 × 10 15 1.290117 × 10 75
LK10 | | f ( x ( k + 1 ) ) | | 2.244598 × 10 3 1.990170 × 10 15 2.732454 × 10 75 1.406960 × 10 154
ACOC--4.6206465.151436
Table 6. Convergence results of Test Example 2.
Table 6. Convergence results of Test Example 2.
γ MTConvk
1234
| | x ( k + 1 ) x ( k ) | | 1.693893 × 10 1 1.149094 × 10 3 3.549501 × 10 18 3.493166 × 10 105
EM1 | | f ( x ( k + 1 ) ) | | 2.056367 × 10 3 6.362117 × 10 18 6.267361 × 10 105 1.065264 × 10 153
ACOC--6.6912555.996265
| | x ( k + 1 ) x ( k ) | | 1.692982 × 10 1 1.428659 × 10 3 2.746116 × 10 17 1.538143 × 10 99
EM2 | | f ( x ( k + 1 ) ) | | 2.583199 × 10 3 4.955568 × 10 17 2.773474 × 10 99 1.214129 × 10 153
ACOC--6.6142895.996680
| | x ( k + 1 ) x ( k ) | | 1.690899 × 10 1 2.112975 × 10 3 6.934014 × 10 16 9.702901 × 10 91
EM3 | | f ( x ( k + 1 ) ) | | 3.845543 × 10 3 1.256302 × 10 15 1.755204 × 10 90 3.803022 × 10 154
ACOC--6.5593515.996044
| | x ( k + 1 ) x ( k ) | | 1.691984 × 10 1 1.784083 × 10 3 1.665966 × 10 16 1.251394 × 10 94
EM4 | | f ( x ( k + 1 ) ) | | 3.226540 × 10 3 3.003719 × 10 16 2.255795 × 10 94 1.718658 × 10 153
ACOC--6.5907295.995838
| | x ( k + 1 ) x ( k ) | | 1.696836 × 10 1 5.056251 × 10 4 5.617116 × 10 21 1.067837 × 10 122
2 / 3 LK1 | | f ( x ( k + 1 ) ) | | 9.769234 × 10 4 1.081373 × 10 20 2.049789 × 10 122 1.022636 × 10 153
ACOC--6.7124235.999711
| | x ( k + 1 ) x ( k ) | | 1.690923 × 10 1 2.505552 × 10 3 1.573370 × 10 15 1.203280 × 10 88
LK2 | | f ( x ( k + 1 ) ) | | 4.438132 × 10 3 2.754313 × 10 15 2.113408 × 10 88 1.214129 × 10 153
ACOC--6.6706405.992134
| | x ( k + 1 ) x ( k ) | | 1.692117 × 10 1 1.691344 × 10 3 1.130119 × 10 16 1.105887 × 10 95
LK3 | | f ( x ( k + 1 ) ) | | 3.075005 × 10 3 2.047404 × 10 16 1.999450 × 10 95 5.273843 × 10 154
ACOC--6.5869005.996870
| | x ( k + 1 ) x ( k ) | | 1.689594 × 10 1 2.958354 × 10 3 6.486550 × 10 15 8.973645 × 10 85
LK4 | | f ( x ( k + 1 ) ) | | 5.258560 × 10 3 1.141154 × 10 14 1.584002 × 10 84 1.162640 × 10 153
ACOC--6.6367745.991836
| | x ( k + 1 ) x ( k ) | | 1.685364 × 10 1 6.414593 × 10 3 1.328419 × 10 12 1.350065 × 10 70
LK5 | | f ( x ( k + 1 ) ) | | 1.156979 × 10 2 2.297856 × 10 12 2.303747 × 10 70 7.963328 × 10 154
ACOC--6.8218855.988638
| | x ( k + 1 ) x ( k ) | | 1.692289 × 10 1 7.728531 × 10 4 1.465253 × 10 18 8.451602 × 10 110
EM5 | | f ( x ( k + 1 ) ) | | 2.749360 × 10 3 3.407724 × 10 18 5.285841 × 10 109 6.712506 × 10 154
ACOC--6.2905166.197380
| | x ( k + 1 ) x ( k ) | | 1.694268 × 10 1 1.299228 × 10 3 3.070758 × 10 17 9.845686 × 10 88
EM6 | | f ( x ( k + 1 ) ) | | 2.233122 × 10 3 1.638869 × 10 16 1.972892 × 10 88 1.103733 × 10 153
ACOC--6.4418575.173325
| | x ( k + 1 ) x ( k ) | | 1.694807 × 10 1 8.695880 × 10 4 3.551305 × 10 21 1.212608 × 10 122
EM7 | | f ( x ( k + 1 ) ) | | 1.531744 × 10 3 2.119838 × 10 20 2.456366 × 10 122 9.164962 × 10 154
ACOC--7.5940585.835131
| | x ( k + 1 ) x ( k ) | | 1.702183 × 10 1 1.263786 × 10 3 8.905640 × 10 17 7.142410 × 10 87
1LK6 | | f ( x ( k + 1 ) ) | | 2.510417 × 10 3 3.912749 × 10 16 1.286062 × 10 87 1.022636 × 10 153
ACOC--6.1765875.329667
| | x ( k + 1 ) x ( k ) | | 1.694697 × 10 1 1.047665 × 10 3 1.817552 × 10 17 1.353168 × 10 88
LK7 | | f ( x ( k + 1 ) ) | | 1.782079 × 10 3 8.463587 × 10 17 2.405609 × 10 88 1.179267 × 10 153
ACOC--6.2297645.168920
| | x ( k + 1 ) x ( k ) | | 1.695517 × 10 1 7.818713 × 10 4 2.349799 × 10 21 2.572134 × 10 124
LK8 | | f ( x ( k + 1 ) ) | | 1.341414 × 10 3 1.199165 × 10 20 4.782437 × 10 124 1.118751 × 10 153
ACOC--7.5003645.876048
| | x ( k + 1 ) x ( k ) | | 1.692138 × 10 1 1.711949 × 10 3 7.627769 × 10 17 2.261434 × 10 88
LK9 | | f ( x ( k + 1 ) ) | | 3.071484 × 10 3 5.523109 × 10 16 6.890397 × 10 88 1.214129 × 10 153
ACOC--6.6924615.357465
| | x ( k + 1 ) x ( k ) | | 1.694361 × 10 1 1.263022 × 10 3 3.859268 × 10 17 4.312641 × 10 88
LK10 | | f ( x ( k + 1 ) ) | | 2.163494 × 10 3 1.989944 × 10 16 8.416084 × 10 88 5.273843 × 10 154
ACOC--6.3521985.249889
Table 7. Convergence results of Test Example 3.
Table 7. Convergence results of Test Example 3.
γ MT k
123
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 6.736388 × 10 13 4.536504 × 10 82
EM1 | | f ( x ( k + 1 ) ) | | 9.957324 × 10 13 7.141684 × 10 82 1.826912 × 10 154
ACOC--6.004539
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 1.267438 × 10 12 4.035285 × 10 80
EM2 | | f ( x ( k + 1 ) ) | | 1.767457 × 10 12 6.182344 × 10 80 3.498271 × 10 154
ACOC--6.002191
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 2.766996 × 10 12 1.050932 × 10 77
EM3 | | f ( x ( k + 1 ) ) | | 3.768017 × 10 12 1.581912 × 10 77 2.473651 × 10 154
ACOC--5.998396
| | x ( k + 1 ) x ( k ) | | 7.495256 × 10 3 1.558196 × 10 21 2.153197 × 10 133
EM4 | | f ( x ( k + 1 ) ) | | 2.256546 × 10 21 3.262086 × 10 133 1.826912 × 10 154
ACOC--5.987504
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 3.364258 × 10 13 4.136278 × 10 84
2 / 3 LK1 | | f ( x ( k + 1 ) ) | | 4.578823 × 10 13 5.464893 × 10 84 1.826912 × 10 154
ACOC--5.998447
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 3.625463 × 10 12 6.626067 × 10 77
LK2 | | f ( x ( k + 1 ) ) | | 4.879379 × 10 12 9.868594 × 10 77 3.498271 × 10 154
ACOC--6.000400
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 1.690241 × 10 12 3.277540 × 10 79
LK3 | | f ( x ( k + 1 ) ) | | 2.350565 × 10 12 5.011579 × 10 79 2.473651 × 10 154
ACOC--5.999111
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 4.876737 × 10 12 5.664655 × 10 76
LK4 | | f ( x ( k + 1 ) ) | | 6.554988 × 10 12 8.425112 × 10 76 4.775668 × 10 154
ACOC--5.997538
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 1.633145 × 10 11 2.761720 × 10 72
LK5 | | f ( x ( k + 1 ) ) | | 2.189664 × 10 11 4.097352 × 10 72 2.790655 × 10 154
ACOC--5.996055
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 1.263880 × 10 13 2.738954 × 10 87
EM5 | | f ( x ( k + 1 ) ) | | 3.328247 × 10 13 5.839804 × 10 87 2.473651 × 10 154
ACOC--6.015056
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 1.440429 × 10 12 3.815658 × 10 75
EM6 | | f ( x ( k + 1 ) ) | | 4.931396 × 10 12 5.169287 × 10 75 2.790655 × 10 154
ACOC--5.592300
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 9.195040 × 10 14 9.072776 × 10 88
EM7 | | f ( x ( k + 1 ) ) | | 3.009661 × 10 13 1.315517 × 10 87 2.473651 × 10 154
ACOC--5.975547
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 1.540472 × 10 12 4.642753 × 10 74
1LK6 | | f ( x ( k + 1 ) ) | | 5.029171 × 10 12 6.674639 × 10 74 2.790655 × 10 154
ACOC--5.503817
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 2.772391 × 10 13 2.131801 × 10 78
LK7 | | f ( x ( k + 1 ) ) | | 9.686043 × 10 13 2.828895 × 10 78 4.775668 × 10 154
ACOC--5.461376
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 8.798181 × 10 14 7.810164 × 10 88
LK8 | | f ( x ( k + 1 ) ) | | 2.973773 × 10 13 1.085267 × 10 87 2.473651 × 10 154
ACOC--5.970019
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 1.521759 × 10 11 4.464500 × 10 75
LK9 | | f ( x ( k + 1 ) ) | | 5.00431 × 10 11 6.355336 × 10 75 4.284490 × 10 154
ACOC--6.238981
| | x ( k + 1 ) x ( k ) | | 2.230120 × 10 1 3.140371 × 10 12 4.176009 × 10 74
LK10 | | f ( x ( k + 1 ) ) | | 7.346393 × 10 12 5.641973 × 10 74 2.473651 × 10 154
ACOC--5.693155
Table 8. Convergence results of Test Example 4.
Table 8. Convergence results of Test Example 4.
γ MT k
123
| | x ( k + 1 ) x ( k ) | | 7.433286 × 10 1 1.760138 × 10 5 1.064249 × 10 32
EM1 | | f ( x ( k + 1 ) ) | | 4.777837 × 10 5 2.888862 × 10 32 1.826912 × 10 154
ACOC--5.884275
| | x ( k + 1 ) x ( k ) | | 7.433279 × 10 1 1.832336 × 10 5 1.699619 × 10 32
EM2 | | f ( x ( k + 1 ) ) | | 4.973815 × 10 5 4.613549 × 10 32 2.358534 × 10 154
ACOC--5.866238
| | x ( k + 1 ) x ( k ) | | 7.060120 × 10 1 2.066996 × 10 5 1.050932 × 10 32
EM3 | | f ( x ( k + 1 ) ) | | 3.968813 × 10 5 4.518912 × 10 32 2.437651 × 10 154
ACOC--6.098760
| | x ( k + 1 ) x ( k ) | | 7.433270 × 10 1 1.925233 × 10 5 2.905803 × 10 32
EM4 | | f ( x ( k + 1 ) ) | | 5.225983 × 10 5 7.887687 × 10 32 7.458340 × 10 155
ACOC--5.847611
| | x ( k + 1 ) x ( k ) | | 7.433299 × 10 1 1.638128 × 10 5 3.979469 × 10 33
2 / 3 LK1 | | f ( x ( k + 1 ) ) | | 4.446645 × 10 5 1.080211 × 10 32 8.204174 × 10 154
ACOC--5.929892
| | x ( k + 1 ) x ( k ) | | 7.433252 × 10 1 2.100907 × 10 5 6.997564 × 10 32
LK2 | | f ( x ( k + 1 ) ) | | 5.702845 × 10 5 1.899461 × 10 31 7.075603 × 10 154
ACOC--5.820793
| | x ( k + 1 ) x ( k ) | | 7.320120 × 10 1 2.190241 × 10 5 3.727540 × 10 32
LK3 | | f ( x ( k + 1 ) ) | | 7.536565 × 10 5 5.115079 × 10 31 7.743651 × 10 154
ACOC--5.917090
| | x ( k + 1 ) x ( k ) | | 6.901011 × 10 1 3.076737 × 10 5 4.664651 × 10 32
LK4 | | f ( x ( k + 1 ) ) | | 4.454988 × 10 5 1.225112 × 10 31 7.375668 × 10 154
ACOC--6.164184
| | x ( k + 1 ) x ( k ) | | 7.433217 × 10 1 1.951033 × 10 5 3.688441 × 10 32
LK5 | | f ( x ( k + 1 ) ) | | 6.653254 × 10 5 1.001212 × 10 31 1.826912 × 10 154
ACOC--5.833646
| | x ( k + 1 ) x ( k ) | | 6.230120 × 10 1 1.263880 × 10 5 7.738954 × 10 32
EM5 | | f ( x ( k + 1 ) ) | | 3.328247 × 10 5 5.839804 × 10 32 2.473651 × 10 154
ACOC--5.585807
| | x ( k + 1 ) x ( k ) | | 6.621130 × 10 1 1.440429 × 10 5 1.815658 × 10 31
EM6 | | f ( x ( k + 1 ) ) | | 4.931396 × 10 5 5.169287 × 10 31 2.790655 × 10 154
ACOC--5.554915
| | x ( k + 1 ) x ( k ) | | 5.930120 × 10 1 1.395040 × 10 5 1.072776 × 10 31
EM7 | | f ( x ( k + 1 ) ) | | 3.009661 × 10 5 1.315517 × 10 31 2.473651 × 10 154
ACOC--5.642046
| | x ( k + 1 ) x ( k ) | | 6.212301 × 10 1 2.450274 × 10 5 3.642573 × 10 31
1LK6 | | f ( x ( k + 1 ) ) | | 5.029171 × 10 5 6.674639 × 10 32 2.790655 × 10 154
ACOC--5.864574
| | x ( k + 1 ) x ( k ) | | 7.433602 × 10 1 1.398344 × 10 5 7.304552 × 10 34
LK7 | | f ( x ( k + 1 ) ) | | 3.795743 × 10 5 1.982791 × 10 33 7.420955 × 10 154
ACOC--5.984871
| | x ( k + 1 ) x ( k ) | | 6.692451 × 10 1 2.819751 × 10 5 5.180461 × 10 32
LK8 | | f ( x ( k + 1 ) ) | | 2.973773 × 10 5 1.085267 × 10 32 2.473651 × 10 154
ACOC--6.110527
| | x ( k + 1 ) x ( k ) | | 3.931569 × 10 1 1.115921 × 10 5 5.460045 × 10 32
LK9 | | f ( x ( k + 1 ) ) | | 5.00431 × 10 5 6.355336 × 10 32 4.284490 × 10 154
ACOC--5.786415
| | x ( k + 1 ) x ( k ) | | 7.433626 × 10 1 1.633627 × 10 5 1.381626 × 10 34
LK10 | | f ( x ( k + 1 ) ) | | 4.43440 × 10 5 3.750370 × 10 34 6.459112 × 10 154
ACOC--6.241405
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, M.-Y.; Kim, Y.I. Development of a Family of Jarratt-Like Sixth-Order Iterative Methods for Solving Nonlinear Systems with Their Basins of Attraction. Algorithms 2020, 13, 303. https://doi.org/10.3390/a13110303

AMA Style

Lee M-Y, Kim YI. Development of a Family of Jarratt-Like Sixth-Order Iterative Methods for Solving Nonlinear Systems with Their Basins of Attraction. Algorithms. 2020; 13(11):303. https://doi.org/10.3390/a13110303

Chicago/Turabian Style

Lee, Min-Young, and Young Ik Kim. 2020. "Development of a Family of Jarratt-Like Sixth-Order Iterative Methods for Solving Nonlinear Systems with Their Basins of Attraction" Algorithms 13, no. 11: 303. https://doi.org/10.3390/a13110303

APA Style

Lee, M. -Y., & Kim, Y. I. (2020). Development of a Family of Jarratt-Like Sixth-Order Iterative Methods for Solving Nonlinear Systems with Their Basins of Attraction. Algorithms, 13(11), 303. https://doi.org/10.3390/a13110303

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