Next Article in Journal
Electrothermal Monte Carlo Simulation of a GaAs Resonant Tunneling Diode
Previous Article in Journal
Parameter Estimation Analysis in a Model of Honey Production
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Analysis of Simple Root Seeker for Nonlinear Equation

School of Mathematical Sciences, Bohai University, Jinzhou 121013, China
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(2), 215; https://doi.org/10.3390/axioms12020215
Submission received: 28 December 2022 / Revised: 15 February 2023 / Accepted: 16 February 2023 / Published: 18 February 2023

Abstract

:
In this paper, the stability of a class of Liu–Wang’s optimal eighth-order single-parameter iterative methods for solving simple roots of nonlinear equations was studied by applying them to arbitrary quadratic polynomials. Under the Riemann sphere and scaling theorem, the complex dynamic behavior of the iterative method was analyzed by fractals. We discuss the stability of all fixed points and the parameter spaces starting from the critical points with the Mathematica software. The dynamical planes of the elements with good and bad dynamical behavior are given, and the optimal parameter element with stable behavior was obtained. Finally, a numerical experiment and practical application were carried out to prove the conclusion.

1. Introduction

Nonlinear problems comprise an extraordinarily active field in modern mathematics. A huge number of nonlinear problems, such as nonlinear finite element problems, economic and nonlinear programming problems, and plentiful central problems in physics, chemistry, and fluid mechanics, are attributed to solving some specific nonlinear equations in many cases. In nonlinear mechanics, there are myriad types of nonlinearity, such as material nonlinearity, geometric nonlinearity, contact nonlinearity, and so on. Solving nonlinear equations is also a universal and crucial problem in scientific and engineering calculation [1,2].
However, except for some rather special cases, the direct method has difficulty solving nonlinear equations. For practical problems, in many cases, it is not required to gain the precise answer to the equation; only an approximate value is needed. Of course, the error between the approximate value and the exact solution should be controlled within the allowable range of the actual problem. The approximate solution can be acquired by numerical methods; hence, the study of the numerical solution of nonlinear equations has pivotal academic significance and realistic application worth. There are several numerical algorithms for nonlinear equations, and more numerical solutions are used, such as the multi-point iterative method [3,4], semi-explicit method [5,6], splitting method [7,8], and analytical numerical method [9].
Among these techniques, the essence of the iterative method is a method of successive approximation. In an effort to resolve the approximate solutions of nonlinear equations as accurately as possible, a large number of researchers have come up with distinctive iterative methods [10,11,12,13,14]. Generally speaking, in these iterative methods, the most-typical is Newton’s method. Newton’s iterative scheme is the earliest and most-classical numerical method for solving nonlinear equations [15]. Its iterative format is
y n = x n f ( x n ) f ( x n ) .
On the basis of Newton’s iterative method, numerous well-known lower-order iteration methods have been proposed [16,17]. In particular, the influential iterative scheme devised by Ostrowski, which is fourth-order, is
y n = x n f ( x n ) f ( x n ) , x n + 1 = y n f ( y n ) f ( x n ) f ( x n ) f ( x n ) 2 f ( y n ) .
Although there are many iterative methods, in practical applications, we will also selectively find iterative methods with higher information efficiency, a more stable convergence process, and higher convergence accuracy. Informational efficiency is usually defined according to the convergence order and computational cost of the iterative method. Based on the concept of the optimal order proposed by Kung [18], we know that the iterative method with the optimal order is the optimal iterative method, and its computational efficiency is the optimal efficiency. For the iterative method without memory, the information efficiency of the multi-point iterative method with the optimal order is higher than that of the non-optimal multi-point iterative method. For the study of iterative stability, its dynamical behavior can be explored. The convergence and stability of the multi-point iterative method are analyzed by drawing fractal diagrams, which are based on fractal theory [19,20].
Based on (1), Liu and Wang proposed an optimal eighth-order scheme, L W M , defined as [21]
y n = x n f ( x n ) f ( x n ) , z n = y n f ( x n ) f ( x n ) 2 f ( y n ) f ( y n ) f ( x n ) , x n + 1 = z n f ( z n ) f ( x n ) [ ( f ( x n ) f ( y n ) f ( x n ) 2 f ( y n ) ) 2 + f ( z n ) f ( y n ) + G ( μ n ) ] ,
where μ n = f ( z n ) f ( x n ) and G ( t ) represents a real-valued function and satisfies G ( 0 ) = 0 , G ( 0 ) = 4 .
In this paper, we chose G ( t ) = λ t 2 + 4 t , λ C , then the corresponding iterative format is L W M ( λ )
y n = x n f ( x n ) f ( x n ) , z n = y n f ( x n ) f ( x n ) 2 f ( y n ) f ( y n ) f ( x n ) , x n + 1 = z n f ( z n ) f ( x n ) [ ( f ( x n ) f ( y n ) f ( x n ) 2 f ( y n ) ) 2 + f ( z n ) f ( y n ) + λ ( f ( z n ) f ( x n ) ) 2 + 4 ( f ( z n ) f ( x n ) ) ] ,
where λ C .
This paper discusses the stability of the iterative method (3) under fractal study. Considering the complexity of the calculation process, we used the mathematical symbol calculation Mathematica software to achieve the calculation. Firstly, the rational operator corresponding to (3) was obtained under the Möbius conjugate map on the Riemann sphere. Secondly, on the basis of the rational operator, we studied its fixed points and critical points. Finally, the relevant parameter spaces and dynamical planes starting from the critical points were analyzed, and the most-stable member in (3) was obtained, which will be introduced in detail in Section 2. In Section 3, we perform numerical experiments to verify the results in the previous section. In the last section, some conclusions are given.

2. Stability of Iterative Method under Fractal Study

Using fractals to study the stability of iterative method is to study the complex dynamic behavior of the iterative method. The complex dynamic behavior refers to the related properties of the rational operator related to the iterative method under the Möbius conjugate map on the Riemann sphere. From a numerical point of view, the study of the dynamic properties of rational operators allows us to draw important conclusions about the stability and reliability of fixed points and critical points. From the perspective of parameter selection, by researching the parameter spaces of the method constructed from the critical points, the performance of different members can be understood, which is conducive to the selection of parameters. The dynamic planes show the stability of these special methods, and we obtain the elements of the most-stable parameter members.
Next, the complex dynamics of the L W M ( λ ) family described in (3) were explored. To start with, rational operators linked to L W M ( λ ) were constructed on nonlinear polynomials. On this basis, the stability of the corresponding fixed point and critical point was studied. Furthermore, the parameter spaces starting from the free critical points were constructed. Finally, the relevant dynamical planes were analyzed, and the influence of the parameter selection on stability is discussed.

