Three-Step Iterative Algorithm for the Extended Cayley–Yosida Inclusion Problem in 2-Uniformly Smooth Banach Spaces: Convergence and Stability Analysis

Ali, Imran; Wang, Yuanheng; Ahmad, Rais

doi:10.3390/math12131977

Open AccessArticle

Three-Step Iterative Algorithm for the Extended Cayley–Yosida Inclusion Problem in 2-Uniformly Smooth Banach Spaces: Convergence and Stability Analysis

by

Imran Ali

^1,†

,

Yuanheng Wang

^2,3,*,†

and

Rais Ahmad

^4,†

¹

Department of Engineering Mathematics, College of Engineering, Koneru Lakshmaiah Education Foundation, Vaddeswaram 522302, Andhra Pradesh, India

²

School of Mathematical Sciences, Zhejiang Normal University, Jinhua 321004, China

³

Mathematics Department of Humanities College, Zhejiang Guangsha Vocational and Technical University of Construction, Jinhua 321004, China

⁴

Department of Mathematics, Aligarh Muslim University, Aligarh 202002, Uttar Pradesh, India

^*

Author to whom correspondence should be addressed.

^†

These authors contributed equally to this work.

Mathematics 2024, 12(13), 1977; https://doi.org/10.3390/math12131977

Submission received: 8 June 2024 / Revised: 22 June 2024 / Accepted: 24 June 2024 / Published: 26 June 2024

(This article belongs to the Special Issue Advanced Research in Functional Analysis and Operator Theory)

Download

Browse Figures

Versions Notes

Abstract

:

In this article, we investigate and study an extended Cayley–Yosida inclusion problem. We show that our problem is equivalent to a fixed-point equation. Based on the fixed-point equation, we develop a three-step iterative algorithm to solve our problem. Finally, we illustrate the convergence of the proposed algorithm with an example, computational table, and convergence graph by using MATLAB 2018b.

Keywords:

three-step iterative algorithm; Cayley–Yosida inclusion; resolvent operator; convergence and stability; smooth Banach space

MSC:

47H09; 49J40

1. Introduction

Variational inclusion represents a valuable and significant generalization of the variational inequality, playing a crucial role in various fields such as engineering science, economics, transportation equilibrium, and optimization and control. For more detail, we direct readers to [1,2,3] and the references therein. The multi-valued inclusion problem is to find

\hat{a} \in \hat{H}

such that

0 \in M (\hat{a}),

(1)

where

M : \hat{H} \to 2^{\hat{H}}

is the multi-valued mapping. This problem serves as a versatile model for a broad spectrum of challenges in both pure and applied sciences. For detail, we refer to [4,5,6,7] and the references therein.

Resolvent and Yosida approximation operators have important roles in convex analysis, partial differential equations, and other analytical problems, as illustrated in [8,9,10,11]. In 1994, Hassouni and Moudafi [12] explored variational inclusions, a type of variational inequality with single-valued mappings, using the resolvent operator approach for maximal monotone mappings. The initial method employed to address problem (1) was the forward–backward splitting algorithm. Subsequently, various iterative methods, including the Douglas–Rachford splitting method [13], the Peaceman–Rachford splitting method [14], and others, have been developed to tackle mixed inclusion problems (MIPs).

Tseng [15] introduced the modified forward–backward approach in 2000, demonstrated weak convergence. In 2018, Gibali and Thong [16] obtained a modified version of Tseng’s splitting method, showcasing its strong convergence. In 2019, Ahmad et al. [17] developed an XOR operation-based, three-step iterative technique for generalized mixed ordered quasi-variational inclusion and independently demonstrated the existence, convergence, and stability results. In 1989, Glowinski et al. [18] demonstrated that three-step iterative methods outperformed and gave better numerical results than two- and one-step iterative schemes.

Due to the Cayley operators’ [19] widespread applicability in fundamental and allied sciences, encompassing fields such as financial modeling, computer programming, economics, and engineering, among others, we investigate a generalized Cayley operator. This operator incorporates a generalized resolvent operator initially introduced by Fang and Huang [20].

Motivated by the above-mentioned work, our aim is to introduce and study an extended Cayley–Yosida inclusion problem. We formulate the proposed problem into a fixed-point equation by using a generalized Yosida approximation operator and a generalized resolvent operator. Based on the fixed-point equation, we develop a three-step iterative algorithm. We provide proof of stability analysis results and strong convergence for the sequence generated by the three-step iterative algorithm. Additionally, convergence graphs and a numerical example are provided, implemented through MATLAB programming.

2. Background and Problem Formulation

The majority of splitting methods are founded on the resolvent operator of type

{[I + λ M (f, g)]}^{- 1}

, where

M

is a multi-valued monotone mapping,

λ > 0

constant, and

I

is an identity operator. Thus, the resolvent operator,

R_{I, λ}^{M (f, g)} : \hat{H} \to \hat{H}

, where

\hat{H}

is a Hilbert space, is defined as

R_{I, λ}^{M (f, g)} (\hat{a}) = {[I + λ M (f, g)]}^{- 1} (\hat{a}), f o r a l l \hat{a} \in \hat{H} .

(2)

The Yosida approximation operator,

Y_{I, λ}^{M (f, g)}

, is defined as

Y_{I, λ}^{M (f, g)} (\hat{a}) = \frac{1}{λ} [I - R_{I, λ}^{M (f, g)}] (\hat{a}), f o r a l l \hat{a} \in \hat{H} .

(3)

The Cayley operator,

C_{I, λ}^{M (f, g)}

, is defined as

C_{I, λ}^{M (f, g)} (\hat{a}) = (2 R_{I, λ}^{M (f, g)} - I) (\hat{a}), f o r a l l \hat{a} \in \hat{H} .

(4)

We consider

R_{I, λ}^{M (f, g)}

,

Y_{I, λ}^{M (f, g)}

, and

C_{I, λ}^{M (f, g)}

defined by (2), (3), and (4), respectively, in the context of 2-Uniformly Smooth Banach Spaces in this article.

Throughout this work, we suppose

\hat{E}

to be a real 2-Uniformly Smooth Banach Space (in short, 2-USBS) whose norm is denoted by

∥ \cdot ∥

, where

{\hat{E}}^{*}

is the topological dual of

\hat{E}

,

〈 \cdot, \cdot 〉

is the duality pairing between

\hat{E}

and

{\hat{E}}^{*}

, and

J : \hat{E} \to 2^{{\hat{E}}^{*}}

is the normalized duality mapping, which is defined as

J (\hat{a}) = \{\hat{f} \in {\hat{E}}^{*} : 〈 \hat{a}, \hat{f} 〉 = ∥ \hat{a} ∥^{2}, ∥ \hat{f} ∥ = ∥ \hat{a} ∥ and \hat{a} \in \hat{E}\} .

The function

τ_{\hat{E}} : [0, \infty) \to [0, \infty)

measures the modulus of smoothness of the Banach space

E

and is defined as

τ_{\hat{E}} (t) = sup \{\frac{∥ \hat{a} + \hat{b} ∥ - ∥ \hat{a} - \hat{b} ∥}{2} - 1 : \hat{a}, \hat{b} \in \hat{E}, ∥ \hat{a} ∥ \leq 1, ∥ \hat{b} ∥ \leq t\} .

The Banach space

\hat{E}

is called the following:

(i): Uniformly Smooth if $lim_{t \to 0} \frac{τ_{\hat{E}} (t)}{t} = 0$ ;
(ii): 2-uniformly smooth if there exists a positive real constant $k$ such that $τ_{\hat{E}} (t) \leq k t^{2} .$

The Lemma stated below plays a significant role in the proof of our main result.

Lemma 1

([21,22]). Let

\hat{E}

be a 2-Uniformly Smooth Banach Space and

J : \hat{E} \to 2^{{\hat{E}}^{*}}

be the normalized duality mapping. Then, for any

\hat{a}, \hat{b}, \in \hat{E}

,

(i): ${∥ \hat{a} + \hat{b} ∥}^{2} \leq ∥ \hat{a} ∥^{2} + 2 〈 \hat{b}, j (\hat{a} + \hat{b}) 〉 + k {∥ \hat{b} ∥}^{2}$ , for all $j (\hat{a} + \hat{b}) \in J (\hat{a} + \hat{b})$ ;
(ii): $〈 \hat{a} - \hat{b}, j (\hat{a}) - j (\hat{b}) 〉 \leq 2 d^{2} τ_{\hat{E}} (\frac{4 ∥ \hat{a} - \hat{b} ∥}{d})$ , where $d = \sqrt{\frac{(∥ \hat{a} ∥^{2} + ∥ \hat{b} ∥^{2})}{2}}, j (\hat{a}) \in J (\hat{a})$ and $j (\hat{b}) \in J (\hat{b}) .$

Definition 1

([21,22]). A single-valued mapping

F : \hat{E} \times \hat{E} \to \hat{E}

is said to be mixed

(\land_{1}, \land_{2})

-Lipschitz continuous if

F (\hat{a}, .)

is

\land_{1}

-Lipschitz continuous in the first component and

F (., \hat{a})

is

\land_{2}

-Lipschitz continuous in the second component. That is,

∥ F (\hat{a}, \hat{z}) - F (\hat{b}, \hat{z}) ∥ \leq \land_{1} ∥ \hat{a} - \hat{b} ∥,

and

