Next Article in Journal
Neutrosophic Triplets in Neutrosophic Rings
Next Article in Special Issue
A Seventh-Order Scheme for Computing the Generalized Drazin Inverse
Previous Article in Journal
Shift Scheduling with the Goal Programming Method: A Case Study in the Glass Industry
Previous Article in Special Issue
The Modified Inertial Iterative Algorithm for Solving Split Variational Inclusion Problem for Multi-Valued Quasi Nonexpansive Mappings with Some Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Generic Family of Optimal Sixteenth-Order Multiple-Root Finders and Their Dynamics Underlying Purely Imaginary Extraneous Fixed Points

1
Department of Applied Mathematics, Dankook University, Cheonan 330-714, Korea
2
Naval Postgraduate School, Department of Applied Mathematics, Monterey, CA 93943, USA
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(6), 562; https://doi.org/10.3390/math7060562
Submission received: 26 April 2019 / Revised: 5 June 2019 / Accepted: 18 June 2019 / Published: 20 June 2019
(This article belongs to the Special Issue Iterative Methods for Solving Nonlinear Equations and Systems)

Abstract

:
A generic family of optimal sixteenth-order multiple-root finders are theoretically developed from general settings of weight functions under the known multiplicity. Special cases of rational weight functions are considered and relevant coefficient relations are derived in such a way that all the extraneous fixed points are purely imaginary. A number of schemes are constructed based on the selection of desired free parameters among the coefficient relations. Numerical and dynamical aspects on the convergence of such schemes are explored with tabulated computational results and illustrated attractor basins. Overall conclusion is drawn along with future work on a different family of optimal root-finders.

1. Introduction

Many nonlinear equations governing real-world natural phenomena cannot be solved exactly by virtue of their intrinsic complexities. It would be certainly an important matter to discuss methods for approximating such solutions of the nonlinear equations. The most widely accepted method under general circumstances is Newton’s method, which has quadratic convergence for a simple-root and linear convergence for a multiple-root. Other higher-order root-finders have been developed by many researchers [1,2,3,4,5,6,7,8,9] with optimal convergence satisfying Kung–Traub’s conjecture [10]. Several authors [10,11,12,13,14] have proposed optimal sixteenth-order simple-root finders, although their applications to real-life problems are limited due to the high degree of their algebraic complexities. Optimal sixteenth-order multiple-root finders are hardly found in the literature to the best of our knowledge at the time of writing this paper. It is not too much to emphasize the theoretical importance of developing optimal sixteenth-order multiple root-finders as well as to apply them to numerically solve real-world nonlinear problems.
In order to develop an optimal sixteenth-order multiple-root finders, we pursue a family of iterative methods equipped with generic weight functions of the form:
y n = x n m f ( x n ) f ( x n ) , z n = y n m Q f ( s ) f ( y n ) f ( x n ) = x n m 1 + s Q f ( s ) f ( x n ) f ( x n ) , w n = z n m K f ( s , u ) f ( z n ) f ( x n ) = x n m 1 + s Q f ( s ) + s u K f ( s , u ) f ( x n ) f ( x n ) , x n + 1 = w n m J f ( s , u , v ) f ( w n ) f ( x n ) = x n m 1 + s Q f ( s ) + s u K f ( s , u ) + s u v J f ( s , u , v ) f ( x n ) f ( x n ) ,
where s = f ( y n ) f ( x n ) 1 / m , u = f ( z n ) f ( y n ) 1 / m , v = f ( w n ) f ( z n ) 1 / m ; Q f : C C is analytic [15] in a neighborhood of 0, K f : C 2 C holomorphic [16,17] in a neighborhood of ( 0 , 0 ) , and J f : C 3 C holomorphic in a neighborhood of ( 0 , 0 , 0 ) . Since s , u and v are respectively one-to-m multiple-valued functions, their principal analytic branches [15] are considered. Hence, for instance, it is convenient to treat s as a principal root given by s = exp [ 1 m Log ( f ( y n ) f ( x n ) ) ] , with Log ( f ( y n ) f ( x n ) ) = Log | f ( y n ) f ( x n ) | + i Arg ( f ( y n ) f ( x n ) ) for π < Arg ( f ( y n ) f ( x n ) ) π ; this convention of Arg(z) for z C agrees with that of Log[z] command of Mathematica [18] to be employed later in numerical experiments.
The case for m = 1 has been recently developed by Geum–Kim–Neta [19]. Many other existing cases for m = 1 are special cases of (1) with appropriate forms of weight functions Q f , K f , and J f ; for example, the case developed in [10] uses the following weight functions:
Q f ( s ) = 1 ( 1 s ) 2 , K f ( s , u ) = 1 + ( 1 u ) s 2 ( 1 s ) 2 ( 1 u ) ( 1 s u ) 2 , J f ( s , u , v ) = 1 + 2 s u 2 ( v 1 ) + s 4 ( u 1 ) u 2 ( v 1 ) ( u v 1 ) + s 2 [ u v 1 u 3 ( v 2 1 ) ] ( 1 s ) 2 ( u 1 ) ( s u 1 ) 2 ( v 1 ) ( u v 1 ) ( s u v 1 ) 2 .
One goal of this paper is to construct a family of optimal sixteenth-order multiple-root finders by characterizing the generic forms of weight functions Q f ( s ) , K f ( s , u ) , and J f ( s , u , v ) . The other goal is to investigate the convergence behavior by exploring their numerical behavior and dynamics through basins of attractions [20] underlying the extraneous fixed points [21] when f ( z ) = ( z a ) m ( z b ) m is applied. In view of the right side of final substep of (1), we can conveniently locate extraneous fixed points from the roots of the weight function m [ 1 + s Q f ( s ) + s u K f ( s , u ) + s u v J f ( s , u , v ) ] .
A motivation undertaking this research is to investigate the local and global characters on the convergence of proposed family of methods (1). The local convergence of an iterative method for solving nonlinear equations is usually guaranteed with an initial guess taken in a sufficiently close neighborhood of the sought zero. On the other hand, effective information on its global convergence is hardly achieved under general circumstances. We can obtain useful information on the global convergence from attractor basins through which relevant dynamics is worth exploring. Especially the dynamics underlying the extraneous fixed points (to be described in Section 3) would influence the dynamical behavior of the iterative methods by the presence of possible attractive, indifferent, repulsive, and other chaotic orbits. One way of reducing such influence is to control the location of the extraneous fixed points. We prefer the location to be the imaginary axis that divides the entire complex plane into two symmetrical half-planes. The dynamics underlying the extraneous fixed points on the imaginary axis would be less influenced by the presence of the possible periodic or chaotic attractors.
The main theorem is presented in Section 2 with required constraints on weight functions, Q f , K f , and J f to achieve the convergence order of 16. Section 2 discusses special cases of rational weight functions. Section 3 extensively investigates the purely imaginary extraneous fixed points and investigates their stabilities. Section 4 presents numerical experiments as well as the relevant dynamics, while Section 5 states the overall conclusions along with the short description of future work.

2. Methods and Special Cases

