In this paper, we construct a new iterative algorithm and show that the newly introduced iterative algorithm converges faster than a number of existing iterative algorithms. We present a numerical example followed by graphs to validate our claim. We prove strong and weak convergence results for approximating fixed points of Suzuki generalized nonexpansive mappings. Again we reconfirm our results by example and table. Further, we utilize our proposed algorithm to solve delay differential equation.
1.
Introduction
Numerical reckoning fixed points for nonlinear operators is nowadays an active research problem of nonlinear analysis owing to its applications to: variational inequalities, equilibrium problems, computer simulation, image encoding and much more. Mann [19], Ishikawa [17] and Halpern [12] are the three basic iterative algorithms to approximate fixed points of nonexpansive mappings. Following these study, several authors constructed numerous algorithms to approximate the fixed points of different classes of nonlinear mappings mainly Noor iteration [20], Agarwal et al. iteration [4], SP iteration [21], Normal-S iteration [23], Abbas and Nazir iteration [1], Thakur et al. iterations [28,29], Karakaya et al. iteration [18] and many others.
In 2008, Suzuki [26] introduced a new generalization of nonexpansive mappings and called the defining condition as Condition (C) which is also referred as Suzuki generalized nonexpansive mappings. A mapping T:K→K defined on a nonempty subset K of a Banach space E is said to satisfy the Condition (C) if
for all x and y∈K.
Suzuki proved that the mappings satisfying the Condition (C) is weaker than nonexpansive and also obtained few results regarding the existence of fixed points for such mappings. In 2011, Phuengrattana [22] used Ishikawa iteration to obtain some convergence results for mappings satisfying Condition (C) in uniformly convex Banach spaces. In the last few years, many authors have studied this particular class of mappings in various domain and have obtained many convergence results (e.g. [2,3,9,10,28,30,31,35]).
Recently, Ullah and Arshad [31] introduced a new algorithm namely M-iteration algorithm as follows:
where {αn} is a sequence in (0,1). They proved some convergence results for Suzuki generalized nonexpansive mappings and showed that M-iteration converges faster than Picard-S [11] and S-iteration [4].
To achieve a better rate of convergence, we introduce a new iterative algorithm for approximating fixed points of Suzuki generalized nonexpansive mappings as follows:
where {αn} is a sequence in (0,1).
The aim of this paper is to prove that newly defined iterative algorithm (1.2) converges faster than algorithm (1.1) for contractive-like mappings. Also, we prove some convergence results involving algorithm (1.2) for Suzuki generalized nonexpansive mappings. Further, we provide a numerical example to show that our iteration (1.2) converges faster than a number of existing iterative algorithms in respect of Suzuki generalized nonexpansive mappings. In the last section, we use our algorithm to find solution of a delay differential equation.
2.
Preliminaries
For making our paper self contained, we collect some basic definitions and needed results.
Definition 2.1. A Banach space E is said to be uniformly convex if for each ϵ∈(0,2] there is a δ>0 such that for x,y∈E with ‖x‖≤1, ‖y‖≤1 and ‖x−y‖>ϵ, we have
Definition 2.2. A Banach space E is said to satisfy the Opial's condition if for any sequence {xn} in E which converges weakly to x∈E i.e. xn⇀x implies that
for all y∈E with y≠x.
Examples of Banach spaces satisfying this condition are Hilbert spaces and all lp spaces (1<p<∞). On the other hand, Lp[0,2π] with 1<p≠2 fail to satisfy Opial's condition.
A mapping T:K→E is demiclosed at y∈E if for each sequence {xn} in K and each x∈E, xn⇀x and Txn→y imply that x∈K and Tx=y.
Let K be a nonempty closed convex subset of a Banach E, and let {xn} be a bounded sequence in E. For x∈E write:
The asymptotic radius of {xn} relative to K is given by
and the asymptotic center A(K,{xn}) of {xn} is defined as:
It is known that, in a uniformly convex Banach space, A(K,{xn}) consists of exactly one point.
The following definitions about the rate of convergence were given by Berinde [7].
Definition 2.3. Let {an} and {bn} be two real sequences converging to a and b respectively. Then, {an} converges faster then {bn} if limn→∞‖an−a‖‖bn−b‖=0.
Definition 2.4. Let {un} and {vn} be two fixed point iteration processes converging to the same fixed point p. If {an} and {bn} are two sequences of positive numbers converging to zero such that ‖un−p‖≤an and ‖vn−p‖≤bn for all n≥1, then we say that {un} converges faster than {vn} to p if {an} converges faster then {bn}.
The following lemma due to Schu [24] is very useful in our subsequent discussion.
Lemma 2.1. Let E be a uniformly convex Banach space and {tn} be any sequence such that 0<p≤tn≤q<1 for some p,r∈R and for all n≥1. Let {xn} and {yn} be any two sequences of E such that lim supn→∞‖xn‖≤r, lim supn→∞‖yn‖≤r and lim supn→∞‖tnxn+(1−tn)yn‖=r for some r≥0. Then, limn→∞‖xn−yn‖=0.
Now, we list few lemmas involving Suzuki generalized nonexpansive mappings.
Lemma 2.2. ([26]) Let K be a nonempty subset of a Banach space E and T:K→K be any mapping. Then,
(ⅰ) If T is nonexpansive then T is Suzuki generalized nonexpansive mapping.
(ⅱ) If T is Suzuki generalized nonexpanisve mapping such that F(T)≠∅, then T is a quasi-nonexpansive mapping.
(ⅲ) If T is a Suzuki generalized nonexpansive mapping, then ‖x−Ty‖≤3‖x−Tx‖+‖x−y‖ for all x and y∈K.
Lemma 2.3. ([27]) Let T be a Suzuki generalized nonexpansive mapping defined on a subset K of a Banach space E with the Opial property. If a sequence {xn} converges weakly to z and limn→∞‖Txn−xn‖=0, then I−T is demiclosed at zero.
Lemma 2.4. ([26]) If T is a Suzuki generalized nonexpansive mapping defined on a compact convex subset K of a uniformly convex Banach space E then, T has a fixed point.
In 1972, Zamfirescu [34] introduced Zamfirescu mappings which serves as an important generalization for Banach contraction principle [5]. In 2004, Berinde [6] gave a more general class of mappings known as quasi-contractive mappings. Following this, Imoru and Olantiwo [16] gave the following definition:
Definition 2.5. A mapping T:K→K is known as contractive-like mapping if there exists a strictly increasing and continuous function φ:[0,∞)→[0,∞) with φ(0)=0 and a constant δ∈[0,1) such that for all x,y∈K, we have
Clearly, the class of contractive-like mappings is wider than the class of quasi-contractive mappings.
3.
Rate of convergence
In this section, first we show that our algorithm (1.2) converges faster than the M-iteration (1.1) for contractive-like mappings.
Theorem 3.1. Let T be a contractive-like mapping defined on a nonempty closed convex subset K of a Banach space E with F(T)≠∅. If {xn} is a sequence defined by (1.2), then {xn} converges faster than the iterative algorithm (1.1).
Proof. From (1.1), for any p∈F(T), we have
and
As, {αn} is a sequence in (0,1), we can find a constant α∈R such that αn≤α<1 for all n∈N. So,
Now, from (1.2) we get
and
So,
Let bn=δ3n(1−(1−δ)α)n‖x1−p‖ and an=δ2n(1−(1−δ)α)n‖d1−p‖, then
Hence, {xn} converges faster than {dn}.
Now, we present a example of a contractive-like mapping which is not a contraction.
Example 1: Let E=R and K=[0,6]. Let T:K→K be a mapping defined as
Proof: Clearly x=0 is the fixed point of T. First, we prove that T is a contractive-like mapping but not a contraction. Since T is not continuous at x=3∈[0,6], so T is not a contraction. We show that T is a contractive-like mapping. For this, define φ:[0,∞)→[0,∞) as φ(x)=x8. Then, φ is a strictly increasing as well as continuous function. Also, φ(0)=0.
We need to show that
for all x,y∈[0,6] and δ is a constant in [0,1).
Before going ahead, let us note the following. When x∈[0,3), then
and
Similarly, when x∈[3,6], then
and
Consider the following cases:
Case A: Let x,y∈[0,3), then using (3.1) we get
So (A) is satisfied with δ=15.
Case B: Let x∈[0,3) and y∈[3,6] then using (3.1) we get
So (A) is satisfied with δ=15.
Case C: Let x∈[3,6] and y∈[0,3) then using (3.2) we get
So (A) is satisfied with δ=15.
Case D: Let x,y∈[3,6] then using (3.2) we get
So (A) is satisfied with δ=15.
Consequently, (A) is satisfied for δ=15 and φ(x)=x8 in all the possible cases. Thus, T is a contractive-like mapping.
Now, using T, we show that our iterative algorithm (1.2) has a better rate of convergence. Set αn=βn=γn=nn+1 for each n∈N. Then, we get the following Table 1, Table 2, Figure 1 and Figure 2 with the initial value 4.5.
Clearly, our algorithm (1.2) converges at a faster rate for contractive-like mappings.
4.
Convergence results
First, we prove few lemmas which will be useful in obtaining convergence results.
Lemma 4.1. Let T be a Suzuki generalized nonexpansive mapping defined on a nonempty closed convex subset K of a Banach space E with F(T)≠∅. Let {xn} be the iterative sequence defined by the algorithm (1.2). Then, limn→∞‖xn−p‖ exists for all p∈F(T).
Proof. Let p∈F(T) and z∈K. Since T is a Suzuki generalized nonexpansive mapping, 12‖p−Tp‖=0≤‖p−z‖ implies that ‖Tp−Tz‖≤‖p−z‖.
Now we have,
and
Using (4.1) and (4.2), we get
Thus, {‖xn−p‖} is bounded and decreasing sequence of reals and hence limn→∞‖xn−p‖ exists.
Lemma 4.2. Let T be a Suzuki generalized nonexpansive mapping defined on a nonempty closed convex subset K of a uniformly convex Banach space E. Let {xn} be the iterative sequence defined by the algorithm (1.2). Then, F(T)≠∅ if and only if {xn} is bounded and limn→∞‖Txn−xn‖=0.
Proof. Suppose F(T)≠∅ and let p∈F(T). Then, by Lemma 4.1, limn→∞‖xn−p‖ exists. Let limn→∞‖xn−p‖=c.
From Eqs (4.1) and (4.2), we have
and
Now,
and
So,
which along with Eq (4.3) implies
Since T is a Suzuki generalized nonexpansive mapping, we get
From Eq (4.4), we obtain
Consider,
Using Lemma 2.3, from Eqs (4.4), (4.5) and (4.6), we get
Now, consider
which on using Eq (4.7) gives
Since,
this together with Eqs (4.7) and (4.8) yields that
Now, using Eqs (4.8) and (4.9), we have
Hence,
Conversely, suppose that {xn} is bounded and limn→∞‖xn−Txn‖=0. Let p∈A(K,{xn}), we have
This implies that Tp∈A(K,{xn}). Since E is uniformly convex, A(K,{xn}) is singleton, therefore we get Tp=p.
Theorem 4.1. Let T be a Suzuki generalized nonexpansive mapping defined on a nonempty closed convex subset K of a Banach space E which satisfies the Opial's condition with F(T)≠∅. If {xn} is the iterative sequence defined by the iterative algorithm (1.2), then {xn} converges weakly to a fixed point of T.
Proof. Let p∈F(T). Then, from Lemma 4.1 limn→∞‖xn−p‖ exists. In order to show the weak convergence of the algorithm (1.2) to a fixed point of T, we will prove that {xn} has a unique weak subsequential limit in F(T). For this, let {xnj} and {xnk} be two subsequences of {xn} which converges weakly to u and v respectively. By Lemma 4.1, we have limn→∞‖Txn−xn‖=0 and using the Lemma 2.3, we have I−T is demiclosed at zero. So u,v∈F(T).
Next, we show the uniqueness. Since u,v∈F(T), so limn→∞‖xn−u‖ and limn→∞‖xn−v‖ exists. Let u≠v. Then, by Opial's condition, we obtain
which is a contradiction, so u=v. Thus, {xn} converges weakly to a fixed point of T.
Next, we establish some strong convergence results for iterative algorithm (1.2).
Theorem 4.2. Let T be a Suzuki generalized nonexpansive mapping defined on a nonempty compact convex subset K of a uniformly convex Banach space E. If {xn} is the iterative sequence defined by the iterative algorithm (1.2), then {xn} converges strongly to a fixed point of T.
Proof. Using Lemma 2.4, we get F(T)≠∅. So, by Lemma 4.2, we have limn→∞‖Txn−xn‖=0. Since K is compact, there exists a subsequence {xnk} of {xn} such that {xnk} converges strongly to p for some p∈K. From Lemma 2.2(ⅲ), we have
for all n≥1. Letting k→∞, we get that {xnk} converges to Tp. This implies that Tp=p, i.e., p∈F(T). Further, limn→∞‖xn−p‖ exists by Lemma 4.1. So, p is the strong limit of the sequence {xn}.
A mapping T:K→K is said to satisfy the Condition (A) ([25]) if there exists a nondecreasing function f:[0,∞)→[0,∞) with f(0)=0 and f(r)>0 for all r∈(0,∞) such that ‖x−Tx‖≥f(d(x,F(T))) for all x∈K, where d(x,F(T))=inf{‖x−p‖:p∈F(T)}.
Theorem 4.3. Let T be a Suzuki generalized nonexpansive mapping defined on a nonempty closed convex subset K of a uniformly convex Banach space E such that F(T)≠∅ and {xn} be the sequence defined by (1.2). If T satisfies Condition (A), then {xn} converges strongly to a fixed point of T.
Proof. By Lemma 4.1, limn→∞‖xn−p‖ exists and ‖xn+1−p‖≤‖xn−p‖ for all p∈F(T).
We get
which yields
This shows that the sequence {d(xn,F(T))} is decreasing and bounded below, so limn→∞d(xn,F(T)) exists.
Let limn→∞‖xn−p‖=r for some r≥0. If r=0 then the result follows. Assume r>0. Also, by Lemma 4.2 we have limn→∞‖xn−Txn‖=0.
It follows from Condition (A) that
so that limn→∞f(d(xn,F(T)))=0.
Since f is a non decreasing function satisfying f(0)=0 and f(r)>0 for all r∈(0,∞), therefore limn→∞d(xn,F(T))=0. So, we have a subsequence {xnk} of {xn} and a sequence {yk}⊂F(T) such that
for all k∈N. Using (4.4), we obtain
Therefore,
This implies that {yk} is a cauchy sequence in F(T). Since F(T) is closed, so {yk} converges to a point p∈F(T). Then, {xnk} converges strongly to p. Since limn→∞‖xn−p‖ exists, we get xn→p∈F(T). This completes the proof.
5.
Numerical example
In this section, first we will construct an example of a Suzuki generalized nonexpansive mapping which is not a nonexpansive mapping. Then, using that example, we will show that our iteration scheme (1.2) has a better speed of convergence than number of existing iteration schemes.
Example 2: Define a mapping T:[0,1]→[0,1] by
First we show that T is not a nonexpansive map. For this, take x=8100 and y=112. Then,
and
Clearly, ‖Tx−Ty‖>‖x−y‖ which proves that T is not a nonexpansive mapping.
Now, we show that T satisfies the condition K. For this, consider the following cases:
Case-Ⅰ: Let x∈[0,112), then 12‖x−Tx‖=12|2x−1|=12(1−2x). For 12‖x−Tx‖≤‖x−y‖, we must have 12(1−2x)≤‖x−y‖, i.e., 12(1−2x)≤|x−y|. Here note that the case y<x is not possible. So, we are left with only one case when y>x, which gives 12(1−2x)≤y−x, which yields y≥12. So, y∈[12,1]. Now, we have x∈[0,112) and y∈[12,1]. So,
and
Hence,
Case-Ⅱ: Let x∈[112,1], then 12‖x−Tx‖=12|x−x+1112|=11−11x24. For 12‖x−Tx‖≤‖x−y‖, we must have 11−11x24≤‖x−y‖, i.e., 11−11x24≤|x−y|. Here we have two possibilities.
A: When x<y, we get 11−11x24≤y−x, i.e., y≥11+13x24. So, y∈[145288,1]⊂[112,1], which gives ‖Tx−Ty‖=112‖x−y‖≤‖x−y‖. Hence,
B: When x>y, then 11−11x24≤x−y, i.e. y≤35x−1124 which gives y∈[0,1]. Also, 24y+1135≤x which yields x∈[1135,1]. Here, for x∈[1135,1] and y∈[112,1] Case IIA can be used. So, we only need to verify when x∈[1135,1] and y∈[0,112). For this,
and
So, ‖Tx−Ty‖≤‖x−y‖. Thus, mapping T satisfies the Condition (C) for all the possible cases.
Now, using above example, we will show that iteration algorithm (1.2) converges faster than Thakur New, M and M∗ iteration. Let αn=βn=nn+10 for all n∈N and x1=0.02, then we get the following Table 3 of iteration values and Figure 3.
It is evident from above table and graph that our algorithm (1.2) converges at a better speed than the above mentioned schemes.
6.
Application
In this section, we show that our iterative algorithm can be used to find a solution of a delay differential equation.
Many physical problems arising in various fields can be easily modeled with the help of ordinary differential equations. Later, it was recognized that a phenomena may have a delayed effect in a differential equation, leading to the development of concept of delay differential equations. Following this, numerous methods have been obtained to solve various kinds of delay differential equations (e.g. [13,14,15,32,33]).
In this paper, we consider the following delay differential equation
with initial condition
Now, we will show that the sequence generated by our iteration scheme (1.2) converges strongly to the solution of (6.1).
It is well known that (C([a,b]),||.||∞) is a Banach space where C([a,b]) denotes the space of all continuous real valued functions on a closed interval [a,b] and ||.||∞ is a Chebyshev norm ||x−y||∞=maxt∈[a,b]|x(t)−y(t)|.
Assume that the following conditions are satisfied
(A1) t0,b∈R,τ>0;
(A2) f∈C([t0,b]×R2,R);
(A3) ψ∈C([t0−τ,b],R);
(A4) there exists Lf>0 such that
|f(t,u1,u2)−f(t,v1,v2)|≤Lf2∑i=1|ui−vi|,∀ui,vi∈R,i=1,2,t∈[t0,b];
(A5) 2Lf(b−t0)<1.
We notice that the solution of (6.1)-(6.2) if it exists is of the following form
Here, x∈C([t0−τ,b],R)∩C1([t0,b],R).
Coman et al. [8] established the following results.
Theorem 6.1. Assume that conditions (A1)−(A5) are satisfied. Then Problem (6.1)−(6.2) has a unique solution, say p∈C([t0−τ,b],R)∩C1([t0,b],R) and
Now, we prove the following result using our iterative process (1.2).
Theorem 6.2. Suppose that conditions (A1)−(A5) are satisfied. Then the problem (6.1)−(6.2) has a unique solution say p∈C([t0−τ,b],R)∩C1([t0,b],R) and sequence generated by the algorithm (1.2) converges to p.
Proof. Let {xn} be a iterative sequence generated by (1.2) for the following operator:
where αn∈(0,1) for all n∈N such that ∑∞n=0αn=∞. Denote by p the fixed point of T. We will show that xn→p as n→∞.
For t∈[t0−τ,t0], it is easy to see that xn→p as n→∞.
For t∈[t0,b], we have
Using (6.3), (6.4), (6.5) and (6.6) we get
On using assumption (A_5), we have
Therefore, inductively we get
Since \alpha_n \in [0, 1], for all n \in \mathbb{N} , assumption (A_5) yields
Using the fact that e^{-x} \geq 1 - x for all x \in [0, 1], we have
which gives \lim\limits_{n\to\infty}\|x_{n} - p\|_ {\infty} = 0.
From the above theorem, we can say that our method will definitely converge to the unique solution of (6.1) which is a main advantage over the other methods available for the same.
7.
Conclusion
In this study a new fixed iteration process (1.2) has been obtained which is utilized to approximate fixed point of Suzuki generalized nonexpansive mappings. Further, We show that our iteration process (1.2) converges faster than the recent M-iteration process (1.1) for contractive-like mappings. It must be noted here that Ullah and Arshad [31] did not give the rate of convergence of their process analytically. They claimed just by an example. However, we not only give the proof analytically but also validate with an example. Further, we performed convergence analysis and a non trivial example has been given to illustrate the convergence behaviour. In the last section, we applied our iteration process to find the solution of delay differential equation.
Acknowledgements
We wish to pay our sincere thanks to learned referees for pointing out many omission and motivating us to study deeply for numerical aspects.
Conflict of interest
The authors declare that they have no competing interests.