∥ F (\hat{w}, \hat{a}) - F (\hat{w}, \hat{b}) ∥ \leq \land_{2} ∥ \hat{a} - \hat{b} ∥, \forall \hat{a}, \hat{b}, \hat{w}, \hat{z} \in \hat{E} a n d \land_{1}, \land_{2} > 0 .

For

\land_{1} = 1

and

\land_{2} = 1

,

F

is non-expansive in both arguments, respectively.

Definition 2

([23]). A multi-valued mapping

M : \hat{E} \times \hat{E} \to 2^{\hat{E}}

with regard to single-valued mappings

f, g : \hat{E} \to \hat{E}

is said to be as follows:

(i): α-strongly accretive with regard to $f$ if there is a constant $α > 0$ such that

$〈 \hat{u} - \hat{v}, j (\hat{a} - \hat{b}) 〉 \geq α {∥ \hat{a} - \hat{b} ∥}^{2}, \forall \hat{a}, \hat{b} \in \hat{E}, \forall \hat{u} \in M (f (\hat{a}), \hat{z}), \hat{v} \in M (f (\hat{b}), \hat{z});$
(ii): β-relaxed strongly accretive with regard to $g$ if there is a constant $β > 0$ such that

$〈 \hat{u} - \hat{v}, j (\hat{a} - \hat{b}) 〉 \geq - β {∥ \hat{a} - \hat{b} ∥}^{2}, \forall \hat{a}, \hat{b} \in \hat{E}, \forall \hat{u} \in M (\hat{z}, g (\hat{a})), \hat{v} \in M (\hat{z}, g (\hat{b}));$
(iii): $α β$ -symmetric accretive with regard to $f$ , $g$ if $M (f, .)$ is α-strongly accretive, and $M (., g)$ is β-relaxed accretive with $α \geq β$ if and only if $\hat{a} = \hat{b}$ .

Definition 3

([23]). A multi-valued mapping

M : \hat{E} \times \hat{E} \to 2^{\hat{E}}

is said to be

α β

–I-accretive with regard to

f

and

g

if

M (f, g)

is

α β

-symmetric accretive and for all

λ > 0

,

[I + λ M (f, g)] (\hat{E}) = \hat{E},

where

f, g : \hat{E} \to \hat{E}

are the single-valued mappings and

I : \hat{E} \to \hat{E}

is an identity mapping.

Definition 4.

The generalized resolvent operator

R_{I, λ}^{M (f, g)} : \hat{E} \to \hat{E}

in terms of I and

M

is defined as

R_{I, λ}^{M (f, g)} (\hat{a}) = {[I + λ M (f, g)]}^{- 1} (\hat{a}), for all \hat{a} \in \hat{E},

(5)

where

f, g : \hat{E} \to \hat{E}

are single-valued mappings,

I : \hat{E} \to \hat{E}

is an identity mapping, and

M : \hat{E} \times \hat{E} \to 2^{\hat{E}}

is

α β

–I-accretive with regard to

f

and

g

.

Definition 5

([24]). For

{\hat{a}}_{0} \in \hat{H}

, let

{\hat{a}}_{n + 1} = G (T, {\hat{p}}_{n})

, where

T : \hat{H} \to \hat{H}

is a single-valued mapping and

G : \hat{H} \times \hat{H} \to \hat{H}

is the bimapping. Assume that the fixed-point set

S (T) \neq \emptyset

and

{{\hat{a}}_{n}} \to {\hat{a}}^{*} ({\hat{a}}^{*} \in S (T))

. Let

{{\hat{u}}_{n}} \subset \hat{H}

and

{\hat{ϑ}}_{n} = ∥ {\hat{u}}_{n + 1} - G (T, {\hat{u}}_{n}) ∥ .

Define an iterative scheme that generates a sequence of points

{{\hat{a}}_{n}}

in

\hat{H} .

Assume that

lim_{n \to \infty} {\hat{ϑ}}_{n} = 0

implies

{\hat{u}}_{n} \to {\hat{a}}^{*}

; then, the sequence

{{\hat{a}}_{n}}

is said to be stable with respect to

T .

Lemma 2

([25]). Let

{ξ_{n}} \subseteq R^{+}

be a sequence in

R^{+}

and

{φ_{n}} \subseteq [0, 1]

be a sequence in

[0, 1]

such that

\sum_{n = 0}^{\infty} φ_{n} = \infty .

If there is a natural number m such that

ξ_{n} \leq (1 - φ_{n}) ξ_{n} + φ_{n} η_{n}, \forall n \geq m,

(6)

where

η_{n} \geq 0,

for all

n \geq 0

and

η_{n} \to 0 (n \to 0),

then

lim_{n \to \infty} ξ_{n} = 0 .

Theorem 1.

Let

f, g : \hat{E} \to \hat{E}

be the single-valued mappings and

I : \hat{E} \to \hat{E}

be an identity mapping. Let

M : \hat{E} \times \hat{E} \to 2^{\hat{E}}

be

α β

–I-accretive with regard to

f

and

g

. Then,

R_{I, λ}^{M (f, g)} : \hat{E} \to \hat{E}

is θ-Lipschitz continuous. That is,

∥R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b})∥ \leq θ ∥ \hat{a} - \hat{b} ∥, f o r a l l \hat{a}, \hat{b} \in \hat{E},

where

θ = [\frac{1}{1 + λ (α - β)}]

, and

λ, α,

and

β

are non-negative real numbers with

α \geq β

.

Proof.

Let

\hat{a} \in \hat{E} .

It follows that

R_{I, λ}^{M (f, g)} (\hat{a}) = {[I + λ M (f, g)]}^{- 1} (\hat{a})

and

R_{I, λ}^{M (f, g)} (\hat{b}) = {[I + λ M (f, g)]}^{- 1} (\hat{b}) .

This implies that

\frac{1}{λ} (\hat{a} - R_{I, λ}^{M (f, g)} (\hat{a})) \in M (f (R_{I, λ}^{M (f, g)} (\hat{a}), g (R_{I, λ}^{M (f, g)} (\hat{a})),

\frac{1}{λ} (\hat{b} - R_{I, λ}^{M (f, g)} (\hat{b})) \in M (f (R_{I, λ}^{M (f, g)} (\hat{b}), g (R_{I, λ}^{M (f, g)} (\hat{b})) .

Since

M (., .)

is

α β

-symmetric with respect to

f

and

g

, we have

\begin{matrix} (α - β) {∥R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b})∥}^{2} \\ \leq 〈\frac{1}{λ} (\hat{a} - R_{I, λ}^{M (f, g)} (\hat{a})) - (\hat{b} - R_{I, λ}^{M (f, g)} (\hat{b})), j (R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b}))〉 \\ \leq 〈\frac{1}{λ} (\hat{a} - \hat{b}), j (R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b}))〉 \\ + 〈\frac{1}{λ} (R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b})), j (R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b}))〉 . \end{matrix}

(7)

On the other hand, we calculate

\begin{matrix} ∥ \hat{a} - \hat{b} ∥ ∥R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b})∥ & \geq 〈\frac{1}{λ} (\hat{a} - \hat{b}), j (R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b}))〉 \\ \geq λ (α - β) {∥R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b})∥}^{2} \\ + 〈 R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b}), \\ j (R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b})) 〉 \\ \geq λ (α - β) {∥R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b})∥}^{2} \\ + {∥R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b})∥}^{2} . \end{matrix}

(8)

Hence,

∥R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b})∥ \leq \frac{1}{[1 + λ (α - β)]} ∥ \hat{a} - \hat{b} ∥ .

That is,

R_{I, λ}^{M (f, g)}

is

θ = \frac{1}{[1 + λ (α - β)]}

-Lipschitz continuous. □

We now demonstrate certain properties of the generalized Yosida approximation and the Cayley operators.

Proposition 1

([26]). The Yosida approximation operator defined by (3) is as follows:

(i): $θ_{Y}$ -Lipschitz continuous, where $θ_{Y} = \frac{1 + θ}{λ}$ ;
(ii): $m$ -strongly accretive, where $m = \frac{1 - θ}{λ}$ .

Proof.

The definition of the generalized Yosida approximation operator leads directly to the proof. □

Proposition 2

([26]). The Cayley operator defined by (4) is

θ_{C}

-Lipschitz continuous, where

θ_{C} = 2 θ + 1

.

Proof.

The definition of the generalized Cayley operator leads directly to the proof. □

Problem Formulation

Let

F : \hat{E} \times \hat{E} \to \hat{E}

be the single-valued bimapping,

f, g : \hat{E} \to \hat{E}

be single-valued mappings, and

M : \hat{E} \times \hat{E} \to 2^{\hat{E}}

be a multi-valued mapping. Let

C_{I, λ}^{M (f, g)}

and

Y_{I, λ}^{M (f, g)}

be generalized Cayley and Yosida approximation operators, respectively. We consider the following problem:

Find

\hat{a} \in \hat{E}

such that

0 \in F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) + M (f (\hat{a}), g (\hat{a})) .

(9)

Problem (9) is called an extended Cayley–Yosida inclusion problem.

Additionally, we summarize a few specific examples of our problem below to demonstrate the generality of problem (9).

(i): If $M (f (\hat{a}), g (\hat{a})) = M (\hat{a})$ and $F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) = C_{I, λ}^{M (f, g)} (\hat{a})$ , then problem (9) reduces to the problem of finding $\hat{a} \in \hat{E}$ such that

$0 \in C_{I, λ}^{M} (\hat{a}) + M (\hat{a}) .$