A main theorem on the convergence of (1) is established here with the error equation and relationships among generic weight functions Q f ( s ) , K f ( s , u ) , and J f ( s , u , v ) :
Theorem 1.
Suppose that f : C C has a multiple root α of multiplicity m 1 and is analytic in a neighborhood of α. Let c j = m ! ( m 1 + j ) ! f ( m 1 + j ) ( α ) f ( m ) ( α ) for j = 2 , 3 , . Let x 0 be an initial guess selected in a sufficiently small region containing α. Assume L f : C C is analytic in a neighborhood of 0. Let Q i = 1 i ! d i d s i Q f ( s ) | ( s = 0 ) for 0 i 6 . Let K f : C 2 C be holomorphic in a neighborhood of ( 0 , 0 ) . Let J f : C 3 C be holomorphic in a neighborhood of ( 0 , 0 , 0 ) . Let K i j = 1 i ! j ! i + j s i u j K f ( s , u ) | ( s = 0 , u = 0 ) for 0 i 12 and 0 j 6 . Let J i j k = 1 i ! j ! k ! i + j + k s i u j v j J f ( s , u , v ) | ( s = 0 , u = 0 , v = 0 ) for 0 i 8 , 0 j 4 and 0 k 2 . If Q 0 = 1 , Q 1 = 2 , K 00 = 1 , K 10 = 2 , K 01 = 1 , K 20 = 1 + Q 2 , K 11 = 4 , K 30 = 4 + 2 Q 2 + Q 3 , J 000 = 1 , J 100 = 2 , J 200 = 1 + Q 2 , J 010 = 1 , J 110 = 4 , J 300 = 4 + 2 Q 2 + Q 3 , J 001 = 1 , J 020 = K 02 , J 210 = 1 + K 21 , J 400 = K 40 , J 101 = 2 , J 120 = 2 + K 12 , J 310 = 4 + K 31 + 2 Q 2 , J 500 = K 50 , J 011 = 2 , J 201 = 1 + Q 2 , J 030 = 1 + K 02 + K 03 , J 220 = 1 + K 21 + K 22 Q 2 , J 410 = 3 + K 40 + K 41 + Q 2 Q 4 , J 600 = K 60 , J 111 = 8 , J 301 = 4 + 2 Q 2 + Q 3 , J 130 = 4 + 2 K 02 + K 12 + K 13 , J 320 = 6 + 2 K 21 + K 31 + K 32 2 Q 2 Q 3 , J 510 = 6 + 2 K 40 + K 50 + K 51 3 Q 3 2 Q 4 Q 5 , J 700 = K 70 are fulfilled, then Scheme (1) leads to an optimal class of sixteenth-order multiple-root finders possessing the following error equation: with e n = x n α for n = 0 , 1 , 2 , ,
e n + 1 = 1 3456 m 15 c 2 ( ρ c 2 2 2 m c 3 ) β 0 c 2 4 + β 1 c 2 2 c 3 + 12 m 2 ( K 02 1 ) c 3 2 12 m 2 c 2 c 4 Ψ e n 16 + O ( e n 17 ) ,
where ρ = 9 + m 2 Q 2 , β 0 = ( 431 + 12 K 40 7 m 2 + 6 m ( 17 + Q 2 ) + 102 Q 2 24 Q 3 12 Q 4 + 6 K 21 ρ + 3 K 02 ρ 2 ) , β 1 = 12 m ( 17 + K 21 2 m + Q 2 + K 02 ρ ) , Ψ = Δ 1 c 2 8 + Δ 2 c 2 6 c 3 + Δ 3 c 2 5 c 4 + Δ 4 c 2 3 c 3 c 4 + Δ 5 c 2 4 + Δ 6 c 2 2 + Δ 7 c 3 4 + Δ 8 c 2 c 3 2 c 4 ,
Δ 1 = ( 255124 + 144 J 800 144 K 80 122577 m 23941 m 2 2199 m 3 79 m 4 + 24 K 40 ( 472 + 93 m + 5 m 2 ) 72 ( 17 + m ) Q 2 3 576 Q 3 2 + Q 3 ( 48 ( 566 + 6 K 40 117 m 7 m 2 ) 576 Q 4 ) + 24 ( 485 + 6 K 40 108 m 7 m 2 ) Q 4 144 Q 4 2 + Q 2 2 ( 36 ( 87 + 14 m + m 2 ) + 288 Q 3 + 144 Q 4 ) + Q 2 ( 18 ( 5300 + 1529 m + 172 m 2 + 7 m 3 8 K 40 ( 18 + m ) ) + 144 ( 35 + m ) Q 3 + 72 ( 29 + m ) Q 4 ) + 18 ρ 3 σ + 6 ρ ( 12 J 610 12 ( 2 K 50 + K 60 + K 61 2 Q 5 Q 6 ) + J 211 ( 12 K 40 + σ 2 ) + 2 K 21 ( J 002 σ 2 + σ 3 + 6 η 0 ) ) + ρ 2 ( 36 J 420 36 J 211 K 21 + 36 J 002 K 21 2 72 K 31 36 J 021 K 40 36 K 41 36 K 42 + 3 J 021 σ 2 + 6 K 02 ( 6 J 401 + 12 J 002 K 40 J 002 σ 2 + σ 3 ) ) + 12 J 401 σ 7 + J 002 σ 7 2 + 9 ρ 4 τ ) ,
Δ 2 = m ( 144 ( Q 2 3 2 Q 5 Q 6 ) + 288 Q 3 ( K 21 + ( 39 + 4 m ) K 02 ρ ) + 144 Q 2 2 ( 2 K 21 ( 7 + m ) K 02 ρ ) + 144 Q 4 ( K 21 + 4 ( 9 + m ) K 02 ρ ) + Q 2 ( 144 K 21 ( 58 + 5 m ) 36 ( 1529 8 K 40 + m ( 302 + 17 m ) ) 288 Q 3 144 Q 4 + 144 K 02 ( 38 + 3 m ) ρ ) 108 ρ 2 σ + 6 ( 40859 24 J 610 + 48 K 50 + 24 K 60 + 24 K 61 + 24 K 40 ( 31 + J 211 3 m ) + m ( 14864 + m ( 1933 + 88 m ) ) 2 J 211 σ 8 ) 72 ρ 3 τ 24 ρ 2 ( J 002 σ 2 6 η 0 ) 24 K 21 ( 1309 + m ( 267 + 14 m ) J 002 σ 8 + 6 η 0 ) + ρ ( 144 J 211 K 21 144 J 002 K 21 2 + 12 ( 12 J 420 + 12 ( 2 K 31 + J 021 K 40 + K 41 + K 42 ) J 021 σ 4 ) 12 K 02 ( 1781 + m ( 360 + 19 m ) 2 J 002 σ 4 + 12 η 0 ) ) ) ,
Δ 3 = 12 m 2 ( 6 ( 63 + 5 m ) Q 2 48 Q 3 24 Q 4 + 6 ( J 211 + K 21 2 J 002 K 21 ) ρ ( 1645 12 J 401 + 12 ( 1 + 2 J 002 ) K 40 + 372 m + 23 m 2 2 J 002 σ 2 ) + 3 ρ 2 σ 5 ) ,
Δ 4 = 144 m 3 ( 3 Q 2 + ( 53 J 211 + ( 1 + 2 J 002 ) K 21 + 6 m 2 J 002 ρ 2 ) ρ σ 5 ) ,
Δ 5 = 72 ρ m 3 c 5 + 12 m 2 c 3 2 ( ( 12 K 21 ( J 211 4 ( 10 + m ) ) + K 02 ( 4778 12 J 401 + m ( 990 + 52 m ) ) + 3 ( 1929 + 4 J 401 + 4 J 420 8 K 31 + 8 K 40 4 K 41 4 K 42 476 m 31 m 2 + 2 J 211 ( 43 + 5 m ) + J 021 ( 115 + 22 m + m 2 ) ) ) 6 ( 88 + 121 J 021 + 4 J 211 + 232 K 02 + 8 K 21 + ( 10 + 13 J 021 + 22 K 02 ) m + 2 J 002 ( 51 4 K 21 + 5 m ) ) Q 2 + 12 ( 1 + J 002 + 3 J 021 + 6 K 02 ) Q 2 2 + 18 ρ σ + 18 ρ 2 τ + Q 4 ( 12 ( 2 + K 02 ) + 12 η 1 ) + Q 3 ( 24 ( 2 + K 02 ) + 24 η 1 ) J 021 η 2 + η 1 ( ( 2165 12 K 40 + m ( 510 + 31 m ) ) + η 2 ) ) ,
Δ 6 = 72 m 3 ( 2 ( 1 + J 002 ) m c 4 2 2 m c 3 c 5 + c 3 3 ( 2 Q 2 ( 1 + 6 K 02 ) 2 K 02 ( 49 J 211 + 5 m + 2 J 002 ρ 3 ) + ( 85 2 J 211 2 J 230 + 4 K 12 4 K 21 + 4 J 002 ( K 21 ρ 2 ) + 2 K 22 + 2 K 23 + 11 m + 2 J 021 ρ 3 ) 4 ρ τ ) ) ,
Δ 7 = 144 m 4 ( 1 + J 002 + J 021 + J 040 K 03 K 04 + ( 2 η 1 ) K 02 + J 002 K 02 2 ) , Δ 8 = 144 m 4 ( 3 + η 1 + ( 1 2 J 002 ) K 02 ) ,
τ = J 040 J 021 K 02 + J 002 K 02 2 K 03 K 04 , ρ 2 = 17 + 2 m Q 2 , ρ 3 = 26 + K 21 3 m + 3 Q 2 , σ = J 230 J 211 K 02 2 K 12 J 021 K 21 + 2 J 002 K 02 K 21 K 22 K 23 , σ 2 = 431 + 7 m 2 6 m ( 17 + Q 2 ) 102 Q 2 + 24 Q 3 + 12 Q 4 , σ 3 = 472 + 5 m 2 + m ( 93 6 Q 2 ) 108 Q 2 + 12 Q 3 + 6 Q 4 , σ 4 = 890 + 13 m 2 231 Q 2 + 6 Q 2 2 3 m ( 69 + 7 Q 2 ) + 24 Q 3 + 12 Q 4 , σ 5 = J 021 + K 02 2 J 002 K 02 , σ 6 = 1255 + 6 K 40 288 m 17 m 2 + 363 Q 2 + 39 m Q 2 18 Q 2 2 12 Q 3 6 Q 4 , σ 7 = 431 12 K 40 + 7 m 2 6 m ( 17 + Q 2 ) 102 Q 2 + 24 Q 3 + 12 Q 4 , σ 8 = 1349 + 19 m 2 + m ( 312 36 Q 2 ) 360 Q 2 + 12 Q 2 2 + 24 Q 3 + 12 Q 4 , η 0 = J 401 + 2 J 002 K 40 , η 1 = 2 J 002 + J 021 , η 2 = 6 K 21 2 6 K 21 ( 43 + 5 m ) + 2 σ 6 K 02 .
Proof. 
Since Scheme (1) employs five functional evaluations, namely, f ( x n ) , f ( x n ) , f ( y n ) , f ( z n ) , and f ( w n ) , optimality can be achieved if the corresponding convergence order is 16. In order to induce the desired order of convergence, we begin by the 16th-order Taylor series expansion of f ( x n ) about α :
f ( x n ) = f ( α ) m ! e n m { 1 + i = 2 17 c i e n i 1 + O ( e n 17 ) } .
It follows that
f ( x n ) = f ( α ) ( m 1 ) ! e n m 1 { 1 + i = 2 16 i m + i 1 m c i e n i 1 + O ( e n 16 ) } .
For brevity of notation, we abbreviate e n as e. Using Mathematica [18], we find:
y n = x n m f ( x n ) f ( x n ) = α + c 2 m e 2 + ( ( m + 1 ) c 2 2 + 2 m c 3 ) m 2 e 3 + Y 4 m 3 e 4 + i = 5 16 Y i m i 1 e n i + O ( e 17 ) ,
where Y 4 = ( 1 + m ) 2 c 2 3 m ( 4 + 3 m ) c 2 c 3 + 3 m 2 c 4 and Y i = Y i ( c 2 , c 3 , , c 16 ) for 5 i 16 .
After a lengthy computation using the fact that f ( y n ) = f ( x n ) | e n ( y n α ) , we get:
s = f ( y n ) f ( x n ) 1 / m = c 2 m e + ( ( m + 2 ) c 2 2 + 2 m c 3 ) m 2 e 2 + γ 3 2 m 3 e 3 + i = 4 15 E i e i + O ( e 16 ) ,
where γ 3 = ( 7 + 7 m + 2 m 2 ) c 2 3 2 m ( 7 + 3 m ) c 2 c 3 + 6 m 2 c 4 , E i = E i ( c 2 , c 3 , , c 16 ) for 4 i 15 .
In the third substep of Scheme (1), w n = O ( e 8 ) can be achieved based on Kung–Traub’s conjecture. To reflect the effect on w n from z n in the second substep, we need to expand z n up to eighth-order terms; hence, we carry out a sixth-order Taylor expansion of Q f ( s ) about 0 by noting that s = O ( e ) and f ( y n ) f ( x n ) = O ( e 2 ) :
Q f ( s ) = Q 0 + Q 1 s + Q 2 s 2 + Q 3 s 3 + Q 4 s 4 + Q 5 s 5 + Q 6 s 6 + O ( e 7 ) ,
where Q j = 1 j ! d j d s j Q f ( s ) for 0 j 6 . As a result, we come up with:
z n = x n m Q f ( s ) f ( y n ) f ( x n ) = α + ( 1 Q 0 ) m e 2 + μ 3 m 2 e 3 + i = 4 16 W i e i + O ( e 17 ) ,
where μ 3 = ( 1 + m ( Q 0 1 ) + 3 Q 0 Q 1 ) c 2 2 2 m ( Q 0 1 ) c 3 and W i = W i ( c 2 , c 3 , , c 16 , Q 0 , , Q 6 ) for 4 i 16 . Selecting Q 0 = 1 and Q 1 = 2 leads us to an expression:
z n = α + c 2 ( ρ c 2 2 2 m c 3 ) m 2 e 4 + i = 5 16 W i e i + O ( e 17 ) .
By a lengthy computation using the fact that f ( z n ) = f ( x n ) | e n ( z n α ) , we deduce:
u = f ( z n ) f ( y n ) 1 / m = ( ρ c 2 2 2 m c 3 ) 2 m 2 e 2 + δ 3 3 m 3 e 3 + i = 4 16 G i e i + O ( e 17 ) ,
where δ 3 = ( 49 + 2 m 2 + m ( 27 6 Q 2 ) 18 Q 2 + 3 Q 3 ) c 2 3 6 m ρ c 2 c 3 + 6 m 2 c 4 and G i = G i ( c 2 , c 3 , , c 16 , Q 2 , , Q 6 ) for 4 i 16 .
In the last substep of Scheme (1), x n + 1 = O ( e 16 ) can be achieved based on Kung-Traub’s conjecture. To reflect the effect on x n + 1 from w n in the third substep, we need to expand w n up to sixteenth-order terms; hence, we carry out a 12th-order Taylor expansion of K f ( s , u ) about ( 0 , 0 ) by noting that: s = O ( e ) , u = O ( e 2 ) and f ( z n ) f ( x n ) = O ( e 4 ) with K i j = 0 satisfying i + 2 j > 12 for all 0 i 12 , 0 j 6 :
K f ( s , u ) = K 00 + K 10 s + K 20 s 2 + K 30 s 3 + K 40 s 4 + K 50 s 5 + K 60 s 6 + K 70 s 7 + K 80 s 8 + K 90 s 9 + K 100 s 10 + K 110 s 11 + K 120 s 12 + ( K 01 + K 11 s + K 21 s 2 + K 31 s 3 + K 41 s 4 + K 51 s 5 + K 61 s 6 + K 71 s 7 + K 81 s 8 + K 91 s 9 + K 101 s 10 ) u + ( K 02 + K 12 s + K 22 s 2 + K 32 s 3 + K 42 s 4 + K 52 s 5 + K 62 s 6 + K 72 s 7 + K 82 s 8 ) u 2 + ( K 03 + K 13 s + K 23 s 2 + K 33 s 3 + K 43 s 4 + K 53 s 5 + K 63 s 6 ) u 3 + ( K 04 + K 14 s + K 24 s 2 + K 34 s 3 + K 44 s 4 ) u 4 + ( K 05 + K 15 s + K 25 s 2 ) u 5 + K 06 u 6 + O ( e 13 ) .
Substituting z n , f ( x n ) , f ( y n ) , f ( z n ) , f ( x n ) , and K f ( s , u ) into the third substep of (1) leads us to:
w n = z n m K f ( s , u ) · f ( z n ) f ( x n ) = α + ( 1 K 00 ) c 2 ( ρ c 2 2 2 m c 3 ) 2 m 3 e 4 + i = 5 16 Γ i e i + O ( e 17 ) ,
where Γ i = Γ i ( c 2 , c 3 , , c 16 , Q 2 , , Q 6 , K j ) , for 5 i 16 , 0 j 12 and 0 6 . Thus K 00 = 1 immediately annihilates the fourth-order term. Substituting K 00 = 1 into Γ 5 = 0 and solving for K 10 , we find:
K 10 = 2 .
Continuing the algebraic operations in this manner at the i-th ( 6 i 7 ) stage with known values of K j , we solve Γ i = 0 for remaining K j to find:
K 20 = 1 + Q 2 , K 01 = 1 .
Substituting K 00 = 1 , K 10 = 2 , K 20 = 1 + Q 2 , K 01 = 1 into (12) and simplifying we find:
v = f ( w n ) f ( z n ) 1 / m = β 0 c 2 4 + β 1 c 2 2 c 3 + 12 m 2 ( K 02 1 ) c 3 2 12 m 2 c 2 c 4 12 m 4 e 4 + i = 5 16 T i e i + O ( e 17 ) ,
where β 0 and β 1 are described in (3) and T i = T i ( c 2 , c 3 , , c 16 , Q 2 , , Q 6 ) for 5 i 16 .
To compute the last substep of Scheme (1), it is necessary to have an eighth-order Taylor expansion of J f ( s , u , v ) about ( 0 , 0 , 0 ) due to the fact that f ( w n ) f ( x n ) = O ( e 8 ) . It suffices to expand J f up to eighth-, fourth-, and second-order terms in s , u , v in order, by noting that s = O ( e ) , u = O ( e 2 ) , v = O ( e 4 ) with J i j k = 0 satisfying i + 2 j + 4 k > 8 for all 0 i 8 , 0 j 4 , 0 k 2 :
J f ( s , u , v ) = J 000 + J 100 s + J 200 s 2 + J 300 s 3 + J 400 s 4 + J 500 s 5 + J 600 s 6 + J 700 s 7 + J 800 s 8 + ( J 010 + J 110 s + J 210 s 2 + J 310 s 3 + J 410 s 4 + J 510 s 5 + J 610 s 6 ) u + ( J 020 + J 120 s + J 220 s 2 + J 320 s 3 + J 420 s 4 ) u 2 + ( J 030 + J 130 s + J 230 s 2 ) u 3 + J 040 u 4 + ( J 001 + J 101 s + J 201 s 2 + J 301 s 3 + J 401 s 4 + ( J 011 + J 111 s + J 211 s 2 ) u + J 021 u 2 ) v + J 002 v 2 .
Substituting w n , f ( x n ) , f ( y n ) , f ( z n ) , f ( w n ) , f ( x n ) and J f ( s , u , v ) in (1), we arrive at:
x n + 1 = w n m J f ( s , u , v ) · f ( w n ) f ( x n ) = α + ϕ e 8 + i = 9 16 Ω i e i + O ( e 17 ) ,
where ϕ = 1 24 m 7 ( 1 J 000 ) c 2 ( ρ c 2 2 2 m c 3 ) β 0 c 2 4 + β 1 c 2 2 c 3 + 12 m 2 ( K 02 1 ) c 3 2 12 m 2 c 2 c 4 and Ω i = Ω i ( c 2 , c 3 , , c 16 , Q 2 , , Q 6 , K δ θ , J j k ), for 9 i 16 , 0 δ 12 , 0 θ 6 , 0 j 8 , 0 k 4 , 0 2 .
Since J 000 = 1 makes ϕ = 0 , we substitute J 000 = 1 into Ω 9 = 0 and solve for J 100 to find:
J 100 = 2 .
Continuing the algebraic operations in the same manner at the i-th ( 10 i 15 ) stage with known values of J j k , we solve Ω i = 0 for remaining J j k to find:
J 200 = 1 + Q 2 , J 010 = 1 , J 110 = 4 , J 300 = 4 + 2 Q 2 + Q 3 , J 001 = 1 , J 020 = K 02 , J 210 = 1 + K 21 , J 400 = K 40 , J 101 = 2 , J 120 = 2 + K 12 , J 310 = 4 + K 31 + 2 Q 2 , J 500 = K 50 , J 011 = 2 , J 201 = 1 + Q 2 , J 111 = 8 , J 030 = 1 + K 02 + K 03 , J 220 = 1 + K 21 + K 22 Q 2 , J 410 = 3 + K 40 + K 41 + Q 2 Q 4 , J 301 = 4 + 2 Q 2 + Q 3 , J 130 = 4 + 2 K 02 + K 12 + K 13 , J 320 = 6 + 2 K 21 + K 31 + K 32 2 Q 2 Q 3 , J 600 = K 60 , J 510 = 6 + 2 K 40 + K 50 + K 51 3 Q 3 2 Q 4 Q 5 , J 700 = K 70 .
Upon substituting Relation (19) into Ω 16 , we finally obtain:
Ω 16 = 1 3456 m 15 c 2 ( ρ c 2 2 2 m c 3 ) β 0 c 2 4 + β 1 c 2 2 c 3 + 12 m 2 ( K 02 1 ) c 3 2 12 m 2 c 2 c 4 Ψ ,
where ρ , β 0 , β 1 , and Ψ as described in (3). This completes the proof. ☐
Remark 1.
Theorem 1 clearly reflects the case for m = 1 with the same constraints on weight functions Q f , K f , J f studied in [19].