2.1. Rational Operator

According to the consequences of Riemann sphere dynamics [22,23] and the scaling theorem [24], we can construct rational operators on any nonlinear function. Given that the quadratic polynomial stability criterion can be extended to other nonlinear functions, we constructed a rational operator on a quadratic polynomial.
Proposition 1.
Let g ( x ) = ( x a ) ( x b ) be any quadratic polynomial, where a , b R are roots. Consequently, the rational operator R λ ( x ) corresponding to the L W M ( λ ) family applied to g ( x ) given in (3) is
R λ ( x ) = x 8 · γ λ ( x ) θ λ ( x ) ,
where γ λ ( x ) = 13 + 46 x ( 100 + λ ) x 2 + 170 x 3 + 223 x 4 + 240 x 5 + 216 x 6 + 160 x 7 + 99 x 8 + 50 x 9 + 20 x 10 + 6 x 11 + x 12 and θ λ ( x ) = 1 + 6 x + 20 x 2 + 50 x 3 + 99 x 4 + 160 x 5 + 216 x 6 + 240 x 7 + 223 x 8 + 170 x 9 ( 100 + λ ) x 10 + 46 x 11 + 13 x 12 , λ C .
Proof. 
Let g ( x ) = ( x a ) ( x b ) be any quadratic polynomial, where a , b R are roots. By applying the iterative scheme given in (3) on g ( x ) , we gain a rational function F g , λ ( x ) , which is only related to a , b , λ . Subsequently, we utilized the Möbius transformation in F g , λ ( x ) with
r ( ε ) = ε a ε b ,
which satisfies r ( a ) = 0 , r ( b ) = , and r ( ) = 1 , and we obtain
R λ ( x ) = ( r F g , λ r 1 ) ( x ) = x 8 · γ λ ( x ) θ λ ( x ) ,
where γ λ ( x ) = 13 + 46 x ( 100 + λ ) x 2 + 170 x 3 + 223 x 4 + 240 x 5 + 216 x 6 + 160 x 7 + 99 x 8 + 50 x 9 + 20 x 10 + 6 x 11 + x 12 and θ λ ( x ) = 1 + 6 x + 20 x 2 + 50 x 3 + 99 x 4 + 160 x 5 + 216 x 6 + 240 x 7 + 223 x 8 + 170 x 9 ( 100 + λ ) x 10 + 46 x 11 + 13 x 12 , which only depends on λ C . □
It can be seen from Proposition 1 that the analysis of the rational operator (4) is equivalent to the analysis of the iterative scheme (3). After that, we study the fixed point and its critical point of (4) by using Mathematica, a mathematical symbolic computing software, in Section 2.2.

2.2. Fixed Points and Critical Points

Here, we analyzed the fixed points and stability of R λ ( x ) . Firstly, it can be obtained from (4):
R λ ( x ) x = x ( x 1 ) · ρ λ ( x ) θ λ ( x ) ,
where ρ λ ( x ) = 1 + 7 x + 27 x 2 + 77 x 3 + 176 x 4 + 336 x 5 + 552 x 6 + 779 x 7 + 956 x 8 + ( 1026 + λ ) x 9 + 956 x 10 + 779 x 11 + 552 x 12 + 336 x 13 + 176 x 14 + 77 x 15 + 27 x 16 + 7 x 17 + x 18 .
Proposition 2.
For ρ λ ( x ) and θ λ ( x ) , there are the following conclusions:
  • When λ = 0 , ρ λ ( x ) and θ λ ( x ) have common factors ( 1 + x ) 2 ( 1 + x 2 ) 2 .
  • When λ = 6848 , ρ λ ( x ) has a factor ( x 1 ) 2 .
  • When λ = 1344 , θ λ ( x ) has a factor ( x 1 ) .
Proof. 
By solving the equations ρ λ ( x ) = 0 and θ λ ( x ) = 0 simultaneously, we have the following:
When λ = 0 , ρ λ ( x ) and θ λ ( x ) have common factors ( 1 + x ) 2 ( 1 + x 2 ) 2 . At this point, ρ λ ( x ) = ( 1 + x ) 2 ( 1 + x 2 ) 2 ( 1 + 5 x + 14 x 2 + 30 x 3 + 51 x 4 + 71 x 5 + 84 x 6 + 71 x 7 + 51 x 8 + 30 x 9 + 14 x 10 + 5 x 11 + x 12 ) and θ λ ( x ) = ( 1 + x ) 2 ( 1 + x 2 ) 2 ( 1 + 4 x + 9 x 2 + 16 x 3 + 21 x 4 + 20 x 5 + 13 x 6 ) .
Bring x = 1 into ρ λ ( x ) and θ λ ( x ) , respectively. We obtain: ρ λ ( 1 ) = 6848 + λ and θ λ ( 1 ) = 1344 λ . As a consequence, ρ λ ( x ) has a factor ( x 1 ) 2 when λ = 6848 and θ λ ( x ) has a factor ( x 1 ) when λ = 1344 . □
Proposition 3.
The fixed points of R λ ( x ) are x = 0 , x = and the following strange fixed points:
  • e x 1 = 1 (when λ 1344 ) and e x i ( λ ) , which correspond to the 18 roots of polynomial 1 + 7 x + 27 x 2 + 77 x 3 + 176 x 4 + 336 x 5 + 552 x 6 + 779 x 7 + 956 x 8 + ( 1026 + λ ) x 9 + 956 x 10 + 779 x 11 + 552 x 12 + 336 x 13 + 176 x 14 + 77 x 15 + 27 x 16 + 7 x 17 + x 18 , where i = 2 , 3 , , 19 .