(10)

Problem (10) was investigated and studied by Ahmad et al. [26], and the results were published in Applicable Analysis.
(ii): If $M (f (\hat{a}), g (\hat{a})) = M (\hat{a})$ and $F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) = Y_{I, λ}^{M (f, g)} (\hat{a})$ , then problem (9) reduces to the problem of finding $\hat{a} \in \hat{E}$ such that

$0 \in Y_{I, λ}^{M} (\hat{a}) + M (\hat{a}) .$

(11)

Problem (11) was investigated and studied by Arvind et al. [27] in 2023, and the results were published in Mathematics.

3. Fixed-Point Equation and Three-Step Iterative Algorithm

The subsequent lemma presents a fixed-point formulation for the extended Cayley–Yosida inclusion problem, (9), incorporating the generalized resolvent and Yosida approximation operators defined by (3) and (4), respectively.

Lemma 3.

The extended Cayley–Yosida inclusion problem in (9) has a solution

\hat{a} \in \hat{E}

if and only if the following equation holds:

\hat{a} = R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} (\hat{a}) - F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a}))) + R_{I, λ}^{M (f, g)} (\hat{a})] .

(12)

Proof.

Let

\hat{a} \in \hat{E}

satisfy Equation (12); then,

\begin{matrix} \hat{a} & = & R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} (\hat{a}) - F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a}))) + R_{I, λ}^{M (f, g)} (\hat{a})] \\ \hat{a} & = & {(I + λ M (f, g))}^{- 1} [λ (Y_{I, λ}^{M (f, g)} (\hat{a}) - F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a}))) \\ + R_{I, λ}^{M (f, g)} (\hat{a})], \end{matrix}

which implies that

\begin{matrix} I (\hat{a}) + λ M (f (\hat{a}), g (\hat{a})) & = & λ \frac{1}{λ} [I (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{a})] - λ F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) \\ + R_{I, λ}^{M (f, g)} (\hat{a}) \\ 0 & \in & F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) + M (f (\hat{a}), g (\hat{a})) . \end{matrix}

Conversely, suppose that

\hat{a}

is the solution of problem (9). That is,

0 \in F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) + M (f (\hat{a}), g (\hat{a})) .

From above, we have

\begin{matrix} - F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) & \in M (f (\hat{a}), g (\hat{a})) \\ I (\hat{a}) - λ F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) & \in I (\hat{a}) + λ M (f (\hat{a}), g (\hat{a})), \\ = [I + λ M (f, g)] (\hat{a}), \end{matrix}

which implies that

{[I + λ M (f, g)]}^{- 1} [I (\hat{a}) - λ F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a}))] = \hat{a}

or

\begin{matrix} \hat{a} & = R_{I, λ}^{M (f, g)} [I (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{a}) + R_{I, λ}^{M (f, g)} (\hat{a}) - λ F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a}))] \\ \hat{a} & = R_{I, λ}^{M (f, g)} [λ \frac{1}{λ} {I - R_{I, λ}^{M (f, g)}} (\hat{a}) + R_{I, λ}^{M (f, g)} (\hat{a}) - λ F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a}))] \\ \hat{a} & = R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} (\hat{a}) - F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a}))) + R_{I, λ}^{M (f, g)} (\hat{a})] . \end{matrix}

Thus, Equation (12) is fulfilled. □

3.1. Existence Result

The following theorem gives the guarantee for the existence of the solution of problem (9).

Theorem 2.

Let

\hat{E}

be a real 2-USBS with a modulus of smoothness

τ_{\hat{E}} (t) \leq k t^{2}

for some constant

k > 0

. Let

f, g : \hat{E} \to \hat{E}

be the single-valued mappings and

I : \hat{E} \to \hat{E}

be an identity mapping. Let

F : \hat{E} \times \hat{E} \to \hat{E}

be a

(\land_{1}, \land_{2})

-Lipschitz continuous mapping with non-negative constants

\land_{1}, \land_{2} > 0

. Let

M : \hat{E} \times \hat{E} \to 2^{\hat{E}}

be

α β

–I-accretive with regard to

f

and

g

. Then, the extended Cayley–Yosida inclusion problem has a unique solution

\hat{a}

in

\hat{E}

if

|θ - (- \frac{£}{2})| < \frac{\sqrt{£^{2} + 4}}{2},

(13)

where

£ = λ [▵ + \land_{1} θ_{C} + \land_{2} θ_{Y} + 1], ▵ = \sqrt{1 - 2 m + 64 k θ_{Y}^{2} + k θ_{Y}}, θ_{C} = 2 θ + 1,

θ_{Y} = \frac{1 + θ}{λ}, m = \frac{1 - θ}{λ}, θ = \frac{1}{1 + λ (α - β)}, and λ, α, a n d β a r e n o n - n e g a t i v e c o n s t a n t s

with

α \geq β .

Proof.

Let us define

P (\hat{a}) = Y_{I, λ}^{M (f, g)} (\hat{a}) - F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a}))

. For any

\hat{a}, \hat{b} \in \hat{E}

, we will prove that

R_{I, λ}^{M (f, g)} \{λ P + R_{I, λ}^{M (f, g)}\}

is a contraction mapping.

\begin{matrix} ∥R_{I, λ}^{M (f, g)} \{λ P (\hat{a}) + R_{I, λ}^{M (f, g)} (\hat{a})\} - R_{I, λ}^{M (f, g)} \{λ P (\hat{b}) + R_{I, λ}^{M (f, g)} (\hat{b})\}∥ \\ \leq θ ∥λ (P (\hat{a}) - P (\hat{b})) + R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b})∥ \\ \leq θ \{λ ∥ P (\hat{a}) - P (\hat{b}) ∥ + ∥R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b})∥\} \\ \leq θ \{λ ∥ P (\hat{a}) - P (\hat{b}) ∥ + θ ∥ \hat{a} - \hat{b} ∥\} . \end{matrix}

(14)

On the other hand, we have

\begin{matrix} ∥ P (\hat{a}) - P (\hat{b}) ∥ & \leq ∥Y_{I, λ}^{M (f, g)} (\hat{a}) - Y_{I, λ}^{M (f, g)} (\hat{b}) - (\hat{a} - \hat{b})∥ + ∥ F (C_{I, λ}^{M (f, g)} (\hat{b}), Y_{I, λ}^{M (f, g)} (\hat{b})) \\ - F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) + (\hat{a} - \hat{b}) ∥ . \end{matrix}

(15)

Using the strong accretivity of the Yosida approximation operator, we obtain

\begin{matrix} {∥(\hat{a} - \hat{a}) - (Y_{I, λ}^{M (f, g)} (\hat{a}) - Y_{I, λ}^{M (f, g)} (\hat{a}))∥}^{2} \\ \leq & ∥ \hat{b} - \hat{a} ∥^{2} + 2 〈 - (Y_{I, λ}^{M (f, g)} (\hat{a}) - Y_{I, λ}^{M (f, g)} (\hat{b})), \\ j ((\hat{a} - \hat{b}) - (Y_{I, λ}^{M (f, g)} (\hat{a}) - Y_{I, λ}^{M (f, g)} (\hat{b}))) 〉 + k {∥ Y_{I, λ}^{M (f, g)} (\hat{a}) - Y_{I, λ}^{M (f, g)} (\hat{b}) ∥}^{2} \\ = & ∥ \hat{a} - \hat{b} ∥^{2} + 2 〈- (Y_{I, λ}^{M (f, g)} (\hat{a}) - Y_{I, λ}^{M (f, g)} (\hat{a})), j (\hat{a} - \hat{b})〉 + 2 〈 - (Y_{I, λ}^{M (f, g)} (\hat{a}) - Y_{I, λ}^{M (f, g)} (\hat{b})), \\ j ((\hat{a} - \hat{b}) - (Y_{I, λ}^{M (f, g)} (\hat{a}) - Y_{I, λ}^{M (f, g)} (\hat{b}))) - j (\hat{a} - \hat{b}) 〉 + k {∥ Y_{I, λ}^{M (f, g)} (\hat{a}) - Y_{I, λ}^{M (f, g)} (\hat{b}) ∥}^{2} \\ \leq & ∥ \hat{a} - \hat{b} ∥^{2} - 2 m ∥ \hat{a} - \hat{b} ∥^{2} + 4 d^{2} τ_{\hat{E}} (\frac{4 ∥Y_{I, λ}^{M (f, g)} (\hat{a}) - Y_{I, λ}^{M (f, g)} (\hat{b})∥}{d}) + k θ_{Y} {∥ \hat{a} - \hat{b} ∥}^{2} \\ \leq & ∥ \hat{a} - \hat{b} ∥^{2} - 2 m ∥ \hat{a} - \hat{b} ∥^{2} + 64 k {∥Y_{I, λ}^{M (f, g)} (\hat{a}) - Y_{I, λ}^{M (f, g)} (\hat{b})∥}^{2} + k θ_{Y} {∥ \hat{a} - \hat{b} ∥}^{2} \\ \leq & [1 - 2 m + 64 k θ_{Y}^{2} + k θ_{Y}] {∥ \hat{a} - \hat{b} ∥}^{2} . \end{matrix}

(16)

Since

F

,

C_{I, λ}^{M (f, g)}

and

Y_{I, λ}^{M (f, g)}

are

(\land_{1}, \land_{2})

-,

θ_{C}

- and

θ_{Y}

