Next Article in Journal
Fractional Maclaurin-Type Inequalities for Multiplicatively Convex Functions
Previous Article in Journal
Improved Results on Delay-Dependent and Order-Dependent Criteria of Fractional-Order Neural Networks with Time Delay Based on Sampled-Data Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Constructing a Family of Sixth-Order Methods for Multiple Roots

Department of Applied Mathematics, Dankook University, Cheonan 31116, Republic of Korea
Fractal Fract. 2023, 7(12), 878; https://doi.org/10.3390/fractalfract7120878
Submission received: 2 November 2023 / Revised: 5 December 2023 / Accepted: 9 December 2023 / Published: 11 December 2023

Abstract

:
A family of three-point, sixth-order, multiple-zero solvers is developed, and special cases of weight functions are investigated based on polynomials and low-order rational functions. The chosen cases of the proposed iterative method are compared with existing methods. The experiments show the superiority of the proposed schemes in terms of the number of divergent points and the average number of function evaluations per point. The dynamical characteristics of the developed methods, along with their illustrations, are represented with detailed analyses, comparisons, and comments.

1. Introduction

Nonlinear equations [1,2,3,4] occur frequently in various fields of science, artificial intelligence, and engineering. Finding the roots of a nonlinear equation [5,6,7,8] involves determining the value of a variable that satisfies a given equation. A zero α of h ( x ) = 0 is called a multiple root with multiplicity m if h ( i ) ( α ) = 0 , i = 0 , 1 , 2 , , m 1 and h ( m ) ( α ) 0 .
It is known that the Newton method is the most used method for solving the equations, given by
x n + 1 = x n m f ( x n ) f ( x n ) , n = 0 , 1 , 2 , .
Researchers [9,10,11,12,13,14,15] are interested in finding multiple roots of nonlinear equations [16,17,18,19,20] and investigating the dynamics by exploring the relevant basins of attraction [21,22,23,24,25,26].
Geum-Kim-Neta [27,28] constructed a class of two-point, sixth-order, multiple-root solvers. Since the following two methods are faster and require fewer iterations per point on average, we have chosen two methods as follows. These two schemes, called here (X1) and (X2), are given by
y n = x n m f ( x ) f ( x ) , x n + 1 = y n m [ 2 ( m 1 ) ( u s ) u 2 2 u s + 1 ] h ,
and
y n = x n m f ( x n ) f ( x n ) , x n + 1 = y n a u + m + c u 2 + 1 + b u · 1 1 + d s ,
where
a = 13 + 2 m ( 4 m 4 16 m 3 + 31 m 2 30 m + 13 ) ( m 1 ) ( 7 + 4 m 2 8 m ) , b = 4 ( 2 m 2 4 m + 3 ) ( m 1 ) ( 4 m 2 8 m + 7 ) , c = 3 + 4 m 2 8 m 7 + 4 m 2 8 m , d = 2 ( m 1 ) , u = ( f ( y n ) f ( x n ) ) 1 m , s = ( f ( y n ) f ( x n ) ) 1 m 1 .
In this paper, we propose a three-point, sixth-order method by adding a third step, as follows:
y n = x n m f ( x n ) f ( x n ) , w n = x n m A f ( k ) f ( x n ) f ( x n ) , k = ( f ( y n ) f ( x n ) ) 1 m 1 , x n = x n m B f ( k , v ) f ( x n ) f ( x n ) , v = ( f ( w n ) f ( x n ) ) 1 m ,
where A f : C C is analytic in a small neighborhood of 0, and B f : C 2 C is holomorphic in a neighborhood of ( 0 , 0 ) .
Definition 1.
Let a 0 , a 1 , a n be a sequence converging to a zero α and e n = a n α be the n-th iterative error. If there exist real numbers β R and δ R { 0 } , such that the following error equations satisfy
e n + 1 = δ e n β + O ( e n β + 1 ) ,
then | δ | or δ is called the asymptotic error constant, and β is defined as the order of convergence [1].
Definition 2.
Suppose that the theoretical asymptotic error constant ω = lim n | e n | | e n 1 | β and the convergence order β 1 are known [2,3]. The computational convergence order β n = log | e n / ω | log | e n 1 | is defined. Then lim n β n = β .
Definition 3.
If a fixed point ρ of a rational map J is attracting, then all nearby points of ρ are attracted to ρ under a rational map J. The collection of all points whose iterates under J converge to ρ is called the basin of attraction of ρ.
We investigate the convergence behavior of the proposed scheme in Section 2. The numerical results and the basins of attraction for the typical examples are shown in Section 3. In addition, the basins of attraction for the equation in the blood rheology model are shown, and the convergent regions for selected values are plotted in Section 3. Finally, Section 4 is discussed in the last section.

2. Convergence Analysis

