7.9
CiteScore
 
3.6
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Correspondence
Corrigendum
Editorial
Full Length Article
Invited review
Letter to the Editor
Original Article
Research Article
Retraction notice
REVIEW
Review Article
SHORT COMMUNICATION
Short review
7.2
CiteScore
3.7
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Correspondence
Corrigendum
Editorial
Full Length Article
Invited review
Letter to the Editor
Original Article
Research Article
Retraction notice
REVIEW
Review Article
SHORT COMMUNICATION
Short review
View/Download PDF

Translate this page into:

Research Article
2025
:37;
2872025
doi:
10.25259/JKSUS_287_2025

A novel AA* method for exploring the interplay between fractals, polynomiographs, and fractional calculus

Department of Sciences and Humanities, National University of Computer and Emerging Sciences, 852-B Milaad St, Block B Faisal Town, Lahore, 54770, Pakistan
Department of Mathematics, Division of Science and Technology University of Education, Township, Lahore, 54770, Pakistan
Department of Solids and Structures, School of Engineering, The University of Manchester, Oxford Road, Manchester, M13 9PL, Pakistan
Department of Mathematics, College of Science, King Saud University, P.O. Box 2455, Riyadh, 11451, Saudi Arabia

*Corresponding author E-mail address: muhammad.azeem@riphah.edu.pk (M Azeem)

Licence
This is an open-access article distributed under the terms of the Creative Commons Attribution-Non Commercial-Share Alike 4.0 License, which allows others to remix, transform, and build upon the work non-commercially, as long as the author is credited and the new creations are licensed under the identical terms.

Abstract

We introduce an innovative numerical approach for estimating fixed points of symmetric generalized nonexpansive mappings within uniformly convex Banach spaces. This approach demonstrates a notably quicker convergence rate relative to existing iterative methods. Its effectiveness is confirmed through comprehensive numerical tests, comprising CPU time comparisons and visual representations created in MATLAB, including polynomiographs. Additionally, using the proposed method, we show that solutions to a type of delay fractional differential equations exist and are unique. Our investigation is performed within the scope of Garcia-Falset type mappings, a wider category that includes nonexpansive and Suzuki mappings, thus broadening the relevance of our method.

Keywords

28A80
34K37
47H09
47H10
Banach space
Convergence mathematical classification
Fixed point
Fractals
Garcia-falset mapping

1. Introduction

The development of efficient algorithms is a central focus in research on iterative methods for approximating fixed-point solutions. To assess the performance of an iterative scheme, the notion of rate of convergence is widely utilized. However, since the use of technology is increasing in the approximation process, one cannot deny the importance of the time efficiency of the algorithms. Therefore, in addition to the rate of convergence, it is of significance to know the efficiency of the algorithm via the amount of time it takes to converge to the solution. Furthermore, fractals, such as polynomiographs, can be used to visualize the convergence zones of certain polynomials with complex coefficients. As a result, an empirical comparison of iteration procedures is useful in determining efficiency. Picard’s iterative technique is recognized as one of the pioneering algorithms in fixed-point estimation. However, this approach is not effective in approximating the fixed points of nonexpansive mappings (Banach 1922; Picard, 1890). As a result, several scholars proposed various iterative strategies for determining fixed points of nonexpansive mappings after the failure or sluggish convergence of previous schemes. Many recent studies have contributed to the creation of precise and stable iterative methods as well as fractional techniques for addressing nonlinear models and systems of equations (Qureshi et al., 2025; Qureshi et al., 2024a; Chang et al., 2025; Naseem et al., 2025; Chang et al., 2024; Argyros et al., 2024; Qureshi et al., 2024b).

Let χ:AA be a mapping, where A is a subset of a Banach space B. Mann introduced a one-step iterative scheme. To initiate the iteration process, we first choose an initial point u1 A, then generate a sequence {un} using the following iterative scheme:

(1.1)
un+1 =(1αn)un+αnTun,n+,

where {αn}(0,1) is a sequence of real numbers. One disadvantage of this method was that it failed to approximate the fixed points of pseudocontractive mappings and had a poor rate of convergence for nonexpansive mappings. To address this issue, Ishikawa, 1974 devised the following two-step iterative technique. To begin the iterative process, we select an initial point u1 A and generate sequences { bn} and {un} using the following iterative scheme.

(1.2)
  u1 A bn= 1βn un+βnχun,   un+1 = 1αn un+αnχbn,      n+,

where {αn}, {βn}(0,1) are sequences of real numbers.

Numerous scholars, including Hacıoğlu,2021; Thakur et al., 2016; Agarwal et al., 2007; Abbas and Nazir, 2014; Noor, 2000, have introduced various iterative methods. Recently, Beg et al., 2022 have proposed a AA-iterative method (1.3) that converges quickly when compared to previously existing iterative techniques in the literature. The algorithm is shown below:

(1.3)
u1 A, dn= 1rn un+rnχun, cn=χ 1qn dn+qnχdn , bn=χ 1sn χdn+snχcn , un+1 =χbn,                                        n+,

where u1 A is an initial point and {sn},{qn} and {rn}(0,1) are sequences of real numbers. The sequence un is generated recursively as follows: for each n, we first compute dn from the current iterate un , then use dn to compute cn . Next, bn is determined using both dn and cn . Finally, un+1 is generated from bn .​

(Browder, 1965) demonstrated the presence of fixed points of nonexpansive mappings under some conditions. (Garcia-Falset et al., 2011) presented a generalization of nonexpansive mappings known as mappings satisfying property E. These mappings generate a larger class than Suzuki mappings, which satisfy condition C. A mapping χ:AB is said to be equipped with property (Eλ), also known as Garcia-Falset mapping, if it satisfies the inequality

(1.4)
aχıλaχa+aı,a,ıA,

where λ1 . We note that χ satisfied condition (E) on A, whenever λ1 such that χ fulfills condition (Eλ). Suzuki, 2007 and Garcia-Falset et al., 2011 have demonstrated convergence results for such mappings under the effects of a uniformly convex Banach space.

Pandey et al., 2019 proved that the class of mappings defined in (1.4) contains numerous classes of generalized nonexpansive mappings. Usurelu et al., 2022 worked for the approximation of common fixed points for mappings satisfying Garcia-Falset and (α,β)-generalized hybrid property (Bejenaru and Ciobanescu, 2022; Chalarux and Chaichana, 2021). Fractional calculus, which involves integrals and derivatives of arbitrary order, has seen a rise in interest over the last thirty years. It has emerged as a captivating field that has attracted the attention of researchers from various fields (Guariglia, 2021; Guariglia, 2018; Hacıoğlu, 2021; Hezenci, 2023; Hezenci and Budak, 2023a; Hezenci and Budak, 2023b; Jiao et al., 2016; Kaur and Goyal, 2019; Kumar and Suat Erturk, 2023). It provides powerful techniques for addressing differential and integral equations, as well as applications in fields such as mathematical biology, chemical processes, and engineering challenges. It has become a powerful tool for modeling complex systems with memory and hereditary properties, wave propagation, fractal-fractional systems, tumor-immune interactions, and nonlinear wave equations (Atangana and Alkahtani, 2021; El-Sayed and Kılıçman, 2021; Singh et al., 2017; Jarad et al., 2024; Inc and Ali, 2013)).

As previously noted, various iterative techniques have been proposed in recent years to achieve a faster rate of convergence for nonexpansive mappings and their extensions. However, there is still a gap in the literature concerning the convergence analysis of many broader classes of mappings, including Garcia-Falset mappings.

In view of the preceding discussion, we alter the AAiterative scheme and demonstrate convergence results by using the AA* iterative strategy for Garcia-Falset mappings in a uniformly convex Banach space. We shall use the following approaches to demonstrate the novel scheme’s efficiency:

  • 1.

    Analytical rate of convergence;

  • 2.

    Comparative Numerical experiments: Convergence behavior, CPU time analysis, and polynomiographs.

In contrast to earlier reported results, we take another strategy to establish weak convergence by adopting the Fréchet differential norm rather than Opial’s property. To further illustrate the utility of our findings, we apply the proposed scheme to delay differential equations. The AA* iterative scheme is defined as:

(1.5)
u1 A, dn=χ 1rn un+rnχun , cn=χ 1qn dn+qnχdn ,                             bn=χ 1sn χdn+snχcn ,      un+1 =χbn,                                        n+,

where u1 A is an initial point and {sn},{qn} and {rn}(0,1) are sequences of real numbers. The sequence un is generated recursively as follows: for each n, we first compute dn from the current iterate un , then use dn to compute cn . Next, bn is determined using both dn and cn . Finally, un+1 is generated from bn .

We have used the following symbols and notations throughout the paper:

+ : Set of positive integers

: Set of real numbers

B: Uniformly convex Banach space

A: Non-empty closed convex subset of B

F(χ): Set of all fixed points in A

2. Preliminaries

This section caters to certain preliminary details and terminology that will be necessary in the sequel. Throughout, we denote the uniformly convex Banach space by B.

Suppose A is a closed convex subset of B denoted as and χ:AB be a mapping. Denote F(χ)={aA:a=χa}. For any a,ıA, a mapping χ is said to be:

  • 1.

    a contraction if for any β(0,1), we have

    χaχıβaı,

  • 2.

    a nonexpansive mapping if

    χaχıaı,

  • 3.

    quasi-nonexpansive if

    χalal,lFχ

  • 4.

    Suzuki generalized nonexpansive mapping if

    1 2 aχaaı then χaχıaı.

Definition 2.1 (Ali et al., 2019)

A mapping χ is said to be demiclosed if for every sequence un Awith unu (weakly) and χun u (strongly), it follows that χu=d.

Definition 2.2 (Opial, 1967) If every sequence u n weakly converges to uB and the following holds:

liminf nunu< liminf nund,

for all dB and ud, then B satisfies Opial’s condition.

Definition 2.3 (Goebel and Kirk, 1990) The existence of the following limit indicates that a B is uniformly smooth.

(2.1)
lim t0 a+tbat,a,bSB,

where SB represents the unit sphere of B.

If for every aB the limit (2.1) is uniformly obtained for bSB, then the norm of B is said to be Fréchet differentiable. Then, we have

(2.2)
1 2 y2 + a,K(y) 1 2 y+a2 1 2 y2 + a,K(y) +f a IIa,yB,

where f: 0, assumed to be an increasing function with lim t0 ft t=0 and Ky stands for the Fréchet derivative of the functional 1 2 .2 at yB.

Definition 2.4 (Senter and Dotson, 1974) The function χ:AA satisfies Condition (I) if there is a non-decreasing function s: 0, 0, such that

d u,χu s d u,Fχ ,uA,where following holds:

  • 1.

    s 0 =0 ;

  • 2.

    s x >0,x>0 ;

  • 3.

    d u,Fχ =inf u,l :lFχ.

The concepts of asymptotic radius and center are defined as follows, (Edelstein, 1974).

Definition 2.5 The mapping  ρ ., un :B+ is defined by:

ρ d, un = limsup n und,

ρ d, un , for each dB, is referred to as the asymptotic radius of un at d. Asymptotic radius of un relative to a subset AB is defined as:

ρ A, un =infρ d, un ,

with dA. Asymptotic center of un corresponds to A is defined as:

C A, un =ρ A, un =ρ(d, un ,

where dA.

In B, the asymptotic center consists of exactly one element (Clarkson, 1936).

Definition 2.6 Let an and bn be two iterative methods that converge to same fixed point l for a mapping χ and {kn} n=0 and {ln} n=0 be two positive sequences approaching to k and l. The error estimates are given as follows:

anlkn,bnlln,n, and R= lim n knk lnl .

  • 1.

    If R=0 , then an iterative algorithm an converges rapidly as compared to bn.

  • 2.

    If 0<R<, then convergence rate of both algorithms is identical.

Proposition 2.1 (Garcia-Falset et al., 2011) If a mapping χ:AA fulfills (1.4) then χ is quasi-nonexpansive provided it ensures the existence of fixed point.

Lemma 2.1 (Schu, 1991) Let pn and vn be two sequences meeting limsup n pnβ, limsup n vnβ and limsup n 1sn pn+sn vn=β for any β0, and 0<sn<1 for all n+, then lim n pn vn=0.

Lemma 2.2 (Soltuz and Otrocol, 2007) Consider a sequence βn on + that satisfies the inequality: βn+1 βn 1αn . If αn 0,1 and n=0 αn=, then lim n βn=0.

3. Results on Convergence

In this section, we prove convergence results by applying the AA* -iterative algorithm on generalized nonexpansive operators satisfying (1.4). Throughout this section, χ:AA satisfies (1.4).

Lemma 3.1. If Fχ and un is a sequence generated by (1.5) then lim n unl exists where lFχ.

Proof. Let lFχ. Employing proposition (2.1), equations (1.4) and (1.5), we have

(3.1)
   dnl = χ 1rn un+rnχ un l   1rn un+rnχ un l 1rn unl+rnχ un l  1rn unl+rn λχll+unl   1rn unl+rnunl   = unl. 

Consider,

cnl = χ 1qn dn+qnχ dn l            1qn dn+qnχ dn l                   1qn dnl+qnχ dn l    dnl. 

Utilizing (3.1), we get

(3.2)
cnlunl.

Also,

bnl = χ 1sn χdn+snχ cn l      1sn χdn+snχ cn l            1sn χdnl+snχcnl  1sn dnl+sncnl.      

Using (3.1) and (3.2)

(3.3)
bnlunl.

Consequently,

un+1 l=χbnl bnl.

Equation (3.3) implies

un+1 lunl,

Thus, unl is monotone decreasing and bounded for every lFχ, which guarantees the existence of its limit.

Lemma 3.2. If un is a sequence generated by (1.5), then Fχ iff un is bounded and lim n unχun=0.

Proof. From Lemma 3.1 we already know that lim n unl exists and un is bounded. Let

(3.4)
lim nunl=l.

From (3.1), (3.2), (3.3) and (3.4), we have

(3.5)
limsup ndnl limsup nunll, limsup ncnl limsup nunll, limsup nbnl limsup nunll.

Since χ is quasi-nonexpansive mapping, so

(3.6)
χunl unl limsup nχunl l.

So,

un+1 l=χbnlbnl.

Applying liminf n , we get

(3.7)
l liminf nbnl.

From (3.5) and (3.7)

lim nbnl=l.

In the same way,

l liminf ncnl lim ncnl=l l liminf ndnl lim ndnl=l.

Furthermore,

l lim ndnl = lim nχ 1rn un+rnχun l l.

Accordingly,

(3.8)
lim ndnl=l.

Using Lemma (3.1), (3.4), (3.6), and (3.8), yields

lim nunχun=0.

The converse holds if un is bounded, unχun=0 and lC A, un . Since χ is a Garcia-Falset mapping, as a result

R χl, un = limsup nunχl limsup nunl =R l, un =R A, un ,

which implies χlC A, un .

As the asymptotic center is unique, C A, un has one element. Moreover, as the underlying Banach space is uniformly convex, χl=l.

Theorem 3.1. Let B be endowed with Opial’s condition. If un is a sequence approximated by (1.5), then un weakly converges to a fixed point of χ.

Proof. Consider lFχ. Applying Lemma 3.1, lim n unl exists. Suppose u ni and u nj are subsequences converge weakly to p 1 and p 2 , respectively. According to Lemma 3.2, lim n unχun=0, Iχ at zero yields.

Iχ p 1 =0, i.e., χ p 1 = p 1 . Similarly, χ p 2 = p 2 . This implies p 1 and p 2 are fixed points.To show uniqueness, when p 1 p 2 . By Opial’s condition, we get

lim nun p 1 = lim n u ni p 1 < lim niu ni p 2 = lim nun p 2 = lim nju nj p 2 < lim nju nj p 1 = lim nun p 1 ,

This leads to a contradiction. Hence p 1 = p 2 . This implies, un converges weakly to a fixed point of χ.

It is known that Hilbert spaces, all finite-dimensional Banach spaces, and the sequence spaces lp for 1<p< possess the Opial property. However, the function spaces Lp[0,2π] fail to satisfy the Opial property when p2 . Instead of relying on Opial’s property, we assume that the Fréchet differential norm exists (Goebel and Kirk, 1990).

Theorem 3.2. Assume the B is equipped with the Fréchet differential norm and let limn tun+ 1t p1 p2 exists for each p1 ,p2 Fχ. In addition, consider demiclosedness of Iχ at zero. If un is a sequence generated by (1.5) then un converges weakly to a fixed point of χ with Fχ.

Proof. Suppose l1 and l2 are weak limits of u ni and u nj respectively. From Theorem 3.2 we have, lim n unχun=0. Since Iχ is demiclosed at zero, it follows that l1 ,l2 Fχ. By replacing p and q with p 1 p 2 and t un p 1 with respectively in (2.2), then

1 2 p 1 p 2 2 +t un p 1 ,K( p 1 p 2 ) 1 2 tun+ 1t p 1 p 2 2 1 2 p 1 p 2 2 +t un p 1 ,K( p 1 p 2 ) +f tun p 1 .

From the given constraint, we have

1 2 p1 p2 2 +t limsup n unp1 ,K(p1 p2 ) 1 2 lim ntun+(1t)p1 p2 2 1 2 p1 p2 2 +tliminf n unp1 ,K(p1 p2 ) +Ot.

Therefore,

limsup n unp1 ,K(p1 p2 ) liminf n unp1 ,K(p1 p2 ) + O(t)t.

Employing t 0+, we ensure the existence of lim n unp1 ,K(p1 p2 ) . Let l1 p1 ,K(p1 p2 ) =e and l2 p1 ,K(p1 p2 ) =e. Hence, l1 l2 ,K(p1 p2 ) =0 p1 ,p2   F(v). This yields

l1 l2 2 = l1 l2 ,K(l1 l2 ) =0,

which cannot be true unless l1 =l2 . Consequently, un converges weakly to a fixed point of χ. ◻

Theorem 3.3. If un is a sequence generated by AA*iterative scheme then un converges to a fixed point of χ iff  liminfn d un,Fχ =0 or limsupn d un,Fχ =0 , with d un,Fχ =inf unl and lFχ.

Proof. When lim n un=l, where lFχ, then the prof is complete.

For the converse, consider liminf n d un,Fχ =0. By Lemma 3.1, lim n unl exists for each lFχ . According to our assumption,

liminf nd un,Fχ =0.

To show un is a Cauchy sequence in A. Since liminf n d un,Fχ =0, for η>0, there exists n0 +, such that

d un,Fχ <η2 , infd un,l <η2 ,where lFχ

Consequently, we can find lFχ such that

u n0 l<η2 .

Choose m,nn0 ,

um+n unum+n l+unl 2u n0 l <η.

Hence, the sequence un is Cauchy sequence in A. As A is closed, there exists a point sA, such that lim n un=s. In addition, liminf n d un,Fχ =0. It follows, sFχ.

Theorem 3.4. Let AB be a nonempty, compact and convex. If un is a sequence generated by (1.5) and Fχ then un converges strongly to lFχ.

Proof. Assume Fχ. Lemma 3.2 implies that B lim n unχun=0. Compactness of A guarantees the existence of a subsequence u nk of un for which u nk l with lA. Garcia-Falset property implies,

u nk χpνu nk χu nk +u nk l.

Letting k, we obtain u nk χl=0. So, we get u nk χl. Hence, lFχ.It follow from Lemma 3.1, lim n unl exists. We obtain un l. ◻

Theorem 3.5. If un is a sequence generated by AA*iterative scheme then un converges strongly to lFχ, where χ satisfies the condition I.

Proof. From Lemma 3.2, we get

lim nunχun=0.

According to the Condition (I), we deduce

lim nunχun lim nh d un,Fχ 0,

This implies,

lim nh d un,Fχ =0.

Considering the properties of h, we yield

lim nd un,Fχ =0.

Since all conditions of Theorem 3.3 hold, we conclude un converges strongly to lFχ. ◻

4. Analytical Rate of Convergence

In this section, we examine the analytical convergence comparison of the AAiterative scheme and the AA* iterative scheme.

Lemma 4.1. Let χ:AA be a contraction mapping with α 0,1 a constant of contraction and un be the sequence generated by AAiterative and AA*iterative scheme. If 1sn<sn,1qn<qn and 1rn<rn, where sequences sn , qn and rn 0,1 , n, then the AA*iterative scheme converges faster than the AAiterative scheme.

Proof. For the AA* iterative scheme, we have

dnl=χ 1rn un+rnχun l α 1rn+αrn unl.

Also,

cnl=χ 1qn dn+qnχdn l α2 1qn+αqn 1rn +αrn unl.

Moreover,

bnl=χ 1sn χdn+snχcn l α3 1rn+αrn 1+sn ααqn+α2 qn sn unl.

Thus,

un+1 l=χbnl α4 1rn+αrn 1+sn ααqn+α2 qn sn unl.

Since 1sn<sn,1qn<qn and 1rn<rn,n, we have 1α rn< 1 2 1α , 1α qn< 1 2 1α and 1α sn< 1 2 1α .

So,

α4 1rn+αrn 1+sn ααqn+α2 qn sn <α4 1 1 2 1α 1+αsn(1 1 2 1α sn .

Set en=α 4n 1 1 2 1α n 1+αsn(1 1 2 1α sn nu1 l.

For the AA-iterative scheme,

dnl= 1rn un+rnχunl 1rn+αrn unl.

Also, we have

cnl = χ((1qn)dn+qnχdn)l (ααqn+α2 qn)(1rn+αrn) unl

Moreover,

bnl=χ 1sn χdn+snχcn l α2 1rn+αrn 1+sn ααqn+α2 qn sn unl.

Thus,

un+1 l=χbnl α3 1rn+αrn 1+sn ααqn+α2 qn sn unl.

Since 1sn<sn,1qn<qn and 1rn<rn,n, we have 1α rn< 1 2 1α , 1α qn< 1 2 1α and 1α sn< 1 2 1α .

So,

α3 1rn+αrn 1+sn ααqn+α2 qn sn <α3 1 1 2 1α 1+snα 1 1 2 1α sn

Set hn=α 3n 1 1 2 1α n {1+αsn 1 1 2 1α sn}nu1 l. Then,

lim n en hn = lim n α 4n 1 1 2 1α n {1+αsn 1 1 2 1α sn}nu1 l α 3n 1 1 2 1α n {1+αsn 1 1 2 1α sn}nu1 l =0.

This implies that the AA* iterative scheme converges faster than the AAiterative scheme.

5. Discussion

5.1 Graphical analysis of convergence and time efficiency

To illustrate the convergence results, we consider the following numerical example.

Example 5.1 Consider A= 0, equipped with the usual norm, and let a mapping χ: 0, 0, be defined as:

χa= 2a 3 ,a3 0,a 0,3

We can see χ is not satisfying Condition C and also χ is Garcia-Falset mapping satifying (1.4).

Suppose a=3.5 and ı=1.5

1 2 aχa = 1 2 3.5 2 3.5 3 =0.58.

Also, we have

|aı|=2.

Hence,

1 2 |aχa|<|aı|,

but |χaχ1|>|a1|. So χ is not satisfying Condition (C).

Let ν=1 , then the following cases are considered:

Case I

If a>3 and ı 0,3 ,

then

|aχ1|ν|aχa|+|a1| |a0||a1|+ a 2a 3 |a||a1|+ a3 ,

which holds.

Case II

If a,ı>3

then

|aχ1|ν|aχa|+|a1| a 2ι 3 a 2a 3 +|a1| a 2ι 3 |a1|+ a3 ,

which is true.

Case III

If a 0,3 and ı>3

then

|aχ1|ν|aχa|+|a1| a 1 2 |a1|+|a0| 2ι 3 |a1|+|a|,

which is true.

Case IV

If a,ı 0,3

then

|aχ1|ν|aχa|+|a1| |a0||a1|+|a|,

which is true.

This implies that χ satisfies (1.4).

A numerical simulation is presented to visualize the convergence behaviors using Garcia Falset mapping, the AA* iterative scheme, and well-known iterative schemes such as AA, Mann, Ishikawa, Thakur, etc. Considering the values of the relevant parameters, sn=n n2 +1 , qn=n n2 +4n+2 , and rn=n n+4 1 3 with initial point a1 =5 . and sn= 1 10n+1 , qn= 2n+1 1 4 9n+10 , and rn= n 2n+1 with initial point a1 =10 as illustrated in Figs. 1(a,b) respectively. It indicates that Mann, Ishikawa, Agarwal, and Thakur processes need more iterations as compared to AA* and AA to converge the fixed point i.e. l=0 indicating that these processes are either not convergent at all or converge very slowly. While AA* and AA take fewer iterations to converge to the fixed point of χ.

(a-b) Graphs of an−l for different parameters.
Fig. 1.
(a-b) Graphs of anl for different parameters.

Moreover, in Fig. 2, we have shown the graphs of errors an+1 l versus the number of iterations for error analysis of AA* and AA iterative scheme, where l is fixed point. We have set stopping criteria an+1 l <5× 10 4 . It is clear that AA* meets stopping criteria before AA. This implies that AA* converges fast comparative to AAiterative scheme.

(a-b) Error analysis of AA* and AA iterative scheme.
Fig. 2.
(a-b) Error analysis of AA* and AA iterative scheme.

In Table 1, we have measured the time for every simulation and also noted the impact of initial points on CPU time for different iterative schemes by choosing different parameters. It is obvious that AA* has a computation time less than other algorithms. width=0.7

Table 1. CPU time comparison of AA* - scheme with other schemes.
Initial point AA(sec) AA* (sec) Mann (sec) Agarwal (sec) Thakur (sec)
3.1 0.0781 0.0469 0.1875 0.1563 0.0781
3.5 0.0938 0.0469 0.2344 0.1406 0.0781
5 0.0625 0.0312 0.2188 0.1563 0.0938
10 0.0938 0.0781 0.25 0.1875 0.1406
40 0.0781 0.0625 1.1406 0.1094 0.1094

6. Polynomiographs

Evaluating the stability and convergence speed is a critical component of assessing any iterative scheme. This evaluation can be performed using visual methods known as polynomiography (Naseem et al., 2025; Kalantari, 2009). The overall process for creating a polynomiograph is outlined in Algorithm 1. By applying different coloring techniques in the final step of Algorithm 1, we can illustrate diverse characteristics of the iterative approach.

In this section, we offer a visual comparison of various iteration processes with AA* -scheme for fixed point approximation of Newton’s iteration operator using polynomiographs of some complex polynomials. Polynomiographs can help show the convergence zones of polynomials with complex values. Polynomiographs twist symmetrically when all parameters are real, but asymmetrically when some are imaginary. The renowned Newton method for locating roots of a complex polynomial P is given by the formula:

τ1 τn+1 = τn P(τn) P (τn) ; n=1,2,3,

The above equation may be rewritten as a fixed point iteration scheme

τn+1 =R(τn),

if we consider R(τ)=τ P(τ) P(τ) .

Note: τ denotes the current root approximation computed by the iterative method R.

Assume that this iteration converges to a fixed point τ of R, then τ=R(τ)=τ P(τ) P(τ) . Therefore, τ is a root of P as P(τn) P(τn) =0 if and only if P(τ)=0 .

We analyze Mann, Ishikawa, Agarwal, Thakur, AA, and AA* iteration processes using Newton’s operator to locate the roots of the following complex polynomials:

  • 1.

    P(τ)=τ3 1 ;

  • 2.

    P(τ)=τ5 1 .

Polynomiographs are developed for two different iteration parameter values: rn=qn=sn=0.01 and rn=qn=sn=0.5 whereas other parameters are set as: A= 2,2 × 2,2 , N=15 , ϵ=0.001 . Algorithm 1 provides pseudocode for generating polynomiographs. The algorithm uses iteration coloring to color points. This coloring method assigns a color to each starting point based on the number of iterations executed, representing the speed of convergence. Therefore, this form of polynomiograph demonstrates the speed of convergence for root finding. The color variation has been given in Fig. 3.

Colorbar.
Fig. 3.
Colorbar.

Algorithm 1: Generation of a polynomiograph

Input: Polynomial Pz with deg P2 , iteration process O, area A, maximum number of iterations N, accuracy ϵ, color map for colors.

Output: Polynomiograph for the complex-valued polynomial P within the area A.

Result: Polynomiograph for P within A.

for z0 A do

n = 0

While n < N and |P( τn )| > ϵ do

τn+1 =O zn,P

n = n + 1

end while

Map n to a color from the color map and

end for

We have the following observations:

  • 1.

    For parameters rn=sn=qn=0.01 : In Figs. 4 and 5, observe that Mann and Ishikawa’s iterations did not reach any of the roots of the equation, resulting in a uniform red color. This corresponds to the maximum of 15 iterations. The speed of convergence varies throughout iterations. Visual analysis shows that the proposed AA* achieves the fastest convergence rate, followed by AA, Thakur, and Agarwal.

  • 2.

    For parameters rn=sn=qn=0.5 : In Figs. 6 and 7, the Mann and Ishikawa iterations yield the slowest convergence rate since the polynomiograph shows an abundance of reddish colors, indicating a high number of iterations. The polynomiograph for the AA iteration shows faster convergence than the Agarwal and Thakur iteration. We may see that the polynomiographs for the AA* scheme exhibit darker blue colors compared to other iteration techniques, indicating fewer iterations. Thus, the AA* iteration is the quickest of the analyzed iterations.

(a-f) Polynomial: P(τ)=τ3 −1 , Parameters rn=qn=sn=0.01 .
Fig. 4.
(a-f) Polynomial: P(τ)=τ3 1 , Parameters rn=qn=sn=0.01 .
(a-f) Polynomial: P(τ)=τ5 −1 , Parameters rn=qn=sn=0.01 .
Fig. 5.
(a-f) Polynomial: P(τ)=τ5 1 , Parameters rn=qn=sn=0.01 .
(a-f) Polynomial: P(τ)=τ3 −1 , Parameters: rn=qn=sn=0.5 .
Fig. 6.
(a-f) Polynomial: P(τ)=τ3 1 , Parameters: rn=qn=sn=0.5 .
(a-f) Polynomial: P(τ)=τ5 −1 , Parameters: rn=qn=sn=0.5.
Fig. 7.
(a-f) Polynomial: P(τ)=τ5 1 , Parameters: rn=qn=sn=0.5.

7. Application to the solution of Delay Differential Equation

Fractional differential equations (FDEs) are becoming a very important and popular topic. It is possible to specify the so-called FDEs by extending the normal integer order derivative to any order. FDEs have been widely employed to explain numerous physical processes, such as the flow of seepage in media with pores and in fluid dynamic traffic models. Because fractional calculus has so many applications in the applied sciences and engineering, many authors have focused their attention on FDEs as a significant area of fractional calculus research. Several methods have been proposed to solve FDEs numerically. We cite the references (Ansari et al., 2023; Kaushik et al., 2023; Kumar and Suat Erturk, 2023; Ünal et al., 2023; Xu et al., 2019; Yu and Li, 2021) therein for interesting theory results and scientific applications of FDEs.

Caputo’s fractional derivatives, which he introduced in 1967, are a novel method of fractional differentiation that is defined as:

c Dsag(s)= 1 Γ(an) bs g n (r) (sr) na1 dr, (n1<a<n).

In this section, we employ the AA* iterative approach (1.5) with a 0,1 to estimate the solution of a DFDEs (Cong and Tuan, 2017). Let g: 0,N be a function such that f Cn 0,N . Consider ωC kh,k : Rn a continuous mapping and h>0 any constant. sConsider DFDE

(7.1)
c Daus=g s,us,u sh , s k,N

with initial conditions

(7.2)
us=ωs,s kr,k ,

where u Rn,g: k,N × Rn× Rn Rn is continuous, r>0 and N>0. Assume that the following assumptions hold.

  • 1.

    g fulfills the Lipschitz condition w.r.t second and third variables, and there is a constant Lf>0 such that

    g s,u,t g s,u^,t^ Lf uu^+tt^

    sR+ and u,u^,t,t^ Rn.

  • 2.

    There is a constant αL that depends upon L such that αL>2L which means 2L αL <1.

u* C kh,N : Rn C1 k,N : Rn is known as solution of IVP (7.2) if it satisfies (7.1) and (7.2). We know that approximating the solution of (7.1) and (7.2) is equivalent to solving the following integral equation:

us=ωs+ 1 Γa ks sr β1 g r,ur,u rh dr,s k,N ,

with u(s)=ω(s),s[kh,N]. We define norm . αL on C([k ‒ h, k] : R n) as

ω αL = supωs Ea αL sa , for any ωC kh,k : Rn ,

where Ea:RR denotes Mittag–Leffler function, given by:

Ea s= n=0 sn Γ βn+1 ,sR.

Observe that C kh,k : Rn ,. αL ) is a Banach space. Now, we present the result for approximation of the solution of (7.1) and (7.2) using (1.5).

Theorem 7.1 Suppose ω and g are functions as mentioned above and assumptions (i) and (ii) are fulfilled, then IVP (7.1) and (7.2) possesses a unique solution u* C kh,N : Rn C1 k,N : Rn and the sequence un given by (1.5) converges to u* .

Proof. For existence of unique solution u* , see (Cong and Tuan, 2017). Let un is defined by (1.5). Consider a map χ on C kh,N : Rn C1 k,N : Rn by

χus= ωk+ 1 Γa ks sr β1 g r,ur,u rh dr,s k,N , ωs, s kr,k .

Let pn= 1rn un+rnχun, using (1.5) and conditions (i) and (ii), we get

pnu* = 1rn un+rnχunu* 1rn unu* +rnχunu* .

We apply supremum over kh,N in the above inequality on both sides, we have

(7.3)
sup s kh,N pnu* 1rn sup s kh,N unu* +rn sup s kh,N χunχu* = 1rn sup s kh,N unu* +rn sup s kh,N 1 Γa ks sr β1 g r,ur,u rh dr 1 Γa ks sr β1 g r,u* r,u* rh dr

Dividing both sides by Ea αL sa , we obtain

sup s kh,N pnu* Ea αL sa 1rn sup s kh,N unu* Ea αL sa +rn Lf Γa ks sr β1 dr sup s kh,N un ru* r Ea αL sa + sup s kh,N un rh u* rh Ea αL sa . pnu* αL 1rn unu* αL + rn Γa ks sr β1 dr Lf un ru* r αL +un rh u* rh αL . = 1rn unu* αL +rn 2 Lf unu* αL 1 Γa ks sr β1 dr = 1rn unu* αL + rn 2 Lf Ea αL sa unu* αL 1 Γa

ks sr β1 Ea αL sa dr = 1rn unu* αL + rn 2 Lf Ea αL sa unu* αL .c Ia cDa Ea αL sa αL = 1rn unu* αL + rn 2 Lf Ea αL sa . Ea αL sa αL unu* αL = 1rn unu* αL + rn 2 Lf αL unu* αL .

Since 2 Lf αL <1, we have

(7.4)
pnu* αL unu* αL .

Therefore,

dnu* =χpnχu* 1 Γa ks sr β1 g r,pr,p rh drωk 1 Γa ks sr β1 g r,u* r,u* rh dr.

Applying supremum over kh,N in above inequality on both sides, we get

sup s kh,N dnu* = sup s kh,N ωk + 1 Γa ks sr β1 g r,pr,p rh drωk 1 Γa ks sr β1 g r,u* r,u* rh dr. Lf Γa ks sr β1 dr sup s kh,N pn ru* r+ sup s kh,N pn rh u* rh .

Dividing both sides by Ea αL sa yields,

dnu* αL 1 Γa ks sr β1 dr Lf pn ru* r αL +pn rh u* rh αL . = 2 Lf pnu* αL 1 Γa ks sr β1 dr = 2 Lf Ea αL sa pnu* αL 1 Γa ks sr β1 Ea αL sa dr = 2 Lf Ea αL sa pnu* αL .c Ia c Da Ea αL sa αL = 2 Lf Ea αL sa . Ea αL sa αL pnu* αL = 2 Lf αL pnu* αL .

As 2 Lf αL <1, this implies

dnu* αL pnu* αL .

From (7.4), we have

dnu* αL unu* αL .

Similarly, we can show that

cnu* αL unu* αL , bnu* αL unu* αL , un+1 u* αL unu* αL .

By setting, unu* αL =mn, we get

mn+1 mn n.

Hence, mn is a monotonic decreasing sequence in R+. Moreover, mn is bounded sequence, we have

lim nmn=inf mn =0.

Hence,

lim nun=inf un =0.

Example 7.2

Consider the DFDE:

(7.5)
C Dt 0.9 u s =u s +0.5 u s0.4 ,s 0,2 ,

with initial condition:

u s =ϕ s =1, s 0.4,0 .

C Ds 0.9 represents the Caputo fractional derivative of order a=0.9 , and τ=0.4 is the constant delay.

Define g s,u,y =u+0.5y. Then: g satisfies Lipschitz condition with constants L1 =1 , L2 =0.5 , g s,1,1 =1+0.5=0.5M=0.5 .

This example satisfies assumptions (i) and (ii).

The DFDE (7.5) is equivalent to the Volterra-type integral equation:

u s =ϕ 0 + 1 Γ 0.9 0t sr 0.1 g r,u r ,u t0.4 dr.

To approximate the Caputo derivative numerically, we use the classical L1 scheme (Xiao et al., 2024). We apply the AA*-iterative method where χ denotes the integral operator from the L1 approximation, and sn , qn , rn 0,1 are control sequences. The numerical solution obtained using the AA* -iterative scheme has been shown in Fig. 8. We set sn=0.5,qn=0.4,rn=0.3 and use h=0.01 .

Solution of DFDE using -iterative scheme.
Fig. 8.
Solution of DFDE using -iterative scheme.

8. Conclusions

The main purpose of developing new schemes in approximation of fixed points via iteration is to provide efficient algorithms, that is, algorithms that not only have a faster rate of convergence but are also time efficient. This study has investigated the fixed points of Garcia-Falset mappings using a novel algorithm AA* -iterative approach that outperforms previously existing schemes. Some convergence results were provided, as well as an analytical examination of the rate of convergence demonstrated that the AA* technique surpasses existing algorithms. We performed some numerical experiments with graphical visualizations and CPU time tables to show the applicability and efficiency of our suggested method in comparison to certain existing techniques in the literature. Various parameters and starting points were chosen to discuss the scheme’s speed and time efficiency. We also see the faster convergence to the root of complex polynomials through symmetric polynomiographs. An application to the solution of a delay fractional differential equation was discussed to illustrate the usefulness of our results.

For future work, we recommend exploring the generation of fractals on the Mandelbrot and Julia sets using the AA* iterative method. This approach could enhance the visualization of complex dynamics and potentially reveal deeper insights into the behavior of these sets.

Acknowledgement

This project was supported by the Deanship of Scientific Research at King Saud University under the research funding program (ORF-2025-158).

CRediT authorship contribution statement

All authors contributed equally to this manuscript.

Declaration of competing interest

Patient's consent not required as there are no patients in this study.

Data availability

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Declaration of Generative AI and AI-assisted technologies in the writing process

The authors confirm that there was no use of Artificial Intelligence (AI)-Assisted Technology for assisting in the writing or editing of the manuscript and no images were manipulated using AI.

References

  1. , . Some new faster iteration process applied to constrained minimization and feasibility problems. 2014
    [Google Scholar]
  2. , , . Iterative construction of fixed points of nearly asymptotically nonexpansive mappings. J Nonlinear Convex Anal. 2007;8:61.
    [Google Scholar]
  3. , , . Approximation of fixed points for Suzuki's generalized non-expansive mappings. Mathematics. 2019;7:522.
    [Google Scholar]
  4. , , , , , . On new updated concept for delay differential equations with piecewise Caputo fractional-order derivative. Waves in Random and Complex Media 2023:1-20. https://doi.org/10.1080/17455030.2023.2187241
    [Google Scholar]
  5. , , , , , . Local and semi-local convergence and dynamic analysis of a time-efficient nonlinear technique. Appl Numerical Mathematics. 2024;201:446-464. https://doi.org/10.1016/j.apnum.2024.03.022
    [Google Scholar]
  6. , , , , , . Fractional optimal control of COVID-19 pandemic model with generalized Mittag-Leffler function. Adv Differ Equ. 2021;2021:387. https://doi.org/10.1186/s13662-021-03546-y
    [Google Scholar]
  7. . Sur les opérations dans les ensembles abstraits et leur application aux équations intégrales. Fund Math. 1922;3:133-181.
    [Google Scholar]
  8. , , . Convergence of iterative algorithm for generalized nonexpansive mappings with an application. Mathematics. 2022;10:4375.
    [Google Scholar]
  9. , . Common fixed points of operators with property (E) in CAT(0) spaces. Mathematics. 2022;10:433. https://doi.org/10.3390/math10030433
    [Google Scholar]
  10. . Nonexpansive nonlinear operators in a Banach space. Proc Natl Acad Sci USA. 1965;54:1041-1044.
    [Google Scholar]
  11. , . Approximation of fixed points for Garcia-Falset mappings in a uniformly convex Banach space. Results Nonlinear Anal. 2021;4:200-206.
    [Google Scholar]
  12. , , , , . A modified two-step optimal iterative method for solving nonlinear models in one and higher dimensions. Math Comput Simul. 2025;229:448-467. https://doi.org/10.1016/j.matcom.2024.09.021
    [Google Scholar]
  13. , , , , . A modified fractional Newton's solver. Axioms. 2024;13:689.
    [Google Scholar]
  14. . Uniformly convex spaces. Trans Amer Math Soc. 1936;40:396-414. https://doi.org/10.1090/s0002-9947-1936-1501880-4
    [Google Scholar]
  15. , . Existence, uniqueness, and exponential boundedness of global solutions to delay fractional differential equations. Mediterr J Math. 2017;14 https://doi.org/10.1007/s00009-017-0997-4
    [Google Scholar]
  16. . Fixed point theorems in uniformly convex Banach spaces. Proc Amer Math Soc. 1974;44:369-374. https://doi.org/10.1090/s0002-9939-1974-0358451-4
    [Google Scholar]
  17. , . Analysis of the fractional tumor-immune-vitamins model with Mittag–Leffler kernel. Math Methods Appl Sci. 2021;44:7652-7665.
    [Google Scholar]
  18. , , . Fixed point theory for a class of generalized nonexpansive mappings. J Math Anal Appl. 2011;375:185-195. https://doi.org/10.1016/j.jmaa.2010.08.069
    [Google Scholar]
  19. , . Topics in metric fixed point theory. Cambridge University Press; . No. 28
  20. . Harmonic sierpinski gasket and applications. Entropy (Basel). 2018;20:714. https://doi.org/10.3390/e20090714
    [Google Scholar]
  21. . Fractional calculus, zeta functions, and Shannon entropy. Open Math. 2021;19:87-100. http://doi.org/10.1515/math-2021-0010
    [Google Scholar]
  22. . A comparative study on iterative algorithms of almost contractions in the context of convergence, stability and data dependency. Comp Appl Math. 2021;40 https://doi.org/10.1007/s40314-021-01671-8
    [Google Scholar]
  23. . A note on fractional simpson type inequalities for twice differentiable functions. Math Slovaca. 2023;73:675-686. https://doi.org/10.1515/ms-2023-0049
    [Google Scholar]
  24. , . An extensive study on parameterized inequalities for conformable fractional integrals. Anal Math Phys. 2023;13 https://doi.org/10.1007/s13324-023-00846-2
    [Google Scholar]
  25. , . Novel results on trapezoid-type inequalities for conformable fractional integrals. Turk J Math. 2023;47:425-438. https://doi.org/10.55730/1300-0098.3371
    [Google Scholar]
  26. , . Improved (G′/G)-Expansion method for the space and time fractional foam drainage and KdV equation. Abstr Appl Anal. 2013;2013 Article ID 414353. https://doi.org/10.1155/2013/414353
    [Google Scholar]
  27. . Fixed points by a new iteration method. Proc Amer Math Soc. 1974;44:147-150. https://doi.org/10.1090/s0002-9939-1974-0336469-5
    [Google Scholar]
  28. , , . On solutions of fractal fractional differential equations. Int J Numer Model. 2024;37:e3283.
    [Google Scholar]
  29. , , . Box dimension of Weyl fractional integral of continuous functions with bounded variation. Anal Theory Appl. 2016;32:174-180.
    [Google Scholar]
  30. . Polynomial root-finding and polynomiography. Singapore: World Scientific; .
  31. , . Riemann zeta fractional derivative: Functional equation and link with primes. Adv Differ Equ. 2019;261
    [Google Scholar]
  32. , , , . Existence of solutions by fixed point theorem of general delay fractional differential equation with p-Laplacian operator. AIMS Math. 2023;8:10160-10176. https://doi.org/10.3934/math.2023514
    [Google Scholar]
  33. , . The analysis of a time delay fractional COVID‐19 model via Caputo type fractional derivative. Math Methods in App Sciences. 2023;46:7618-7631. https://doi.org/10.1002/mma.6935
    [Google Scholar]
  34. , , , , , , . Memory based approaches to one-dimensional nonlinear models. Acta Appl Math. 2025;195 https://doi.org/10.1007/s10440-024-00703-9
    [Google Scholar]
  35. . New approximation schemes for general variational inequalities. J Math Anal Appl. 2000;251:217-229. https://doi.org/10.1006/jmaa.2000.7042
    [Google Scholar]
  36. . Weak convergence of the sequence of successive approximations for nonexpansive mappings. Bull Amer Math Soc. 1967;73:591-597. https://doi.org/10.1090/s0002-9904-1967-11761-0
    [Google Scholar]
  37. , , , . Approximating fixed points of a general class of nonexpansive mappings in Banach spaces with applications. Results Math. 2019;74:1-24. https://doi.org/10.1007/s00025-018-0930-6
    [Google Scholar]
  38. . Memoire sur la theorie des equations aux derivees partielles et la methode des approximations successives. J Math Pures Appl. 1890;6:145-210.
    [Google Scholar]
  39. , , , , . A highly accurate family of stable and convergent numerical solvers based on Daftardar–Gejji and Jafari decomposition technique for systems of nonlinear equations. MethodsX. 2024;13:102865. https://doi.org/10.1016/j.mex.2024.102865
    [Google Scholar]
  40. , , , , , . A new optimal numerical root-solver for solving systems of nonlinear equations using local, semi-local, and stability analysis. Axioms. 2024;13:341. https://doi.org/10.3390/axioms13060341
    [Google Scholar]
  41. , , , , , . Use of fractional calculus to avoid divergence in Newton-like solver for solving one-dimensional nonlinear polynomial-based models. Commun Nonlinear Sci Numer Simul. 2025;143:108631. https://doi.org/10.1016/j.cnsns.2025.108631
    [Google Scholar]
  42. . Weak and strong convergence to fixed points of asymptotically nonexpansive mappings. Bull Austral Math Soc. 1991;43:153-159. https://doi.org/10.1017/s0004972700028884
    [Google Scholar]
  43. , . Approximating fixed points of nonexpansive mappings. Proc Am Math Soc. 1974;44:375. https://doi.org/10.2307/2040440
    [Google Scholar]
  44. , , . Numerical solution of fractional telegraph differential equations by theta-method. Eur Phys J Spec Top. 2017;226:3547-3562.
    [Google Scholar]
  45. , . Classical results via Mann-Ishikawa iteration. J Numer Anal Approx Theory. 2007;36:193-197. https://doi.org/10.33993/jnaat362-868
    [Google Scholar]
  46. . Fixed point theorems and convergence theorems for some generalized nonexpansive mappings. J Math Anal Appl. 2008;340:1088-1095. https://doi.org/10.1016/j.jmaa.2007.09.023
    [Google Scholar]
  47. , , . A new iteration scheme for approximating fixed points of nonexpansive mappings. Filomat. 2016;30:2711-2720. https://doi.org/10.2298/fil1610711t
    [Google Scholar]
  48. , , . Conformable fractional Newton-type inequalities with respect to differentiable convex functions. J Inequal Appl. 2023;2023 https://doi.org/10.1186/s13660-023-02996-0
    [Google Scholar]
  49. , , . Algorithm for two generalized nonexpansive mappings in uniformly convex spaces. Mathematics. 2022;10:318. https://doi.org/10.3390/math10030318
    [Google Scholar]
  50. , , , . L1 Schemes for time-fractional differential equations: A brief survey and new development. 2024
    [Google Scholar]
  51. . Multiscale modeling of porous media. In: Heat and mass transfer in drying of porous media Heat and mass transfer in drying of porous media. Boca Raton: Taylor & Francis, a CRC title, part of the Taylor & Francis imprint, a member of the Taylor & Francis Group, the academic division of T&F Informa, plc, [2020] | Series: Advances in drying science and technology: CRC Press; p. :55-82. https://doi.org/10.1201/9781351019224-3
    [Google Scholar]
  52. , . Some fractal characters of porous media. Fractals. 2021;9:365-372.
    [Google Scholar]
Show Sections