Choose different λ-values; the number of fixed points is also different:
  • R λ ( x ) has 21 fixed points when λ C and λ { 0 , 1344 } .
  • R λ ( x ) has 20 fixed points excluding e x 1 = 1 when λ = 1344 .
  • R λ ( x ) has 15 fixed points when λ = 0 .
  • R λ ( x ) has 21 fixed points, and e x 1 = 1 is a triple root when λ = 6848 .
  • The strange fixed points of R λ ( x ) satisfy e x m = 1 e x n for m n ; each pair is conjugate to each other, being e x 2 and e x 3 , e x 4 and e x 5 , e x 6 and e x 7 , e x 8 and e x 9 , e x 10 and e x 11 , e x 12 and e x 13 , e x 14 and e x 15 , e x 16 and e x 17 , and e x 18 and e x 19 .
Proof. 
From (6), we can obtain
R λ ( x ) x = x ( x 1 ) · ρ λ ( x ) θ λ ( x ) = 0 .
When λ = 1344 , θ λ ( x ) = ( x 1 ) ( 1 + 7 x + 27 x 2 + 77 x 3 + 176 x 4 + 336 x 5 + 552 x 6 + 792 x 7 + 1015 x 8 + 1185 x 9 + 59 x 10 + 13 x 11 ) , and
R 0 ( x ) x = x ( x 1 ) ( 1 + 5 x + 14 x 2 + 30 x 3 + 51 x 4 + 71 x 5 + 84 x 6 + 71 x 7 + 51 x 8 + 30 x 9 + 14 x 10 + 5 x 11 + x 12 ) 1 + 4 x + 9 x 2 + 16 x 3 + 21 x 4 + 20 x 5 + 13 x 6 ,
R 1344 ( x ) x = x · ρ 1 ( x ) θ 1 ( x ) ,
where ρ 1 ( x ) = 1 + 7 x + 27 x 2 + 77 x 3 + 176 x 4 + 336 x 5 + 552 x 6 + 779 x 7 + 956 x 8 + 2370 x 9 + 956 x 10 + 779 x 11 + 552 x 12 + 336 x 13 + 176 x 14 + 77 x 15 + 27 x 16 + 7 x 17 + x 18 and θ 1 ( x ) = 1 + 7 x + 27 x 2 + 77 x 3 + 176 x 4 + 336 x 5 + 552 x 6 + 792 x 7 + 1015 x 8 + 1185 x 9 59 x 10 13 x 11 .
R 6848 ( x ) x = x ( x 1 ) 3 · ρ 2 ( x ) θ 2 ( x ) ,
where ρ 2 ( x ) = 1 + 9 x + 44 x 2 + 156 x 3 + 444 x 4 + 1068 x 5 + 2244 x 6 + 4199 x 7 + 7110 x 8 + 4199 x 9 + 2244 x 10 + 1068 x 11 + 444 x 12 + 156 x 13 + 44 x 14 + 9 x 15 + x 16 and θ 2 ( x ) = 1 + 6 x + 20 x 2 + 50 x 3 + 99 x 4 + 160 x 5 + 216 x 6 + 240 x 7 + 223 x 8 + 170 x 9 + 6948 x 10 + 46 x 11 + 13 x 12 . □
According to the above proposition, we determined that there are a maximum of 21 and a minimum of 15 fixed points, where 0, correspond to a , b and e x 1 = 1 (when λ 1344 ) corresponds to the divergence of the original method.
Proposition 4.
For the stability of e x 1 = 1 , λ C 1344 , verify:
(1)
e x 1 is a repulsive point when | 8192 1344 + λ | > 1 , that is | λ 1344 | > 8192 ;
(2)
e x 1 is an attractive point when | 8192 1344 + λ | < 1 , that is | λ 1344 | < 8192 ;
(3)
e x 1 is a parabolic point when | 8192 1344 + λ | = 1 , that is | λ 1344 | = 8192 ;
(4)
e x 1 is never a superattracting point because | 8192 1344 + λ | 0 .
Proof. 
By (4), we easily obtain
R λ ( x ) = 2 x 7 ( 1 + x ) 6 ( 1 + x 2 ) 5 · ω λ ( x ) θ λ ( x ) 2 ,
where ω λ ( x ) = 52 + 168 x + ( 336 5 λ ) x 2 + 3 ( 168 + λ ) x 3 + 2 ( 284 + λ ) x 4 + 3 ( 168 + λ ) x 5 + ( 336 5 λ ) x 6 + 168 x 7 + 52 x 8 .
Substituting e x 1 = 1 into (11), we acquire
| R λ ( 1 ) | = | 8192 1344 + λ | .
It is easy to obtain
| 8192 1344 + λ | = 1 | 8192 | = | 1344 + λ | . Let λ = z + i y , λ C , then the following formula holds:
( z 1344 ) 2 + y 2 = 8192 2 .
Therefore,
| R λ ( 1 ) | > 1 | λ 1344 | > 8192 .
Obviously, no matter what the λ -value is, 0 and are superattracting fixed points, while the stability of other fixed points varies with λ C . Table 1 summarizes the stability results of strange fixed points e x i corresponding to the special λ -values related to Proposition 2.
Figure 1 shows the stability surface of e x 1 = 1 . In Figure 1b, the gray surface represents the repulsion area, while the gold surface represents the attraction area. The attraction area is significantly smaller than the exclusion area. For the λ -values inside the disk, e x 1 is a repulsion, while for the λ -values outside the disk, e x 1 is an attractor. We are always interested in values within the disk.
According to Proposition 2, the stability of each pair of conjugate strange fixed points’ representations is the same, so the study of their stability can be reduced from eighteen strange fixed points to nine pairs of strange fixed points. Figure 2 shows the stability surfaces of these nine pairs of strange fixed points. The value of λ determines whether these fixed points attract or not. It can be seen from Figure 2 that there are few points that converge to the strange fixed points.
According to (11), the critical points of R λ ( x ) can be obtained.
Proposition 5.
According to the definition of the critical points, the critical points of R λ ( x ) are the roots of R λ ( x ) = 0 , that is to say the process of finding the critical points of R λ ( x ) is to solve R λ ( x ) = 0 . That is, x = 0 , x = , and the following free critical points:
  • c r 1 = 1 ;
  • c r 2 = i ;
  • c r 3 = i ;
  • c r 4 = 1 2 ( 21 26 + M ( λ ) + 1 2 P ( λ ) N ( λ ) 4 + ( 21 26 M ( λ ) 1 2 P ( λ ) N ( λ ) ) 2 ) ;
  • c r 5 = 1 2 ( 21 26 + M ( λ ) + 1 2 P ( λ ) N ( λ ) + 4 + ( 21 26 M ( λ ) 1 2 P ( λ ) N ( λ ) ) 2 ) ;
  • c r 6 = 1 2 ( 21 26 + M ( λ ) 1 2 P ( λ ) N ( λ ) 4 + ( 21 26 M ( λ ) + 1 2 P ( λ ) N ( λ ) ) 2 ) ;
  • c r 7 = 1 2 ( 21 26 + M ( λ ) 1 2 P ( λ ) N ( λ ) + 4 + ( 21 26 M ( λ ) + 1 2 P ( λ ) N ( λ ) ) 2 ) ;
  • c r 8 = 1 2 ( 21 26 M ( λ ) + 1 2 P ( λ ) + N ( λ ) 4 + ( 21 26 + M ( λ ) 1 2 P ( λ ) + N ( λ ) ) 2 ) ;
  • c r 9 = 1 2 ( 21 26 M ( λ ) + 1 2 P ( λ ) + N ( λ ) + 4 + ( 21 26 + M ( λ ) 1 2 P ( λ ) + N ( λ ) ) 2 ) ;
  • c r 10 = 1 2 ( 21 26 M ( λ ) 1 2 P ( λ ) + N ( λ ) 4 + ( 21 26 + M ( λ ) + 1 2 P ( λ ) + N ( λ ) ) 2 ) ;
  • c r 11 = 1 2 ( 21 26 M ( λ ) 1 2 P ( λ ) + N ( λ ) + 4 + ( 21 26 + M ( λ ) + 1 2 P ( λ ) + N ( λ ) ) 2 ) ;
