Next Article in Journal
Fractal Interpolation Using Harmonic Functions on the Koch Curve
Next Article in Special Issue
A New Approach for Dynamic Stochastic Fractal Search with Fuzzy Logic for Parameter Adaptation
Previous Article in Journal
Uncertainty Quantification of Random Microbial Growth in a Competitive Environment via Probability Density Functions
Previous Article in Special Issue
A Characterization of the Dynamics of Schröder’s Method for Polynomials with Two Roots
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Local Convergence and Dynamical Analysis of a Third and Fourth Order Class of Equation Solvers

by
Debasis Sharma
1,†,
Ioannis K. Argyros
2,*,†,
Sanjaya Kumar Parhi
3,† and
Shanta Kumari Sunanda
1,†
1
Department of Mathematics, IIIT Bhubaneswar, Odisha 751003, India
2
Department of Mathematical Sciences, Cameron University, Lawton, OK 73505, USA
3
Department of Mathematics, Fakir Mohan University, Odisha 756020, India
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 17 March 2021 / Revised: 31 March 2021 / Accepted: 2 April 2021 / Published: 5 April 2021
(This article belongs to the Special Issue Convergence and Dynamics of Iterative Methods: Chaos and Fractals)

Abstract

:
In this article, we suggest the local analysis of a uni-parametric third and fourth order class of iterative algorithms for addressing nonlinear equations in Banach spaces. The proposed local convergence is established using an ω -continuity condition on the first Fréchet derivative. In this way, the utility of the discussed schemes is extended and the application of Taylor expansion in convergence analysis is removed. Furthermore, this study provides radii of convergence balls and the uniqueness of the solution along with the calculable error distances. The dynamical analysis of the discussed family is also presented. Finally, we provide numerical explanations that show the suggested analysis performs well in the situation where the earlier approach cannot be implemented.

1. Introduction

The primary objective of this research is to provide an approximate solution α of:
F ( s ) = 0 ,
where F : O S 1 S 2 is a Fréchet derivable operator and O S 1 is nonempty, open and convex. S 1 and S 2 are Banach spaces. The requirement for the solutions of nonlinear equations in the form (1) is present in many problems of applied sciences and engineering as a foundation to solve other complicated ones. It is not always possible to determine the exact solutions to this kind of equations, so scientists and researchers frequently use iterative algorithms to approximate the required solution. The most regularly used iterative procedure for addressing (1) is Newton’s scheme that converges quadratically. Over the past several years, numerous authors have developed and are currently improving higher order iterative approaches [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17] for solving nonlinear equations.
Local convergence results for different iterative processes have been studied by many authors [18,19,20,21,22,23,24,25,26,27,28], and many important findings have been obtained. In these studies, essential results including convergence radii, measurements on error distances and the extended utility of efficient iterative procedures have been discussed.
In addition, analysis of dynamical characteristics of a class of iterative schemes applied on complex polynomials provides significant results about the stability and reliability of its elements. Researchers such as Amat et al. [29,30], Argyros and Magreñán [19,20], Cordero et al. [31,32,33,34] and others [27,35,36,37] have described complex dynamical behaviors of some famous classes of methods including Jarratt, King, Kim, Chebyshev–Halley, etc.
Recently, Argyros and Magreñán [19] studied the convergence analysis and dynamics of an optimal fourth order family given by:
t n = s n 2 3 F ( s n ) F ( s n ) , s n + 1 = s n [ ( β 2 22 β 27 ) F ( s n ) + 3 ( β 2 + 10 β + 5 ) F ( t n ) ] F ( s n ) 2 ( β F ( s n ) + 3 F ( t n ) ) ( 3 ( β + 1 ) F ( t n ) ( β + 5 ) F ( s n ) ) , n = 0 , 1 , 2 ,
where β C . Cordero et al. [34] developed a third and fourth order parametric class of nonlinear equation solvers and studied its stability analysis using a complex dynamics technique.
In this report, our primary aim is to analyze the local convergence and complex dynamical results for a third and fourth order convergent class of iterative procedures. Kou et al. [8] derived a family of iterative algorithms for addressing (1). This family is expressed as follows:
t n = s n 2 3 F ( s n ) 1 F ( s n ) , s n + 1 = s n I 3 4 A B F ( s n ) 1 F ( s n ) , n = 0 , 1 , 2 ,
where B = F ( s n ) 1 [ F ( t n ) F ( s n ) ] , A = I 3 4 θ B and θ R . The authors [8] also proved that the convergence rate of (3) is at least three, and for θ = 2 the convergence order increases to four. Observe that only F is involved in the iterative structure of (3). However, the convergence results of (3) were shown in [8] using conditions on the fourth order derivative. Such convergence criteria [8] restrict the usefulness of these solvers. As an example, we take a function F defined on Ω = [ 1 2 , 3 2 ] by:
F ( s ) = s 3 ln ( s 2 ) + s 5 s 4 , if s 0 0 , if s = 0 .
For the above function F , the third derivative F is not bounded on its domain O . The convergence theorem for the methods (3) derived in [8] does not work for this problem. Additionally, the authors proved the convergence theorem using the assumption that “if the initial estimation s 0 is sufficiently close to the solution s , then the sequence of iterates { s n } converges to s ”. But how close to the solution s the initial guess s 0 should be? The earlier work [8] offers no answer to this question. In other words, no information about the distance between s 0 and s was provided for ensuring the convergence of (3). The earlier convergence analysis was established to solve equations only in R . Additionally, complex dynamical analysis was not studied by the authors [8] to extract the stable schemes of the family (3). In this study, we establish the local convergence for the class of methods (3) using a set of conditions only on F in a general Banach space setting. More specifically, the ω -continuity of the first derivative is applied to eliminate the use of higher order derivatives in convergence study. The proposed analysis is useful to address the above discussed problem (4) as well as the problems where higher order derivatives of the involved operator do not exist. This analysis also enables us to apply the methods (3) for solving Banach space valued nonlinear operator equations. Hence, the utility of the family of methods (3) is expanded. This study also provides radii of convergence balls and the uniqueness of the solution along with the calculable error distances. Our technique is so general that it can be used to extend the applicability of other solvers [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37]. These are our motivations for this document. Furthermore, we analyze the dynamical properties of the parametric class (3). This analysis provides an idea to determine such values of the parameter θ whose corresponding iterative members are numerically stable. Additionally, various chaotic behaviors of the class, such as convergence to n-periodic orbits or attracting strange fixed points and divergence to , are presented using the techniques described in [20,37]. These anomalies are demonstrated with the help of a parameter space and some dynamical planes.
The rest of this document is structured as follows. Section 2 offers the local analysis of the discussed family of methods (3). The dynamical study of this family is explored in Section 3. Numerical illustrations are presented in the final section.