In this section, we determine the maximal convergence order of the developed schemes. We investigate the main theorem describing the convergent behavior of the proposed schemes and address how to construct weight functions A f and B f for sixth-order convergence. It is sufficient to consider weight functions A f and B f up to the fifth-order terms in e n in terms of O ( f ( x n ) f ( x n ) ) = O ( e n ) .
Let a function f : C C have a root α of multiplicity m > 1 and be analytic in a small neighborhood of α .
y n = x n m f ( x n ) f ( x n ) , w n = x n m · A f ( k ) · f ( x n ) f ( x n ) , k = ( f ( y n ) f ( x n ) ) 1 m 1 , x n = x n m · B f ( k , v ) · f ( x n ) f ( x n ) , v = ( f ( w n ) f ( x n ) ) 1 m ,
where A f : C C is analytic in a small neighborhood of the origin 0 and B f : C 2 C is holomorphic in a small neighborhood of ( 0 , 0 ) .
Using the Taylor expansion about α , we have
f ( x n ) = e m Γ ( 1 + θ 2 e n + θ 3 e n 2 + θ 4 e n 3 + θ 5 e n 4 + θ 6 e n 5 + θ 7 e n 6 + O ( e n 7 ) ) ,
f ( x n ) = e m 1 Γ ( m + ( 1 + m ) θ 2 e n + ( 2 + m ) θ 3 e n 2 + ( 3 + m ) θ 4 e n 3 + ( 4 + m ) θ 5 e n 4 + ( 5 + m ) θ 6 e n 5 + ( 6 + m ) θ 7 e n 6 + O ( e n 7 ) ) ,
where Γ = f ( m ) ( α ) m ! , θ j = m ! ( m + j 1 ) ! f ( m + j 1 ) ( α ) f ( m ) ( α ) for j N { 1 } .
For convenience, we denote e n by e without subscript n. Dividing (5) by (6), we have
f ( x n ) f ( x n ) = e m θ 2 e 2 m 2 + t 3 e 3 m 3 + t 4 e 4 m 4 + t 5 e 5 m 5 + t 6 e 6 m 6 + O ( e 7 ) ,
where
t 3 = ( ( 1 + m ) θ 2 2 2 m θ 3 ) , t 4 = ( ( 1 + m ) 2 θ 2 3 + m ( 4 + 3 m ) θ 2 θ 3 3 m 2 θ 4 ) , t 5 = ( 1 + m ) 3 θ 2 4 2 m ( 3 + 5 m + 2 m 2 ) θ 2 2 θ 3 + 2 m 2 ( 3 + 2 m ) θ 2 θ 4 + 2 m 2 ( ( 2 + m ) θ 3 2 2 m θ 5 ) ,
and
t 6 = ( 1 + m ) 4 θ 2 5 + m ( 1 + m ) 2 ( 8 + 5 m ) θ 2 3 θ 3 m 2 ( 9 + 14 m + 5 m 2 ) θ 2 2 θ 4 + m 2 θ 2 ( ( ( 12 + 16 m + 5 m 2 ) θ 3 2 ) + m ( 8 + 5 m ) θ 5 ) + m 3 ( ( 12 + 5 m ) θ 3 θ 4 5 m θ 6 ) .
By Taylor’s expansion, we have expression k, as follows:
k = ( f ( y n ) f ( x n ) ) 1 m 1 , = θ 2 m e + m ( 1 + m ) θ 2 2 2 m ( 1 m ) θ 3 m ( 1 + m ) e 2 + w 3 2 ( m 1 ) 2 e 3 + w 4 6 ( 1 + m ) 3 e 4 + w 5 24 ( 1 + m ) 4 e 5 + w 6 120 ( 1 + m ) 5 e 6 + O ( e 7 ) ) ,
where w i = w i ( θ 2 , θ 3 , θ 6 ) for 3 i 6 .
Using k in (8), and expanding the Taylor series of A f ( k ) about 0 up to the fifth-order term, we find
A f ( k ) = A 0 + A 1 k + A 2 k 2 + A 3 k 3 + A 4 k 4 + A 5 k 5 + O ( e 6 ) ,
where A j = A f ( j ) ( 0 ) j ! for 0 j 5 .
Substituting (5)–(9) into w n , we have
w n = α + ( 1 A 0 ) e + ( A 0 A 1 ) m θ 2 e 2 + z 3 e 3 + z 4 e 4 + z 5 e 5 + z 6 e 6 + O ( e 7 ) ,
where
z 3 = ( A 0 A 1 + A 2 + 2 A 1 m A 2 m A 0 m 2 + A 1 m 2 ) θ 2 2 + ( 2 A 0 m + 2 A 1 m + 2 A 0 m 2 2 A 2 m 2 ) θ 3 m 2 ( 1 + m ) ,
z i = z i ( θ 2 , θ n , A 0 , A 1 , A 4 ) for 4 i 6 .
Choosing A 0 = A 1 = 1 , A 2 = 2 m 1 + m , we have
w n = α + ( ( 2 + m + m 3 2 m 2 ( 4 + A 3 ) 2 A 3 + 4 m A 3 ) θ 2 3 2 ( 1 + m ) m 2 θ 2 θ 3 ) 2 m 3 ( 1 + m ) 2 e 4 + z 5 e 5 + z 6 e 6 + O ( e 7 ) .
Then we obtain
f ( w n ) = e 4 m Γ ( 1 + m ) 3 m ( 1 2 ( 2 + m + m 3 2 m 2 ( 4 + A 3 ) 2 A 3 + 4 m A 3 ) θ 2 3 ( 1 + m ) m 2 θ 2 θ 3 m + 1 3 ( ( 2 + m + m 3 2 m 2 ( 4 + A 3 ) 2 A 3 + 4 m A 3 ) θ 2 3 2 ( 1 + m ) m 2 θ 2 θ 3 ) 1 + m ( ( 7 m 5 + 2 m 4 ( 35 + 9 A 3 ) + m ( 17 + 36 A 3 18 A 4 ) + 6 ( 2 A 3 + A 4 ) 6 m 3 ( 9 + 2 A 3 + A 4 ) + m 2 ( 14 36 A 3 + 18 A 4 ) ) θ 2 4 + 12 ( 1 + m ) m ( 2 + m + 2 m 3 3 m 2 ( 4 + A 3 ) 3 A 3 + 6 m A 3 ) θ 2 2 θ 3 12 ( 1 + m ) 2 m 3 θ 3 2 12 ( 1 + m ) 2 m 3 θ 2 θ 4 ) e + w 2 e 2 + w 3 e 3 + O ( e 4 ) ) ,
where w i = w i ( θ 2 , θ n , A 3 , A 4 ) , i = 2 , 3 .
With the use of (5) and (12), we have v as follows:
v = ( ( 2 + m + m 3 2 m 2 ( 4 + A 3 ) 2 A 3 + 4 m A 3 ) θ 2 3 2 ( 1 + m ) m 2 θ 2 θ 3 2 ( 1 + m ) 2 m 3 e 3 + 1 6 ( 1 + m ) 3 m 4 ( ( 7 m 5 + m 4 ( 73 + 18 A 3 ) + 2 m ( 7 + 27 A 3 9 A 4 ) + 6 ( 3 2 A 3 + A 4 ) 3 m 3 ( 25 + 2 A 3 + 2 A 4 ) + m 2 ( 35 54 A 3 + 18 A 4 ) ) θ 2 4 + 6 ( 1 + m ) m ( 4 + m + 4 m 3 + m 2 ( 25 6 A 3 ) 6 A 3 + 12 m A 3 ) θ 2 2 θ 3 12 ( 1 + m ) 2 m 3 θ 3 2 12 ( 1 + m ) 2 m 3 θ 2 θ 4 ) e 4 + v 5 e 5 + O ( e 6 ) ,
where v 5 = v 5 ( θ 2 , θ n , A 1 , A 4 ) .
Expanding the Taylor series of B f ( k , v ) about (0,0) up to the fifth-order term, we find
B f ( k , v ) = B 00 + B 30 k 3 + B 40 k 4 + B 50 k 5 + B 01 v + k ( B 10 + B 11 v ) + k 2 ( B 20 + B 21 v ) + O ( e 6 ) ,
where B i j = 1 i ! j ! i + j k i v j B f ( k , v ) | ( k = 0 , v = 0 ) for 0 i 5 , 0 j 1 .
By Substituting (5)–(14) into the proposed scheme (4), we have
x n + 1 α = x n α B f ( k , v ) f ( x n ) f ( x n )
= C 1 e + C 2 e 2 + C 3 e 3 + C 4 e 4 + C 5 e 5 + C 6 e 6 + O ( e 7 ) ,
where C 1 = 1 B 00 and the coefficients C i ( 2 i 6 ) depend of m, A j ( j = 0 , 1 , 5 ) , and θ i ( i = 1 , 2 , ) .
Solving C 1 for B 00 , we have
B 00 = 1 .
Using B 00 = 1 into C 2 , we have
B 10 = 1 .
Substituting B 00 = B 10 = 1 into C 3 , we obtain ( B 20 + 2 m B 20 m ) θ 2 2 m ( m 1 ) = 0 .
From which we find
B 20 = 2 m 1 + m .
Substituting B 00 = B 10 = 1 and B 20 = 2 m 1 + m into C 4 = 0 , we have
C 4 = 1 2 ( 1 + m ) 2 m 3 m 1 + m ( 2 ( 4 + ( 4 + A 3 ) B 01 B 30 ) m 2 + 3 1 + m 2 ( 1 + B 01 A 3 B 01 + B 30 ) m 3 1 + m + ( 1 B 01 ) m 3 m 1 + m + ( 1 B 01 4 A 4 B 01 + 4 B 30 ) m 2 + m 1 + m ) θ 2 3 + m 3 m 1 + m ( 2 m 2 + 3 1 + m 2 B 01 m 2 + 3 1 + m 2 ( 1 B 01 ) m 3 m 1 + m ) θ 2 θ 3 2 ( 1 + m ) 2 .
and
B 01 = 1 , B 30 = A 3 .
Substituting B 00 = B 10 = 1 , B 20 = 2 m 1 + m B 01 = 1 , B 30 = A 3 into C 5 = 0 , we have
C 5 = 1 2 ( 1 + m ) 3 m 4 m 1 + m ( ( 14 6 A 4 6 A 3 ( 2 + B 11 ) + 7 B 11 + 6 B 40 ) ) m 2 + 4 1 + m + ( 14 + 2 A 4 + 2 A 3 ( 2 + B 11 ) 7 B 11 2 B 40 ) m 3 + 4 1 + m + 2 ( 2 A 4 A 3 ( 2 + B 11 ) + B 11 + B 40 ) m 4 1 + m ) + ( 2 B 11 ) m 4 m 1 + m + ( 2 + 6 A 4 + 6 A 3 ( 2 + B 11 ) B 11 6 B 40 ) m 3 + m 1 + m θ 2 4 + m 4 m 1 + m ( 2 ( 2 + B 11 ) m 2 + 4 1 + m 4 ( 2 + B 11 ) m 3 + 4 1 + m 2 ( 2 B 11 ) m 4 m 1 + m ) θ 2 2 θ 3 2 ( 1 + m ) 3 ,
and
B 40 = A 4 , B 11 = 2 .
Putting (17)–(21) into C 6 , we obtain
C 6 = p 1 θ 2 5 + p 2 θ 2 3 θ 3 + p 3 θ 2 θ 3 2 4 ( 1 + m ) 3 m 5 ,
where
p 1 = 4 A 5 ( 1 + m ) 3 + 4 A 3 B 21 ( 1 + m ) 3 4 B 50 ( 1 + m ) 3 2 B 21 ( 1 + m ) ( 2 + m + 8 m 2 + m 3 ) + ( 2 2 A 3 ( 1 + m ) 2 + m + m 2 ( 8 + m ) ) ( 3 + m ( 10 + m ) ) , p 2 = 4 ( 1 + m ) m ( 1 + m 3 A 3 m 2 ( 9 + B 21 + A 3 ) + m ( 1 + B 21 + 2 A 3 ) ) , p 3 = 4 ( 1 + m ) 2 m 3 .
Then the proposed method, (4), is of the sixth order and possesses the following error equation:
e n + 1 = p 1 θ 2 5 + p 2 θ 2 3 θ 3 + p 3 θ 2 θ 3 2 4 ( 1 + m ) 3 m 5 e n 6 + O ( e n 7 ) ,
where
p 1 = 4 A 5 ( 1 + m ) 3 + 4 A 3 B 21 ( 1 + m ) 3 4 B 50 ( 1 + m ) 3 2 B 21 ( 1 + m ) ( 2 + m + 8 m 2 + m 3 ) + ( 2 2 A 3 ( 1 + m ) 2 + m + m 2 ( 8 + m ) ) ( 3 + m ( 10 + m ) ) , p 2 = 4 ( 1 + m ) m ( 1 + m 3 A 3 m 2 ( 9 + B 21 + A 3 ) + m ( 1 + B 21 + 2 A 3 ) ) , p 3 = 4 ( 1 + m ) 2 m 3 .
Theorem 1.
Let m N { 1 } . Let a function f : C C have a multiple root α with multiplicity m. Assume that f is analytic in a neighborhood of α. Let θ j = m ! ( m + j 1 ) ! f ( n + j 1 ) ( α ) f ( m ) ( α ) for j N { 1 } . Let x 0 be an initial value in a neighborhood of zero α. Let A j ( 0 j 5 ) be defined in (9) and let B i j ( 0 j 5 , 0 j 1 ) be defined in (14). Suppose A 0 = A 1 = 1 , A 2 = 2 m 1 + m , | A 3 |   < , | A 4 |   < , | A 5 |   < and B 00 = B 10 = B 01 = 1 , B 20 = 2 m 1 + m , B 30 = A 3 , B 40 = A 4 , B 11 = 2 , | B 50 | < , | B 21 | < hold. Then the iterative scheme, (4), is of the sixth-order and has the following error equation:
e n + 1 = p 1 θ 2 5 + p 2 θ 2 3 θ 3 + p 3 θ 2 θ 3 2 4 ( 1 + m ) 3 m 5 e n 6 + O ( e n 7 ) ,
where p i ( i = 1 , 2 , 3 ) is defined in (22).