where:
  • M ( λ ) = 1 2 a ( λ ) + b ( λ ) + c ( λ ) ,
  • P ( λ ) = 882 169 + 1 39 ( 128 + 5 λ ) b ( λ ) c ( λ ) N ( λ ) ;
  • N ( λ ) = g ( λ ) h ( λ ) ;
  • a ( λ ) = 441 169 + 1 156 ( 128 5 λ ) + 1 52 ( 128 + 5 λ ) ;
  • b ( λ ) = d ( λ ) 156 ( e ( λ ) + 18 f ( λ ) ) 1 / 3 ;
  • c ( λ ) = 1 156 ( e ( λ ) + 18 f ( λ ) ) 1 / 3 ;
  • d ( λ ) = 16384 + 4696 λ + 25 λ 2 ;
  • e ( λ ) = 2097152 + 1160832 λ + 139578 λ 2 125 λ 3 ;
  • f ( λ ) = 3355443200 λ + 2558365696 λ 2 + 643302144 λ 3 + 54034441 λ 4 134875 λ 5 ;
  • g ( λ ) = 74088 2197 42 169 ( 128 5 λ ) + 6 λ 13 ;
  • h ( λ ) = 4 441 169 + 1 156 ( 128 5 λ ) + 1 52 ( 128 + 5 λ ) + b ( λ ) + c ( λ ) .
Choose different λ-values; the number of critical points is also different:
  • R λ ( x ) has 13 critical points, when λ { 0 , 1344 } and λ C ;
  • R λ ( x ) has 7 critical points, when λ = 0 ;
  • R λ ( x ) has 11 critical points, when λ = 1344 ;
  • The free critical points of R λ ( x ) satisfy c r m = 1 c r n for m n ; each pair is conjugate to each other, being c r 2 and c r 3 , c r 4 and c r 5 , c r 6 and c r 7 , c r 8 and c r 9 , and c r 10 and c r 11 .
Proof. 
From (11), let R λ ( x ) = 0 :
2 x 7 ( 1 + x ) 6 ( 1 + x 2 ) 5 · ω λ ( x ) θ λ ( x ) 2 = 0 ,
where ω λ ( x ) = 52 + 168 x + ( 336 5 λ ) x 2 + 3 ( 168 + λ ) x 3 + 2 ( 284 + λ ) x 4 + 3 ( 168 + λ ) x 5 + ( 336 5 λ ) x 6 + 168 x 7 + 52 x 8 .
From (12), the free critical points of R λ ( x ) are 1 , ± i and the eight roots of polynomial ω λ ( x ) = 0 . Next is the process of solving the expression of the roots of the polynomial ω λ ( x ) = 0 .
First, we can express ω λ ( x ) as the product of four quadratic polynomials ( 1 + d i x + x 2 ) , 1 i 4 . We have the following formulas:
ω λ ( x ) = 52 i = 1 4 ( 1 + d i x + x 2 ) .
Next, by the corresponding coefficients are equal, we can obtain
52 ( d 1 + d 2 + d 3 + d 4 ) = 168 , 52 ( 4 + d 1 d 2 + d 1 d 3 + d 2 d 3 + d 1 d 4 + d 2 d 4 + d 3 d 4 ) = 336 5 λ , 52 ( 3 d 1 + 3 d 2 + 3 d 3 + 3 d 4 + d 1 d 2 d 3 + d 1 d 2 d 4 + d 1 d 3 d 4 + d 2 d 3 d 4 ) = 504 + 3 λ , 52 ( 6 + 2 d 1 d 2 + 2 d 1 d 3 + 2 d 2 d 3 + 2 d 1 d 4 + 2 d 2 d 4 + 2 d 3 d 4 + d 1 d 2 d 3 d 4 ) = 568 + 2 λ .
By eliminating d 2 , d 3 , and d 4 , we can obtain a quartic polynomial about d 1 :
52 d 1 4 168 d 1 3 + ( 128 5 λ ) d 1 2 + 3 λ d 1 + 12 λ = 0 .
After that, the specific expression of d i , 1 i 4 can be obtained by (15). Finally, by taking d i into (13), we can obtain the specific expression of the roots of ω λ ( x ) = 0 , that is c r 4 c r 11 .
Additionally,
R 0 ( x ) = 8 x 7 ( 1 + x ) 4 ( 1 + x 2 ) 3 ( 13 + 16 x + 13 x 2 ) ( 1 + 4 x + 9 x 2 + 16 x 3 + 21 x 4 + 20 x 5 + 13 x 6 ) 2 ,
R 1344 ( x ) = 8 x 7 ( 1 + x ) 6 ( 1 + x 2 ) 5 ( 13 + 68 x 1473 x 2 1880 x 3 1473 x 4 + 68 x 5 + 13 x 6 ) θ 1344 ( x ) 2 ,
where θ 1344 ( x ) = 1 + 7 x + 27 x 2 + 77 x 3 + 176 x 4 + 336 x 5 + 552 x 6 + 792 x 7 + 1015 x 8 + 1185 x 9 59 x 10 13 x 11 .
R 6848 ( x ) = 8 x 7 ( 1 + x ) 6 ( 1 + x 2 ) 5 ( 13 + 42 x + 8644 x 2 5010 x 3 3282 x 4 5010 x 5 + 8644 x 6 + 42 x 7 + 13 x 8 ) θ 6848 ( x ) 2 ,
where θ 6848 ( x ) = 1 + 6 x + 20 x 2 + 50 x 3 + 99 x 4 + 160 x 5 + 216 x 6 + 240 x 7 + 223 x 8 + 170 x 9 + 6948 x 10 + 46 x 11 + 13 x 12 . □
Table 2 shows the number of free critical points corresponding to R λ ( x ) , where λ is given a particular value. According to the above proposition, 0 and correspond to the roots of g ( x ) = 0 , and the number of free critical points of R λ ( x ) is at most 11 and at least 5. Since the free critical points c r 1 , c r 2 , and c r 3 are the pre-images of the strange fixed point e x 1 , the stability of c r 1 , c r 2 , and c r 3 corresponds to the stability of e x 1 . In addition, as with the strange fixed point in this section, each pair of conjugate free critical points has the same stability properties, so we only need to study the dynamics of half of the free critical points, which will be discussed in Section 2.3.

2.3. Parameter Spaces and Dynamical Planes

The dynamic behavior of the rational operator R λ ( x ) (4) varies with the parameter selection. In this section, we discuss the progressive behavior of the free critical points of L W M ( λ ) through parameter spaces. At first, the parameter space is represented by a 400 × 400 grid on the complex plane, and the points on the grid correspond to different λ -values. The parameter space shows the convergence analysis of the λ -related iterative method (3), where the free critical points c r i in Proposition 5 is used as the initial estimate value and the maximum number of iterations is 25. If it converges to 0, it is expressed as orange. If it converges to , it is expressed as faint yellow, and other cases are expressed as black.
Our goal was to find a relatively stable region in the parameter space, that is the orange and faint yellow regions, because the λ -values in these regions are the best parameter members in terms of numerical stability.
The L W M ( λ ) family has a maximum of eleven free critical points. Proposition 5 states that it is sufficient to study four different parameter spaces. These parameter spaces are called P 1 ( f o r x = c r 4 , c r 5 ) , P 2 ( f o r x = c r 6 , c r 7 ) , P 3 ( f o r x = c r 8 , c r 9 ) , and P 4 ( f o r x = c r 10 , c r 11 ) , as described in Figure 3. Among them, the large area of Figure 3 is orange and faint yellow, that is, in most cases, it converges to the roots. However, by carefully observing the details of the parameter spaces, we can see that the rare black area near the imaginary part of the parameter λ is 0. For example, the points corresponding to λ = 7 , λ = 11.68 , β = 9.69 , and β = 4.999 in P 1 P 4 are in the black regions. Combining Figure 1 and Figure 2, we know that only a quite small part of the points in the iterative method (3) converge to the strange fixed points, which corresponds to the results shown in P 1 P 4 , indicating that the members in L W M ( λ ) are numerically stable.
Then, we explored the stability of the corresponding iterative method by selecting various parameter values to generate a dynamical plane.
Firstly, a grid generated by 400 × 400 points is defined on the complex plane, where each point on the grid corresponds to a different value of the initial estimate x 0 . Its graphical representation indicates that this method has a maximum of 25 iterations to any root starting from x 0 . Among them, the white asterisk “∗” represents the attractive point. Fixed points are illustrated with a red circle “∘”. Periodic points are represented as blue squares “⋄”. In addition, we represent the orbit of a periodic point in blue and the orbit of a fixed point in red.
Figure 4, Figure 5, Figure 6 and Figure 7 show the dynamical planes of L W M ( λ ) corresponding to a given λ value. Different colors represent different basins of attraction: gray represents convergence to 0, khaki convergence to , green convergence to fixed point e x 1 = 1 , and black non-convergence.
Next, the stability of some L W M ( λ ) was studied by generating a dynamical plane. We discuss the choice of parameter λ in four parts. In the first part, we chose the parameter values in the black regions of the parameter spaces, that is λ = 7 , 11.68 , 9.69 , 4.999 ; see Figure 4.
In the second part, we selected the parameter values that satisfy | λ 1344 | < 8192 , which is divided into two parts: one is the smaller value, λ = 5 , 50 , 100 , and the other part is the larger value, λ = 1300 , 1340 , 1500 . At this time, e x 1 = 1 is the attractive point, which corresponds to the conclusion of Proposition 4 ( b ) ; see Figure 5.
In the third part, we drew λ = 7000 , 10,000, which correspond to the conclusion of the proposition in Figure 4a; see Figure 6.
Finally, we drew the dynamical plane of λ = 0 , 6848 , 1344 , where λ are the special values obtained in Proposition 2; see Figure 7.
In Figure 4, Figure 5, Figure 6 and Figure 7, the dynamical planes of most parameters converge to 0 and , but there are also parameter values that make the iterative method not converge, such as the parameter λ = 4.999 in Figure 4 and λ = 5 , 100 in Figure 5. Note from Figure 6 that the iterative method corresponding to the parameters satisfying | λ 1344 | > 8192 converges to two extra fixed points in addition to 0 , , e x 1 = 1 . This indicates that the iterative method corresponding to these parameters is unstable. What we want to study more is stable parameters, that is the dynamical planes only show gray and khaki.
Comparing Figure 4, Figure 5, and Figure 7, we can see that the larger the absolute value of the parameter values, the more complex the dynamical plane displayed. Carefully observing the black area in Figure 4d and Figure 5a,c, we can also observe that the larger the absolute value of the parameter values, the larger the black area is. In summary, the iterative method corresponding to λ = 0 is more stable.

3. Numerical Experiment and Practical Application