Special Cases of Weight Functions

Theorem 1 enables us to obtain Q f ( s ) , K f ( s , u ) , and J f ( s , u , v ) by means of Taylor polynomials:
Q f ( s ) = 1 + 2 s + Q 2 s 2 + Q 3 s 3 + Q 4 s 4 + Q 5 s 5 + Q 6 s 6 + O ( e 7 ) , K f ( s , u ) = 1 + 2 s + ( 1 + Q 2 ) s 2 + ( 2 Q 2 + Q 3 4 ) s 3 + K 40 s 4 + K 50 s 5 + K 60 s 6 + K 70 s 7 + K 80 s 8 + K 90 s 9 + K 100 s 10 + K 110 s 11 + K 120 s 12 + ( 1 + 4 s + K 21 s 2 + K 31 s 3 + K 41 s 4 + K 51 s 5 + K 61 s 6 + K 71 s 7 + K 81 s 8 + K 91 s 9 + K 101 s 10 ) u + ( K 02 + K 12 s + K 22 s 2 + K 32 s 3 + K 42 s 4 + K 52 s 5 + K 6 s 6 + K 72 s 7 + K 82 s 8 ) u 2 + ( K 03 + K 13 s + K 23 s 2 + K 33 s 3 + K 43 s 4 + K 53 s 5 + K 63 s 6 ) u 3 + ( K 04 + K 14 s + K 24 s 2 + K 34 s 3 + K 44 s 4 ) u 4 + ( K 05 + K 15 s + K 25 s 2 ) u 5 + K 06 u 6 + O ( e 13 ) , J f ( s , u , v ) = 1 + 2 s + ( 1 + Q 2 ) s 2 + ( 2 Q 2 + Q 3 4 ) s 3 + K 40 s 4 + K 50 s 5 + K 60 s 6 + K 70 s 7 + J 800 s 8 + ( 1 + 4 s + ( 1 + K 21 ) s 2 + ( K 31 + 2 Q 2 4 ) s 3 + ( K 40 + K 41 3 + Q 2 Q 4 ) s 4 + ( 2 K 40 + K 50 + K 51 + 6 3 Q 3 2 Q 4 Q 5 ) s 5 + J 610 s 6 ) u + ( K 02 + ( 2 + K 12 ) s + ( K 21 + K 22 Q 2 + 1 ) s 2 + ( 2 K 21 + K 31 + K 32 6 2 Q 2 Q 3 ) s 3 + J 420 s 4 ) u 2 + ( K 02 + K 03 1 + ( 2 K 02 + K 12 + K 13 4 ) s + J 230 s 2 ) u 3 + J 040 u 4 + ( 1 + 2 s + ( 1 + Q 2 ) s 2 + ( 2 Q 2 + Q 3 4 ) s 3 + J 401 s 4 + ( 2 + 8 s + J 211 s 2 ) u + J 021 u 2 ) v + J 002 v 2 + O ( e 9 ) ,
where parameters Q 2 Q 6 , K 40 , K 50 , K 60 , K 70 , K 80 , K 90 , K 100 , K 110 , K 120 , K 21 , K 31 , K 41 , K 51 , K 61 , K 71 , K 81 , K 91 , K 101 , K 02 , K 12 , K 22 , K 32 , K 42 , K 52 , K 62 , K 72 , K 82 , K 03 , K 13 , K 23 , K 33 , K 43 , K 53 , K 63 , K 04 , K 14 , K 24 , K 34 , K 44 , K 05 , K 15 , K 25 , K 06 and J 040 , J 002 , J 021 , J 211 , J 230 , J 401 , J 420 , J 610 , J 800 may be free.
Although various forms of weight functions Q f ( s ) , K f ( s , u ) and J f ( s , u , v ) are available, in the current study we limit ourselves to all three weight functions in the form of rational functions, leading us to possible purely imaginary extraneous fixed points when f ( z ) = ( z 2 1 ) m is employed. In the current study, we will consider two special cases described below:
The first case below will represent the best scheme, W3G7, studied in [19] only for m = 1 .
Case 1:
Q f ( s ) = 1 1 2 s , K f ( s , u ) = Q f ( s ) · ( s 1 ) 2 1 2 s u + 2 s 2 u , J f ( s , u , v ) = K f ( s , u ) · 1 + i = 1 3 q i s i + u i = 4 8 q i s i 4 + u 2 i = 9 14 q i s i 9 + u 3 i = 15 21 q i s i 15 A 1 ( s , u ) + v · ( i = 22 25 r i s i 22 + u ( r 26 + r 27 s + λ s 2 ) ,
where
q 9 = q 12 = q 13 = q 17 = q 18 = q 19 = q 20 = r 9 = r 10 = r 19 = r 20 = 0 , q 1 = 3055820263252 76497245 λ 142682111242 , q 2 = 56884034112404 + 44614515451 λ 285364222484 , q 3 = 45802209949332 44308526471 λ 142682111242 , q 4 = 3 ( 17778426888128 + 67929066997 λ ) 1426821112420 , q 5 = 2 ( 21034820227211 + 132665343294 λ ) 356705278105 , q 6 = 1589080655012451 + 134087681464 λ 142682111242 , q 7 = 2 ( 780300304419180 + 71852971399 λ ) 71341055621 , q 8 = 12288 ( 727219117761 + 128167952 λ ) 71341055621 , q 10 = 2 , q 11 = 2 ( 741727036224277 + 126275739062 λ ) 71341055621 , q 14 = 8192 ( 3964538065856 + 615849113 λ ) 71341055621 , q 15 = 8 ( 226231159891830 + 34083208621 λ ) 71341055621 , q 16 = 24 ( 908116719056544 + 136634733499 λ ) 356705278105 , q 21 = 131072 ( 918470889768 + 136352293 λ ) 356705278105 , r 1 = q 1 , r 2 = q 2 , r 3 = q 3 , r 4 = q 4 , r 5 = q 5 , r 6 = q 6 1 , r 7 = q 7 q 1 2 , r 8 = q 8 + q 3 2 , r 11 = 29558910226378916 + 5256346708371 λ 1426821112420 , r 12 = 55018830261476 109759858153 λ 142682111242 , r 13 = 25 ( 75694849962572 + 11301475999 λ ) 71341055621 , r 14 = 4096 ( 1500792372416 + 228734011 λ ) 15508925135 , r 15 = q 15 , r 16 = 43641510974266076 6354680006961 λ 713410556210 , r 17 = 2 ( 1060205894022116 + 202907726307 λ ) 71341055621 , r 18 = 2 ( 2870055173156756 + 475573395275 λ ) 71341055621 , r 21 = q 21 2 , r 22 = 1 , r 23 = q 1 , r 24 = q 2 , r 25 = q 3 , r 26 = 1 q 4 , r 27 = 2 q 1 q 5 , λ = 1353974063793787 212746858830 ,
and A 1 ( s , u ) = 1 + i = 1 3 r i s i + u i = 4 8 r i s i 4 + u 2 i = 9 14 r i s i 9 + u 3 i = 15 21 r i s i 15 .
As a second case, we will consider the following set of weight functions:
Case 2:
Q f ( s ) = 1 1 2 s , K f ( s , u ) = Q f ( s ) · ( s 1 ) 2 1 2 s u + 2 s 2 u , J f ( s , u , v ) = 1 + i = 1 3 q i s i + u i = 4 8 q i s i 4 + u 2 i = 9 14 q i s i 9 + u 3 i = 15 19 q i s i 15 A 0 ( s , u ) + v · ( i = 20 23 r i s i 20 + u i = 24 28 r i s i 24 + r 29 u 2 ) ,
where A 0 ( s , u ) = 1 + i = 1 3 r i s i + u i = 4 8 r i s i 4 + u 2 i = 9 14 r i s i 9 + u 3 i = 15 19 r i s i 15 and determination of the 48 coefficients q i , r i of J f is described below. Relationships were sought among all free parameters of J f ( s , u , v ) , giving us a simple governing equation for extraneous fixed points of the proposed family of methods (1).
To this end, we first express s , u and v for f ( z ) = ( z 2 1 ) m as follows with t = z 2 :
s = 1 4 ( 1 1 t ) , u = 1 4 · ( t 1 ) 2 ( t + 1 ) 2 , v = ( t 1 ) 4 4 ( 1 + 6 t + t 2 ) 2 .
In order to obtain a simple form of J f ( s , u , v ) , we needed to closely inspect how it is connected with K f ( s , u ) . When applying to f ( z ) = ( z 2 1 ) m , we find K f ( s , u ) with t = z 2 as shown below:
K f ( s , u ) = 4 t ( 1 + t ) t 2 + 6 t + 1 .
Using the two selected weight functions Q f , K f , we continue to determine coefficients q i , r i of J f yielding a simple governing equation for extraneous fixed points of the proposed methods when f ( z ) = ( z 2 1 ) m is applied. As a result of tedious algebraic operations reflecting the 25 constraints (with possible rank deficiency) given by (18) and (19), we find only 23 effective relations, as follows:
q 1 = 1 4 ( 8 r 23 ) , q 2 = 3 2 q 1 , q 3 = 2 + q 1 , q 4 = r 24 , q 5 = 2 q 4 r 25 , q 6 = 5 + 7 q 3 4 + 13 q 4 4 + 9 q 5 8 5 q 9 4 25 q 10 8 + 5 q 12 8 + r 8 4 5 r 11 4 5 r 12 8 , q 7 = 4 q 3 5 + 4 q 4 5 + q 5 5 2 q 6 5 2 r 8 5 , q 8 = q 4 + q 5 2 q 7 2 , q 9 = q 15 r 15 , q 10 = 2 q 4 2 q 15 + q 16 r 16 , q 11 = 6 + q 5 2 5 q 7 2 + q 9 + 2 q 10 r 8 + r 11 , r 1 = 2 + q 1 , r 2 = 2 ( 1 + q 2 ) , r 3 = 4 q 3 , r 4 = 1 + q 4 , r 5 = 2 q 3 2 q 4 + q 5 , r 6 = 1 + 2 q 3 q 4 2 q 5 + q 6 , r 7 = 2 + 5 q 3 4 q 4 2 q 5 + 6 q 7 + 2 r 8 , r 9 = q 4 + q 9 , r 10 = 2 q 5 2 q 9 + q 10 , r 20 = 1 , r 21 = 4 q 3 , r 22 = 4 ( 1 q 3 ) .
The three relations, J 500 = K 50 , J 600 = K 60 , and J 700 = K 70 give one relation r 22 = 4 ( 1 q 3 ) .
Due to 23 constraints in Relation (27), we find that 18 free parameters among 48 coefficients of J f in (24) are available. We seek relationships among the free parameters yielding purely imaginary extraneous fixed points of the proposed family of methods when f ( z ) = ( z 2 1 ) m is applied.
To this end, after substituting the 23 effective relations given by (27) into J f in (24) and by applying to f ( z ) = ( z 2 1 ) m , we can construct H ( z ) = 1 + s Q f ( s ) + s u K f ( s , u ) + s u v J f ( s , u , v ) in (1) and seek its roots for extraneous fixed points with t = z 2 :
H ( z ) = A · G ( t ) t ( 1 + t ) 2 ( 1 + 6 t + t 2 ) · W ( t ) ,
where A is a constant factor, G ( t ) = i = 0 20 g i t i , with g 0 = q 14 , g 1 = 16 2 q 12 4 q 13 8 q 14 16 q 15 + 10 q 16 4 r 8 + 4 r 11 + 2 r 12 16 14 4 r 15 10 r 16 + 3 r 20 + 60 r 21 + 10 r 22 , g i = g i ( q 12 , r 12 , , r 25 ) , f o r 2 i 20 and W ( t ) = i = 0 15 w i t i , with w 0 = r 14 , w 1 = 16 r 8 + 4 r 13 5 r 14 + 4 r 25 , w i = w i ( q 12 , r 12 , , r 25 ) , for 2 i 15 . The coefficients of both polynomials, G ( t ) and W ( t ) , contain at most 18 free parameters.
We first observe that partial expressions of H ( z ) with t = z 2 , namely, 1 + s Q f ( s ) = 1 + 3 t 2 ( 1 + t ) , 1 + s Q f ( s ) + s u K f ( s , u ) = 1 + 21 t + 35 t 2 + 7 t 3 4 ( 1 + t ) ( 1 + 6 t + t 2 ) and the denominator of (28) contain factors t , ( 1 + 3 t ) , ( 1 + t ) , ( 1 + 6 t + t 2 ) , ( 1 + 21 t + 35 t 2 + 7 t 3 ) when f ( z ) = ( z 2 1 ) m is applied. With an observation of presence of such factors, we seek a special subcase in which G ( t ) may contain all the interested factors as follows:
G ( t ) = t ( 1 + 3 t ) ( 1 + t ) λ ( 1 + 6 t + t 2 ) β ( 1 + 21 t + 35 t 2 + 7 t 3 ) · ( 1 + 10 t + 5 t 2 ) ( 1 + 92 t + 134 t 2 + 28 t 3 + t 4 ) · Φ ( t ) ,
where Φ ( t ) is a polynomial of degree ( 9 ( λ + 2 β ) ) , with λ { 0 , 1 , 2 } , β { 1 , 2 , 3 } and 1 λ + β 3 ; two polynomial factors ( 1 + 10 t + 5 t 2 ) and ( 1 + 92 t + 134 t 2 + 28 t 3 + t 4 ) were found in Case 3G of the previous study done by Geum–Kim–Neta [19]. Notice that factors ( 1 + 6 t + t 2 ) , ( 1 + 21 t + 35 t 2 + 7 t 3 ) , ( 1 + 10 t + 5 t 2 ) and ( 1 + 92 t + 134 t 2 + 28 t 3 + t 4 ) of G ( t ) are all negative, i.e., the corresponding extraneous fixed points are all purely imaginary.
In fact, the degree of Φ ( t ) will be decreased by annihilating the relevant coefficients containing free parameters to make all its roots negative. We take the 6 pairs of ( λ , β ) { ( 0 , 1 ) , ( 0 , 2 ) , ( 0 , 3 ) , ( 1 , 1 ) , ( 1 , 2 ) , ( 2 , 1 ) } to form 6 subcases named as Case 2A–2F in order. The lengthy algebraic process eventually leads us to additional constraints to each subcase described below:
Case 2A: ( λ , β ) = ( 0 , 1 )
q 12 q 15 r 8 r 11 r 12 r 13 r 14 r 15 r 16 r 18 r 23 r 24 = 3 9 9 13 2 7 2 13 0 0 12 23 2 13 134 0 1 2 1 4 1 8 1 16 0 0 0 0 0 0 0 7 4 17 4 5 59 16 2 1 4 1 16 181 28 699 112 99 16 1807 224 1163 16 11 4 34 171 4 125 4 159 8 3 3 4 291 7 274 7 75 2 629 14 454 11 2 14 23 2 35 4 9 2 1 1 4 255 14 425 28 14 821 56 693 4 7 30 33 49 2 29 2 1 1 2 297 7 537 14 37 1199 28 937 2 0 1 1 3 4 1 2 0 0 9 7 15 14 1 27 28 27 2 1 2 21 4 23 4 67 16 21 8 1 4 1 16 165 28 599 112 79 16 1291 224 967 16 3 2 16 35 2 51 4 8 1 1 4 241 14 439 28 29 2 947 56 707 4 2 21 23 67 4 21 2 0 1 152 7 279 14 37 2 603 28 447 2 0 0 0 0 0 0 0 2 7 4 7 1 12 7 0 0 0 0 0 0 0 0 2 7 1 14 0 1 28 1 2 q 13 q 16 q 17 q 18 q 19 r 17 r 19 r 25 r 26 r 27 r 28 r 29 + 74 0 2207 56 1900 7 1355 14 1929 7 61 7 1987 56 1453 14 919 7 16 7 2 7
and q 14 = 0 . These 12 additional constraints q 12 , q 15 , r 8 , r 11 , r 12 , r 13 , r 14 , r 15 , r 16 , r 18 , r 23 , r 24 are expressed in terms of 12 parameters q 13 , q 16 , q 17 , q 18 , q 19 , r 17 , r 19 , r 25 , r 26 , r 27 , r 28 , r 29 that are arbitrarily free for the purely imaginary extraneous fixed points. Those 12 free parameters are chosen at our disposal. Then, using Relations (27) and (30), the desired form of J f ( s , u , v ) in (24) can be constructed.
Case 2B: ( λ , β ) = ( 0 , 2 )
q 12 q 15 r 8 r 11 r 12 r 13 r 14 r 15 r 16 r 18 r 23 r 24 r 25 r 26 = 3 9 9 13 2 7 2 0 0 1 6 1 3 178 3 0 1 2 1 4 1 8 1 16 0 0 0 0 0 7 4 17 4 5 59 16 2 1 4 1 16 1 48 161 96 917 24 11 4 34 171 4 125 4 159 8 3 3 4 5 6 25 6 668 3 11 2 14 23 2 35 4 9 2 1 1 4 1 6 61 24 317 6 7 30 33 49 2 29 2 1 1 2 0 7 4 217 0 1 1 3 4 1 2 0 0 0 1 4 5 1 2 21 4 23 4 67 16 21 8 1 4 1 16 11 48 5 96 617 24 3 2 16 35 2 51 4 8 1 1 4 2 3 5 24 455 6 2 21 23 67 4 21 2 0 1 5 6 5 12 293 3 0 0 0 0 0 0 0 1 3 4 3 4 3 0 0 0 0 0 0 0 0 1 4 3 0 0 0 0 0 0 0 1 3 2 3 44 3 0 0 0 0 0 0 0 4 3 1 3 29 3 q 13 q 16 q 17 q 18 q 19 r 17 r 19 r 27 r 28 r 29 + 10 0 87 8 76 25 2 57 1 43 8 33 2 23 0 4 16 12
and q 14 = 0 . These 14 additional constraints are expressed in terms of 10 parameters q 13 , q 16 , q 17 , q 18 , q 19 , r 17 , r 19 , r 27 , r 28 , r 29 that are arbitrarily free for the purely imaginary extraneous fixed points. Those 10 free parameters are chosen at our disposal. Then, using Relations (27) and (31), the desired form of J f ( s , u , v ) in (24) can be constructed.
Case 2C: ( λ , β ) = ( 0 , 3 )
q 12 q 15 q 16 r 8 r 11 r 12 r 13 r 14 r 15 r 16 r 18 r 23 r 24 r 25 r 26 r 27 = 3 0 1 4 1 0 0 23 16 29 4 0 1 4 1 4 3 16 0 0 1 32 21 8 0 1 3 4 1 2 0 0 1 16 21 4 7 4 3 4 1 2 1 8 1 4 1 16 15 8 33 2 11 4 35 4 23 4 23 8 3 3 4 9 20 11 2 5 2 7 4 5 2 1 1 4 9 8 51 2 7 3 2 1 2 1 1 2 29 8 119 2 0 0 0 0 0 0 3 16 1 4 1 2 1 2 1 4 0 1 4 1 16 9 8 17 2 3 2 3 2 3 4 0 1 1 4 27 8 55 2 2 2 1 0 0 1 71 16 147 4 0 0 0 0 0 0 1 4 11 0 0 0 0 0 0 1 4 3 0 0 0 0 0 0 7 4 5 0 0 0 0 0 0 4 29 0 0 0 0 0 0 13 4 29 q 13 q 17 q 18 q 19 r 17 r 19 r 28 r 29 + 57 4 1 8 1 4 9 100 19 2 99 2 3 4 13 77 2 201 4 13 4 29 64 39
and q 14 = 0 . These 16 additional constraints are expressed in terms of 8 parameters q 17 , q 18 , q 19 , r 17 , r 18 , r 19 , r 28 , r 29 that are arbitrarily free for the purely imaginary extraneous fixed points. Those 8 free parameters are chosen at our disposal. Then, using Relations (27) and (32), the desired form of J f ( s , u , v ) in (24) can be constructed.
Case 2D: ( λ , β ) = ( 1 , 1 ) , being identical with Case 2A.
Case 2E: ( λ , β ) = ( 1 , 2 ) , being identical with Case 2B.
Case 2F: ( λ , β ) = ( 2 , 1 ) ,
q 12 q 15 q 16 r 8 r 11 r 12 r 13 r 14 r 15 r 16 r 18 r 23 r 24 = 3 0 1 4 1 0 0 13 12 23 2 13 134 0 1 4 1 4 3 16 0 0 0 0 0 0 0 0 1 3 4 1 2 0 0 0 0 0 0 0 7 4 3 4 1 2 1 8 1 4 1 16 181 28 699 112 99 16 1807 224 1163 16 11 4 35 4 23 4 23 8 3 3 4 291 7 274 7 75 2 629 14 454 11 2 5 2 7 4 5 2 1 1 4 255 14 425 28 14 821 56 693 4 7 3 2 1 2 1 1 2 297 7 537 14 37 1199 28 937 2 0 0 0 0 0 0 9 7 15 14 1 27 28 27 2 1 2 1 2 1 4 0 1 4 1 16 165 28 599 112 79 16 1291 224 967 16 3 2 3 2 3 4 0 1 1 4 241 14 439 28 29 2 947 56 707 4 2 2 1 0 0 1 152 7 279 14 37 2 603 28 447 2 0 0 0 0 0 0 2 7 4 7 1 12 7 0 0 0 0 0 0 0 2 7 1 14 0 1 28 1 2 q 13 q 17 q 18 q 19 r 17 r 19 r 25 r 26 r 27 r 28 r 29 + 74 0 0 2207 56 1900 7 1355 14 1929 7 61 7 1987 56 1453 14 919 7 16 7 2 7
and q 14 = 0 . These 13 additional constraints are expressed in terms of 11 parameters q 13 , q 17 , q 18 , q 19 , r 17 , r 19 , r 25 , r 26 , r 27 , r 28 , r 29 that are arbitrarily free for the purely imaginary extraneous fixed points. Those 11 free parameters are chosen at our disposal. Then, using Relations (27) and (33), the desired form of J f ( s , u , v ) in (24) can be constructed. After a process of careful factorization, we find the expression for H ( z ) in (28) stated in the following lemma.
Proposition 1.
The expression H ( z ) in (28) is identical in each subcase of 2A–2F and given by a unique relation below:
H ( z ) = ( 1 + 3 t ) ( 1 + 10 t + 5 t 2 ) ( 1 + 92 t + 134 t 2 + 28 t 3 + t 4 ) 8 ( 1 + t ) ( 1 + 6 t + t 2 ) ( 1 + 28 t + 70 t 2 + 28 t 3 + t 4 ) , t = z 2 ,
despite the possibility of different coefficients in each subcase.
Proof. 
Let us write G ( t ) in (28) as G ( t ) = t ( 1 + 3 t ) · ψ 1 ( t ) · ψ 2 ( t ) · Φ ( t ) · ( 1 + t ) λ ( 1 + 6 t + t 2 ) β 1 with ψ 1 ( t ) = ( 1 + 6 t + t 2 ) ( 1 + 21 t + 35 t 2 + 7 t 3 ) and ψ 2 ( t ) = ( 1 + 10 t + 5 t 2 ) ( 1 + 92 t + 134 t 2 + 28 t 3 + t 4 ) . Then after a lengthy process of a series of factorizations with the aid of Mathematica symbolic ability, we find Φ ( t ) and W ( t ) in each subcase as follows.
(1) Case 2A: with λ = 0 and β = 1 , we get
Φ ( t ) = 2 7 ( 1 + t ) · Γ 1 ( t ) , W ( t ) = 16 7 ψ 1 ( t ) ( 1 + 28 t + 70 t 2 + 28 t 3 + t 4 ) Γ 1 ( t ) ,
where Γ 1 ( t ) = 244 + 28 q 16 + 28 q 17 + 21 q 18 + 14 q 19 36 r 25 30 r 26 28 r 27 27 r 28 378 r 29 + 14 t ( 72 + 4 q 16 + 4 q 17 + 3 q 18 + 2 q 19 4 r 25 + 2 r 27 + 4 r 28 72 r 29 ) + t 2 ( 1692 476 q 16 476 q 17 357 q 18 238 q 19 + 548 r 25 + 578 r 26 + 672 r 27 + 957 r 28 + 6006 r 29 ) + 4 t 3 ( 2288 + 196 q 16 + 196 q 17 + 147 q 18 + 98 q 19 148 r 25 100 r 26 42 r 27 6 r 28 1540 r 29 ) 7 t 4 ( 1636 + 68 q 16 + 68 q 17 + 51 q 18 + 34 q 19 + 4 r 25 + 22 r 26 + 36 r 27 + 55 r 28 + 386 r 29 ) + t 5 ( 4176 + 56 q 16 + 56 q 17 + 42 q 18 + 28 q 19 + 648 r 25 + 400 r 26 + 140 r 27 32 r 28 + 7168 r 29 ) + t 6 ( 4332 + 28 q 16 + 28 q 17 + 21 q 18 + 14 q 19 484 r 25 394 r 26 392 r 27 545 r 28 2926 r 29 ) .
(2) Case 2B: with λ = 0 and β = 2 , we get
Φ ( t ) = 2 3 ( 1 + t ) · Γ 2 ( t ) , W ( t ) = 16 3 ( 1 + 6 t + t 2 ) ψ 1 ( t ) ( 1 + 28 t + 70 t 2 + 28 t 3 + t 4 ) Γ 2 ( t ) ,
where Γ 2 ( t ) = ( 1 + 6 t + t 2 ) ( 3 ( 4 + 4 q 16 + 4 q 17 + 3 q 18 + 2 q 19 + r 28 20 r 29 ) + t ( 24 48 q 16 48 q 17 36 q 18 24 q 19 + 4 r 27 + 22 r 28 + 280 r 29 ) + 6 t 2 ( 32 + 12 q 16 + 12 q 17 + 9 q 18 + 6 q 19 + 2 r 27 + 6 r 28 16 r 29 ) 2 t 3 ( 396 + 24 q 16 + 24 q 17 + 18 q 18 + 12 q 19 + 2 r 27 + 11 r 28 + 140 r 29 ) + 3 t 4 ( 188 + 4 q 16 + 4 q 17 + 3 q 18 + 2 q 19 4 r 27 13 r 28 + 52 r 29 ) ) .
(3) Case 2C: with λ = 0 and β = 3 , we get
Φ ( t ) = 1 2 ( 1 + t ) · Γ 3 ( t ) , W ( t ) = 4 ( 1 + 6 t + t 2 ) 2 ψ 1 ( t ) ( 1 + 28 t + 70 t 2 + 28 t 3 + t 4 ) Γ 3 ( t ) ,
where Γ 3 ( t ) = ( 1 + 6 t + t 2 ) 2 ( 12 3 r 28 + 2 t ( 60 + r 28 84 r 29 ) 4 r 29 + t 2 ( 124 + r 28 + 172 r 29 ) ) .
(4) Case 2D: with λ = 1 and β = 1 , we get
Φ ( t ) = 2 7 · Γ 1 ( t ) , W ( t ) = 16 7 ψ 1 ( t ) ( 1 + 28 t + 70 t 2 + 28 t 3 + t 4 ) Γ 1 ( t ) .
(5) Case 2E: with λ = 1 and β = 2 , we get
Φ ( t ) = 2 3 · Γ 2 ( t ) , W ( t ) = 16 3 ( 1 + 6 t + t 2 ) ψ 1 ( t ) ( 1 + 28 t + 70 t 2 + 28 t 3 + t 4 ) Γ 2 ( t ) ,
(6) Case 2F: with λ = 2 and β = 1 , we get
Φ ( t ) = 2 7 · Γ 4 ( t ) , W ( t ) = 2 7 ( 1 + t ) ψ 1 ( t ) ( 1 + 28 t + 70 t 2 + 28 t 3 + t 4 ) Γ 4 ( t ) ,
where Γ 4 ( t ) = 244 + 36 r 25 + 30 r 26 + 28 r 27 + 27 r 28 + 378 r 29 + t ( 764 + 20 r 25 30 r 26 56 r 27 83 r 28 + 630 r 29 ) 2 t 2 ( 1228 + 284 r 25 + 274 r 26 + 308 r 27 + 437 r 28 + 3318 r 29 ) + 2 t 3 ( 5804 + 580 r 25 + 474 r 26 + 392 r 27 + 449 r 28 + 6398 r 29 ) t 4 ( 156 + 1132 r 25 + 794 r 26 + 532 r 27 + 513 r 28 + 10094 r 29 ) + t 5 ( 4332 + 484 r 25 + 394 r 26 + 392 r 27 + 545 r 28 + 2926 r 29 ) .
Substituting each pair of ( Φ ( t ) , W ( t ) ) into (28) yields an identical Relation (34) as desired. ☐
Remark 2.
The factorization process in the above proposition yields the additional constraints given by (30)(33) for subcases 2A–2F, after a lengthy computation. Case 2D and Case 2E are found to be identical with Case 2A and Case 2B, respectively, by direct computation.
In Table 1, we list free parameters selected for typical subcases of 2A–2F. Combining these selected free parameters with Relations (27) and (30)–(33), we can construct special iterative schemes named as W2A1, W2A2, , W2F3, W2F4. Such schemes together with W3G7 for Case 1 shall be used in Section 4 to display results on their numerical and dynamical aspects.

3. The Dynamics behind the Extraneous Fixed Points

The dynamics behind the extraneous fixed points [21] of iterative map (1) have been investigated by Stewart [20], Amat et al. [22], Argyros–Magreñan [23], Chun et al. [24], Chicharro et al. [25], Chun–Neta [26], Cordero et al. [27], Geum et al. [14,19,28,29,30], Rhee at al. [9], Magreñan [31], Neta et al. [32,33], and Scott et al. [34].
We locate a root α of a given function f ( x ) as a fixed point ξ of the iterative map R f :
x n + 1 = R f ( x n ) , n = 0 , 1 , ,
where R f is the iteration function associated with f. Typically, R f is written in the form: R f ( x n ) = x n f ( x n ) f ( x n ) H f ( x n ) , where H f is a weight function whose zeros are other fixed points ξ α called extraneous fixed points of R f . The dynamics of R f might be influenced by presence of possible attractive, indifferent, or repulsive, and other periodic or chaotic orbits underlying the extraneous fixed points. For ease of analysis, we rewrite the iterative map (41) in a more specific form:
x n + 1 = R f ( x n ) = x n m f ( x n ) f ( x n ) H f ( x n ) ,
where H f ( x n ) = 1 + s Q f ( s ) + s u K f ( s , u ) + s u v J f ( s , u , v ) can be regarded as a weight function in the classical modified Newton’s method for a multiple root of integer multiplicity m. Notice that α is a fixed point of R f , while ξ α for which H f ( ξ ) = 0 are extraneous fixed points of R f .
The influence of extraneous fixed points on the convergence behavior was well demonstrated for simple zeros via König functions and Schröder functions [21] applied to a class of functions { f k ( x ) = x k 1 , k 2 } . The basins of attraction may be altered due to the trapped sequence { x n } by the attractive extraneous fixed points of R f . An initial guess x 0 chosen near a desired root may converge to another unwanted remote root when repulsive or indifferent extraneous fixed points are present. These aspects of the Schröder functions were observed when applied to the same class of functions { f k ( x ) = x k 1 , k 2 } .
To simply treat dynamics underlying the extraneous fixed points of iterative map (42), we select a member f ( z ) = ( z 2 1 ) m . By a similar approach made by Chun et al. [35] and Neta et al. [33,36], we construct H f ( x n ) = s · Q f ( s ) + s · u · K f ( s , u ) + s · u · v · J f ( s , u , v ) in (42). Applying f ( z ) = ( z 2 1 ) m to H f , we find a rational function H ( z ) with t = z 2 :
H ( z ) = N ( t ) D ( t ) ,
where both D ( t ) and N ( t ) are co-prime polynomial functions of t. The underlying dynamics of the iterative map (42) can be favorably investigated on the Riemann sphere [37] with possible fixed points “0(zero)” and “∞”. As can be seen in Section 5, the relevant dynamics will be illustrated in a 6 × 6 square region centered at the origin.
Indeed, the roots t of N ( t ) provide the extraneous fixed points ξ of R f in Map (42) by the relation:
ξ = t 1 2 , if t 0 , 0 ( double   root ) , if t = 0 .

Extraneous Fixed Points and their Stability

The following proposition describes the stability of the extraneous fixed points of (42).
Proposition 2.
Let f ( z ) = ( z 2 1 ) m . Then the extraneous fixed points ξ for Case 2 discussed earlier are all found to be repulsive.
Proof. 
By direct computation of R f ( z ) with f ( z ) = ( z 2 1 ) m , we write it as with t = z 2 :
R f ( z ) = Ψ n ( t ) Ψ d ( t ) ,
where Ψ n ( t ) = ( 1 + t ) 15 and Ψ d ( t ) = 16 t ( 1 + t ) 2 ( 1 + 6 t + t 2 ) 2 ( 1 + 28 t + 70 t 2 + 28 t 3 + t 4 ) 2 . With the help of Mathematica, we are able to express Ψ n ( t ) = 1 61509375 ( 1 + 3 t ) ( 1 + 10 t + 5 t 2 ) ( 1 + 92 t + 134 t 2 + 28 t 3 + t 4 ) · Q n ( t ) 2097152 · R n ( t ) and Ψ d ( t ) = 1 61509375 16 ( 1 + 3 t ) ( 1 + 10 t + 5 t 2 ) ( 1 + 92 t + 134 t 2 + 28 t 3 + t 4 ) · Q d ( t ) 131072 · R d ( t ) , with Q n ( t ) and Q d ( t ) as six- and eight-degree polynomials, while R n ( t ) = ( 327 , 923 , 929 , 643 + 34 , 417 , 198 , 067 , 010 t + 446 , 061 , 306 , 116 , 505 t 2 + 1621107643125740 t 3 + 2 , 036 , 953 , 856 , 667 , 405 t 4 + 892 , 731 , 761 , 917 , 554 t 5 + 108 , 873 , 731 , 877 , 775 t 6 ) and R d ( t ) = ( 327 , 923 , 929 , 643 + 34417198067010 t + 446 , 061 , 306 , 116 , 505 t 2 + 1621107643125740 t 3 + 2 , 036 , 953 , 856 , 667 , 405 t 4 + 892 , 731 , 761 , 917 , 554 t 5 + 108 , 873 , 731 , 877 , 775 t 6 ) . Further, we express R n ( t ) = ( 1 + 10 t + 5 t 2 ) Q ν ( t ) + R ν ( t ) and R d ( t ) = ( 1 + 10 t + 5 t 2 ) Q δ ( t ) + R δ ( t ) , with R ν ( t ) = 10 , 077 , 696 25 ( 36 + 341 t ) = R δ ( t ) . Now let t = ξ 2 , then
R f ( ξ ) = 16
using the fact that ( 1 + 3 t ) ( 1 + 10 t + 5 t 2 ) ( 1 + 92 t + 134 t 2 + 28 t 3 + t 4 ) = 0 . Hence ξ for Case 2 are all found to be repulsive. ☐
Remark 3.
Although not described here in detail due to limited space, by means of a similar proof as shown in Proposition 2, extraneous fixed points ξ for Case 1 was found to be indifferent in [19].
If f ( z ) = p ( z ) is a generic polynomial other than ( z 2 1 ) m , then theoretical analysis of the relevant dynamics may not be feasible as a result of the highly increased algebraic complexity. Nevertheless, we explore the dynamics of the iterative map (42) applied to f ( z ) = p ( z ) , which is denoted by R p as follows:
z n + 1 = R p ( z n ) = z n m p ( z n ) p ( z n ) H p ( z n ) .
Basins of attraction for various polynomials are illustrated in Section 5 to observe the complicated dynamics behind the fixed points or the extraneous fixed points. The letter W was conveniently prefixed to each case number in Table 1 to symbolize a way of designating the numerical and dynamical aspects of iterative map (42).

4. Results and Discussion on Numerical and Dynamical Aspects

We first investigate numerical aspects of the local convergence of (1) with schemes W3G7 and W2A1–W2F4 for various test functions; then we explore the dynamical aspects underlying extraneous fixed points based on iterative map (45) applied to f ( z ) = ( z 2 1 ) m , whose attractor basins give useful information on the global convergence/
Results of numerical experiments are tabulated for all selected methods in Table 2, Table 3 and Table 4. Computational experiments on dynamical aspects have been illustrated through attractor basins in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7. Both numerical and dynamical aspects have strongly confirmed the desired convergence.
Throughout the computational experiments with the aid of Mathematica, $ M i n P r e c i s i o n = 400 has been assigned to maintain 400 digits of minimum number of precision. If α is not exact, then it is given by an approximate value with 416 digits of precision higher than $ M i n P r e c i s i o n .
Limited paper space allows us to list x n and α with up to 15 significant digits. We set error bound ϵ = 1 2 × 10 360 to meet | x n α | < ϵ . Due to the high-order of convergence and root multiplicity, close initial guesses have been selected to achieve a moderate number of accurate digits of the asymptotic error constants.
Methods W3G7, W2A1, W2C2 and W2F2 successfully located desired zeros of test functions F 1 F 4 :
W 3 G 7 : F 1 ( x ) = [ cos ( π x 2 ) + 2 x 2 3 π ] 4 , α 2.27312045629419 , m = 4 , W 2 A 1 : F 2 ( x ) = [ cos ( x 2 + 1 ) x log ( x 2 π + 2 ) + 1 ] 4 · ( x 2 + 1 π ) , α = π 1 , m = 5 , W 2 C 2 : F 3 ( x ) = [ sin 1 ( x 2 1 ) + 3 e x 2 x 3 ] 2 , α 0.477696831914490 , m = 2 , W 2 F 2 : F 4 ( x ) = ( x 2 + 1 ) 4 + log [ 1 + ( x 2 + 1 ) 3 ] , α = i , m = 3 , where log z ( z C ) is   a   principal   analytic   branch   with π < I m ( log z ) π .
We find that Table 2 ensures sixteenth-order convergence. The computational asymptotic error constant | e n | / | e n 1 | 16 is in agreement with the theoretical one η = lim n | e n | / | e n 1 | 16 up to 4 significant digits. The computational convergence order p n = log | e n / η | / log | e n 1 | well approaches 16.
Additional test functions in Table 3 confirm the convergence of Scheme (1). The errors | x n α | are listed in Table 4 for comparison among the listed methods W3G7 and W2A1–W2F4. In the current experiments, W3G7 has slightly better convergence for f 5 and slightly poor convergence for all other test functions than the rest of the listed methods. No specific method performs better than the other among methods W2A1–W2F4 of Case 2.
According to the definition of the asymptotic error constant η ( c i , Q f , K f , J f ) = lim n | R f ( x n ) α | / | x n α | 16 , the convergence is dependent on iterative map R f ( x n ) , f ( x ) , x 0 , α and the weight functions Q f , K f and J f . It is clear that no particular method always achieves better convergence than the others for any test functions.
The proposed family of methods (1) has efficiency index E I [38], which is 16 1 / 5 1.741101 and larger than that of Newton’s method. In general, the local convergence of iterative methods (45) is guaranteed with good initial values x 0 that are close to α. Selection of good initial values is a difficult task, depending on precision digits, error bound, and the given function f ( x ) .
The global convergence with appropriate initial values x 0 is effectively described by means of a basin of attraction that is the set of initial values leading to long-time behavior approaching the attractors under the iterative action of R f . Basins of attraction contain information about the region of convergence. A method occupying a larger region of convergence is likely to be a more robust method. A quantitative analysis will play the important role for measuring the region of convergence.
The basins of attraction, as well as the relevant statistical data, are constructed in a similar manner shown in the work of Geum–Kim–Neta [19]. Because of the high order, we take a smaller square [ 1.5 , 1.5 ] 2 and use 601 × 601 initial points uniformly distributed in the domain. Maple software has been used to perform the desired dynamics with convergence stopping criteria satisfying | x n + 1 x n | < 10 6 within the maximum number of 40 iterations. An initial point is painted with a color whose intensity measures the number of iterations converging to a root. The brighter color implies the faster convergence. The black point means that its orbit did not converge within 40 iterations.
Despite the limited space, we will explore the dynamics of all listed maps W3G7 and W2A1–W2F4, with applications to p k ( z ) , ( 1 k 7 ) through the following seven examples. In each example, we have shown dynamical planes for the convergence behavior of iterative map x n + 1 = R f ( x n ) (42) with f ( z ) = p k ( z ) by illustrating the relevant basins of attraction through Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 and displaying relevant statistical data in Table 5, Table 6 and Table 7 with colored fonts indicating best results.
Example 1.
As a first example, we have taken a quadratic polynomial raised to the power of two with all real roots:
p 1 ( z ) = ( z 2 1 ) 2 .
Clearly the roots are ± 1 . Basins of attraction for W3G7, W2A1–W2F4 are given in Figure 1. Consulting Table 5, Table 6 and Table 7, we find that the methodsW2B2andW2F4use the least number (2.71) of iterations per point on average (ANIP) followed by W2F1 with 2.72 ANIP, W2C3 with 2.73 and W2B1 with 2.74. The fastest method is W2A2 with 969.374 s followed closely by W2A3 with 990.341 s. The slowest is W2A4 with 4446.528 s. Method W2C4 has the lowest number of black points (601) and W2A4 has the highest number (78843). We will not includeW2A4in the coming examples.
Example 2.
As a second example, we have taken the same quadratic polynomial now raised to the power of three:
p 2 ( z ) = ( z 2 1 ) 3 .
The basins for the best methods are plotted in Figure 2. This is an example to demonstrate the effect of raising the multiplicity from two to three. In one case, namely W3G7, we also have m = 5 with CPU time of 4128.379 s. Based on the figure we see that W2B4, W2C4 and W2F3 were chaotic. The worst are W2B4, W2C4 and W2F3. In terms of ANIP, the best was W2A2 (3.15) followed by W2F4 (3.19) and the worst was W2B4 (3.91). The fastest was W2B3 using (2397.111 s) followed by W2F1 using 2407.158 s and the slowest was W2C4 (4690.295 s) preceded by W3G7 (2983.035 s). Four methods have the highest number of black points (617). Those were W2A1, W2B4, W2C1 and W2F2. The lowest number was 601 for W2A2, W2C2, W2C4 and W2F1.
Comparing the CPU time for the cases m = 2 and m = 3 of W3G7, we find it is about doubled. But when increasing from three to five, we only needed about 50% more.
Example 3.
In our third example, we have taken a cubic polynomial raised to the power of three:
p 3 ( z ) = ( 3 z 3 + 4 z 2 10 ) 3 .
Basins of attraction are given in Figure 3. It is clear that W2B4, W2C4 and W2F3 were too chaotic and they should be eliminated from further consideration. In terms of ANIP, the best was W2F4 (3.50) followed by W2A2 (3.52), W2C3 (3.53) and W2F1 (3.55) and the worst were W2B4 and W2C4 with 4.99 and 4.96 ANIP, respectively. The fastest was W2C3 using 2768.362 s and the slowest was W2B3 (7193.034 s). There were 13 methods with only one black point and one with two points. The highest number of black points was 101 for W2F2.
Example 4.
As a fourth example, we have taken a different cubic polynomial raised to the power of four:
p 4 ( z ) = ( z 3 z ) 4 .
The basins are given in Figure 4. We now see that W2F3 is the worst. In terms of ANIP, W2A2 and W2F1 were the best (3.84 each) and the worst was W2F3 (4.41). The fastest was W2A2 (2956.702 s) and the slowest was W2F3 (4130.158 s). The lowest number of black points (78) was for method W2F3 and the highest number (614) for W2F2. We did not include W2F3 in the rest of the experiments.
Example 5.
As a fifth example, we have taken a quintic polynomial raised to the power of three:
p 3 ( z ) = ( z 5 1 ) 3 .
The basins for the best methods left are plotted in Figure 5. The worst were W2A1 and W2C2. In terms of ANIP, the best was W2F4 (3.42) followed by W2F1 (3.49) and the worst were W2A1 (6.84) and W2C2 (6.70). The fastest was W2C3 using 2778.518 s followed by W2F1 using 2832.23 s and W2B1 using 2835.755 s. The slowest was W2A1 using 5896.635 s. There were three methods with one black point (W2A2, W2B2 and W2C3) and four others with 10 or less such points, namely W2B3 (3), W2A3 and W2B1 (9) and W2F1 (10). The highest number was for W2A1 (34,396) preceded by W2C2 with 17,843 black points.
Example 6.
As a sixth example, we have taken a quartic polynomial raised to the power of three:
p 6 ( z ) = ( z 4 1 ) 3 .
The basins for the best methods left are plotted in Figure 6. It seems that most of the methods left were good except W2B3 and W2C1. Based on Table 5 we find that W2F4 has the lowest ANIP (3.53) followed by W2F1 (3.57). The fastest method was W2A2 (2891.478 s) followed by W2C3 (2914.941 s). The slowest was W2C1 (4080.019 s) preceded by W3G7 using 3901.679 s. The lowest number of black points was for W2A1, W2A2, W2B1 and W2C3 (1201) and the highest number was for W2C1 with 18,157 black points.
Example 7.
As a seventh example, we have taken a non-polynomial equation having ± i as its triple roots:
p 6 ( z ) = ( z + i ) 3 ( e z i 1 ) 3 , w i t h i = 1 .
The basins for the best methods left are plotted in Figure 7. It seems that most of the methods left have a larger basin for the root i , i.e., the boundary does not match the real line exactly. Based on Table 5 we find that W2A2 has the lowest ANIP (4.84) followed by W2C3 (4.94) and W2A3 (4.98). The fastest method was W2A2 (2981.179 seconds) followed by W2B3 (3139.084 s), W2A3 (3155.307 s) and W2B2 (3155.619 s). The slowest was W2C1 (4802.662 s). The lowest number of black points was for W2B1 (13,946) and the highest number was for W3G7 with 33,072 black points. In general all methods had higher number of black points compared to the polynomial examples.
We now average all these results across the seven examples to try and pick the best method. W2A2 had the lowest ANIP (3.63), followed by W2C3 with 3.66, W2F4 with 3.68 and W2F1 with 3.69. The fastest method was W2A2 (2682.252 seconds), followed by W2C3 (2700.399 s) and W2A3 using 2703.600 s of CPU. W2B1 has the lowest number of black points on average (2379), followed by W2A3 (2490 black points). The highest number of black points was for W2A1.
Based on these seven examples we see that W2F4 has four examples with the lowest ANIP, W2A2 had three examples and W2F1 has one example. On average, though, W2A2 had the lowest ANIP. W2A2 was the fastest in four examples and on average. W2C3 was the fastest in two examples and W2B3 in one example. In terms of black points, W2A2, W2B1 and W2B3 had the lowest number in three examples and W2F1 in two examples. On average W2B1 has the lowest number. Thus, we recommend W2A2, since it is in the top in all categories.

5. Conclusions

Both numerical and dynamical aspects of iterative map (1) support the main theorem well through a number of test equations and examples. The W2C2 and W2B3 methods were observed to occupy relatively slower CPU time. Such dynamical aspects would be greatly strengthened if we could include a study of parameter planes with reference to appropriate parameters in Table 1.
The proposed family of methods (1) employing generic weight functions favorably cover most of optimal sixteenth-order multiple-root finders with a number of feasible weight functions. The dynamics behind the purely imaginary extraneous fixed points will choose best members of the family with improved convergence behavior. However, due to the high order of convergence, the algebraic difficulty might arise resolving its increased complexity. The current work is limited to univariate nonlinear equations; its extension to multivariate ones becomes another task.

Author Contributions

investigation, M.-Y.L.; formal analysis, Y.I.K.; supervision, B.N.

Conflicts of Interest

The authors have no conflict of interest to declare.

References

  1. 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]
  2. 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]
  3. Geum, Y.H.; Kim, Y.I. A uniparametric family of three-step eighth-order multipoint iterative methods for simple roots. Appl. Math. Lett. 2011, 24, 929–935. [Google Scholar] [CrossRef]
  4. Lee, S.D.; Kim, Y.I.; Neta, B. An optimal family of eighth-order simple-root finders with weight functions dependent on function-to-function ratios and their dynamics underlying extraneous fixed points. J. Comput. Appl. Math. 2017, 317, 31–54. [Google Scholar] [CrossRef]
  5. 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]
  6. Petković, M.S.; Neta, B.; Petković, L.D.; Džunić, J. Multipoint Methods for Solving Nonlinear Equations; Elsevier: New York, NY, USA, 2012. [Google Scholar]
  7. Petković, M.S.; Neta, B.; Petković, L.D.; Džunić, J. Multipoint methods for solving nonlinear equations: A survey. Appl. Math. Comput. 2014, 226, 635–660. [Google Scholar] [CrossRef] [Green Version]
  8. Sharma, J.R.; Arora, H. A new family of optimal eighth order methods with dynamics for nonlinear equations. Appl. Math. Comput. 2016, 273, 924–933. [Google Scholar] [CrossRef]
  9. Rhee, M.S.; Kim, Y.I.; Neta, B. An optimal eighth-order class of three-step weighted Newton’s methods and their dynamics behind the purely imaginary extraneous fixed points. Int. J. Comput. Math. 2017, 95, 2174–2211. [Google Scholar] [CrossRef]
  10. Kung, H.T.; Traub, J.F. Optimal order of one-point and multipoint iteration. J. Assoc. Comput. Mach. 1974, 21, 643–651. [Google Scholar] [CrossRef]
  11. Maroju, P.; Behl, R.; Motsa, S.S. Some novel and optimal families of King’s method with eighth and sixteenth-order of convergence. J. Comput. Appl. Math. 2017, 318, 136–148. [Google Scholar] [CrossRef]
  12. Sharma, J.R.; Argyros, I.K.; Kumar, D. On a general class of optimal order multipoint methods for solving nonlinear equations. J. Math. Anal. Appl. 2017, 449, 994–1014. [Google Scholar] [CrossRef]
  13. Neta, B. On a family of Multipoint Methods for Non-linear Equations. Int. J. Comput. Math. 1981, 9, 353–361. [Google Scholar] [CrossRef]
  14. Geum, Y.H.; Kim, Y.I.; Neta, B. Constructing a family of optimal eighth-order modified Newton-type multiple-zero finders along with the dynamics behind their purely imaginary extraneous fixed points. J. Comput. Appl. Math. 2018, 333, 131–156. [Google Scholar] [CrossRef]
  15. Ahlfors, L.V. Complex Analysis; McGraw-Hill Book, Inc.: New York, NY, USA, 1979. [Google Scholar]
  16. Hörmander, L. An Introduction to Complex Analysis in Several Variables; North-Holland Publishing Company: Amsterdam, The Netherlands, 1973. [Google Scholar]
  17. Shabat, B.V. Introduction to Complex Analysis PART II, Functions of several Variables; American Mathematical Society: Providence, RI, USA, 1992. [Google Scholar]
  18. Wolfram, S. The Mathematica Book, 5th ed.; Wolfram Media: Champaign, IL, USA, 2003. [Google Scholar]
  19. Geum, Y.H.; Kim, Y.I.; Neta, B. Developing an Optimal Class of Generic Sixteenth-Order Simple-Root Finders and Investigating Their Dynamics. Mathematics 2019, 7, 8. [Google Scholar] [CrossRef]
  20. Stewart, B.D. Attractor Basins of Various Root-Finding Methods. Master’s Thesis, Naval Postgraduate School, Department of Applied Mathematics, Monterey, CA, USA, June 2001. [Google Scholar]
  21. 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]
  22. Amat, S.; Busquier, S.; Plaza, S. Review of some iterative root-finding methods from a dynamical point of view. Scientia 2004, 10, 3–35. [Google Scholar]
  23. Argyros, I.K.; Magreñá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]
  24. Chun, C.; Lee, M.Y.; Neta, B.; Džunić, J. On optimal fourth-order iterative methods free from second derivative and their dynamics. Appl. Math. Comput. 2012, 218, 6427–6438. [Google Scholar] [CrossRef] [Green Version]
  25. Chicharro, F.; Cordero, A.; Gutiérrez, J.M.; Torregrosa, J.R. Complex dynamics of derivative-free methods for nonlinear equations. Appl. Math. Comput. 2013, 219, 7023–7035. [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]
  28. Geum, Y.H.; Kim, Y.I.; Magreñá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. Geum, Y.H.; Kim, Y.I.; Neta, B. A sixth-order family of three-point modified Newton-like multiple-root finders and the dynamics behind their extraneous fixed points. Appl. Math. Comput. 2016, 283, 120–140. [Google Scholar] [CrossRef] [Green Version]
  31. Magreñan, Á.A. Different anomalies in a Jarratt family of iterative root-finding methods. Appl. Math. Comput. 2014, 233, 29–38. [Google Scholar]
  32. Neta, B.; Scott, M.; Chun, C. Basin attractors for various methods for multiple roots. Appl. Math. Comput. 2012, 218, 5043–5066. [Google Scholar] [CrossRef]
  33. Neta, B.; Chun, C.; Scott, M. Basins of attraction for optimal eighth order methods to find simple roots of nonlinear equations. Appl. Math. Comput. 2014, 227, 567–592. [Google Scholar] [CrossRef]
  34. Scott, M.; Neta, B.; Chun, C. Basin attractors for various methods. Appl. Math. Comput. 2011, 218, 2584–2599. [Google Scholar] [CrossRef]
  35. Chun, C.; Neta, B. Basins of attraction for Zhou-Chen-Song fourth order family of methods for multiple roots. Math. Comput. Simul. 2015, 109, 74–91. [Google Scholar] [CrossRef]
  36. Neta, B.; Chun, C. Basins of attraction for several optimal fourth order methods for multiple roots. Math. Comput. Simul. 2014, 103, 39–59. [Google Scholar] [CrossRef]
  37. Beardon, A.F. Iteration of Rational Functions; Springer: New York, NY, USA, 1991. [Google Scholar]
  38. Traub, J.F. Iterative Methods for the Solution of Equations; Chelsea Publishing Company: Chelsea, VT, USA, 1982. [Google Scholar]