-Lipschitz continuous, respectively, we have

\begin{matrix} ∥F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) - F (C_{I, λ}^{M (f, g)} (\hat{b}), Y_{I, λ}^{M (f, g)} (\hat{b})) + (\hat{a} - \hat{b})∥ \\ \leq & ∥ \hat{a} - \hat{b} ∥ + ∥F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) - F (C_{I, λ}^{M (f, g)} (\hat{b}), Y_{I, λ}^{M (f, g)} (\hat{a}))∥ \\ + ∥F (C_{I, λ}^{M (f, g)} (\hat{b}), Y_{I, λ}^{M (f, g)} (\hat{a})) - F (C_{I, λ}^{M (f, g)} (\hat{b}), Y_{I, λ}^{M (f, g)} (\hat{b}))∥ \\ \leq & ∥ \hat{a} - \hat{b} ∥ + \land_{1} ∥C_{I, λ}^{M (f, g)} (\hat{a}) - C_{I, λ}^{M (f, g)} (\hat{b})∥ + \land_{2} ∥Y_{I, λ}^{M (f, g)} (\hat{a}) - Y_{I, λ}^{M (f, g)} (\hat{b})∥ \\ \leq & ∥ \hat{a} - \hat{b} ∥ + \land_{1} θ_{C} ∥ \hat{a} - \hat{b} ∥ + \land_{2} θ_{Y} ∥ \hat{a} - \hat{b} ∥ \\ = & [1 + \land_{1} θ_{C} + \land_{2} θ_{Y}] ∥ \hat{a} - \hat{b} ∥, \end{matrix}

(17)

Using (16) and (17) in (15), we obtain

∥ P (\hat{a}) - P (\hat{b}) ∥ \leq \{\sqrt{[1 - 2 m + 64 k θ_{Y}^{2} + k θ_{Y}]} + [1 + \land_{1} θ_{C} + \land_{2} θ_{Y}]\} ∥ \hat{a} - \hat{b} ∥ .

(18)

From (18) and (14), we obtain

∥R_{I, λ}^{M (f, g)} \{λ P (\hat{a}) + R_{I, λ}^{M (f, g)} (\hat{a})\} - R_{I, λ}^{M (f, g)} \{λ P (\hat{b}) + R_{I, λ}^{M (f, g)} (\hat{b})\}∥ \leq ¶ ∥ \hat{a} - \hat{b} ∥,

where

¶ = θ [λ \{\sqrt{[1 - 2 m + 64 k θ_{Y}^{2} + k θ_{Y}]} + [1 + \land_{1} θ_{C} + \land_{2} θ_{Y}]\} + θ] .

From (13),

¶ < 1

. Hence,

R_{I, λ}^{M (f, g)} \{λ P + R_{I, λ}^{M (f, g)}\}

is a contraction mapping. Since

\hat{E}

is a 2-Uniformly Smooth Banach Space, from Lemma 3, problem (9) has a unique solution

\hat{a}

which is a fixed point of

R_{I, λ}^{M (f, g)} \{λ P + R_{I, λ}^{M (f, g)}\} .

□

3.2. Three-Step Iterative Algorithm

Drawing upon Lemma 3 as a cornerstone, we formulate a three-step iterative algorithm to ascertain the solution of the extended Cayley–Yosida inclusion problem in (9). This iterative algorithm incorporates both the generalized resolvent and Yosida approximation operators, presenting itself as a novel approach.

The following result is from the convergence and stability analysis of Cayley–Yosida inclusion problem (9).

Theorem 3.

Assume

\hat{E}

,

f, g, F, I

, and

M

are all the same as as in Theorem 2, such that all the conditions specified in Theorem 2 are fulfilled. Moreover, the following condition is satisfied:

|θ - (- \frac{£}{2})| < \frac{\sqrt{£^{2} + 4}}{6},

(19)

where £ is the same as in (13) and

\underset{n \to \infty}{l i m} {\hat{δ_{1}}}_{n} = \underset{n \to \infty}{l i m} {\hat{δ_{2}}}_{n} = \underset{n \to \infty}{l i m} {\hat{δ_{3}}}_{n}

. Then, we have the following:

(I): The sequence ${{\hat{a}}_{n}}$ produced by Algorithm 1 converges strongly to the unique solution of problem (9).
(II): Additionally, the sequence ${{\hat{a}}_{n}}$ produced by Algorithm 1 is stable with respect to $R_{I, λ}^{M (f, g)} .$

Proof.

(I) Convergence of

{\hat{a}}_{n}

:

On the basis of Lemma 3, we can write

\begin{matrix} {\hat{a}}^{*} & = (1 - α_{n}) {\hat{a}}^{*} + α_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*}) - F (C_{I, λ}^{M (f, g)} ({\hat{a}}^{*}), Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*}))) \\ + R_{I, λ}^{M (f, g)} ({\hat{a}}^{*})] . \end{matrix}

Now, from Algorithm 1 and Theorem 1, we have

\begin{matrix} ∥ {\hat{a}}_{n + 1} - {\hat{a}}^{*} ∥ \\ = ∥ (1 - α_{n}) {\hat{a}}_{n} + α_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{b}}_{n}) - F (C_{I, λ}^{M (f, g)} ({\hat{b}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{b}}_{n}))) \\ + R_{I, λ}^{M (f, g)} ({\hat{b}}_{n})] + α_{n} {\hat{δ_{1}}}_{n} - {(1 - α_{n}) {\hat{a}}^{*} + α_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*}) \\ - F (C_{I, λ}^{M (f, g)} ({\hat{a}}^{*}), Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*}))) + R_{I, λ}^{M (f, g)} ({\hat{a}}^{*})]} ∥ \\ \leq | (1 - α_{n}) | ∥ {\hat{a}}_{n} - {\hat{a}}^{*} ∥ + | α_{n} | ∥ R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{b}}_{n}) - F (C_{I, λ}^{M (f, g)} ({\hat{b}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{b}}_{n})) \\ + R_{I, λ}^{M (f, g)} ({\hat{b}}_{n})] - R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*}) - F (C_{I, λ}^{M (f, g)} ({\hat{a}}^{*}), Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*})) \\ + R_{I, λ}^{M (f, g)} ({\hat{a}}^{*})] ∥ + ∥ α_{n} {\hat{δ_{1}}}_{n} ∥ . \end{matrix}

(20)

Using a Lipschitz continuity of

R_{I, λ}^{M (f, g)}

, we have

\begin{matrix} \leq | (1 - α_{n}) | ∥ {\hat{a}}_{n} - {\hat{a}}^{*} ∥ + | α_{n} | θ ∥ ∥ λ (Y_{I, λ}^{M (f, g)} ({\hat{b}}_{n}) - Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*})) \\ - λ F (C_{I, λ}^{M (f, g)} ({\hat{b}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{b}}_{n})) + λ F (C_{I, λ}^{M (f, g)} ({\hat{a}}^{*}), Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*})) \\ + R_{I, λ}^{M (f, g)} ({\hat{b}}_{n}) - R_{I, λ}^{M (f, g)} ({\hat{a}}^{*}) ∥ + ∥ α_{n} {\hat{δ_{1}}}_{n} ∥ \\ \leq | (1 - α_{n}) | ∥ {\hat{a}}_{n} - {\hat{a}}^{*} ∥ + | α_{n} | θ λ ∥ (Y_{I, λ}^{M (f, g)} ({\hat{b}}_{n}) - Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*})) - ({\hat{a}}_{n} - {\hat{a}}^{*}) ∥ \\ + λ ∥ F (C_{I, λ}^{M (f, g)} ({\hat{b}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{b}}_{n})) - F (C_{I, λ}^{M (f, g)} ({\hat{a}}^{*}), Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*})) + ({\hat{b}}_{n} - {\hat{a}}^{*}) ∥ \\ + ∥ R_{I, λ}^{M (f, g)} ({\hat{b}}_{n}) - R_{I, λ}^{M (f, g)} ({\hat{a}}^{*}) ∥ + ∥ α_{n} {\hat{δ_{1}}}_{n} ∥ . \end{matrix}

(21)

Since

F

is

(\land_{1}, \land_{2})

-Lipschitz continuous,

C_{I, λ}^{M (f, g)}

is

θ_{C}

-Lipschitz continuous, and

Y_{I, λ}^{M (f, g)}

is

m

-strongly accretive and

θ_{Y}

-Lipschitz continuous, from Lemma 1, we obtain

\begin{matrix} ∥ {\hat{a}}_{n + 1} - {\hat{a}}^{*} ∥ \\ \leq | (1 - α_{n}) | ∥ {\hat{a}}_{n} - {\hat{a}}^{*} ∥ + | α_{n} | θ λ (\sqrt{[1 - 2 m + 64 k θ_{Y}^{2} + k θ_{Y}]}) ∥ {\hat{b}}_{n} - x^{*} ∥ \\ + | α_{n} | θ λ [1 + \land_{1} θ_{C} + \land_{2} θ_{Y}] ∥ {\hat{b}}_{n} - {\hat{a}}^{*} ∥ + | α_{n} | θ^{2} ∥ {\hat{b}}_{n} - {\hat{a}}^{*} ∥ + ∥ α_{n} {\hat{δ_{1}}}_{n} ∥ \\ = | (1 - α_{n}) | ∥ {\hat{a}}_{n} - {\hat{a}}^{*} ∥ + ¶ | α_{n} | ∥ {\hat{b}}_{n} - {\hat{a}}^{*} ∥ + ∥ α_{n} {\hat{δ_{1}}}_{n} ∥, \end{matrix}