In this section, the more stable parameter elements obtained in Section 2, namely { 0 , 50 , 1300 , 1340 , 1500 } , are verified by a numerical experiment and applied to solve practical problem.
First, we used the iterative method corresponding to these parameters to solve the nonlinear matrix sign function:
X 2 I = 0 , I i s a n i d e n t i t y m a t r i x .
After that, we obtained the following special cases:
  • When λ = 0 , X n + 1 = 128 X n 5 ( I + X n 2 ) 4 [ I 9 X n 2 + 57 X n 4 + 247 X n 6 + 763 X n 8 + 637 X n 10 + 331 X n 12 + 21 X n 14 ] 1 ;
  • When λ = 50 , X n + 1 = 8192 X n 7 ( I + X n 2 ) 6 [ 50 ( I + X n 2 ) 10 + 64 ( X n + X n 3 ) 2 ( I 9 X n 2 + 57 X n 4 + 247 X n 6 + 763 X n 8 + 637 X n 10 + 331 X n 12 + 21 X n 14 ) ] 1 ;
  • When λ = 1300 , X n + 1 = 8192 X n 7 ( I + X n 2 ) 6 [ 1300 ( I + X n 2 ) 10 + 64 ( X n + X n 3 ) 2 ( I 9 X n 2 + 57 X n 4 + 247 X n 6 + 763 X n 8 + 637 X n 10 + 331 X n 12 + 21 X n 14 ) ] 1 ;
  • When λ = 1340 , X n + 1 = 8192 X n 7 ( I + X n 2 ) 6 [ 1340 ( I + X n 2 ) 10 + 64 ( X n + X n 3 ) 2 ( I 9 X n 2 + 57 X n 4 + 247 X n 6 + 763 X n 8 + 637 X n 10 + 331 X n 12 + 21 X n 14 ) ] 1 ;
  • When λ = 1500 , X n + 1 = 8192 X n 7 ( I + X n 2 ) 6 [ 1500 ( I + X n 2 ) 10 + 64 ( X n + X n 3 ) 2 ( I 9 X n 2 + 57 X n 4 + 247 X n 6 + 763 X n 8 + 637 X n 10 + 331 X n 12 + 21 X n 14 ) ] 1 .
Finally, we used the Matlab software to record the iteration times and computer running time of the experimental results of solving the matrix symbol function corresponding to the random square matrix. The calculation stopping standard is X n 2 I 2 10 100 . The experimental results are shown in Table 3.
Through the results in Table 3, we noticed that, in { 0 , 50 , 1300 , 1340 , 1500 } , when λ = 0 , the iterative method (3) is used to solve the matrix function with the least number of iterations and the shortest computer running time.
In 1873, van der Waals revised the equation of state of the ideal gas and proposed the van der Waals equation of state based on the two assumptions of the ideal gas model, namely the non-occupying volume of molecules and the non-interacting force between molecules. The volume occupied by molecules decreases their free space of movement, and the frequency of molecules hitting the container wall increases at the same temperature, so the pressure increases correspondingly. If V m b is used to represent the free movement space per mole of gas molecules, the gas pressure should be P = R T / ( V m b ) by referring to the ideal gas equation of state.
The specific form of Van der Waals equation is:
( P + a n 2 V 2 ) ( V n b ) = n R T
where a and b are the relative constants corrected for gas pressure and volume, respectively, which are called the van der Waals constants. Each gas has a specific value for a and b. n , R , T , V , and P represent the number of moles, the universal gas constant (0.0820578), the absolute temperature, the volume, and the absolute pressure, respectively.
By simplification, the following polynomial of the nonlinear form can be obtained:
H ( V ) = P V 3 ( R T + b P ) V 2 n + a V n 2 a b n 3 .
When n = 1.5 moles of benzene vapor, P = 50 a t m , a = 1.378 , and b = 0.3193 , then there is
H 1 ( V ) = 50 V 3 119.0185191 V 2 + 3.1005 V 1.484984475 .
The exact solution to H 1 ( V ) = 0 is
V = 2.3594236086921194141576714329391903003579204650023
(50 significant digits).
Next, we used the iterative method (3) with λ = 0 and the other three eighth-order iterative methods to solve H 1 ( V ) = 0 and compared the results. According to Theorem 1 in [21], when applying the iterative method L W M (3) to solve nonlinear equations, in order to maintain the convergence rate of the eighth order, we must select the initial point near the exact solution, so that the initial point is close enough to the exact solution. Therefore, we chose the initial estimate of x 0 = 2.36 . The experimental results are shown in Table 4.
The method by Sharma et al., S M (see [25]), is
w n = x n f ( x n ) f ( x n ) , z n = w n f ( w n ) f ( x n ) f ( x n ) + a f ( w n ) f ( x n ) + ( a 2 ) f ( w n ) , x n + 1 = x n P + Q + R P f [ z n , x n ] + Q f ( x n ) + R f [ w n , x n ] f ( x n ) .
where P = ( x n w n ) f ( x n ) f ( w n ) , Q = ( w n z n ) f ( w n ) f ( z n ) and R = ( z n x n ) f ( z n ) f ( x n ) . We chose a = 0 for the numerical experiment part.
The method by Soleymani, S S M (see [26]), is
y n = x n f ( x n ) f ( x n ) , z n = y n f ( y n ) f ( x n ) f ( x n ) + β f ( y n ) f ( x n ) + ( β 2 ) f ( y n ) , x n + 1 = z n f ( z n ) f ( x n ) { f ( x n ) + β f ( y n ) f ( x n ) + ( β 2 ) f ( y n ) } ( 1 + 2 t + μ + μ 2 + ρ 2 + 4 ρ 3 + 12 ρ 4 σ τ 2 ) .
where t = f ( z n ) f ( x n ) , μ = f ( z n ) f ( y n ) , ρ = f ( y n ) f ( x n ) , σ = f ( z n ) f ( x n ) , and τ = f ( y n ) f ( x n ) . We chose β = 0 for the numerical experiment part.
It can be seen from the experimental results in Table 4 that the L W M with λ = 0 is more dominant in convergence accuracy.

4. Conclusions