2. Local Convergence Analysis of the Discussed Family of Iterative Algorithms

We establish, with the help of the ω -continuity of F , the local convergence of the class of algorithms defined by (3). For c S 1 and ρ > 0 , we define the following notations:
B ( c , ρ ) = { s S 1 : | | c s | | < ρ } ,
B ¯ ( c , ρ ) = { s S 1 : | | c s | | ρ }
and
L ( S 1 , S 2 ) = { T : S 1 S 2 is bounded and linear } .
Let I = [ 0 , + ) and suppose there exists a non-decreasing and continuous function ω 0 : I I such that ω 0 ( 0 ) = 0 . Let the smallest positive solution D 0 of ω 0 ( y ) = 1 exist. Set I 0 = [ 0 , D 0 ) . Suppose there exist continuous and non-decreasing functions ω , V : I 0 I such that ω ( 0 ) = 0 and 1 3 V ( 0 ) < 1 . We introduce the following functions on the interval I 0 .
J 1 ( y ) = 0 1 ω ( ( 1 X ) y ) d X + 1 3 0 1 V ( X y ) d X 1 ω 0 ( y ) ,
H 1 ( y ) = J 1 ( y ) 1 ,
J 2 ( y ) = 1 1 ω 0 ( y ) [ 0 1 ω ( ( 1 X ) y ) d X + 3 4 ω ( ( J 1 ( y ) + 1 ) y ) 1 ω 0 ( y ) × 1 + 3 4 | θ | ω ( ( J 1 ( y ) + 1 ) y ) 1 ω 0 ( y ) × 0 1 V ( X y ) d X ]
and
H 2 ( y ) = J 2 ( y ) 1 .
Then, H i ( 0 ) < 0 and H i ( y ) + as y D 0 for i = 1 , 2 . It is confirmed by applying the intermediate value theorem that the smallest zeros D i of the functions H i exist and lie in ( 0 , D 0 ) , where i = 1 , 2 . Let us denote:
R = min { D 1 , D 2 } .
Hence, we have the following inequalities for each y [ 0 , R ) :
0 J 1 ( y ) < 1
and
0 J 2 ( y ) < 1 .
Furthermore, we assume that the following hold true for the Fréchet derivable operator F : O S 1 S 2 :
F ( α ) = 0 , F ( α ) 1 L ( S 2 , S 1 ) ,
| | F ( α ) 1 ( F ( s ) F ( α ) ) | | ω 0 ( | | s α | | ) , s O ,
| | F ( α ) 1 ( F ( s ) F ( t ) ) | | ω ( | | s t | | ) , s , t O 0 : = O B ( α , D 0 )
and
| | F ( α ) 1 F ( s ) | | V ( | | s α | | ) , s O 0 .
Next, we discuss the local convergence theorem for the family (3).
Theorem 1.
Suppose F : O S 1 S 2 is Fréchet derivable and α O . Let F obey the conditions defined in (10)–(13) and:
B ¯ ( α , R ) O .
Then, for any starter s 0 B ( α , R ) , the sequence of iterates { s n } n 0 constructed by the class of algorithms (3) is well defined, { s n } n 0 B ( α , R ) and lim n | | s n α | | = 0 . Moreover, the following items are true for all n 0 :
| | t n α | | J 1 ( | | s n α | | ) | | s n α | | < | | s n α | | < R
and
| | s n + 1 α | | J 2 ( | | s n α | | ) | | s n α | | < | | s n α | | < R ,
where the functions ‘ J ’ are discussed in Equations (5) and (6). In addition, if there exist D 3 R such that:
0 1 ω 0 ( X D 3 ) d X < 1 ,
then α satisfies Equation (1) uniquely in O 1 : = O B ¯ ( α , D 3 ) .
Proof. 
Let a B ( α , R ) . Using Equations (7) and (11), we obtain:
| | F ( α ) 1 ( F ( a ) F ( α ) ) | | ω 0 ( | | a α | | ) < ω 0 ( R ) < 1 .
Equation (18) together with the Banach lemma on invertible operators [1,2,13,16] imply that:
| | F ( a ) 1 F ( α ) | | 1 1 ω 0 ( | | a α | | ) < 1 1 ω 0 ( R ) .
It also follows from Formula (3) for n = 0 that t 0 is well defined. We also have from (3):
t 0 α = s 0 α F ( s 0 ) 1 F ( s 0 ) + 1 3 F ( s 0 ) 1 F ( s 0 ) = F ( s 0 ) 1 F ( α ) 0 1 F ( α ) 1 ( F ( α + X ( s 0 α ) ) F ( s 0 ) ) ( s 0 α ) d X + 1 3 F ( s 0 ) 1 F ( s 0 ) .
Then, we deduce using (3), (7), (8), (11), (19) (for a = s 0 ) and (20),
| | t 0 α | | | | F ( s 0 ) 1 F ( α ) | | | | 0 1 F ( α ) 1 ( F ( α + X ( s 0 α ) ) F ( s 0 ) ) ( s 0 α ) d X | | + 1 3 | | F ( s 0 ) 1 F ( s 0 ) | | | | F ( s 0 ) 1 F ( α ) | | | | 0 1 F ( α ) 1 ( F ( α + X ( s 0 α ) ) F ( s 0 ) ) ( s 0 α ) d X | | + 1 3 | | F ( s 0 ) 1 F ( α ) | | | | F ( α ) 1 F ( s 0 ) | | 0 1 ω ( ( 1 X ) | | s 0 α | | ) d X | | s 0 α | | 1 ω 0 ( | | s 0 α | | ) + 1 3 0 1 V ( X | | s 0 α | | ) d X | | s 0 α | | 1 ω 0 ( | | s 0 α | | ) = 0 1 ω ( ( 1 X ) | | s 0 α | | ) d X + 1 3 0 1 V ( X | | s 0 α | | ) d X 1 ω 0 ( | | s 0 α | | ) | | s 0 α | | = J 1 ( | | s 0 α | | ) | | s 0 α | | < | | s 0 α | | < R .
Hence, (15) holds for n = 0 . Next, we find a bound for | | B | | . We use Equations (12), (19) (for a = s 0 ) and (21) to obtain:
| | B | | = | | F ( s 0 ) 1 ( F ( t 0 ) 1 F ( s 0 ) 1 ) | | | | F ( s 0 ) 1 F ( α ) | | | | F ( α ) 1 ( F ( t 0 ) 1 F ( s 0 ) 1 ) | | ω ( | | t 0 s 0 | | ) 1 ω 0 ( | | s 0 α | | ) ω ( | | t 0 α | | + | | s 0 α | | ) 1 ω 0 ( | | s 0 α | | ) ω ( J 1 ( | | s 0 α | | ) | | s 0 α | | + | | s 0 α | | ) 1 ω 0 ( | | s 0 α | | ) = ω [ ( J 1 ( | | s 0 α | | ) + 1 ) | | s 0 α | | ] 1 ω 0 ( | | s 0 α | | ) .
It follows from the above inequality that:
| | A | | 1 + 3 4 | θ | ω [ ( J 1 ( | | s 0 α | | ) + 1 ) | | s 0 α | | ] 1 ω 0 ( | | s 0 α | | ) .
Additionally, s 1 is well defined by virtue of Equation (19). Next, we use (3), (6), (7), (9), (12), (19) (for a = s 0 ), (21), (22) and (23) to establish:
| | s 1 α | | | | F ( s 0 ) 1 F ( α ) | | | | 0 1 F ( α ) 1 ( F ( α + X ( s 0 α ) ) F ( s 0 ) ) ( s 0 α ) d X | | + 3 4 | | A | | | | B | | | | F ( s 0 ) 1 F ( s 0 ) | | 0 1 ω ( ( 1 X ) | | s 0 α | | ) d X 1 ω 0 ( | | s 0 α | | ) + 3 4 | | A | | | | B | | | | F ( s 0 ) 1 F ( α ) | | | | F ( α ) 1 F ( s 0 ) | | 0 1 ω ( ( 1 X ) | | s 0 α | | ) d X 1 ω 0 ( | | s 0 α | | ) + 3 4 1 + 3 4 | θ | ω [ ( J 1 ( | | s 0 α | | ) + 1 ) | | s 0 α | | ] 1 ω 0 ( | | s 0 α | | ) × ω [ ( J 1 ( | | s 0 α | | ) + 1 ) | | s 0 α | | ] 1 ω 0 ( | | s 0 α | | ) × 0 1 V ( X | | s 0 α | | ) d X 1 ω 0 ( | | s 0 α | | ) 1 1 ω 0 ( | | s 0 α | | ) [ 0 1 ω ( ( 1 X ) | | s 0 α | | ) d X + 3 4 1 + 3 4 | θ | ω [ ( J 1 ( | | s 0 α | | ) + 1 ) | | s 0 α | | ] 1 ω 0 ( | | s 0 α | | ) × ω [ ( J 1 ( | | s 0 α | | ) + 1 ) | | s 0 α | | ] 1 ω 0 ( | | s 0 α | | ) × 0 1 V ( X | | s 0 α | | ) d X ] = J 3 ( | | s 0 α | | ) | | s 0 α | | < | | s 0 α | | < R .
This validates the estimate (16) for n = 0 . The induction for the estimates (15) and (16) is completed for n = 0 . Assuming that (15) and (16) hold true for all j = 1 , 2 , , n 1 and then by replacing the estimations, we show the items (15) and (16) are true for j = n . From the estimation | | s n + 1 α | | J 2 ( R ) | | s n α | | < R , we derive that s n + 1 B ( α , R ) and lim n | | s n α | | = 0 . Let there exist c ( α ) O 1 such that F ( c ) = 0 . Consider W = 0 1 F ( α + X ( c α ) ) d X . From (11) and (17), we obtain:
| | F ( α ) 1 ( W F ( α ) ) | | 0 1 ω 0 ( X | | c α | | ) d X 0 1 ω 0 ( X D 3 ) d X < 1 ,
and this confirms that W 1 L ( S 2 , S 1 ) . We obtain α = c due to the invertability of W and the identity 0 = F ( α ) F ( c ) = W ( α c ) . □