(22)

where

¶ = θ [λ \{\sqrt{[1 - 2 m + 64 k θ_{Y}^{2} + k θ_{Y}]} + [1 + \land_{1} θ_{C} + \land_{2} θ_{Y}]\} + θ] .

From Algorithm 1, we calculate

\begin{matrix} ∥ {\hat{b}}_{n} - {\hat{a}}^{*} ∥ \\ = ∥ (1 - β_{n}) {\hat{a}}_{n} + β_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{z}}_{n}) - F (C_{I, λ}^{M (f, g)} ({\hat{z}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{z}}_{n})) \\ + R_{I, λ}^{M (f, g)} ({\hat{z}}_{n})] + β_{n} {\hat{δ_{2}}}_{n} - {(1 - β_{n}) {\hat{a}}^{*} + β_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*}) \\ - F (C_{I, λ}^{M (f, g)} ({\hat{a}}^{*}), Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*})) + R_{I, λ}^{M (f, g)} ({\hat{a}}^{*})]} ∥ \\ \leq | (1 - β_{n}) | ∥ {\hat{a}}_{n} - {\hat{a}}^{*} ∥ + | β_{n} | θ λ ∥ (Y_{I, λ}^{M (f, g)} ({\hat{z}}_{n}) - Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*})) - (\hat{b} - {\hat{a}}^{*}) ∥ \\ + λ ∥ F (C_{I, λ}^{M (f, g)} ({\hat{z}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{z}}_{n}) - F (C_{I, λ}^{M (f, g)} ({\hat{a}}^{*}), Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*})) + (\hat{b} - {\hat{a}}^{*}) ∥ \\ + | β_{n} | ∥ R_{I, λ}^{M (f, g)} ({\hat{b}}_{n}) - R_{I, λ}^{M (f, g)} ({\hat{a}}^{*}) ∥ + ∥ β_{n} {\hat{δ_{2}}}_{n} ∥ . \end{matrix}

(23)

By using the same argument as above, we obtain

\begin{matrix} ∥ {\hat{b}}_{n} - {\hat{a}}^{*} ∥ \\ \leq | (1 - β_{n}) | ∥ {\hat{a}}_{n} - {\hat{a}}^{*} ∥ + | β_{n} | θ λ (\sqrt{[1 - 2 m + 64 k θ_{Y}^{2} + k θ_{Y}]}) ∥ {\hat{z}}_{n} - {\hat{a}}^{*} ∥ \\ + [1 + \land_{1} θ_{C} + \land_{2} θ_{Y}] ∥ {\hat{b}}_{n} - {\hat{a}}^{*} ∥ + ∥ β_{n} {\hat{δ_{2}}}_{n} ∥ \\ = | (1 - β_{n}) | ∥ {\hat{a}}_{n} - {\hat{a}}^{*} ∥ + ¶ | β_{n} | ∥ {\hat{z}}_{n} - {\hat{a}}^{*} ∥ + ∥ β_{n} {\hat{δ_{2}}}_{n} ∥ . \end{matrix}

(24)

Similarly, we may calculate

\begin{matrix} ∥ {\hat{z}}_{n} - {\hat{a}}^{*} ∥ \leq | (1 - γ_{n}) | ∥ {\hat{a}}_{n} - {\hat{a}}^{*} ∥ + ¶ | γ_{n} | ∥ {\hat{z}}_{n} - {\hat{a}}^{*} ∥ + ∥ γ_{n} {\hat{δ_{3}}}_{n} ∥ . \end{matrix}

(25)

Algorithm 1 Three-Step Iterative Algorithm

Let

{\hat{a}}_{0} \in \hat{E},

λ > 0

,

F : \hat{E} \times \hat{E} \to \hat{E}

be a single-valued mapping, and

M : \hat{E} \times \hat{E} \to 2^{E}

be a multi-valued mapping. Let

C_{I, λ}^{M (f, g)}

and

Y_{I, λ}^{M (f, g)}

be the generalized Cayley and Yosida approximation operators defined in (3) and (4), respectively. Then, we develop the following iterative scheme:

\begin{matrix} {\hat{a}}_{n + 1} & = (1 - α_{n}) {\hat{a}}_{n} + α_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{b}}_{n}) - F (C_{I, λ}^{M (f, g)} ({\hat{b}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{b}}_{n}))) \\ + R_{I, λ}^{M (f, g)} ({\hat{b}}_{n})] + α_{n} {\hat{δ_{1}}}_{n} \\ {\hat{b}}_{n} & = (1 - β_{n}) {\hat{a}}_{n} + β_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{z}}_{n}) - F (C_{I, λ}^{M (f, g)} ({\hat{z}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{z}}_{n}))) \\ + R_{I, λ}^{M (f, g)} ({\hat{z}}_{n})] + β_{n} {\hat{δ_{2}}}_{n} \\ {\hat{z}}_{n} & = (1 - γ_{n}) {\hat{a}}_{n} + γ_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{a}}_{n}) - F (C_{I, λ}^{M (f, g)} ({\hat{a}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{a}}_{n}))) \\ + R_{I, λ}^{M (f, g)} ({\hat{a}}_{n})] + γ_{n} {\hat{δ_{3}}}_{n} . \end{matrix}

Let

{{\hat{u}}_{n}} \subseteq \hat{E}

be a sequence and define

{{\hat{v}}_{n}}

by

\begin{matrix} {\hat{v}}_{n} & = ∥ {\hat{u}}_{n + 1} - {(1 - α_{n}) {\hat{a}}_{n} + α_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{t}}_{n}) - F (C_{I, λ}^{M (f, g)} ({\hat{t}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{t}}_{n}))) \\ + R_{I, λ}^{M (f, g)} ({\hat{t}}_{n})] + α_{n} {\hat{δ_{1}}}_{n}} ∥ \\ {\hat{t}}_{n} & = (1 - β_{n}) {\hat{u}}_{n} + β_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{s}}_{n}) - F (C_{I, λ}^{M (f, g)} ({\hat{s}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{s}}_{n}))) \\ + R_{I, λ}^{M (f, g)} ({\hat{s}}_{n})] + β_{n} {\hat{δ_{2}}}_{n} \\ {\hat{s}}_{n} & = (1 - γ_{n}) {\hat{u}}_{n} + γ_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{u}}_{n}) - F (C_{I, λ}^{M (f, g)} ({\hat{u}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{u}}_{n}))) \\ + R_{I, λ}^{M (f, g)} ({\hat{u}}_{n})] + γ_{n} {\hat{δ_{3}}}_{n}, \end{matrix}

where

α_{n}, β_{n}, γ_{n} \in [0, 1]

,

\sum_{n = 0}^{\infty} α_{n} = \infty, a n d \forall n \geq 0

. Here,

{{\hat{δ_{1}}}_{n}}, {{\hat{δ_{2}}}_{n}}, {{\hat{δ_{3}}}_{n}} \subseteq \hat{E}

are sequences in

\hat{E}

considered to accommodate possible computational inaccuracies.

Using (24) and (25), (22) becomes

\begin{matrix} ∥ {\hat{a}}_{n + 1} - {\hat{a}}^{*} ∥ & \leq & {(1 - α_{n}) + ¶ α_{n} (1 - β_{n}) + ¶^{2} α_{n} β_{n} (1 - γ_{n}) \\ + ¶^{3} α_{n} β_{n} γ_{n} + ¶^{2} α_{n} β_{n} {\hat{δ_{3}}}_{n} γ_{n}} ∥ {\hat{a}}_{n} - {\hat{a}}^{*} ∥ + ¶ α_{n} β_{n} {\hat{δ_{2}}}_{n} + {\hat{δ_{1}}}_{n} α_{n} \\ \leq & (1 - α_{n} (1 - 3 ¶)) ∥ {\hat{a}}_{n} - {\hat{a}}^{*} ∥ + α_{n} (1 - 3 ¶) \{\frac{¶ β_{n} {\hat{δ_{2}}}_{n} + {\hat{δ_{1}}}_{n}}{1 - 3 ¶}\} . \end{matrix}

(26)

On setting

η_{n} : = \{\frac{¶ β_{n} {\hat{δ_{2}}}_{n} + {\hat{δ_{1}}}_{n}}{1 - 3 ¶}\}, ξ_{n} = ∥ {\hat{a}}_{n} - {\hat{a}}^{*} ∥

, and

ζ_{n} : = {1 - α_{n} (1 - 3 ¶)},

then (26) can be written as

ξ_{n} \leq (1 - ζ_{n}) ξ_{n} + ζ_{n} η_{n} .

By Lemma 2, we can conclude that

ξ_{n} \to 0

, as

n \to \infty

, and so

{{\hat{a}}_{n}} \to {\hat{a}}^{*}

strongly.

(II) Stability of

{\hat{a}}_{n}

:

Let

H ({\hat{a}}^{*}) = R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*}) - F (C_{I, λ}^{M (f, g)} ({\hat{a}}^{*}), Y_{I, λ}^{M (f, g)} ({\hat{a}}^{*})) + R_{I, λ}^{M (f, g)} ({\hat{a}}^{*})]

. From Algorithm 1, we have