3. Numerical Experiments

Based on the convergence analysis, Taylor-polynomial forms of A f ( k ) and B f ( k , v ) are given by
A f ( k ) = A 0 + A 1 k + A 2 k 2 + A 3 k 3 + A 4 k 4 + A 5 k 5 , B f ( k , v ) = B 00 + B 10 k + B 20 k 2 + B 30 k 3 + B 40 k 4 + B 50 k 5 + ( B 01 + B 11 k + B 21 k 2 ) v ,
where A 0 = A 1 = 1 , A 2 = 2 m / ( 1 + m ) , A 3 , A 4 , A 5 ( f r e e ) and B 00 = B 01 = B 10 = 1 , B 20 = 2 m / ( 1 + m ) , B 11 = 2 , B 30 = A 3 , B 40 = A 4 .
Even if various forms of A f ( k ) and B f ( k , v ) are possible, we limit ourselves to consider four forms, as follows:
( M 1 ) B 50 = B 21 = 0
A f ( k ) = 1 + k + 2 m 1 + m k 2 , B k ( k , v ) = 1 + k + 2 m 1 + m k 2 + ( 1 + 2 k ) v .
( M 2 ) B 50 = 0 , B 21 = 1
A f ( k ) = 1 + k + 2 m 1 + m k 2 , B k ( k , v ) = 1 + k + 2 m 1 + m k 2 + ( 1 + 2 k + k 2 ) v .
( M 3 ) B 50 = 1 , B 21 = 0
A f ( k ) = 1 + k + 2 m 1 + m k 2 , B k ( k , v ) = 1 + k + 2 m 1 + m + k 5 + ( 1 + 2 k ) v .
( M 4 ) The product of two univariate functions
A f ( k ) = 1 + 2 m 1 + m k 2 1 k , B k ( k , v ) = 1 + 2 m 1 + m k 2 1 k v .
The numerical experiments are carried out using Mathematica programming to confirm the developed theory. For the experiments, we maintain a minimum of 100 digits of precision to achieve the specified accuracy. If the root is not exact, it is approximated by a value with greater precision, possessing more significant digits than the specified number of precision.
Various numerical experiments have been conducted with Mathematica software 13.0 [29] to confirm the developed theory. In Table 1, the computational convergence order approaches 6. They confirm the sixth-order convergence with test functions h 1 ( x ) , h 2 ( x ) , h 3 ( x ) , and h 4 ( x ) , as follows:
h 1 ( x ) = ( cos [ π 2 x ] + x 2 π ) 5 , m = 5 , α 2.0347 . h 2 ( x ) = ( cos [ x 2 1 ] + 1 x log [ x 2 π ] ) ) 2 ( x 2 1 π ) , m = 3 , α = 1 + π . h 3 ( x ) = ( sin 1 [ x 1 ] + x x 2 3 ) 3 , m = 3 , α 1.0414818 . h 4 ( x ) = ( cos 2 x + 9 2 x 2 x 4 ) ( sin 2 x + 5 x x 4 ) m = 2 , α 1.29173329 .
Methods ( M 1 ) , ( M 2 ) , ( M 3 ) , ( M 4 ) , ( X 1 ) , and ( X 2 ) have been applied to the test functions f 1 ( x ) , f 2 ( x ) , f 3 ( x ) below:
f 1 ( x ) = exp [ ( 1 + x 3 ) 2 7 cos ( x 3 + 1 ) + x 5 ] 1 , α = 1 + i 3 2 , m = 2 , x 0 = 0.51 + 0.83 i , i = 1 f 2 ( x ) = ( π + 2 x + cos [ x ] log ( x 2 + 1 ) ) 3 , α = π / 2 , m = 3 , x 0 = 1.5 , f 3 ( x ) = ( x 3 x 3 cos [ π x 6 ] + 1 x 2 + 1 11 5 + 4 3 ) ( x 2 ) 3 , α = 2 , m = 4 , x 0 = 1.93 ,
In Table 2, we compare errors | x n α | of methods ( M 1 ) , ( M 2 ) , ( M 3 ) , ( M 4 ) , ( X 1 ) , and ( X 2 ) . The least errors within the prescribed error bound are marked in boldface. Method ( M 1 ) shows slightly better convergence in this experiment. The local convergence depends on the function, an initial guess, the root, the multiplicity, and the weight functions. So, for a selected set of test functions, one scheme is hardly expected to show better achievement than the others.
We study the dynamics of numerical methods ( M 1 ) , ( M 2 ) , ( M 3 ) , ( M 4 ) , ( X 1 ) , and ( X 2 ) . In Table 3, abbreviations CPU, tcon, avg, and tdiv denote the value of the CPU time for sixth-order convergence, the total number of convergent points, the value of the average iterative number for sixth-order convergence, and the number of divergent points. Table 3 shows the statistical data for the basins of attraction.
In Figure 1, Figure 2, Figure 3 and Figure 4, a six-by-six square region centered at the origin includes all the roots of the test functions. A 600 × 600 uniform grid in this region is used to mark initial values for the iterative schemes to draw the basins of attraction. The grid point of a square is colored differently according to the number of iterations required for convergence. The point is colored in black if the methods do not converge.
Example 1.
We choose a quadratic polynomial raised to the power of two:
p 1 ( z ) = ( 2 z 2 + 1 ) 2
with roots α = ± 0.707107 i with a multiplicity of 2. Basins of attraction for (M1), (M2), (M3), (M4), (X1), and (X2) are illustrated in Figure 1. Each basin is painted in various colors. Basins for α = + 0.707107 i are painted purple and basins for α = 0.707107 i are painted green in Figure 1 (a)–(c). At zero, its color is light (white), while becoming darker for more iterative numbers required for convergence within the set iteration. Method (M1) performs better in terms of CPU, and M4 is better in terms of avg for p 1 ( z ) . There are similar lobes along the boundaries of the basins in (a), (b), and (c). Picture (d) has colored (yellow, magenta, green, blue, and so on) carrot shapes around the boundary. Methods (X1) and (X2) do not have the value of the total convergent point, and pictures (e) and (f) show the black region.
Example 2.
As a second test example, we selected a quadratic polynomial raised to the power of three with all real roots:
p 2 ( z ) = ( z 2 1 ) 3
The roots are α = ± 1 with a multiplicity of 3. Basins of attraction for (M1)-(X2) are illustrated in Figure 2. Method (M4) performs better in terms of avg. Basins for α = 1 are painted red and basins for α = 1 are blue. There are similar lobes along the boundaries of the basins in Figure 2a–d has colored (blue, green, magenta) carrot shapes around the boundary. Methods (X1) and (X2) have considerable black points and circular patterns made up of blue diamond shapes in Figure 2e,f.
Example 3.
We choose a cubic polynomial raised to the power of two:
p 3 ( z ) = ( z 3 1 ) 2 ,
with roots α = 0.5 ± 0.866025 i , 1 . Method (M4) is better regarding CPU and avg. The basins of method (X1) are represented by black regions and the basins of method (X2) are represented by circular light blue patterns in Figure 3.
The local convergence depends on the function, an initial value, the multiplicity, and a root. It cannot be said that one iterative method is always better than other methods. It is essential to choose an initial guess that guarantees the convergence of the numerical method. Deciding how close the initial guess should be to the root is not easy, as it depends on computational precision, the selected function, and the error bound. To choose a stable initial value, we can use basins of attraction in Figure 1, Figure 2, Figure 3 and Figure 4.
We use the equation in the blood rheology model [30]
f ( x ) = ( x 8 441 + 8 x 5 63 2 , 857 , 144 , 357 x 4 50 , 000 , 000 , 000 + 16 x 2 9 906 , 122 , 449 x 250 , 000 , 000 + 3 10 ) 4
to carry out the experiment (M1) for the basins in Figure 4.
Substituting A 0 = A 1 = 1 , A 2 = 2 m / ( 1 + m ) , A 3 = λ , A 4 = 0 , A 5 = 0 and B 00 = B 01 = B 10 = 1 , B 20 = 2 m / ( 1 + m ) , B 11 = 2 , B 21 = 0 , B 30 = A 3 = λ , B 40 = A 4 = 0 into (24) to plot the convergent region, we select the following method, called (M5):
A f ( k ) = 1 + k + ( 2 m / ( 1 + m ) ) k 2 + λ k 3 , B f ( k , v ) = 1 + k + ( 2 m / ( 1 + m ) ) k 2 + λ k 3 + ( 1 + 2 k ) v ,
Figure 5 shows the convergence region [31] of (M5) according to λ = 2.375 and λ = 3.3305 + 0.0712 i .