3. Complex Dynamics of the Discussed Class of Methods

In this section, the complex dynamical analysis of the class (3)is carried out. When analyzing an iterative formula, convergence analysis is not the only required matter to study. The effectiveness of an algorithm also depends on other components, such as verifying how it performs on the basis of the chosen starting points. That is, how large the set of starting points is for which the scheme converges to a solution. For this reason, a wide range of tools are required for a more extensive analysis.
In the research of iterative procedures, the dynamical analysis of a family of iterative algorithms is now gaining increasing interest because it enables separating different iterative approaches in terms of their rate of convergence. It also allows for evaluating their behavior depending on the chosen starting point. This analysis helps one to visualize the set of initial estimates that converge to a root of the equation or other points. Furthermore, it demonstrates the reliability and productiveness of the iterative process.
This research investigates the complex dynamical behavior of the class of methods (3) on a polynomial N ( z ) : C C defined by N ( z ) = ( z v 1 ) ( z v 2 ) . We discuss the detailed study of the fixed points and their stability employing MATHEMATICA software [19,20]. In addition, a parameter space and various dynamical planes were constructed and analyzed to show the chaotic nature of the considered family.
We applied the class of schemes (3) on the polynomial N ( z ) = ( z v 1 ) ( z v 2 ) , where v 1 v 2 . Then, we considered the Möbius map M ( z ) = z v 1 z v 2 with the properties M ( v 1 ) = 0 , M ( v 2 ) = and M ( ) = 1 to deduce the rational operator:
P N ( z , θ ) = z 3 z 3 + 4 z 2 + 5 z + 2 θ ( 2 θ ) z 3 + 5 z 2 + 4 z + 1 , θ C .
related to the class of algorithms (3).