\begin{matrix} ∥ {\hat{u}}_{n + 1} - {\hat{a}}^{*} ∥ & = ∥ {\hat{u}}_{n + 1} - {(1 - α_{n}) {\hat{a}}^{*} + α_{n} H ({\hat{a}}^{*})} ∥ \\ \leq ∥{\hat{u}}_{n + 1} - (1 - α_{n}) {\hat{u}}_{n} - α_{n} H (t_{n}) - {\hat{δ_{1}}}_{n} α_{n}∥ \\ + ∥(1 - α_{n}) {\hat{u}}_{n} + α_{n} H ({\hat{t}}_{n}) + {\hat{δ_{1}}}_{n} α_{n} - (1 - α_{n}) {\hat{a}}^{*} - α_{n} H ({\hat{a}}^{*})∥ \\ \leq ∥{\hat{u}}_{n + 1} - (1 - α_{n}) {\hat{u}}_{n} - α_{n} H ({\hat{t}}_{n}) - {\hat{δ_{1}}}_{n} α_{n}∥ \\ + | 1 - α_{n} | ∥ {\hat{u}}_{n} - {\hat{a}}^{*} ∥ + α_{n} ∥ H ({\hat{t}}_{n}) - H ({\hat{a}}^{*}) ∥ + ∥ {\hat{δ_{1}}}_{n} α_{n} ∥ \\ \leq ∥{\hat{u}}_{n + 1} - (1 - α_{n}) {\hat{u}}_{n} - α_{n} H ({\hat{t}}_{n}) - {\hat{δ_{1}}}_{n} α_{n}∥ \\ + | 1 - α_{n} | ∥ {\hat{δ_{1}}}_{n} - {\hat{a}}^{*} ∥ + ¶ | α_{n} | ∥ {\hat{t}}_{n} - {\hat{a}}^{*} ∥ + ∥ {\hat{δ_{1}}}_{n} α_{n} ∥, \end{matrix}

(27)

where

¶ = θ [λ \{\sqrt{[1 - 2 m + 64 k θ_{Y}^{2} + k θ_{Y}]} + [1 + \land_{1} θ_{C} + \land_{2} θ_{Y}]\} + θ] .

\begin{matrix} ∥ {\hat{t}}_{n} - {\hat{a}}^{*} ∥ & \leq ∥(1 - β_{n}) {\hat{u}}_{n} + β_{n} H ({\hat{s}}_{n}) + {\hat{δ_{2}}}_{n} β_{n} - (1 - β_{n}) {\hat{a}}^{*} - β_{n} H ({\hat{a}}^{*})∥ \\ \leq | 1 - β_{n} | ∥ {\hat{u}}_{n} - {\hat{a}}^{*} ∥ + β_{n} ∥ H ({\hat{s}}_{n}) - H ({\hat{a}}^{*}) ∥ + ∥ {\hat{δ_{2}}}_{n} β_{n} ∥ \\ \leq | 1 - β_{n} | ∥ {\hat{u}}_{n} - {\hat{a}}^{*} ∥ + ¶ | β_{n} | ∥ {\hat{s}}_{n} - {\hat{a}}^{*} ∥ ∥ {\hat{δ_{2}}}_{n} β_{n} ∥ . \end{matrix}

(28)

Similarly, we obtain

\begin{matrix} ∥ {\hat{s}}_{n} - {\hat{a}}^{*} ∥ & \leq ∥(1 - γ_{n}) {\hat{u}}_{n} + γ_{n} H ({\hat{u}}_{n}) + {\hat{δ_{3}}}_{n} γ_{n} - (1 - γ_{n}) {\hat{a}}^{*} - γ_{n} H ({\hat{a}}^{*})∥ \\ \leq | 1 - γ_{n} | ∥ {\hat{u}}_{n} - {\hat{a}}^{*} ∥ + γ_{n} ∥ H ({\hat{u}}_{n}) - H ({\hat{a}}^{*}) ∥ + ∥ {\hat{δ_{3}}}_{n} γ_{n} ∥ \\ \leq | 1 - γ_{n} | ∥ {\hat{u}}_{n} - {\hat{a}}^{*} ∥ + ¶ | γ_{n} | ∥ {\hat{u}}_{n} - {\hat{a}}^{*} ∥ ∥ {\hat{δ_{3}}}_{n} γ_{n} ∥ . \end{matrix}

(29)

Using (28) and (29) in (27), we have

\begin{matrix} ∥ {\hat{u}}_{n + 1} - {\hat{a}}^{*} ∥ & \leq ∥{\hat{u}}_{n + 1} - (1 - α_{n}) {\hat{u}}_{n} - α_{n} H ({\hat{t}}_{n}) - {\hat{δ_{1}}}_{n} α_{n}∥ \\ + | (1 - α_{n}) ∥ {\hat{u}}_{n} - {\hat{a}}^{*} ∥ + (¶ α_{n} (1 - β_{n}) + ¶^{2} α_{n} β_{n} (1 - γ_{n}) \\ + β_{n} α_{n} γ_{n} ¶^{3}) ∥ {\hat{u}}_{n} - {\hat{a}}^{*} ∥ + α_{n} β_{n} ¶ ∥ {\hat{δ_{3}}}_{n} γ_{n} ∥ + ∥ {\hat{δ_{1}}}_{n} α_{n} ∥ \\ ∥ {\hat{u}}_{n + 1} - {\hat{a}}^{*} ∥ & \leq {\hat{v}}_{n} + (1 - α_{n} (1 - 3 ¶)) ∥ {\hat{u}}_{n} - {\hat{a}}^{*} ∥ + α_{n} (1 - 3 ¶) \\ \times \{\frac{¶ β_{n} {\hat{δ_{2}}}_{n} + {\hat{δ_{1}}}_{n}}{1 - 3 ¶}\} \end{matrix}

(30)

Assume that

\underset{n \to \infty}{l i m} {\hat{v}}_{n} = 0

. Hence,

\underset{n \to \infty}{l i m} {\hat{u}}_{n} = {\hat{a}}^{*}

, where

\underset{n \to \infty}{l i m} {\hat{δ_{1}}}_{n} = \underset{n \to \infty}{l i m} {\hat{δ_{2}}}_{n} = \underset{n \to \infty}{l i m} {\hat{δ_{3}}}_{n} = 0 .

Similarly, one can prove that if

\underset{n \to \infty}{l i m} {\hat{u}}_{n} = {\hat{a}}^{*}

and condition (19) holds then

\underset{n \to \infty}{l i m} {\hat{v}}_{n} = 0 .

Hence, the convergence of

{{\hat{a}}_{n}}

is stable with regard to

R_{I, λ}^{M (f, g)} .

□

4. Numerical Example

We now present the following numerical example to illustrate the convergence analysis of the sequence

{{\hat{a}}_{n}}

towards the solution of problem (9).

Example 1.

Assume

\hat{E} = R

,

f, g : \hat{E} \to \hat{E}

are single-valued mappings, and

M : \hat{E} \times \hat{E} \to 2^{\hat{E}}

is a multi-valued mapping defined by

f (\hat{a}) = (\frac{6 \hat{a}}{5}), g (\hat{a}) = (\frac{5 \hat{a}}{3}), M (\hat{a}, \hat{b}) = \{\frac{\hat{a} + \hat{b}}{3}\} .

Now, for

λ = 1 / 2,

the associate resolvent operator

R_{I, λ}^{M (f, g)} (\hat{a}) = {[I + λ M (f, g)]}^{- 1} (\hat{a}) = (\frac{90 \hat{a}}{133}) .

Furthermore,

\begin{matrix} ∥ R_{I, λ}^{M (f, g)} (\hat{a}) - R_{I, λ}^{M (f, g)} (\hat{b}) ∥ & = & ∥\frac{90}{133} \hat{a} - \frac{90}{133} \hat{b}∥ \\ = & \frac{90}{133} ∥\hat{a} - \hat{b}∥ \leq \frac{9}{11} ∥ \hat{a} - \hat{b} ∥, \end{matrix}

i.e.,

R_{I, λ}^{M (f, g)}

is Lipschitz continuous with constant

θ = \frac{9}{11}

.

Now, we define the Yosida approximation operator as

Y_{I, λ}^{M (f, g)} (\hat{a}) = \frac{1}{λ} [I - R_{I, λ}^{M (f, g)}] (\hat{a}) = (\frac{43 \hat{a}}{266}) .

Next, the Cayley operator is defined as

C_{I, λ}^{M (f, g)} (\hat{a}) = 2 R_{I, λ}^{M (f, g)} (\hat{a}) - I (\hat{a}) = (\frac{47 \hat{a}}{133}) .

It is straightforward that

Y_{I, λ}^{M (f, g)}

and

C_{I, λ}^{M (f, g)}

are Lipschitz continuous with constants

θ_{Y} = \frac{40}{11}

and

θ_{C} = \frac{29}{11}

, respectively.

For

λ = \frac{1}{2}

, we calculate

\begin{matrix} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} (\hat{a}) - F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) + R_{I, λ}^{M (f, g)} (\hat{a})] \\ = \frac{90}{133} [\frac{1}{2} (Y_{I, λ}^{M (f, g)} (\hat{a}) - F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) + R_{I, λ}^{M (f, g)} (\hat{a})] \\ = \frac{90}{133} [\frac{1}{2} (\frac{43}{266} (\hat{a}) - \frac{C_{I, λ}^{M (f, g)} (\hat{a}) + Y_{I, λ}^{M (f, g)} (\hat{a})}{3}) + \frac{90}{133} (\hat{a})] \\ = \frac{90}{133} [\frac{1}{2} (\frac{43}{266} (\hat{a}) - \frac{137}{266 \times 3} (\hat{a})) + \frac{90}{133} (\hat{a})] \\ = \frac{8040}{17689} (\hat{a}) . \end{matrix}