4. Conclusions

Even though the proposed methods (M1)–(M4) are of the sixth order, like the existing methods (X1) and (X2), the number of divergent points for (M1)–(M4) is much smaller than for (X1) and (X2). The basins for (X1) and (X2) are black for p 1 ( z ) and the basins for (X1) and (X2) for p 2 ( z ) have some blue patterns but also a lot of black points. In summary, we find that the proposed methods are faster and require fewer iterations per point on average. It is important to choose the type of initial guess for the iterative schemes, ensuring it is chosen near the root for convergence. In future work, the fractal behavior behind the improved methods will be studied in more detail.

Funding

The author is supported by the Basic Science Research Program through the National Research Foundation of Korea, funded by the Ministry of Education under research grant project number NRF-2021R1A2C1012922.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

We warmly thank two anonymous referees for several comments and suggestions, which helped improve this work.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Ahlfors, L.V. Complex Analysis; McGraw-Hill Book, Inc.: New York, NY, USA, 1979. [Google Scholar]
  2. Hörmander, L. An Introduction to Complex Analysis in Several Variables; North-Holland Publishing Company: Amsterdam, The Netherlands, 1973. [Google Scholar]
  3. Shabat, B.V. Introduction to Complex Analysis PART II, Functions of Several Variables; American Mathematical Society: Providence, RI, USA, 1992. [Google Scholar]
  4. Devaney, R.L. An Introduction to Chaotic Dynamical Systems; Addison-Wesley Publishing Company, Inc.: Boston, MA, USA, 1987. [Google Scholar]
  5. 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]
  6. Soleymani, F.; Babajee, D.K.R. Computing multiple zeros using a class of quartically convergent methods. Alex. Eng. J. 2013, 52, 531–541. [Google Scholar] [CrossRef]
  7. Chun, C.; Neta, B. Comparative study of eighth order methods for finding simple roots of nonlinear equations. Numer. Algorithms 2017, 74, 1169–1201. [Google Scholar] [CrossRef]
  8. Cordero, A.; Torresgrosa, 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]
  9. 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]
  10. Liu, L.; Wang, X. Eighth-order methods with high efficiency index for solving nonlinear equations. Appl. Math. Comput. 2010, 15, 3449–3454. [Google Scholar] [CrossRef]
  11. 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]
  12. 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]
  13. Traub, J.F. Iterative Methods for the Solution of Equations; Chelsea Publishing Company: Chelsea, VT, USA, 1982. [Google Scholar]
  14. 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]
  15. 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]
  16. Li, S.G.; Cheng, L.Z.; Neta, B. Some fourth-order nonliear solvers with closed formulae for multiple roots. Comput. Math. Appl. 2010, 59, 126–135. [Google Scholar] [CrossRef]
  17. Liu, B.; Zhou, X. A new family of fourth-order methods for mulltiple roots of nonlinear equations. Nonlinear Anal. Model. Contral 2013, 18, 143–152. [Google Scholar] [CrossRef]
  18. Magreñán, A.A.; Cordero, A.; Gitiérrez, J.M.; Torregrosa, J.R. Real qualitative behavior of a fourth-order family of iterative methods by using the convergence plane. Math. Comput. Simulat. 2014, 105, 49–61. [Google Scholar] [CrossRef]
  19. Magreñán, A.A. Different anomalies in a Jarratt family of iterative root-finding methods. Appl. Math. Comput. 2014, 233, 29–38. [Google Scholar] [CrossRef]
  20. Geum, Y.H. Study on the dynamical analysis of a family of third order multiple-zero finders. Eur. J. Pure Appl. Math. 2023, 16, 2775–2785. [Google Scholar] [CrossRef]
  21. Neta, B.; Scott, N.; Chun, C. Basins of attractors for various methods for multiple root. Appl. Math. Comput. 2012, 218, 5043–5066. [Google Scholar] [CrossRef]
  22. Neta, B.; Chun, C. On a family of Laguerre methods to find multiple roots of nonlinear equations. Appl. Math. Comput. 2013, 219, 10987–11004. [Google Scholar] [CrossRef]
  23. Neta, B.; Johnson, A.N. High order nonliear solver for multiple roots. Comput. Math. Appl. 2008, 55, 2012–2017. [Google Scholar] [CrossRef]
  24. Neta, B. Extension of Murakami’s high order nonlinear solver to multiple roots. Int. J. Comput. Math. 2010, 8, 1023–1031. [Google Scholar] [CrossRef]
  25. 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]
  26. 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]
  27. 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]
  28. 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]
  29. Wolfram, S. The Mathematica Book, 5th ed.; Wolfram Media: Champaign, IL, USA, 2003. [Google Scholar]
  30. Fournier, R.L. Basic Transport Phenomena in Biomedical Engineering; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar]
  31. Deng, J.J.; Chiang, H.D. Convergence Region of Newton Iterative Power Flow Method: Numerical Studies. J. Appl. Math. 2013, 2013, 509496. [Google Scholar] [CrossRef]
Figure 1. Basins of attraction for P 1 ( z ) = ( 2 z 2 + 1 ) 2 . (a) M 1 ; (b) M 2 ; (c) M 3 ; (d) M 4 ; (e) X 1 ; (f) X 2 .
Figure 1. Basins of attraction for P 1 ( z ) = ( 2 z 2 + 1 ) 2 . (a) M 1 ; (b) M 2 ; (c) M 3 ; (d) M 4 ; (e) X 1 ; (f) X 2 .
Fractalfract 07 00878 g001
Figure 2. Basins of attraction for P 2 ( z ) = ( z 2 1 ) 3 . (a) M 1 ; (b) M 2 ; (c) M 3 ; (d) M 4 ; (e) X 1 ; (f) X 2 .
Figure 2. Basins of attraction for P 2 ( z ) = ( z 2 1 ) 3 . (a) M 1 ; (b) M 2 ; (c) M 3 ; (d) M 4 ; (e) X 1 ; (f) X 2 .
Fractalfract 07 00878 g002
Figure 3. Basins of attraction for P 3 ( z ) = ( z 3 1 ) 2 . (a) M 1 ; (b) M 2 ; (c) M 3 ; (d) M 4 ; (e) X 1 ; (f) X 2 .
Figure 3. Basins of attraction for P 3 ( z ) = ( z 3 1 ) 2 . (a) M 1 ; (b) M 2 ; (c) M 3 ; (d) M 4 ; (e) X 1 ; (f) X 2 .
Fractalfract 07 00878 g003
Figure 4. Basins of attraction for the blood rheology model.
Figure 4. Basins of attraction for the blood rheology model.
Fractalfract 07 00878 g004
Figure 5. Convergent region for (M5). (a) λ = 2.375 ; (b) λ = 3.3305 + 0.0712 i .
Figure 5. Convergent region for (M5). (a) λ = 2.375 ; (b) λ = 3.3305 + 0.0712 i .
Fractalfract 07 00878 g005
Table 1. Convergent process for g i ( x ) , ( i = 1 , 2 , 3 , 4 ) .
Table 1. Convergent process for g i ( x ) , ( i = 1 , 2 , 3 , 4 ) .
MethodFunctionn x n | x n α | | e n / e n 1 6 | p n
M 1 h 1 0 2.1 0.0652751
1 2.0372492017726 4.913 × 10 8 0.30894310956.07933
2 2.0372476627913 5.378 × 10 45 0.4282207006.00000
3 2.0372476627913 0.0 × 10 98
M 2 h 2 0 2.0 0.0350903
1 2.0350902813353 2.389 × 10 7 26.317219536.10370
2 2.0350903306632 7.978 × 10 45 38.0171676.00000
3 2.0350903306632 0.0 × 10 99
M 3 h 3 0 1.084 0.0425181
1 1.04148199694193 1.263 × 10 7 21.387333546.06612
2 1.041481808433 1.076 × 10 40 26.442054496.00000
3 1.041481808433 0.0 × 10 99
M 4 h 4 0 1.35 0.0582667
1 1.29173359504767 3.07 × 10 7 7.73306876.23175
2 1.29173329265770 9.550 × 10 40 12.734657936.00000
3 1.29173329265770 0.0 × 10 99
Table 2. Comparison of | x n α | for f 1 ( x ) , f 2 ( x ) and f 3 ( x ) .
Table 2. Comparison of | x n α | for f 1 ( x ) , f 2 ( x ) and f 3 ( x ) .
Function | x n α | M 1 M 2 M 3 M 4 X 1 X 2
f 1 | x 1 α | 2.31 × 10 11 2.61 × 10 7 7.76 × 10 10 1.21 × 10 8 3.11 × 10 7 2.63 × 10 6
| x 2 α | 3 . 71 × 10 62 2.35 × 10 43 2.42 × 10 55 3.51 × 10 47 8.52 × 10 40 7.16 × 10 41
f 2 | x 1 α | 5.16 × 10 7 5.57 × 10 8 3.12 × 10 7 2.12 × 10 7 5.32 × 10 6 2.22 × 10 5
| x 2 α | 2 . 31 × 10 41 1.66 × 10 40 4.16 × 10 39 3.51 × 10 38 2.12 × 10 37 5.31 × 10 39
f 3 | x 1 α | 2.71 × 10 6 3.46 × 10 6 1.67 × 10 5 1.88 × 10 6 7.62 × 10 5 4.42 × 10 8
| x 2 α | 3 . 71 × 10 42 5.23 × 10 42 5.11 × 10 41 6.75 × 10 40 6.51 × 10 35 4.12 × 10 41
Table 3. Comparison of CPU time, tcon, avg, and tdiv.
Table 3. Comparison of CPU time, tcon, avg, and tdiv.
p m MethodCPUtconavgtdiv
p 1 ( z ) M1120.922360,0008.739610
M2165.969360,0008.967830
M3177.515360,0008.967830
M4183.984354,4408.151075560
X11234.780-360,000
X2333.5780-360,000
p 2 ( z ) M1274.985360,0009.091360
M2340.906360,0009.197660
M3266.187360,0009.079980
M4306.485360,0007.921030
X110.82220,880142.582339,120
X21575.34611837.9915353,882
p 3 ( z ) M13104.25360,0009.30217107,540
M23584.16249,0228.94829110,978
M33660.17249,0228.94829110,978
M4497.906287,6545.4836172,346
X1519.5630-360,000
X23443.61105426.759358,946
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Geum, Y.H. On Constructing a Family of Sixth-Order Methods for Multiple Roots. Fractal Fract. 2023, 7, 878. https://doi.org/10.3390/fractalfract7120878

AMA Style

Geum YH. On Constructing a Family of Sixth-Order Methods for Multiple Roots. Fractal and Fractional. 2023; 7(12):878. https://doi.org/10.3390/fractalfract7120878

Chicago/Turabian Style

Geum, Young Hee. 2023. "On Constructing a Family of Sixth-Order Methods for Multiple Roots" Fractal and Fractional 7, no. 12: 878. https://doi.org/10.3390/fractalfract7120878

APA Style

Geum, Y. H. (2023). On Constructing a Family of Sixth-Order Methods for Multiple Roots. Fractal and Fractional, 7(12), 878. https://doi.org/10.3390/fractalfract7120878

Article Metrics

Back to TopTop