Figure 1. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2A4 (left), W2B1 (center left), W2B2 (center right) and W2B3 (right). The third row for W2B4 (left), W2C1 (center left), W2C2 (center right) and W2C3 (right). The third row for W2C4 (left), W2F1 (center left), W2F2 (center right), and W2F3 (right). The bottom row for W2F4 (center), for the roots of the polynomial equation ( z 2 - 1 ) 2 .
Figure 1. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2A4 (left), W2B1 (center left), W2B2 (center right) and W2B3 (right). The third row for W2B4 (left), W2C1 (center left), W2C2 (center right) and W2C3 (right). The third row for W2C4 (left), W2F1 (center left), W2F2 (center right), and W2F3 (right). The bottom row for W2F4 (center), for the roots of the polynomial equation ( z 2 - 1 ) 2 .
Mathematics 07 00562 g001
Figure 2. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2B1 (left), W2B2 (center left), W2B3 (center right) and W2B4 (right). The third row for W2C1 (left), W2C2 (center left), W2C3 (center right) and W2C4 (right). The fourth row for W2F1 (left), W2F2 (center left), W2F3 (center right), and W2F4 (right). The bottom row for W3G7m5 (center), for the roots of the polynomial equation ( z 2 - 1 ) k , k { 3 , 5 } .
Figure 2. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2B1 (left), W2B2 (center left), W2B3 (center right) and W2B4 (right). The third row for W2C1 (left), W2C2 (center left), W2C3 (center right) and W2C4 (right). The fourth row for W2F1 (left), W2F2 (center left), W2F3 (center right), and W2F4 (right). The bottom row for W3G7m5 (center), for the roots of the polynomial equation ( z 2 - 1 ) k , k { 3 , 5 } .
Mathematics 07 00562 g002
Figure 3. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2B1 (left), W2B2 (center left), W2B3 (center right) and W2B4 (right). The third row for W2C1 (left), W2C2 (center left), W2C3 (center right) and W2C4 (right). The bottom row for W2F1 (left), W2F2 (center left), W2F3 (center right), and W2F4 (right), for the roots of the polynomial equation ( 3 z 3 + 4 z 2 - 10 ) 3 .
Figure 3. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2B1 (left), W2B2 (center left), W2B3 (center right) and W2B4 (right). The third row for W2C1 (left), W2C2 (center left), W2C3 (center right) and W2C4 (right). The bottom row for W2F1 (left), W2F2 (center left), W2F3 (center right), and W2F4 (right), for the roots of the polynomial equation ( 3 z 3 + 4 z 2 - 10 ) 3 .
Mathematics 07 00562 g003
Figure 4. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2B1 (left), W2B2 (center left), W2B3 (center right) and W2C1 (right). The third row for W2C2 (left), W2C3 (center left), W2F1 (center right) and W2F2 (right). The bottom row for W2F3 (left) and W2F4 (right), for the roots of the polynomial equation ( z 3 - z ) 4 .
Figure 4. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2B1 (left), W2B2 (center left), W2B3 (center right) and W2C1 (right). The third row for W2C2 (left), W2C3 (center left), W2F1 (center right) and W2F2 (right). The bottom row for W2F3 (left) and W2F4 (right), for the roots of the polynomial equation ( z 3 - z ) 4 .
Mathematics 07 00562 g004
Figure 5. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2B1 (left), W2B2 (center left), W2B3 (center right) and W2C1 (right). The third row for W2C2 (left), W2C3 (center left), W2F1 (center right) and W2F2 (right). The bottom row for W2F4 (center), for the roots of the polynomial equation ( z 5 - 1 ) 3 .
Figure 5. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2B1 (left), W2B2 (center left), W2B3 (center right) and W2C1 (right). The third row for W2C2 (left), W2C3 (center left), W2F1 (center right) and W2F2 (right). The bottom row for W2F4 (center), for the roots of the polynomial equation ( z 5 - 1 ) 3 .
Mathematics 07 00562 g005
Figure 6. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2B1 (left), W2B2 (center left), W2B3 (center right) and W2C1 (right). The third row for W2C2 (left), W2C3 (center left), W2F1 (center right) and W2F2 (right). The bottom row for W2F4 (center), for the roots of the polynomial equation ( z 4 - 1 ) 3 .
Figure 6. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2B1 (left), W2B2 (center left), W2B3 (center right) and W2C1 (right). The third row for W2C2 (left), W2C3 (center left), W2F1 (center right) and W2F2 (right). The bottom row for W2F4 (center), for the roots of the polynomial equation ( z 4 - 1 ) 3 .
Mathematics 07 00562 g006
Figure 7. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2B1 (left), W2B2 (center left), W2B3 (center right) and W2C1 (right). The third row for W2C2 (left), W2C3 (center left), W2F1 (center right) and W2F2 (right). The bottom row for W2F4 (center), for the roots of the non-polynomial equation ( z + i ) 3 ( e z - i - 1 ) 3 .
Figure 7. The top row for W3G7 (left), W2A1 (center left), W2A2 (center right) and W2A3 (right). The second row for W2B1 (left), W2B2 (center left), W2B3 (center right) and W2C1 (right). The third row for W2C2 (left), W2C3 (center left), W2F1 (center right) and W2F2 (right). The bottom row for W2F4 (center), for the roots of the non-polynomial equation ( z + i ) 3 ( e z - i - 1 ) 3 .
Mathematics 07 00562 g007
Table 1. Free parameters selected for typical subcases of 2A1–2F4.
Table 1. Free parameters selected for typical subcases of 2A1–2F4.
SCN q 13 q 16 q 17 q 18 q 19 r 17 r 19 r 25 r 26 r 27 r 28 r 29
2A1000000000000
2A22000000101242-800
2A3 89 26 000000 149 26 0000
2A400 711 26 622 13 222 13 00 149 26 0000
2B10000000--000
2B200000096--−31−1 1 4
2B3−19000000--−1800
2B40045−52−1200--−1800
2C10-00000---00
2C2−34-0001740---40
2C30-0000280---40
2C4 39 4 - 375 4 627 4 000---00
2F10-000−400 122 3 290 3 184 3 00
2F216-0000010−1000
2F3 138 7 - 356 7 3963 14 00000000
2F40-000−320−3378−46−40
Table 2. Convergence of methods W3G7, W2A1, W2C2, W2F2 for test functions F 1 ( x ) F 4 ( x ) .
Table 2. Convergence of methods W3G7, W2A1, W2C2, W2F2 for test functions F 1 ( x ) F 4 ( x ) .
Method Fn x n | F ( x n ) | | x n α | | e n / e n 1 16 | η p n
W3G7 F 1 02.2735 1.873 × 10 10 0.000379544
12.27312045629419 1.927 × 10 233 6.798 × 10 60 0.00003666355445 0.00003666729357
22.27312045629419 0.0 × 10 400 1.004 × 10 237 16.00000
W2A1 F 2 01.4634 1.93 × 10 21 0.0000181404
11.46341814037882 3.487 × 10 366 2.040 × 10 74 148.4148965 148.4575003
21.46341814037882 0.0 × 10 400 0.0 × 10 399 16.00000
W2C2 F 3 0 0.4777 1.890 × 10 10 3.168 × 10 6
1 0.477696831914490 6.116 × 10 183 1.802 × 10 92 0.00017500020630.0001749999826
2 0.477696831914490 0.0 × 10 400 8.522 × 10 367 16.00000
W2F2 F 4 0 0.99995 i 1.000 × 10 12 0.00005
1 1.00000000000000 i 4.820 × 10 215 1.391 × 10 72 0.0010378384360.001041219259
2 1.00000000000000 i 0.0 × 10 400 0.0 × 10 400 16.00030
i = 1 , η = lim n | e n | | e n 1 | 16 , p n = log | e n / η | log | e n 1 | .
Table 3. Additional test functions f i ( x ) with zeros α and initial values x 0 and multiplicities.
Table 3. Additional test functions f i ( x ) with zeros α and initial values x 0 and multiplicities.
i f i ( x ) α x 0 m
1 [ 4 + 3 sin x 2 x 2 ] 4 1.85471014256339 1.86 4
2 [ 2 x P i + x cos x ] 5 π 2 1.5707 5
3 [ 2 x 3 + 3 e x + 4 sin ( x 2 ) 5 ] 2 0.402282449584416 0.403 2
4 [ 3 x 2 · cos π x 6 + 1 x 3 + 1 1 28 ] · ( x 3 ) 3 3 3.0005 4
5 ( x 1 ) 2 + 1 12 log [ 25 12 2 x + x 2 ] 1 i 3 6 0.99995 0.28 i 2
6 [ x log x x + x 3 ] 3 1 1.0001 3
Here, log z ( z C ) represents a principal analytic branch with π I m ( log z ) < π .
Table 4. Comparison of | x n α | among selected methods applied to various test functions.
Table 4. Comparison of | x n α | among selected methods applied to various test functions.
Method | x n α | f ( x )
f 1 f 2 f 3 f 4 f 5 f 6
W3G7 | x 1 α | 1.77 × 10 40 *1.62 × 10 57 4.89 × 10 51 1.50 × 10 61 5.76 × 10 7 1.19 × 10 62
| x 2 α | 1.02 × 10 159 1.13 × 10 225 1.24 × 10 201 3.27 × 10 245 1.08 × 10 95 2.40 × 10 247
W2A1 | x 1 α | 2.83 × 10 42 1.05 × 10 58 1.92 × 10 52 1.23 × 10 62 1.29 × 10 6 1.11 × 10 63
| x 2 α | 0.0 × 10 399 4.24 × 10 230 0.0 × 10 400 0.0 × 10 399 6.62 × 10 90 3.61 × 10 251
W2A2 | x 1 α | 1.63 × 10 41 2.33 × 10 58 1.45 × 10 51 1.05 × 10 61 1.32 × 10 6 2.34 × 10 63
| x 2 α | 0.0 × 10 399 1.11 × 10 228 0.0 × 10 400 0.0 × 10 399 2.53 × 10 89 7.39 × 10 250
W2A3 | x 1 α | 2.53 × 10 43 1.82 × 10 60 4.56 × 10 54 1.20 × 10 63 4.43 × 10 6 3.85 × 10 65
| x 2 α | 0.0 × 10 399 1.40 × 10 236 8.40 × 10 213 0.0 × 10 399 3.03 × 10 83 1.53 × 10 256
W2A4 | x 1 α | 1.24 × 10 42 1.35 × 10 59 1.53 × 10 52 8.38 × 10 63 1.58 × 10 4 1.12 × 10 64
| x 2 α | 1.70 × 10 125 5.96 × 10 424 5.22 × 10 155 8.16 × 10 187 3.34 × 10 57 2.79 × 10 424
W2B1 | x 1 α | 2.39 × 10 42 2.28 × 10 59 1.30 × 10 52 7.80 × 10 63 2.14 × 10 6 2.81 × 10 64
| x 2 α | 0.0 × 10 399 1.57 × 10 232 0.0 × 10 400 0.0 × 10 399 6.04 × 10 87 2.23 × 10 253
W2B2 | x 1 α | 4.44 × 10 42 2.73 × 10 59 3.03 × 10 52 1.79 × 10 62 4.30 × 10 6 3.29 × 10 64
| x 2 α | 0.0 × 10 399 6.69 × 10 232 0.0 × 10 400 0.0 × 10 399 6.01 × 10 82 7.78 × 10 e 253
W2B3 | x 1 α | 9.85 × 10 43 3.11 × 10 61 4.26 × 10 53 3.01 × 10 63 4.46 × 10 6 3.26 × 10 65
| x 2 α | 0.0 × 10 399 1.17 × 10 239 0.0 × 10 400 0.0 × 10 399 3.06 × 10 83 7.91 × 10 257
W2B4 | x 1 α | 1.04 × 10 42 1.92 × 10 59 1.77 × 10 52 1.12 × 10 62 1.77 × 10 4 1.53 × 10 64
| x 2 α | 1.12 × 10 125 4.68 × 10 405 9.06 × 10 155 2.23 × 10 186 2.32 × 10 56 0.0 × 10 400
W2C1 | x 1 α | 4.87 × 10 42 4.27 × 10 59 2.95 × 10 52 1.50 × 10 62 1.08 × 10 6 5.14 × 10 64
| x 2 α | 0.0 × 10 399 0.0 × 10 399 0.0 × 10 400 0.0 × 10 399 2.41 × 10 91 2.22 × 10 191
W2C2 | x 1 α | 9.31 × 10 42 1.01 × 10 58 5.94 × 10 52 2.61 × 10 62 1.47 × 10 6 1.18 × 10 63
| x 2 α | 0.0 × 10 399 4.23 × 10 230 0.0 × 10 400 0.0 × 10 399 7.11 × 10 89 4.95 × 10 251
W2C3 | x 1 α | 9.17 × 10 42 8.88 × 10 59 6.85 × 10 52 4.30 × 10 62 4.36 × 10 6 9.66 × 10 64
| x 2 α | 0.0 × 10 399 7.84 × 10 230 0.0 × 10 400 0.0 × 10 399 2.11 × 10 81 6.05 × 10 251
W2C4 | x 1 α | 5.36 × 10 42 6.72 × 10 59 6.55 × 10 52 4.38 × 10 62 1.02 × 10 4 6.13 × 10 64
| x 2 α | 9.92 × 10 124 5.96 × 10 424 2.97 × 10 153 8.55 × 10 185 1.30 × 10 59 0.0 × 10 400
W2F1 | x 1 α | 4.36 × 10 42 8.67 × 10 60 2.55 × 10 52 1.29 × 10 62 4.34 × 10 6 1.60 × 10 64
| x 2 α | 0.0 × 10 399 7.12 × 10 234 0.0 × 10 400 0.0 × 10 399 4.57 × 10 82 4.61 × 10 254
W2F2 | x 1 α | 1.20 × 10 42 1.74 × 10 60 5.52 × 10 53 3.67 × 10 63 4.08 × 10 6 5.06 × 10 65
| x 2 α | 0.0 × 10 399 1.16 × 10 236 0.0 × 10 400 0.0 × 10 399 2.33 × 10 83 4.56 × 10 256
W2F3 | x 1 α | 1.08 × 10 41 1.54 × 10 58 1.55 × 10 51 1.14 × 10 61 8.67 × 10 5 1.396 × 10 63
| x 2 α | 1.41 × 10 424 5.27 × 10 172 8.47 × 10 424 0.0 × 10 399 1.66 × 10 60 5.22 × 10 188
W2F4 | x 1 α | 3.80 × 10 42 5.01 × 10 60 2.15 × 10 52 1.07 × 10 62 4.35 × 10 6 1.18 × 10 64
| x 2 α | 0.0 × 10 399 1.16 × 10 235 0.0 × 10 400 0.0 × 10 399 3.65 × 10 82 1.38 × 10 254
Table 5. Average number of iterations per point for each example (1–7).
Table 5. Average number of iterations per point for each example (1–7).
Map Example
1: m = 22: m = 33: m = 34: m = 45: m = 36: m = 37: m = 3Average
W3G72.943.483.833.933.953.976.774.12
W2A12.843.503.704.046.843.745.494.31
W2A22.763.153.523.843.623.664.843.63
W2A32.783.213.613.893.703.744.983.70
W2A411.40-------
W2B12.743.253.703.883.673.725.013.71
W2B22.953.423.664.013.753.775.153.82
W2B32.783.283.643.893.694.655.133.86
W2B43.293.914.99-----
W2C12.883.663.874.083.895.456.254.30
W2C22.933.683.954.156.704.675.754.55
W2C32.733.223.533.983.603.624.943.66
W2C43.143.814.96-----
W2F12.723.243.553.843.493.575.413.69
W2F22.813.283.804.065.024.505.294.10
W2F32.913.544.364.41----
W2F42.713.193.503.863.423.535.523.68
Table 6. CPU time (in seconds) required for each example (1–7) using a Dell Multiplex-990.
Table 6. CPU time (in seconds) required for each example (1–7) using a Dell Multiplex-990.
MapExample
1: m = 22: m = 33: m = 34: m = 45: m = 36: m = 37: m = 3Average
W3G71254.0772983.0353677.8483720.6703944.9373901.6794087.1023367.050
W2A11079.6672694.5373528.1493119.9115896.6352938.7473526.8403254.927
W2A2969.3742471.7273287.0812956.7023218.2232891.4782981.1792682.252
W2A3990.3412843.7892859.0932999.7123002.1463074.8113155.3072703.600
W2A44446.528-------
W2B11084.7522634.8263295.1623051.9412835.7553238.3633272.6672773.352
W2B21075.3932429.9963130.2233051.1922929.1063581.4563155.6192764.712
W2B31180.3662397.1117193.0343000.3832970.7113739.7663139.0843374.351
W2B41274.6532932.0084872.972-----
W2C11132.0692685.3553242.6373287.0663147.6634080.0194802.6623196.782
W2C21112.1622881.6973189.7063873.0375211.6193665.7733950.8963412.127
W2C31052.5702421.0262768.3623014.0332778.5182914.9413953.3462700.399
W2C42051.7104690.2957193.034-----
W2F11071.5402407.1582909.9653472.3172832.2303490.8963246.5842775.813
W2F21015.0512438.4833031.8023061.2703703.1523737.3943324.5372901.670
W2F31272.1882596.2003603.6554130.158----
W2F41216.8392620.0523589.1773233.1683534.3123521.6603934.8453092.865
Table 7. Number of points requiring 40 iterations for each example (1–7).
Table 7. Number of points requiring 40 iterations for each example (1–7).
MapExample
1: m = 22: m = 33: m = 34: m = 45: m = 36: m = 37: m = 3Average
W3G7677605125040126533,0725130.000
W2A1657617116634,396120118,9397996.714
W2A269760111621120115,3852578.286
W2A3675605551529122114,7112489.714
W2A478,843-------
W2B167961312049120113,9462379.000
W2B263560911161121715,9952653.429
W2B36796131146310,64516,3424061.286
W2B46596172-----
W2C168961714002018,15724,2396303.286
W2C2669601117417,843126518,3825562.143
W2C365960911841120114,8632502.571
W2C46016011-----
W2F1681601112610122518,7723059.429
W2F26796171016143515159317,4693512.571
W2F3663605178----
W2F4645605113012121720,0203232.857

Share and Cite

MDPI and ACS Style

Lee, M.-Y.; Kim, Y.I.; Neta, B. A Generic Family of Optimal Sixteenth-Order Multiple-Root Finders and Their Dynamics Underlying Purely Imaginary Extraneous Fixed Points. Mathematics 2019, 7, 562. https://doi.org/10.3390/math7060562

AMA Style

Lee M-Y, Kim YI, Neta B. A Generic Family of Optimal Sixteenth-Order Multiple-Root Finders and Their Dynamics Underlying Purely Imaginary Extraneous Fixed Points. Mathematics. 2019; 7(6):562. https://doi.org/10.3390/math7060562

Chicago/Turabian Style

Lee, Min-Young, Young Ik Kim, and Beny Neta. 2019. "A Generic Family of Optimal Sixteenth-Order Multiple-Root Finders and Their Dynamics Underlying Purely Imaginary Extraneous Fixed Points" Mathematics 7, no. 6: 562. https://doi.org/10.3390/math7060562

APA Style

Lee, M. -Y., Kim, Y. I., & Neta, B. (2019). A Generic Family of Optimal Sixteenth-Order Multiple-Root Finders and Their Dynamics Underlying Purely Imaginary Extraneous Fixed Points. Mathematics, 7(6), 562. https://doi.org/10.3390/math7060562

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