(31)

Clearly, 0 is the fixed point of

R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} (\hat{a}) - F (C_{I, λ}^{M (f, g)} (\hat{a}), Y_{I, λ}^{M (f, g)} (\hat{a})) + R_{I, λ}^{M (f, g)} (\hat{a})] .

Let

{\hat{δ_{1}}}_{n} = \frac{1}{2 n}, {\hat{δ_{2}}}_{n} = \frac{n + 2}{n + 3}, {\hat{δ_{3}}}_{n} = \frac{2 n}{2 n + 1}, α_{n} = \frac{n}{2 n + 1}, β_{n} = \frac{1}{n + 1}, and γ_{n} = \frac{n + 1}{n (n + 2)} .

It is clear that the sequences

{\hat{δ_{1}}}_{n}, {\hat{δ_{2}}}_{n}, {\hat{δ_{3}}}_{n}, α_{n}, β_{n}, a n d γ_{n}

satisfy condition (19). We compute the sequence

{{\hat{p}}_{n}}

using the following iterative scheme:

\begin{matrix} {\hat{a}}_{n + 1} & = & (1 - α_{n}) {\hat{a}}_{n} + α_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{b}}_{n}) - F (C_{I, λ}^{M (f, g)} ({\hat{b}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{b}}_{n})) \\ + R_{I, λ}^{M (f, g)} ({\hat{b}}_{n})] + α_{n} {\hat{δ_{1}}}_{n} \\ {\hat{a}}_{n + 1} & = & (\frac{n + 1}{2 n + 1}) {\hat{a}}_{n} + (\frac{8040 n}{17689 (2 n + 1)}) {\hat{b}}_{n} + (\frac{1}{2 (2 n + 1)}), \end{matrix}

\begin{matrix} {\hat{b}}_{n} & = & (1 - β_{n}) {\hat{a}}_{n} + β_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{z}}_{n}) - F (C_{I, λ}^{M (f, g)} ({\hat{z}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{z}}_{n})) \\ + R_{I, λ}^{M (f, g)} ({\hat{z}}_{n})] + β_{n} {\hat{δ_{2}}}_{n} \\ {\hat{b}}_{n} & = & (\frac{n}{n + 1}) {\hat{a}}_{n} + (\frac{8040}{17689 (n + 1)}) {\hat{z}}_{n} + (\frac{n + 2}{(n + 1) (n + 3)}), \end{matrix}

and

\begin{matrix} {\hat{z}}_{n} & = & (1 - γ_{n}) {\hat{a}}_{n} + γ_{n} R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)} ({\hat{a}}_{n}) - F (C_{I, λ}^{M (f, g)} ({\hat{a}}_{n}), Y_{I, λ}^{M (f, g)} ({\hat{a}}_{n})) \\ + R_{I, λ}^{M (f, g)} ({\hat{a}}_{n})] + γ_{n} {\hat{δ_{3}}}_{n} . \\ {\hat{z}}_{n} & = & (\frac{n^{2} + n - 1}{n (n + 2)}) {\hat{a}}_{n} + (\frac{8040 (n + 1)}{17689 n (n + 2)}) {\hat{a}}_{n} + (\frac{2 (n + 1)}{(n + 2) (2 n + 1)}) . \end{matrix}

All codes were written in MATLAB 2018b, for the following different initial values:

{\hat{a}}_{0} = - 0.5

,

{\hat{a}}_{0} = 0.2

, and

{\hat{a}}_{0} = 1.0

. The sequence

{{\hat{a}}_{n}} \to 0

, which is a fixed point of

R_{I, λ}^{M (f, g)} [λ (Y_{I, λ}^{M (f, g)}

- F (C_{I, λ}^{M (f, g)}, Y_{I, λ}^{M (f, g)}) + R_{I, λ}^{M (f, g)}]

. That is, the solution of problem (9),

{\hat{a}}^{*} = 0

, exists. In this regard, a computational table (Table 1) and convergence graph (Figure 1) are shown below.

5. Discussion

We selected the same operators as those in Example 1 and will now conduct a comparison between our proposed Algorithm 1 and the Ishikawa- and-Mann type algorithms. By setting

γ_{n} = 0

, we may compute

{{\hat{a}}_{n}}

and

{{\hat{b}}_{n}}

using the Ishikawa-type schemes outlined below:

\begin{matrix} {\hat{a}}_{n + 1} & = (\frac{n + 1}{2 n + 1}) {\hat{a}}_{n} + (\frac{8040 n}{17689 (2 n + 1)}) {\hat{b}}_{n} + (\frac{1}{2 (2 n + 1)}), \\ {\hat{b}}_{n} & = (\frac{17689 n + 8040}{17689 (n + 1)}) {\hat{a}}_{n} + (\frac{n + 2}{(n + 1) (n + 3)}) . \end{matrix}

(32)

Furthermore, on selecting

β_{n} = γ_{n} = 0

, we may approximate

{{\hat{a}}_{n}}

by the Mann-type iterative scheme:

{\hat{a}}_{n + 1} = (\frac{25729 n + 17689}{17689 (2 n + 1)}) {\hat{a}}_{n} + (\frac{1}{2 (2 n + 1)}) .

(33)

When the stopping criterion

∥ {\hat{a}}_{n + 1} - {\hat{a}}_{n} ∥ \leq 10^{- 10}

is achieved, the iteration methods will cease. Table 2 and Figure 2 compare Algorithm 1 with the Ishikawa-type algorithm (32) and the Mann-type algorithm (33), all initiated with the initial value

{\hat{a}}_{0} = 1 .

6. Conclusions

We examined and demonstrated the existence of a solution to an extended Cayley–Yosida inclusion problem in a real 2-USBS. For the extended Cayley–Yosida inclusion problem, we also developed a three-step iterative method. This algorithm is much more comprehensive than the Mann- and Ishikawa-type iterative schemes, as well as numerous other iterative schemes that have been examined by different writers (see, e.g., [28,29,30]). We proved that our suggested algorithm converges to a unique solution of the problem under some suitable consideration and that this convergence is stable with respect to

R_{I, λ}^{M (f, g)}

. Lastly, we supported our main result with a numerical example.

Author Contributions

Conceptualization, I.A.; Methodology, I.A. and Y.W.; Validation, Y.W. and R.A.; Writing—original draft, I.A.; Funding acquisition, Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant number 12171435).

Data Availability Statement

No new data were created or analyzed in this study. Data is not applicable for this article.

Acknowledgments

The authors are grateful to the anonymous reviewers for their valuable remarks which improved the results and presentation of this paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