3.1. Study of Fixed Points and Their Stability

The fixed points of the operator P N ( z , θ ) are obtained by solving the equation P N ( z , θ ) = z . The points z = 0 and z = are superattracting fixed points of P N ( z , θ ) as they are related to v 1 and v 2 , respectively. According to the definition, z = 1 is a strange fixed point. Moreover, depending on θ , there exist another four strange fixed points. These points can be expressed as follows:
1.
e x f 1 ( θ ) = 5 4 + 1 4 1 4 θ + 1 4 10 10 1 4 θ 4 θ
2.
e x f 2 ( θ ) = 5 4 + 1 4 1 4 θ 1 4 10 10 1 4 θ 4 θ
3.
e x f 3 ( θ ) = 5 4 1 4 1 4 θ + 1 4 10 + 10 1 4 θ 4 θ
4.
e x f 4 ( θ ) = 5 4 1 4 1 4 θ 1 4 10 + 10 1 4 θ 4 θ
These strange fixed points are presented in Figure 1 for e x f j R , j = 1 , 2 , 3 , 4 .
The relations among e x f j ( θ ) ( j = 1 , 2 , 3 , 4 ) are discussed below.
Lemma 1.
The fixed point operator P N ( z , θ ) has five simple strange fixed points except in the following situations:
(i) 
If θ = 0 , then P N ( z , θ ) has three strange fixed points.
(ii) 
If θ = 20 , then e x f 1 ( θ ) = e x f 2 ( θ ) = 1 and the number of strange fixed points is three.
(iii) 
If θ = 1 4 , then e x f 1 ( θ ) = e x f 3 ( θ ) = 1 2 and e x f 2 ( θ ) = e x f 4 ( θ ) = 2 . Hence the number of strange fixed points of P N ( z , θ ) is three.
Furthermore, e x f 1 ( θ ) = 1 e x f 2 ( θ ) and e x f 3 ( θ ) = 1 e x f 4 ( θ ) for all values of the parameter θ.
To analyze the stability of these strange fixed points and obtain the critical points, we require the first order derivative of P N ( z , θ ) . The simplified expression of the differentiated operator is:
P N ( z , θ ) = 3 z 2 ( z + 1 ) 4 ( ( θ 2 ) z 2 + ( 4 3 θ 4 ) z + θ 2 ) ( 1 + ( 2 θ ) z 3 + 5 z 2 + 4 z ) 2 .
It is easy to observe from Equation (25) that the superattracting fixed points of P N ( z , θ ) are z = 0 and z = . The stability analysis of z = 1 and z = e x f j ( θ ) , j = 1 , 2 , 3 , 4 will be provided in a different approach. We begin with investigating the stability of z = 1 . This fixed point is associated with the convergence to . We denote the stability function | P N ( 1 , θ ) | of the fixed point z = 1 by S t 1 ( θ ) . Using Equation (25), we get:
S t 1 ( θ ) = 32 θ 12
Theorem 2.
The strange fixed point z = 1 is categorized as follows:
1. 
If | θ 12 | > 32 , then z = 1 is attracting. In addition, it cannot be superattracting for any value of θ.
2. 
If | θ 12 | = 32 , then z = 1 is a parabolic fixed point.
3. 
z = 1 is a repulsor when | θ 12 | < 32 .
Therefore, the stability zone of z = 1 is the region on C where | θ 12 | 32 . This stability zone is presented graphically in Figure 2.
It is incredibly challenging to determine the stability of other fixed points ( e x f j , j = 1 , 2 , 3 , 4 ) analytically. However, we have used the graphical tool MATHEMATICA to display the stability zones for these points. It is noticed that the stability of e x f 1 ( θ ) and e x f 2 ( θ ) is the same. The stability regions of the strange fixed points e x f 1 ( θ ) and e x f 2 ( θ ) are presented in the Figure 3a,b, respectively. For the fixed points e x f 3 ( θ ) and e x f 4 ( θ ) , corresponding stability functions are demonstrated in Figure 4a,b, respectively.
Considering the stability regions presented in the Figure 3 and Figure 4, the following results are established.
Theorem 3.
In terms of stability, the strange fixed points e x f j ( θ ) , j = 1 , 2 , 3 , 4 are characterized as follows:
1. 
The stability of e x f 1 ( θ ) and e x f 2 ( θ ) is the same. In detail,
   * If θ { 10.69380 , 0.24936 } , then both fixed points are superattractive fixed points.
    * Additionally, these fixed points are attractors if θ lies in the oval or cardioid presented in Figure 3a,b, respectively.
2. 
The fixed points e x f 3 ( θ ) and e x f 4 ( θ ) are always repulsive in nature for any value of θ, and this can be observed in Figure 4a,b, respectively.

3.2. Study of Critical Points and Parameter Spaces

We next discuss the critical points of the class (3) with the parameter space. We obtained the critical points of P N ( z , θ ) by solving P N ( z , θ ) = 0 . It is clear from (25) that 0 and are critical points of P N ( z , θ ) that are associated with v 1 and v 2 , respectively. The point z = 1 is a free critical point of P N ( z , θ ) . The other free critical points are represented by f c r j ( θ ) , j = 1 , 2 . The expressions of these free critical points f c r j ( θ ) are as follows:
1.
f c r 1 ( θ ) = 1 3 2 θ + 6 + 5 θ 2 + 60 θ θ 2
2.
f c r 2 ( θ ) = 1 3 2 θ + 6 5 θ 2 + 60 θ θ 2
The number of free critical points of P N ( z , θ ) depends on θ , as summarized below.
Lemma 2.
(i) 
If θ { 0 , 2 , 12 } , then the number of distinct free critical points of P N ( z , θ ) is one.
(ii) 
For the rest of scenarios, the operator P N ( z , θ ) has three distinct free critical points.
We observe that f c r 1 ( θ ) = 1 f c r 2 ( θ ) . Hence, f c r 1 ( θ ) and f c r 2 ( θ ) are not independent. This is why we present the parameter plane related to f c r 1 ( θ ) . The behavior of z = 1 and f c r j ( θ ) (when f c r j R , j = 1 , 2 ) is provided in Figure 5.
The dynamical properties of the considered family (3) were analyzed using the procedure discussed in [19,20]. The parameter space associated to f c r 1 ( θ ) is given in Figure 6. The detailed views of this parameter space are displayed in Figure 7a,b. The free critical point z 0 = f c r 1 ( θ ) was considered as a stating estimation for the iterative processes of the class. We applied different colors on the initial choice z 0 according to the convergence of the corresponding sequence of iterates. Cyan is used to show the convergence of the iteration sequence to 0 (related to v 1 ) or (related to the v 2 ). We assigned the yellow color to z 0 to show the convergence to z = 1 (related to ). We executed a maximum of 2000 iterations with 10 6 as the error tolerance. Moreover, the magenta color indicates the convergence to e x f j ( θ ) , j = 1 , 2 , 3 , 4 . Convergences to n-periodic orbits for n = 2 , 3 , 4 , 8 are displayed in other colors, including light green, orange, blue, dark orange, dark green, dark red and white, respectively. For n 9 the n-periodic orbits are shown in black. In the parameter plane, the points that are not painted with cyan colors are not the appropriate choices for θ from the numerical point of view. In these zones, the iteration converges to any of the e x f j ( θ ) , j = 1 , 2 , 3 , 4 or to attracting n-periodic orbits or even to . In addition, one can find broad regions where the points are displayed in cyan. This verifies that some iterative procedures of the discussed family are numerically stable.
Next, we provide dynamical planes to display some of the detected anomalies. In these figures, we use cyan and red to demonstrate the convergence to 0 and , respectively. Black indicates that the iterative member of the family (3) does not converge to either 0 or with the tolerance error 10 6 and at most 1000 iterations. In the dynamical planes shown in Figure 8a,b, we have used yellow to present the convergence of the iterative member to the strange fixed point z = 1 . In Figure 8c,d, attracting strange fixed points e x f 1 and e x f 2 appear, which are displayed in magenta and green, respectively.
In Figure 9a, the dynamical plane related to a scheme of the discussed family (for θ = 20 ) shows the appearance of an attracting 2-cycle { 0.6143 + 0.7891 i , 0.6143 0.7891 i } . Moreover, various attracting 2-periodic orbits are given in Figure 9b,c. In Figure 9b, the existence of two 2-periodic orbits { 0.2674 0.9636 i , 0.0854 0.9963 i } and { 0.2674 + 0.9636 i , 0.0854 + 0.9963 i } is seen. Additionally, the 2-cycles { 0.0015 1.0000 i , 0.9001 + 0.4357 i } and { 0.0015 + 1.0000 i , 0.9001 0.4357 i } are presented in Figure 9c. The occurrence of an attracting 3-cycle is provided in Figure 9d. The convergence to n-cycles is presented in black, since this convergence is not related to v 1 and v 2 .
At the end, we provide dynamical planes for θ = 0.5 , θ = 0 , θ = 0.2 and θ = 2 in Figure 10a–d, respectively. The attraction basins in these planes are related to v 1 or v 2 only; as a consequence, the corresponding elements of the considered family are numerically stable. Thus, these methods are preferable over other elements of the family in most numerical applications.

4. Numerical Examples

In this section, we discuss numerical problems to demonstrate the efficacy of our theoretical findings. In these examples, we apply the proposed results to calculate the convergence radius of four methods derived from the discussed family (3) by using θ = 0.5 , θ = 0 , θ = 0.2 and θ = 2 . Notice that no iterative method is used to find the radius R . However, only Formula (7) is related to the solutions of scalar equations. Iterative method (3) is not used, since the solutions are known. The obtained radii show the degree of difficulty in choosing initial points.
Example 1.
[22] Let F be defined on R 3 for ( s 1 , s 2 , s 3 ) t by:
F ( s ) = ( 10 s 1 + sin ( s 1 + s 2 ) 1 , 8 s 2 cos 2 ( s 3 s 2 ) 1 , 12 s 3 + sin ( s 3 ) 1 ) t
We have α = ( 0.0689783491 , 0.2464424186 , 0.0769289119 ) t , ω 0 ( y ) = ω ( y ) = 0.269812 y and V ( y ) = 2 . The values of R are presented in Table 1.
Example 2.
[22] We chose the nonlinear Hammerstein type integral equation given by:
F ( s ) ( x ) = s ( x ) 5 0 1 s u s ( u ) 3 d u ,
where s ( x ) C [ 0 , 1 ] . We have α ( x ) = 0 , ω 0 ( y ) = 7.5 y , ω ( y ) = 15 y and V ( y ) = 2 . The values of R are provided in Table 2.
Example 3.
[25] Finally, we address the motivational problem discussed in the introduction part. Then α = 1 , ω 0 ( y ) = ω ( y ) = 96.662907 y and V ( y ) = 2 . The values of R are given in Table 3.

Author Contributions

Conceptualization, D.S. and I.K.A.; methodology, D.S., I.K.A., S.K.P. and S.K.S.; software, D.S. and I.K.A.; validation, D.S., I.K.A., S.K.P. and S.K.S.; formal analysis, D.S., I.K.A., S.K.P. and S.K.S.; investigation, D.S., I.K.A., S.K.P. and S.K.S.; resources, D.S., I.K.A., S.K.P. and S.K.S.; data curation, S.K.P. and S.K.S.; writing—original draft preparation, D.S., I.K.A., S.K.P. and S.K.S.; writing—review and editing, D.S. and I.K.A.; visualization, D.S., I.K.A., S.K.P. and S.K.S.; supervision, I.K.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by University Grants Commission of India grant number NOV2017-402662.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Argyros, I.K. Convergence and Application of Newton-Type Iterations; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  2. Argyros, I.K.; Cho, Y.J.; Hilout, S. Numerical Methods for Equations and Its Applications; Taylor & Francis, CRC Press: New York, NY, USA, 2012. [Google Scholar]
  3. Behl, R.; Cordero, A.; Motsa, S.S.; Torregrosa, J.R.; Kanwar, V. An optimal fourth-order family of methods for multiple roots and its dynamics. Appl. Math. Comput. 2016, 71, 775–796. [Google Scholar] [CrossRef]
  4. Candela, V.; Marquina, A. Recurrence relations for rational cubic methods I: The Halley method. Computing 1990, 44, 169–184. [Google Scholar] [CrossRef]
  5. Ezquerro, J.; Hernández, M.A. On Halley-type iteration with free second derivative. J. Comput. Appl. Math. 2004, 170, 455–459. [Google Scholar] [CrossRef] [Green Version]
  6. Ezquerro, J.A.; González, D.; Hernández, M.A. Majorizing sequences for Newton’s method from initial value problems. J. Comput. Appl. Math. 2012, 2012 236, 2246–2258. [Google Scholar] [CrossRef] [Green Version]
  7. Grau, M.; Diaz-Barrero, J.L. An improvement of the Euler-Chebyshev iterative method. J. Math. Anal. Appl. 2006, 315, 1–7. [Google Scholar] [CrossRef] [Green Version]
  8. Kou, J.; Li, Y.; Wang, X. Fourth-order iterative methods free from second derivative. Appl. Math. Comput. 2007, 184, 880–885. [Google Scholar] [CrossRef]
  9. Maroju, P.; Magreñán, Á.A.; Motsa, S.S.; Sarría, Í. Second derivative free sixth order continuation method for solving nonlinear equations with applications. J. Math. Chem. 2018, 56, 2099–2116. [Google Scholar] [CrossRef]
  10. Neta, B.; Scott, M.; Chun, C. Basins of attraction for several methods to find simple roots of nonlinear equations. Appl. Math. Comput. 2012, 218, 10548–10556. [Google Scholar] [CrossRef]
  11. Özban, A.Y. Some new variants of Newton’s method. Appl. Math. Lett. 2004, 17, 677–682. [Google Scholar] [CrossRef] [Green Version]
  12. Petković, M.S.; Neta, B.; Petković, L.; Dz̃unić, D. Multipoint Methods for Solving Nonlinear Equations; Elsevier: Amsterdam, The Netherlands, 2013. [Google Scholar]
  13. Rall, L.B. Computational Solution of Nonlinear Operator Equations; Robert E. Krieger: New York, NY, USA, 1979. [Google Scholar]
  14. Ren, H.; Wu, Q.; Bi, W. New variants of Jarratt method with sixth-order convergence. Numer. Algor. 2009, 52, 585–603. [Google Scholar] [CrossRef]
  15. Sharma, J.R.; Guna, R.K.; Sharma, R. Efficient Jarratt-like methods for solving systems of nonlinear equations. Calcolo 2014, 51, 193–210. [Google Scholar] [CrossRef]
  16. Traub, J.F. Iterative Methods for Solution of Equations; Prentice-Hall: Englewood Cliffs, NJ, USA, 1964. [Google Scholar]
  17. Weerakoon, S.; Fernando, T.G.I. A variant of Newton’s method with accelerated third-order convergence. Appl. Math. Lett. 2000, 13, 87–93. [Google Scholar] [CrossRef]
  18. Amat, S.; Argyros, I.K.; Busquier, S.; Hernández-Verón, M.A.; Martínez, E. On the local convergence study for an efficient k-step iterative method. J. Comput. Appl. Math. 2018, 343, 753–761. [Google Scholar] [CrossRef]
  19. 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]
  20. Argyros, I.K.; Magreñán, Á.A. A study on the local convergence and the dynamics of Chebyshev-Halley-type methods free from second derivative. Numer. Algor. 2015, 71, 1–23. [Google Scholar] [CrossRef]
  21. Argyros, I.K.; Cho, Y.J.; George, S. Local convergence for some third order iterative methods under weak conditions. J. Korean Math. Soc. 2016, 53, 781–793. [Google Scholar] [CrossRef] [Green Version]
  22. Argyros, I.K.; George, S. Local convergence for an almost sixth order method for solving equations under weak conditions. SeMA J. 2017, 75, 163–171. [Google Scholar] [CrossRef]
  23. Argyros, I.K.; Sharma, D.; Parhi, S.K. On the local convergence of Weerakoon-Fernando method with ω continuity condition in Banach spaces. SeMA J. 2020, 77, 291–304. [Google Scholar] [CrossRef]
  24. Argyros, I.K.; George, S. On the complexity of extending the convergence region for Traub’s method. J. Complex. 2020, 56, 101423. [Google Scholar] [CrossRef]
  25. Argyros, I.K.; Behl, R.; González, D.; Motsa, S.S. Ball convergence for combined three-step methods under generalized conditions in Banach space. Stud. Univ. Babes-Bolyai Math. 2020, 65, 127–137. [Google Scholar] [CrossRef]
  26. Maroju, P.; Magreñán, Á.A.; Sarría, Í.; Kumar, A. Local convergence of fourth and fifth order parametric family of iterative methods in Banach spaces. J. Math. Chem. 2020, 58, 686–705. [Google Scholar] [CrossRef]
  27. Sharma, D.; Parhi, S.K. Local Convergence and Complex Dynamics of a Uni-parametric Family of Iterative Schemes. Int. J. Appl. Comput. Math. 2020, 6, 1–16. [Google Scholar] [CrossRef]
  28. Singh, S.; Gupta, D.K.; Badoni, R.P.; Martínez, E.; Hueso, J.L. Local convergence of a parameter based iteration with Hölder continuous derivative in Banach spaces. Calcolo 2017, 54, 527–539. [Google Scholar] [CrossRef]
  29. Amat, S.; Busquier, S.; Plaza, S. Review of some iterative root-finding methods from a dynamical point of view. Sci. Ser. A Math. Sci. 2004, 10, 3–35. [Google Scholar]
  30. Amat, S.; Busquier, S.; Plaza, S. Dynamics of the King and Jarratt iterations. Aequationes Math. 2005, 69, 212–223. [Google Scholar] [CrossRef]
  31. Cordero, A.; García-Maimó, J.; Torregrosa, J.R.; Vassileva, M.P.; Vindel, P. Chaos in King’s iterative family. Appl. Math. Lett. 2013, 26, 842–848. [Google Scholar] [CrossRef] [Green Version]
  32. Cordero, A.; Torregrosa, J.R.; Vindel, P. Dynamics of a family of Chebyshev-Halley type methods. Appl. Math. Comput. 2013, 219, 8568–8583. [Google Scholar] [CrossRef] [Green Version]
  33. Cordero, A.; Guasp, L.; Torregrosa, J.R. Choosing the most stable members of Kou’s family of iterative methods. J. Comput. Appl. Math. 2018, 330, 759–769. [Google Scholar] [CrossRef]
  34. Cordero, A.; Villalba, E.G.; Torregrosa, J.R.; Triguero-Navarro, P. Convergence and Stability of a Parametric Class of Iterative Schemes for Solving Nonlinear Systems. Mathematics 2021, 9, 86. [Google Scholar] [CrossRef]
  35. Chicharro, F.; Cordero, A.; Gutiérrez, J.M.; Torregrosa, J.R. Complex dynamics of derivative-free methods for nonlinear equations. Appl. Math. Comput. 2013, 2013 219, 7023–7035. [Google Scholar] [CrossRef]
  36. Chicharro, F.; Cordero, A.; Torregrosa, J.R. Drawing dynamical and parameters planes of iterative families and methods. Sci. World J. 2013, 2013, 780153. [Google Scholar] [CrossRef] [PubMed]
  37. Magreñán, Á.A. Different anomalies in a Jarratt family of iterative root-finding methods. Appl. Math. Comput. 2014, 233, 29–38. [Google Scholar]
Figure 1. Behavior of the strange fixed points.
Figure 1. Behavior of the strange fixed points.
Fractalfract 05 00027 g001
Figure 2. Stability region of z = 1 .
Figure 2. Stability region of z = 1 .
Fractalfract 05 00027 g002
Figure 3. Stability regions. (a) Stability region of e x f 1 ( θ ) and e x f 2 ( θ ) , (b) Stability region of e x f 1 ( θ ) and e x f 2 ( θ ) .
Figure 3. Stability regions. (a) Stability region of e x f 1 ( θ ) and e x f 2 ( θ ) , (b) Stability region of e x f 1 ( θ ) and e x f 2 ( θ ) .
Fractalfract 05 00027 g003
Figure 4. Stability functions. (a) Stability function for e x f 3 ( θ ) , (b) Stability function for e x f 4 ( θ ) .
Figure 4. Stability functions. (a) Stability function for e x f 3 ( θ ) , (b) Stability function for e x f 4 ( θ ) .
Fractalfract 05 00027 g004
Figure 5. Behavior of free critical points.
Figure 5. Behavior of free critical points.
Fractalfract 05 00027 g005
Figure 6. Parameter plane associated with the free critical point f c r 1 ( θ ) .
Figure 6. Parameter plane associated with the free critical point f c r 1 ( θ ) .
Fractalfract 05 00027 g006
Figure 7. (a,b) Detailed views of the parameter plane associated with f c r 1 ( θ ) .
Figure 7. (a,b) Detailed views of the parameter plane associated with f c r 1 ( θ ) .
Fractalfract 05 00027 g007
Figure 8. Dynamical Planes. (a) θ = 44.5, (b) θ = −22, (c) θ = −10, (d) θ = 0.249.
Figure 8. Dynamical Planes. (a) θ = 44.5, (b) θ = −22, (c) θ = −10, (d) θ = 0.249.
Fractalfract 05 00027 g008
Figure 9. Dynamical Planes. (a) θ = 20, (b) θ = 15.5, (c) θ = −6.5, (d) θ = 25 − 26i.
Figure 9. Dynamical Planes. (a) θ = 20, (b) θ = 15.5, (c) θ = −6.5, (d) θ = 25 − 26i.
Fractalfract 05 00027 g009
Figure 10. Dynamical planes. (a) θ = 0.5 , (b) θ = 0 , (c) θ = 0.2 , (d) θ = 2 .
Figure 10. Dynamical planes. (a) θ = 0.5 , (b) θ = 0 , (c) θ = 0.2 , (d) θ = 2 .
Fractalfract 05 00027 g010
Table 1. Radii of convergence balls for Example 1.
Table 1. Radii of convergence balls for Example 1.
θ = 0.5 θ = 0 θ = 0.2 θ = 2
D 1 = 0.823619 D 1 = 0.823619 D 1 = 0.823619 D 1 = 0.823619
D 2 = 0.665576 D 2 = 0.722293 D 2 = 0.696849 D 2 = 0.568151
R = 0.665576 R = 0.722293 R = 0.696849 R = 0.568151
Table 2. Radii of convergence balls for Example 2.
Table 2. Radii of convergence balls for Example 2.
θ = 0.5 θ = 0 θ = 0.2 θ = 2
D 1 = 0.022222 D 1 = 0.022222 D 1 = 0.022222 D 1 = 0.022222
D 2 = 0.014295 D 2 = 0.015806 D 2 = 0.015118 D 2 = 0.011853
R = 0.014295 R = 0.015806 R = 0.015118 R = 0.011853
Table 3. Radii of convergence balls for Example 3.
Table 3. Radii of convergence balls for Example 3.
θ = 0.5 θ = 0 θ = 0.2 θ = 2
D 1 = 0.002299 D 1 = 0.002299 D 1 = 0.002299 D 1 = 0.002299
D 2 = 0.001945 D 2 = 0.002175 D 2 = 0.002066 D 2 = 0.001613
R = 0.001945 R = 0.002175 R = 0.002066 R = 0.001613
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sharma, D.; Argyros, I.K.; Parhi, S.K.; Sunanda, S.K. Local Convergence and Dynamical Analysis of a Third and Fourth Order Class of Equation Solvers. Fractal Fract. 2021, 5, 27. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5020027

AMA Style

Sharma D, Argyros IK, Parhi SK, Sunanda SK. Local Convergence and Dynamical Analysis of a Third and Fourth Order Class of Equation Solvers. Fractal and Fractional. 2021; 5(2):27. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5020027

Chicago/Turabian Style

Sharma, Debasis, Ioannis K. Argyros, Sanjaya Kumar Parhi, and Shanta Kumari Sunanda. 2021. "Local Convergence and Dynamical Analysis of a Third and Fourth Order Class of Equation Solvers" Fractal and Fractional 5, no. 2: 27. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5020027

Article Metrics

Back to TopTop