In this paper, as stated by fractal theory, we explored the complex dynamic behavior of a kind of iterative method (3) for solving nonlinear equations. Among them, Liu–Wang’s method is based on the well-known Ostrowski method, which adds a step to make it the optimal eighth-order iterative method with faster convergence speed.
Based on the theory of the Möbius conjugate transformation on the Riemann sphere and the scaling theorem, we obtained the corresponding rational operator by applying any quadratic polynomial to the iterative method. The operator can be simplified by giving parameters to facilitate the subsequent research. First, we analyzed its fixed points and the critical points and noticed that the strange fixed points always appeared in pairs, because the two strange points were conjugate, as were the free critical points. Given this relationship, then the number of studies on them could be reduced by half. All of this was performed with the Mathematica software. By drawing the stability surfaces of the strange fixed points, we can observe the range of parameters that can attract the sequence generated by the initial point to the strange fixed point during the iteration. Since we expected that these sequences can converge to the roots rather than the strange fixed points. Based on the above parameter range, we can more effectively select the desired parameters and avoid selecting bad parameters. What is more, the parameter space generated by the iteration with the critical point as the initial point was drawn, so as to obtain a stable parameter selection region. Combining the dynamical planes’ analysis of the given parameters, we finally obtained some parameter families that made the iterative method stable, such as λ = 0 , 50 , 1300 , 1340 , 1500 . When these parameters were selected, the corresponding iterative method had good numerical stability, especially at λ = 0 . Finally, this result was proven by solving the matrix sign function and the van der Waals equation.

Author Contributions

Methodology, X.W.; writing—original draft preparation, W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (No. 61976027), the Natural Science Foundation of Liaoning Province (No. 2022-MS-371), the Educational Commission Foundation of Liaoning Province of China (Nos. LJKMZ20221492 and LJKMZ20221498), and the Key Project of Bohai University (No. 0522xn078).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sidorov, N. Special Issue Editorial “Solvability of Nonlinear Equations with Parameters: Branching, Regularization, Group Symmetry and Solutions Blow-Up”. Symmetry 2022, 14, 226. [Google Scholar] [CrossRef]
  2. Hassan, M.; Muhammed, A.A. Globally convergent diagonal Polak-Ribière-Polyak like algorithm for nonlinear equations. Numer. Algorithms 2022, 91, 1441–1460. [Google Scholar]
  3. Wang, X.; Chen, X. Derivative-Free Kurchatov-Type Accelerating Iterative Method for Solving Nonlinear Systems: Dynamics and Applications. Fractal Fract. 2022, 6, 59. [Google Scholar] [CrossRef]
  4. 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]
  5. Nilay, K.R. A new semi-explicit atomistic molecular dynamics simulation method for membrane proteins. J. Comput. Methods Sci. Eng. 2019, 19, 259–286. [Google Scholar]
  6. Zhang, Z.; Choon, C.F.; Liu, G.R. A Semi-Explicit Finite Element Method for Dynamic Analysis of Dielectric Elastomers. Int. J. Comp. Meth-sing. 2015, 12, 1350108. [Google Scholar] [CrossRef]
  7. Zhao, X.; Liu, M.; Wang, P.; Guan, J.; Wang, X. Highly efficient preparation of multi-angle continuous carbon fibre reinforced hydroxyapatite composites by electrostatic splitting method. J. Eur. Ceram. Soc. 2022, 42, 7631–7647. [Google Scholar] [CrossRef]
  8. Dong, Y. A new splitting method for systems of monotone inclusions in Hilbert spaces. Math. Comput. Simulat. 2023, 203, 518–537. [Google Scholar] [CrossRef]
  9. Dai, L.; Singh, M.C. An analytical and numerical method for solving linear and nonlinear vibration problems. Int. J. Solids Struct. 1997, 34, 2709–2731. [Google Scholar] [CrossRef]
  10. 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]
  11. Cordero, A.; Moscoso-Martínez, M.; Torregrosa, J.R. Chaos and Stability in a New Iterative Family for Solving Nonlinear Equations. Algorithms 2021, 14, 101. [Google Scholar] [CrossRef]
  12. Cordero, A.; Leonardo, S.M.A.; Torregrosa, J.R. Dynamics and Stability on a Family of Optimal Fourth-Order Iterative Methods. Algorithms 2022, 15, 387. [Google Scholar] [CrossRef]
  13. Kanwar, V.; Bhatia, S.; Kansal, M. New optimal class of higher-order methods for multiple roots. College Math. J. 2013, 222, 564–574. [Google Scholar]
  14. Wang, X.; Zhang, T.; Qin, Y. Efficient two-step derivative-free iterative methods with memory and their dynamics. Int. J. Comput. Math. 2016, 93, 1423–1446. [Google Scholar] [CrossRef]
  15. Ortega, J.M.; Rheinbolt, W.C. Iterative Solution of Nonlinear Equations in Several Variables; Academic Press: New York, NY, USA, 1970. [Google Scholar]
  16. Kou, J.; Li, Y.; Wang, X. Modified Halley’s method free from second derivative. Appl. Math. Comput. 2006, 183, 704–708. [Google Scholar] [CrossRef]
  17. Neta, B. On Popovski’s method for nonlinear equations. Appl. Math. Comput. 2008, 201, 710–715. [Google Scholar] [CrossRef] [Green Version]
  18. Kung, H.T.; Traub, J.F. Optimal order of one-point and multipoint iteration. J. Assoc. Comput. Mach. 1974, 21, 634–651. [Google Scholar] [CrossRef]
  19. Husain, A.; Nanda, M.N.; Chowdary, M.S.; Sajid, M. Fractals: An Eclectic Survey, Part-I. Fractal Fract. 2022, 6, 89. [Google Scholar] [CrossRef]
  20. Husain, A.; Nanda, M.N.; Chowdary, M.S.; Sajid, M. Fractals: An Eclectic Survey, Part II. Fractal Fract. 2022, 6, 379. [Google Scholar] [CrossRef]
  21. Liu, L.; Xia, W. Eighth-order methods with high efficiency index for solving nonlinear equations. Appl. Math. Comput. 2010, 215, 3449–3454. [Google Scholar] [CrossRef]
  22. Lee, M.Y.; Kim, Y.I. The dynamical analysis of a uniparametric family of three-point optimal eighth-order multiple-root finders under the Möbius conjugacy map on the Riemann sphere. Numer. Algorithms 2020, 83, 1063–1090. [Google Scholar] [CrossRef]
  23. Wang, X.; Li, W. Choosing the Best Members of the Optimal Eighth-Order Petković’s Family by Its Fractal Behavior. Fractal Fract. 2022, 6, 749. [Google Scholar] [CrossRef]
  24. Amat, S.; Busquier, S.; Plaza, S. Chaotic dynamics of a third-order Newton-type method. J. Math. Anal. Appl. 2010, 366, 24–32. [Google Scholar] [CrossRef] [Green Version]
  25. Sharma, J.R.; Guha, R.K.; Gupta, P. Improved King’s methods with optimal order of convergence based on rational approximations. Appl. Math. Lett. 2013, 26, 473–480. [Google Scholar] [CrossRef]
  26. Soleymani, F.; Vanani, S.K.; Khan, M.; Sharifi, M. Some modifications of King’s family with optimal eighth order of convergence. Math. Comput. Model. 2012, 55, 1373–1380. [Google Scholar] [CrossRef]
Figure 1. Stability surface of e x 1 = 1 .
Figure 1. Stability surface of e x 1 = 1 .
Axioms 12 00215 g001
Figure 2. Stability surfaces of other strange fixed points.
Figure 2. Stability surfaces of other strange fixed points.
Axioms 12 00215 g002
Figure 3. Parameter spaces P i , i = 1 , 2 , 3 , 4 .
Figure 3. Parameter spaces P i , i = 1 , 2 , 3 , 4 .
Axioms 12 00215 g003
Figure 4. Dynamical planes for special λ -values.
Figure 4. Dynamical planes for special λ -values.
Axioms 12 00215 g004
Figure 5. Dynamical planes for special λ -values.
Figure 5. Dynamical planes for special λ -values.
Axioms 12 00215 g005
Figure 6. Dynamical planes for special λ -values.
Figure 6. Dynamical planes for special λ -values.
Axioms 12 00215 g006
Figure 7. Dynamical planes for special λ -values.
Figure 7. Dynamical planes for special λ -values.
Axioms 12 00215 g007
Table 1. Stability of strange fixed points e x i for special λ -values.
Table 1. Stability of strange fixed points e x i for special λ -values.
λ ex i No. of ex i
| R λ ( ex i ) | : t 1
0 1.62291 ± 0.455126 i 0.611529 ± 1.28544 i 0.234438 ± 1.48434 i 6
6.31492 : r 4.91253 : r 7.6094 : r
1344 2.57424 0.085032 0.123846 11
4.9157 : r 1.17175 × 10 6 : a 6.42351 × 10 6 : a
2.16337 ± 1.21147 i 1.09844 ± 2.04782 i 0.10715 ± 2.1662 i
4.67559 : r 4.68214 : r 5.53101:r
0.922374 ± 1.48996 i
5.75625 : r
6848 2.81972 ± 0.753267 i 1.9133 ± 1.99919 i 0.512724 ± 2.5536 i 12
4.07847 : r 3.98574 : r 4.3251 : r
0.00328839 ± 0.0431351 i 0.793402 ± 2.17265 i 1.45563 ± 0.903056 i
7.52394 × 10 9 : a 4.60277 : r 3.6229 : r
1 | R λ ( ex i ) | : t implies that exi is attractive, parabolic, and repulsive, if t = a ( | R | < 1 ) , t = a ( | R | = 1 ) , and t = a ( | R | > 1 ) , respectively.
Table 2. Free critical points c r i from R λ ( c r i ) for special λ -values.
Table 2. Free critical points c r i from R λ ( c r i ) for special λ -values.
λ cr i No. of cr i
0 1 ± i 0.615385 ± 0.788227 i 5
1344 1 ± i 13.0923 9
0.0763809 0.109819 9.10591
0.638919 ± 0.769274 i
6848 1 ± i 1.90256 ± 25.7656 i 11
0.639 ± 0.769207 i 0.00285034 ± 0.038601 i 0.929025 ± 0.370017 i
Table 3. Results of the comparisons.
Table 3. Results of the comparisons.
λ = 0 λ = 50 λ = 1300 λ = 1340 λ = 1500
Matricesn 1 TimenTimenTimenTimenTime
R 4 × 4 3 0.006994 5 0.079566 6 0.074327 6 0.103540 5 0.108107
R 5 × 5 3 0.016522 13 0.066437 13 0.094517 9 0.122004 16 0.129702
R 10 × 10 5 0.069684 10 0.139051 27 0.136659 41 0.157447 15 0.059954
R 15 × 15 6 0.070848 21 0.160829 22 0.187887 13 0.158099 86 0.279862
R 20 × 20 6 0.071416 19 0.180107 19 0.186666 17 0.162618 103 0.718463
R 25 × 25 9 0.111521 14 0.143166 20 0.209045 13 0.158114 97 0.380726
R 30 × 30 8 0.136161 28 0.203730 16 0.177654 15 0.176462 191 0.832457
1n represents the number of iterations.
Table 4. Comparison of iterative methods for H 1 ( V ) = 0 .
Table 4. Comparison of iterative methods for H 1 ( V ) = 0 .
n LWM SM SSM
| x n x n 1 | | f ( x n ) | | x n x n 1 | | f ( x n ) | | x n x n 1 | | f ( x n ) |
1 5.7639 × 10 4 1.818 × 10 21 1.1528 × 10 3 1.593 × 10 1 5.7639 × 10 4 1.4902 × 10 19
2 6.5748 × 10 24 4.1048 × 10 47 1.1528 × 10 3 1.5945 × 10 1 5.3893 × 10 22 1.6823 × 10 45
3 1.4846 × 10 49 9.2685 × 10 73 1.1528 × 10 3 1.593 × 10 1 6.0843 × 10 48 1.8993 × 10 71
4 3.352 × 10 75 2.0928 × 10 98 1.1528 × 10 3 1.5945 × 10 1 6.869 × 10 74 2.1442 × 10 97
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

Wang, X.; Li, W. Stability Analysis of Simple Root Seeker for Nonlinear Equation. Axioms 2023, 12, 215. https://doi.org/10.3390/axioms12020215

AMA Style

Wang X, Li W. Stability Analysis of Simple Root Seeker for Nonlinear Equation. Axioms. 2023; 12(2):215. https://doi.org/10.3390/axioms12020215

Chicago/Turabian Style

Wang, Xiaofeng, and Wenshuo Li. 2023. "Stability Analysis of Simple Root Seeker for Nonlinear Equation" Axioms 12, no. 2: 215. https://doi.org/10.3390/axioms12020215

APA Style

Wang, X., & Li, W. (2023). Stability Analysis of Simple Root Seeker for Nonlinear Equation. Axioms, 12(2), 215. https://doi.org/10.3390/axioms12020215

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