Ansari, Q.H.; Yao, J.C. A fixed point theorem and its applications to system of variational inequalities. Bull. Aust. Math. Soc. 1999, 59, 433–442. [Google Scholar] [CrossRef]
Pang, J.S. Asymmetric variational inequality problems over product of sets: Applications and iterative methods. Math. Program. 1985, 31, 206–219. [Google Scholar] [CrossRef]
Yan, W.Y.; Fang, Y.P.; Huang, N.J. A new system of set-valued variational inclusions with H-monotone operators. Math. Inequal. Appl. 2005, 8, 537–546. [Google Scholar] [CrossRef]
Abubakar, J.; Kumam, P.; Deepho, J. Multistep hybrid viscosity method for split monotone variational inclusion and fixed point problems in Hilbert spaces. AIMS Math. 2020, 5, 5969–5992. [Google Scholar] [CrossRef]
Agarwal, R.P.; Verma, R.U. Generalized system of (A, η)-maximal relaxed monotone variational inclusion problems based on generalized hybrid algorithms. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 238–251. [Google Scholar] [CrossRef]
Rockafellar, R.T. On the maximality of sums ofnonlinear monotone operators. Trans. Am. Math. Soc. 1970, 149, 75–88. [Google Scholar] [CrossRef]
Lan, H.Y.; Kim, J.H.; Cho, Y.J. On a new system of nonlinear A-monotone multivalued variational inclusions. J. Math. Anal. Appl. 2007, 327, 481–493. [Google Scholar] [CrossRef]
De, A. Hille-Yosida Theorem and Some Applications. Ph.D. Thesis, Department of Mathematics and its Applications, Central European University, Budapest, Hungary, 2015. [Google Scholar]
Balooee, J. Resolvent algorithms for system of generalized nonlinear variational inclusions and fxed point problems. Afr. Mat. 2014, 25, 1023–1042. [Google Scholar] [CrossRef]
Cao, H.W. Yosida approximation equations technique for system of generalized set-valued variational inclusions. J. Inequal. Appl. 2013, 2013, 455. [Google Scholar] [CrossRef]
Lan, H.Y. Generalized Yosida approximation based on relatively A-maximal m-relaxed monotonicity frameworks. Abs. Appl. Anal. 2013, 2013, 157190. [Google Scholar]
Hassouni, A.; Moudafi, A. A perturbed algorithm for variational inclusions. J. Math. Anal. Appl. 1994, 183, 706–712. [Google Scholar] [CrossRef]
Douglas, J.; Rachford, H.H. On the numerical solution of heat conduction problems in two and three space variables. Trans. Am. Math. Soc. 1956, 82, 421–439. [Google Scholar] [CrossRef]
Peaceman, D.W.; Rachford, H.H., Jr. The numerical solution of parabolic and elliptic differential equations. J. Appl. Ind. Math. 1955, 3, 28–41. [Google Scholar] [CrossRef]
Tseng, P. A modified forward-backward splitting method for maximal monotone mappings. SIAM J. Control Optim. 2000, 38, 431–446. [Google Scholar] [CrossRef]
Gibali, A.; Thong, D.V. Tseng type methods for solving inclusion problems and its applications. Calcolo 2018, 55, 49. [Google Scholar] [CrossRef]
Ahmad, I.; Rahaman, M.; Ahmad, R.; Ali, I. Convergence analysis and stability of perturbed three-step iterative algorithm for generalized mixed ordered quasi-variational inclusion involving XOR operator. Optimization 2020, 69, 821–845. [Google Scholar] [CrossRef]
Glowinski, G.; Tallec, P.L. Augmented Lagrangian and Operator Spliting Methods in Nonlinear Mechanics; SIAM: Philadelphia, PA, USA, 1989. [Google Scholar]
Helmberg, G. Introduction to Spectral Theory in Hilbert Space: The Cayley Transform; North-Holland Series in Applied Mathematics and Mechanics, Applied Mathematics and Mechanics #6; North-Holland Publishing Company: Amsterdam, The Netherlands, 1969. [Google Scholar]
Fang, Y.P.; Huang, N.J. H-accretive operator and resolvent operator technique for variational inclusions in Banach spaces. Appl. Math. Lett. 2004, 17, 647–653. [Google Scholar] [CrossRef]
Agarwal, R.P.; O’Regan, D.; Sahu, D.R. Fixed Point Theory for Lipschitzian-Type Mappings with Applications; Springer: New York, NY, USA, 2009; Volume 6. [Google Scholar]
Xu, H.K. Inequalities in Banach spaces with applications. Nonlinear Anal. Theory Methods Appl. 1991, 16, 1127–1138. [Google Scholar] [CrossRef]
Kazmi, K.R.; Khan, F.A.; Shahzad, M. A system of generalized variational inclusions involving generalized H(·, ·)-accretive mapping in real q-uniformly smooth Banach spaces. Appl. Math. Comput. 2011, 217, 9679–9688. [Google Scholar] [CrossRef]
Osilike, M.O. Stability results for the Ishikawa fixed point iteration procedure. Ind. J. Pure Appl. Math. 1995, 26, 937–945. [Google Scholar]
Weng, X.L. Fixed point iteration for local srictly pseudocontractive mapping. Proc.Am.Math.Soc. 1991, 113, 727–731. [Google Scholar] [CrossRef]
Ahmad, R.; Ali, I.; Rahaman, M.; Ishtyak, M.; Yao, J.C. Cayley inclusion problem with its corresponding generalized resolvent equation problem in uniformly smooth Banach spaces. Appl. Anal. 2022, 101, 1354–1368. [Google Scholar] [CrossRef]
Rajpoot, A.K.; Ishtyak, M.; Ahmad, R.; Wang, Y.; Yao, J.C. Convergence analysis for Yosida variational inclusion problem with its corresponding Yosida resolvent equation problem through inertial extrapolation scheme. Mathematics 2023, 11, 763. [Google Scholar] [CrossRef]
Akram, M.; Dilshad, M.; Rajpoot, A.K.; Babu, F.; Ahmad, R.; Yao, J.C. Modified Iterative Schemes for a Fixed-Point Problem and a Split Variational Inclusion Problem. Mathematics 2022, 10, 2098. [Google Scholar] [CrossRef]
Zhao, X.; Sahu, D.R.; Wen, C.F. Iterative methods for system of variational inclusions involving accretive operators and applications. Fixed Point Theory 2018, 19, 801–822. [Google Scholar] [CrossRef]
Iqbal, J.; Ali, I.; Arora, P.K.; Mir, W.A. Set-valued variational inclusion problem with fuzzy mappings involving XOR-operation. AIMS Math. 2021, 6, 3288–3304. [Google Scholar] [CrossRef]

Figure 1. Convergence of

{{\hat{a}}_{n}}

starting with values

a_{0} =

−

0.5,

a_{0} = 0.2

, and

a_{0} = 1

.

Figure 1. Convergence of

{{\hat{a}}_{n}}

starting with values

a_{0} =

−

0.5,

a_{0} = 0.2

, and

a_{0} = 1

.

Figure 2. Convergence of

{{\hat{a}}_{n}}

with starting value

{\hat{a}}_{0} = 1

.

Figure 2. Convergence of

{{\hat{a}}_{n}}

with starting value

{\hat{a}}_{0} = 1

.

Table 1. The computational table of

{{\hat{a}}_{n}}

starting with

{\hat{a}}_{0} = - 0.5,

{\hat{a}}_{0} = 0.2

, and

{\hat{a}}_{0} = 1

.

Table 1. The computational table of

{{\hat{a}}_{n}}

starting with

{\hat{a}}_{0} = - 0.5,

{\hat{a}}_{0} = 0.2

, and

{\hat{a}}_{0} = 1

.

No. of Iterations	For ${\hat{a}}_{0} = - 0.55$ ${\hat{a}}_{n}$	For ${\hat{a}}_{0} = 0.2$ . ${\hat{a}}_{n}$	For ${\hat{a}}_{0} = 1$ ${\hat{a}}_{n}$
n = 1	−0.5	0.2	1
n = 2	0.398825035	0.961323961	1.139101739
n = 3	0.269533992	0.498765092	0.557440987
n = 4	0.132297719	0.205800419	0.221155261
n = 5	0.050226267	0.054881135	0.055909079
n = 10	0.022189534	0.022190736	0.022191014
n = 15	0.01479883	0.01479883	0.01479883
n = 20	0.011134277	0.011134277	0.011134277
n = 25	0.008932186	0.008932186	0.008932186
n = 30	0.007222087	0.007222087	0.007222087
n = 35	0.006062878	0.006062878	0.006062878
n = 40	0.005224884	0.005224884	0.005224884
n = 45	0.004499671	0.004412204	0.004412204
n = 55	0.003818503	0.003755349	0.003755349
n = 70	0.003268839	0.003268839	0.003268839

Table 2. The computational table of

{{\hat{a}}_{n}}

with starting value

{\hat{a}}_{0} = 1

.

Table 2. The computational table of

{{\hat{a}}_{n}}

with starting value

{\hat{a}}_{0} = 1

.

No. of Iterations	Algorithm 1 ${\hat{a}}_{n}$	Algorithm (32) ${\hat{a}}_{n}$	Algorithm (33) ${\hat{a}}_{n}$
n = 1	1	1	1
n = 2	0.513406206	0.732158368	0.844808005
n = 3	0.219329442	0.298731566	0.668676662
n = 4	0.081548826	0.197299587	0.515703683
n = 5	0.039620365	0.107807382	0.392252818
n = 10	0.012456521	0.031819764	0.091455223
n = 15	0.007980668	0.020022186	0.020119998
n = 20	0.004324994	0.005890604	0.014881626
n = 25	0.000918114	0.004672308	0.011861738
n = 30	0.000193369	0.00387296	0.009864815
n = 35	4.05 $\times 10^{- 5}$	0.00844483	0.003307718
n = 40	8.45 $\times 10^{- 6}$	0.007382855	0.002886711
n = 45	1.76 $\times 10^{- 6}$	0.006558468	0.002560908
n = 50	3.65 $\times 10^{- 7}$	3.65 $\times 10^{- 2}$	3.65 $\times 10^{- 5}$

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Ali, I.; Wang, Y.; Ahmad, R. Three-Step Iterative Algorithm for the Extended Cayley–Yosida Inclusion Problem in 2-Uniformly Smooth Banach Spaces: Convergence and Stability Analysis. Mathematics 2024, 12, 1977. https://doi.org/10.3390/math12131977

AMA Style

Ali I, Wang Y, Ahmad R. Three-Step Iterative Algorithm for the Extended Cayley–Yosida Inclusion Problem in 2-Uniformly Smooth Banach Spaces: Convergence and Stability Analysis. Mathematics. 2024; 12(13):1977. https://doi.org/10.3390/math12131977

Chicago/Turabian Style

Ali, Imran, Yuanheng Wang, and Rais Ahmad. 2024. "Three-Step Iterative Algorithm for the Extended Cayley–Yosida Inclusion Problem in 2-Uniformly Smooth Banach Spaces: Convergence and Stability Analysis" Mathematics 12, no. 13: 1977. https://doi.org/10.3390/math12131977

APA Style

Ali, I., Wang, Y., & Ahmad, R. (2024). Three-Step Iterative Algorithm for the Extended Cayley–Yosida Inclusion Problem in 2-Uniformly Smooth Banach Spaces: Convergence and Stability Analysis. Mathematics, 12(13), 1977. https://doi.org/10.3390/math12131977

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Three-Step Iterative Algorithm for the Extended Cayley–Yosida Inclusion Problem in 2-Uniformly Smooth Banach Spaces: Convergence and Stability Analysis

Abstract

1. Introduction

2. Background and Problem Formulation

Problem Formulation

3. Fixed-Point Equation and Three-Step Iterative Algorithm

3.1. Existence Result

3.2. Three-Step Iterative Algorithm

4. Numerical Example

5. Discussion

6. Conclusions

Author Contributions

Funding

Data Availability Statement

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI