Next Article in Journal
Research on Corporate Indebtedness Determinants: A Case Study of Visegrad Group Countries
Next Article in Special Issue
Preface to the Special Issue “State-of-the-Art Mathematical Applications in Europe”
Previous Article in Journal
A Review of Combinatorial Optimization Problems in Reverse Logistics and Remanufacturing for End-of-Life Products
Previous Article in Special Issue
Two Families of Continuous Probability Distributions Generated by the Discrete Lindley Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Existence and Stability of Solutions for a Class of Fractional Riemann–Liouville Initial Value Problems

by
Luís P. Castro
*,† and
Anabela S. Silva
CIDMA—Center for Research and Development in Mathematics and Applications, University of Aveiro, 3810-193 Aveiro, Portugal
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2023, 11(2), 297; https://doi.org/10.3390/math11020297
Submission received: 15 December 2022 / Revised: 3 January 2023 / Accepted: 4 January 2023 / Published: 6 January 2023
(This article belongs to the Special Issue State-of-the-Art Mathematical Applications in Europe)

Abstract

:
This article deals with a class of nonlinear fractional differential equations, with initial conditions, involving the Riemann–Liouville fractional derivative of order α ( 1 , 2 ) . The main objectives are to obtain conditions for the existence and uniqueness of solutions (within appropriate spaces), and to analyze the stabilities of Ulam–Hyers and Ulam–Hyers–Rassias types. In fact, different conditions for the existence and uniqueness of solutions are obtained based on the analysis of an associated class of fractional integral equations and distinct fixed-point arguments. Additionally, using a Bielecki-type metric and some additional contractive arguments, conditions are also obtained to guarantee Ulam–Hyers and Ulam–Hyers–Rassias stabilities for the problems under analysis. Examples are also included to illustrate the theory.

1. Introduction

Given the importance that fractional derivatives and integrals [1,2,3,4,5,6,7] have shown to have in the optimization and improvement of mathematical models of real events or even of those associated with other areas of knowledge (namely through making these models more accurate when compared to what they effectively model), we have recently witnessed a large development in the mathematical analysis of classes of fractional order differential and integral equations.
In this context, it is essential to know about the possible existence of several solutions to the problems in question, possible sufficient conditions to obtain a unique solution and even conditions that eventually guarantee distinct forms of stability of the solutions (this being a crucial aspect, in particular, for the study of approximate solutions to the problems in analysis). The most used techniques in these problems involve the consideration and identification of operators that (in a sense) represent the problem (in some “equivalent” way) and usually involve different principles of contraction, as well as different estimates, usually framed, or dependent, on norms (or metrics), within the spaces framework most suited to the problems under study.
For this type of problem, the analysis of their eventual stability is also a study of significant importance. Namely, through the Ulam–Hyers and Ulam–Hyers–Rassias stabilities [1,8,9,10,11,12,13,14,15,16,17] which, with their specific characteristics, make it possible to identify forms of a slight disturbance in the system (that defines the problem) does not have a too disturbing effect on that system.
Having this general framework in mind, we would like to start by emphasizing that in [18], Chai studied the existence of solutions to the boundary value problem
C D 0 + α x ( t ) + r C D 0 + α 1 x ( t ) = f ( t , x ( t ) ) , t ( 0 , 1 ) , x ( 0 ) = x ( 1 ) , x ( ξ ) = η , ξ ( 0 , 1 ) ,
where C D a + α and C D a + α 1 denote the standard Caputo derivatives of order α and α 1 , respectively, in this case with 1 < α 2 , and r 0 . Additionally, more recently, Xu et al. [19] considered the existence of solutions and the Ulam–Hyers stability for the fractional boundary value problem
λ D 0 + α x ( t ) + D 0 + β x ( t ) = f ( t , x ( t ) ) , t ( 0 , T ) , x ( 0 ) = 0 , μ D 0 + γ 1 x ( T ) + I 0 + γ 2 x ( η ) = γ 3 ,
where D 0 + ϑ denotes the Riemann–Liouville fractional derivative operator of order ϑ , 1 < α 2 , 1 β < α , 0 < λ 1 , 0 < μ 1 , 0 γ 1 α β , γ 2 0 , I 0 + γ 2 denotes the Riemann–Liouville fractional integral operator of order γ 2 , and 0 < η < T . Moreover, in [20], Ahmad et al. investigated the existence of solutions and the Ulam–Hyers stability for a fractional initial value problem given by
( C D a + α x ( t ) + λ 1 C D a + α 1 x ( t ) + λ 2 C D a + α 2 x ( t ) = f ( t , x ( t ) ) , t [ a , T ] , x ( k ) ( a ) = b k , k = 0 , 1 , 2 ,
where C D a + α is again the Caputo fractional derivative of order α ( 2 , 3 ) , and λ 1 and λ 2 are nonzero constants. In [21], Alvan et al. investigated the existence of solutions for the fractional boundary value problem
C D 0 + α x ( t ) + 2 r C D 0 + α 1 x ( t ) + r 2 C D 0 + α 2 x ( t ) = f ( t , x ( t ) , D 0 + σ 1 ) , r > 0 , t ( 0 , 1 ) , x ( 0 ) = x ( 1 ) , x ( 0 ) = x ( 1 ) , x ( ξ ) + r x ( ξ ) = η , ξ ( 0 , 1 ) ,
where 2 α < 3 and η is a positive real number. Bilgici and Şan [22] considered the existence and uniqueness of solutions to the problem
λ D 0 + α x ( t ) = f ( t , x ( t ) , D 0 + α 1 x ( t ) ) , t > 0 , x ( 0 ) = 0 , D 0 + α 1 x ( t ) | t = 0 = b ,
where α ( 1 , 2 ) and b 0 .
Motivated by the analysis and the results already achieved for the above-mentioned problems (included in the works [18,19,20,21,22]), we investigate in this paper the stabilities of Ulam–Hyers and Ulam–Hyers–Rassias types [1,8,9,10,11,14,16], and the existence and uniqueness of solutions to the following initial value problem of fractional order (IVPFO)
D a + α x ( t ) + λ ( D a + α 1 x ) ( t ) = f ( t , x ( t ) ) , t [ a , b ] , x ( a ) = x ( a ) = 0 ,
where 1 < α < 2 , λ is a nonzero constant, a , b R (with a < b ) and f : [ a , b ] × R R is a continuous function. Thus, this problem can also be viewed as a class of problems depending on the parameter λ , and with the form of a single-point boundary problem “a” of a two-term fractional differential equation.
The remaining part of the work is organized as follows: Section 2 contains the necessary definitions and the fundamental tools that are used in the sections that follow; in Section 3, we derive different conditions for the existence and uniqueness of solutions for the IVPFO (1); in Section 4, we discuss the Ulam–Hyers and the Ulam–Hyers–Rassias stabilities and obtain conditions for their existence. Finally, some examples are included to describe the obtained results in a more concrete way.

2. Preliminaries and Background Material

We start this section by presenting the known basic definitions of the main objects that we will use.
Definition 1.
The Riemann–Liouville fractional integral of order α R + of a function x (on [ a , b ] ) is defined by
I a + α x ( t ) = 1 Γ ( α ) a t ( t s ) α 1 x ( s ) d s ( a t b )
provided the right-hand side is pointwise defined and where Γ denotes the Euler Gamma function (given by Γ ( α ) = 0 t α 1 e t d t , α > 0 ).
Definition 2.
The Riemann–Liouville fractional derivative of order α > 0 of a function x (on [ a , b ] ) is defined by
D a + α x ( t ) = 1 Γ ( n α ) d n d t n a t ( t s ) n α 1 x ( s ) d s ,
with n = [ α ] + 1 .
In what follows, we denote by L 1 ( [ a , b ] ) the Banach space of Lebesgue integrable functions from [ a , b ] into R with the norm x L 1 = a b | x ( t ) | d t and by C ( [ a , b ] ) the Banach space of all continuous functions g : [ a , b ] R endowed with the norm g = sup t [ a , b ] | g ( t ) | .
Lemma 1
([3]). Assume that x C ( [ a , b ] ) L 1 ( [ a , b ] ) with a fractional derivative of order α > 0 . Then
D a + α I a + α x ( t ) = x ( t )
and
I a + α D a + α x ( t ) = x ( t ) + c 1 ( t a ) α 1 + c 2 ( t a ) α 2 + + c n ( t a ) α n ,
for some c i R , i = 1 , 2 , , n , where n is the smallest integer greater than or equal to α.
For the reader’s convenience, let us recall some classic principles of contraction and inequalities that we will use later.
Theorem 1
(Banach contraction principle). Let ( X , d ) be a generalized complete metric space, and consider a mapping T : X X which is a strictly contractive operator, that is
d ( T x , T y ) L d ( x , y ) , x , y X ,
for some constant 0 L < 1 . Then:
(a) 
the mapping T has a unique fixed point x * = T x * ;
(b) 
the fixed point x * is globally attractive, namely, for any starting point x X , the following identity holds:
lim n T n x = x * ;
(c) 
we have the following inequalities:
d ( T n x , x * ) L n d ( x , x * ) , n 0 , x X ; d ( T n x , x * ) 1 1 L d ( T n x , T n + 1 x ) , n 0 , x X ; d ( x , x * ) 1 1 L d ( x , T x ) , x X .
Theorem 2
(Schauder’s fixed point theorem). If   Ω is a closed, bounded, convex subset of a Banach space X and the mapping T : Ω Ω is completely continuous, then T has a fixed point in Ω .
Keeping in mind some parts of the proofs of the next results, let us recall an important integral inequality that we will actually use later.
Theorem 3
([23], [Theorem 11.2]). Let u ( t ) , b ( t ) , σ ( t ) and k ( t , s ) be nonnegative continuous functions for a s t b and suppose that
u ( t ) c 1 + σ ( t ) c 2 + a t b ( s ) u ( s ) d s + a t a s k ( s , τ ) u ( τ ) d τ d s ) ,
for t [ a , b ] , where c 1 , c 2 0 are constants. Then,
u ( t ) c 2 e a t B ( s ) σ ( s ) d s + a t c 1 B ( s ) e s t B ( τ ) σ ( τ ) d τ d s ,
where B ( s ) = b ( s ) + a s k ( s , τ ) d τ .
We denote by C 2 ( [ a , b ] ) the space of functions x which are 2-times continuously differentiable on [ a , b ] endowed with the norm
x C 2 = k = 0 2 sup t [ a , b ] | x ( k ) ( t ) | .
It is well-known that ( C 2 ( [ a , b ] ) , · C 2 ) is a Banach space.
In our next analysis of the existence and uniqueness of solutions for the IVPFO (1), we will make use of the following auxiliary property (which may be considered as a very natural and expectable property; cf., e.g., [24]).
Lemma 2
(See also [24]). Let α ( 1 , 2 ) and x C 2 ( [ a , b ] ) with x ( a ) = x ( a ) = 0 . Then D a + α x C ( [ a , b ] ) and
( D a + α x ) ( t ) = 1 Γ ( 2 α ) a t ( t s ) 1 α x ( s ) d s .
Moreover,
( D a + α x ) ( t ) = ( D a + α 1 x ) ( t ) .
Proof. 
For the reader’s convenience, we have chosen to include here a proof of this lemma. Within the stated conditions, we simply have to use integration by parts to obtain
a t ( t s ) 1 α x ( s ) d s = 1 2 α a t ( t s ) 2 α x ( s ) d s , a t ( t s ) 2 α x ( s ) d s = 1 3 α a t ( t s ) 3 α x ( s ) d s .
And so, it follows
( D a + α x ) ( t ) = 1 Γ ( 2 α ) d d t 2 a t ( t s ) 1 α x ( s ) d s = 1 Γ ( 3 α ) d d t 2 a t ( t s ) 2 α x ( s ) d s = 1 Γ ( 4 α ) d d t 2 a t ( t s ) 3 α x ( s ) d s = 1 Γ ( 2 α ) a t ( t s ) 1 α x ( s ) d s .
Since under the present conditions a t ( t s ) 1 α x ( s ) d s is continuous on [ a , b ] , we conclude that D a + α x is continuous on [ a , b ] .
Moreover,
( D a + α x ) ( t ) = 1 Γ ( 2 α ) d d t 2 a t ( t s ) 1 α x ( s ) d s = 1 α Γ ( 2 α ) d d t a t ( t s ) α x ( s ) d s .
Integrating by parts, and using the circumstance that x ( a ) = 0 , we obtain
( D a + α x ) ( t ) = 1 Γ ( 2 α ) d d t a t ( t s ) 1 α x ( s ) d s = ( D a + α 1 x ) ( t ) ,
which concludes the proof. □
Remark 1.
Proceeding in a similar way as in the previous lemma, for α ( 1 , 2 ) , x C 2 ( [ a , b ] ) and x ( a ) = x ( a ) = 0 , it follows that D a + α 1 x C ( [ a , b ] ) and
( D a + α 1 x ) ( t ) = 1 Γ ( 3 α ) a t ( t s ) 2 α x ( s ) d s .

3. Different Conditions for the Existence and Uniqueness of Solutions

In the present section, we will analyse conditions to ensure the existence of solutions to the IVPFO (1) and also conditions to guarantee the uniqueness of the solution. In view of this, let us first start to “translate” the IVPFO (1) through a fractional integral equation.
Proposition 1.
As before, let α ( 1 , 2 ) , f : [ a , b ] × R R be a continuous function and λ 0 . A function x C 2 ( [ a , b ] ) is a solution of the IVPFO (1) if and only if x satisfies the integral equation
x ( t ) = e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , x ( s ) ) d s d u .
Proof. 
Let x C 2 ( [ a , b ] ) be the solution of IVPFO (1). By Lemma 2, we have that D a + α x , D a + α 1 x C ( [ a , b ] ) and D a + α x = D a + α 1 x . Thus, we can rewrite our main equation in (1),
( D a + α x ) ( t ) + λ ( D a + α 1 x ) ( t ) = f ( t , x ( t ) ) ,
in the form
( D a + α 1 x ) ( t ) + λ ( D a + α 1 x ) ( t ) = f ( t , x ( t ) ) .
In view of Lemma 1, one has
( I a + α 1 D a + α 1 x ) ( t ) = x ( t ) + c 1 ( t a ) α 2 , ( I a + α 1 D a + α 1 x ) ( t ) = x ( t ) + d 1 ( t a ) α 2 , t [ a , b ] .
Thus, applying I a + α 1 to both members of Equation (4), we obtain
x ( t ) + λ x ( t ) + ( λ c 1 + d 1 ) ( t a ) α 2 = [ I a + α 1 f ( · , x ( · ) ) ] ( t ) .
Since x ( a ) = x ( a ) = 0 , we conclude that
λ c 1 + d 1 = 0 ,
and so it follows
x ( t ) + λ x ( t ) = [ I a + α 1 f ( · , x ( · ) ) ] ( t ) .
Let y ( t ) = e λ t x ( t ) . One has that
x ( t ) = λ e λ t y ( t ) + e λ t y ( t ) .
Substituting the last two identities in (6), we obtain
y ( t ) = e λ t [ I a + α 1 f ( · , e λ · y ( · ) ) ] ( t ) .
Since x C 2 ( [ a , b ] ) , we have that y C 1 ( [ a , b ] ) . Moreover, I a + α 1 f is a continuously differentiable function. Thus, integrating Equation (7) from a to t, we obtain
y ( t ) = y ( a ) + 1 Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , e λ s y ( s ) ) d s d u .
Taking into account that y ( t ) = e λ t x ( t ) , it follows that
x ( t ) = e λ ( t a ) x ( a ) + e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , x ( s ) ) d s d u ,
and using the initial conditions, we conclude that
x ( t ) = e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , x ( s ) ) d s d u .
Conversely, assume that x is given by (3), and thus
e λ t x ( t ) = 1 Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , x ( s ) ) d s d u .
It is clear that x ( a ) = 0 and since x is continuously differentiable on [ a , b ] , differentiating both sides of (8), we get
e λ t x ( t ) + λ e λ t x ( t ) = e λ t Γ ( α 1 ) a t ( t s ) α 2 f ( s , x ( s ) ) d s ,
which is equivalent to
x ( t ) + λ x ( t ) = 1 Γ ( α 1 ) a t ( t s ) α 2 f ( s , x ( s ) ) d s .
Thus x ( a ) = 0 and since x C 2 ( [ a , b ] ) , accordingly to Lemma 2, we have that D a + α x and D a + α 1 x exist. Applying D a + α 1 to both sides of Equation (9), using Lemma 1 and (2), we also obtain
( D a + α x ) ( t ) + λ ( D a + α 1 x ) ( t ) = f ( t , x ( t ) ) ,
which completes the proof. □
Having in mind Proposition 1, we realize that studying the solutions of IVPFO (1) is the same as studying the solutions of
x = T x ,
where T is the fractional integral operator given by
( T x ) ( t ) = e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , x ( s ) ) d s d u ,
for x C 2 ( [ a , b ] ) and λ R { 0 } .
Remark 2.
Another way to discover an integral form of x ( t ) is to consider the integral equation
x ( t ) = λ a t x ( s ) d s + 1 Γ ( α ) a t ( t s ) α 1 f ( s , x ( s ) ) d s .
In fact, applying I a + α to both members of equation ( D a + α x ) ( t ) + λ ( D a + α 1 x ) ( t ) = f ( t , x ( t ) ) , and using Lemma 1, we obtain
x ( t ) + a 1 ( t a ) α 1 + a 2 ( t a ) α 2 + λ a t x ( s ) + b 1 ( s a ) α 2 d s = [ I a + α f ( · , x ( · ) ) ] ( t )
( a 1 , a 2 , b 1 R ), which is equivalent to
x ( t ) = a 1 + λ b 1 α 1 ( t a ) α 1 a 2 ( t a ) α 2 λ a t x ( s ) d s + [ I a + α f ( · , x ( · ) ) ] ( t ) .
Since x ( a ) = 0 , it follows that a 2 = 0 . Observing that
x ( t ) = ( α 1 ) a 1 + λ b 1 ( t a ) α 2 λ x ( t ) + [ I a + α 1 f ( · , x ( · ) ) ] ( t ) ,
and using the initial condition x ( a ) = 0 , we also conclude that a 1 + λ b 1 α 1 = 0 , and thus, Equation (11) is obtained.
Let us fix the following notation
k = ( b a ) α 1 λ Γ ( α ) 1 ( 1 λ + λ 2 ) e λ ( b a ) , k + = ( b a ) α 1 λ Γ ( α ) 1 + 2 λ + 2 λ 2 ( 1 + λ + λ 2 ) e λ ( b a ) ,
and
K = K ( λ ) : = k , λ < 0 k + , λ > 0 .
Theorem 4.
If f : [ a , b ] × R R is continuously differentiable, then the IVPFO (1) has at least one solution in C 2 ( [ a , b ] ) .
Proof. 
We will use the Schauder fixed point theorem for the fractional integral operator T, defined in (10). The continuity of T x follows from the continuity of f. Moreover, we have that
( T x ) ( t ) = λ e λ t Γ ( α 1 ) a t a u e λ u ( u s ) α 2 f ( s , x ( s ) ) d s d u + 1 Γ ( α 1 ) a t ( t s ) α 2 f ( s , x ( s ) ) d s = λ e λ t Γ ( α 1 ) a t a u e λ u ( u s ) α 2 f ( s , x ( s ) ) d s d u + 1 Γ ( α ) ( t a ) α 1 f ( a , 0 ) + a t ( t s ) α 1 f ( s , x ( s ) ) d s ,
and
( T x ) ( t ) = λ 2 e λ t Γ ( α 1 ) a t a u e λ u ( u s ) α 2 f ( s , x ( s ) ) d s d u λ Γ ( α 1 ) a t ( t s ) α 2 f ( s , x ( s ) ) d s + 1 Γ ( α 1 ) ( t a ) α 2 f ( a , 0 ) + a t ( t s ) α 2 f ( s , x ( s ) ) d s .
Since f is continuously differentiable, there exist positive constants A and B such that | f ( t , x ( t ) | A and | f ( t , x ( t ) ) | B , t [ a , b ] . Define Ω = { x C 2 ( [ a , b ] ) : x C 2 R } with R being a positive real number satisfying
R K A + ( b a ) α 2 Γ ( α 1 ) f ( a , 0 ) + ( b a ) α 1 Γ ( α ) B .
It is clear that Ω is a closed, bounded and convex subset of C 2 ( [ a , b ] ) . Moreover, we have that
| ( T x ) ( t ) | = e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , x ( s ) ) d s d u e λ t Γ ( α 1 ) a t e λ u a u ( u s ) α 2 | f ( s , x ( s ) ) | d s d u e λ t A Γ ( α 1 ) a t e λ u a u ( u s ) α 2 d s d u e λ t Γ ( α ) A ( b a ) α 1 a t e λ u d u = ( b a ) α 1 λ Γ ( α ) ( 1 e λ ( t a ) ) A ,
| ( T x ) ( t ) | λ e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , x ( s ) ) d s d u + 1 Γ ( α 1 ) a t ( t s ) α 2 f ( s , x ( s ) ) d s ( b a ) α 1 λ Γ ( α ) | λ | ( 1 e λ ( t a ) ) + λ A ,
and
| ( T x ) ( t ) | λ 2 e λ t Γ ( α 1 ) a t a u e λ u ( u s ) α 2 f ( s , x ( s ) ) d s d u + λ Γ ( α 1 ) a t ( t s ) α 2 f ( s , x ( s ) ) d s + 1 Γ ( α 1 ) ( t a ) α 2 f ( a , 0 ) + a t ( t s ) α 2 f ( s , x ( s ) ) d s ( b a ) α 1 λ Γ ( α ) λ 2 ( 1 e λ ( t a ) ) λ | λ | A + ( b a ) α 2 Γ ( α 1 ) f ( a , 0 ) + ( b a ) α 1 Γ ( α ) B .
Thus, we have that
T x C 2 sup t [ a , b ] ( b a ) α 1 λ Γ ( α ) ( 1 e λ ( t a ) ) A + sup t [ a , b ] ( b a ) α 1 λ Γ ( α ) | λ | ( 1 e λ ( t a ) ) + λ A + sup t [ a , b ] ( b a ) α 1 λ Γ ( α ) λ 2 ( 1 e λ ( t a ) ) + λ | λ | A + ( b a ) α 2 f ( a , 0 ) Γ ( α 1 ) + ( b a ) α 1 Γ ( α ) B .
Thus, if λ < 0 , we have that
T x C 2 ( b a ) α 1 λ Γ ( α ) 1 λ + λ 2 ( 1 λ + λ 2 ) e λ ( b a ) + λ λ 2 A + ( b a ) α 2 Γ ( α 1 ) f ( a , 0 ) + ( b a ) α 1 Γ ( α ) B = k A + ( b a ) α 2 Γ ( α 1 ) f ( a , 0 ) + ( b a ) α 1 Γ ( α ) B R ,
and if λ > 0 , we have
T x C 2 ( b a ) α 1 λ Γ ( α ) 1 + λ + λ 2 ( 1 + λ + λ 2 ) e λ ( b a ) + λ + λ 2 A + ( b a ) α 2 Γ ( α 1 ) f ( a , 0 ) + ( b a ) α 1 Γ ( α ) B = k + A + ( b a ) α 2 Γ ( α 1 ) f ( a , 0 ) + ( b a ) α 1 Γ ( α ) B R .
Consequently, we conclude that T is a bounded operator on Ω C 2 ( [ a , b ] ) .
Let us prove that operator T : Ω Ω is completely continuous. For t 1 , t 2 [ a , b ] , t 1 < t 2 , one has
| ( T x ) ( t 2 ) ( T x ) ( t 1 ) | = e λ t 2 Γ ( α 1 ) a t 2 a u ( u s ) α 2 e λ u f ( s , x ( s ) ) d s d u e λ t 1 Γ ( α 1 ) a t 1 a u ( u s ) α 2 e λ u f ( s , x ( s ) ) d s d u e λ t 2 Γ ( α 1 ) t 1 t 2 a u | ( u s ) α 2 e λ u f ( s , x ( s ) ) | d s d u + | e λ t 2 e λ t 1 | Γ ( α 1 ) a t 1 a u | ( u s ) α 2 e λ u f ( s , x ( s ) ) | d s d u ,
which tends to zero as t 2 t 1 (independently of x and λ ). In the same way, we get
| ( T x ) ( t 2 ) ( T x ) ( t 1 ) | = | λ | | ( T x ) ( t 2 ) ( T x ) ( t 1 ) | + 1 Γ ( α 1 ) t 1 t 2 ( t 2 s ) α 2 | f ( s , x ( s ) ) | d s + + 1 Γ ( α 1 ) a t 1 [ ( t 2 s ) α 2 ( t 1 s ) α 2 ] | f ( s , x ( s ) ) | d s ,
which tends to zero as t 2 t 1 . Finally, we observe that
| ( T x ) ( t 2 ) ( T x ) ( t 1 ) | = | λ | | ( T x ) ( t 2 ) ( T x ) ( t 1 ) | + 1 Γ ( α 1 ) t 1 t 2 ( t 2 s ) α 2 | f ( s , x ( s ) ) | d s + + 1 Γ ( α 1 ) a t 1 [ ( t 2 s ) α 2 ( t 1 s ) α 2 ] | f ( s , x ( s ) ) | d s + ( t 2 a ) α 2 ( t 1 a ) α 2 Γ ( α 1 ) f ( a , 0 )
tends to zero as t 2 t 1 . Thus, we conclude that T Ω is equicontinuous. Following Arzelà-Ascoli Theorem, we obtain that T is completely continuous. Applying Schauder’s fixed point theorem (cf. Theorem 2), we conclude that the operator T has at least one fixed point, which means that the IVPFO (1) has at least one solution and the proof is completed. □
We will now exhibit other conditions under which, besides the existence of solutions, we will also guarantee the uniqueness of the solution to the IVPFO (1).
Theorem 5.
Let f : [ a , b ] × R R be a continuously differentiable function and suppose that there are L 1 and L 2 0 such that, for t [ a , b ] ,
| f ( t , x ( t ) ) f ( t , y ( t ) ) | L 1 | x ( t ) y ( t ) | ,
| f ( t , x ( t ) ) f ( t , y ( t ) ) | L 2 | x ( t ) y ( t ) | + | x ( t ) y ( t ) | .
If
K L 1 + L 2 ( b a ) α 1 Γ ( α ) < 1 ,
then the problem (1) has a unique solution on C 2 ( [ a , b ] ) .
Proof. 
Since f is a continuously differentiable function, according to Theorem 4, the IVPFO (1) admits at least one solution. Let us assume that conditions (13)–(14) hold. Thus, we can obtain that, for x , y C 2 ( [ a , b ] ) ,
| ( T x ) ( t ) ( T y ) ( t ) | e λ t Γ ( α 1 ) a t e λ u a u ( u s ) α 2 | f ( s , x ( s ) ) f ( s , y ( s ) ) | d s d u L 1 e λ t Γ ( α 1 ) a t e λ u a u ( u s ) α 2 | x ( s ) y ( s ) | d s d u L 1 x y C 2 ( b a ) α 1 ( 1 e λ ( t a ) ) λ Γ ( α ) ,
| ( T x ) ( t ) ( T y ) ( t ) | 1 Γ ( α 1 ) a t ( t s ) α 2 | f ( s , x ( s ) ) f ( s , y ( s ) ) | d s + | λ | e λ t Γ ( α 1 ) a t e λ u a u ( u s ) α 2 | f ( s , x ( s ) ) f ( s , y ( s ) ) | d s d u L 1 Γ ( α 1 ) a t ( t s ) α 2 | x ( s ) y ( s ) | d s + L 1 | λ | e λ t Γ ( α 1 ) a t e λ u a u ( u s ) α 2 | x ( s ) y ( s ) | d s d u L 1 x y C 2 ( b a ) α 1 λ + | λ | ( 1 e λ ( t a ) ) λ Γ ( α ) ,
and
| ( T x ) ( t ) ( T y ) ( t ) | λ 2 e λ t Γ ( α 1 ) a t a u e λ u ( u s ) α 2 | f ( s , x ( s ) ) f ( s , y ( s ) ) | d s d u + | λ | Γ ( α 1 ) a t ( t s ) α 2 | f ( s , x ( s ) ) f ( s , y ( s ) ) | d s + 1 Γ ( α 1 ) a t ( t s ) α 2 | f ( s , x ( s ) ) f ( s , y ( s ) ) | d s L 1 λ 2 e λ t Γ ( α 1 ) a t a u e λ u ( u s ) α 2 | x ( s ) y ( s ) | d s d u + | λ | L 1 Γ ( α 1 ) a t ( t s ) α 2 | x ( s ) y ( s ) | d s + L 2 Γ ( α 1 ) a t ( t s ) α 2 | x ( s ) y ( s ) | + | x ( s ) y ( s ) | d s L 1 x y C 2 ( b a ) α 1 λ 2 ( 1 e λ ( t a ) ) + λ | λ | λ Γ ( α ) + L 2 x y C 2 ( b a ) α 1 Γ ( α ) .
Thus, we conclude that, for λ > 0
T x T y C 2 = sup t [ a , b ] | ( T x ) ( t ) ( T y ) ( t ) | + sup t [ a , b ] | ( T x ) ( t ) ( T y ) ( t ) | + sup t [ a , b ] | ( T x ) ( t ) ( T y ) ( t ) | x y C 2 L 1 ( b a ) α 1 [ 1 + 2 λ + 2 λ 2 ( 1 + λ + λ 2 ) e λ ( b a ) ] λ Γ ( α ) + L 2 ( b a ) α 1 Γ ( α ) = k + L 1 + L 2 ( b a ) α 1 Γ ( α ) x y C 2 ,
and for λ < 0 ,
T x T y C 2 x y C 2 L 1 ( b a ) α 1 [ 1 e λ ( b a ) + λ e λ ( b a ) λ 2 e λ ( b a ) ] λ Γ ( α ) + L 2 ( b a ) α 1 Γ ( α ) = k L 1 + L 2 ( b a ) α 1 Γ ( α ) x y C 2 .
Since K L 1 + L 2 ( b a ) α 1 Γ ( α ) < 1 , we have that T is a contractive operator. Thus, by Banach contraction principle (cf. Theorem 1), we conclude that T has a unique fixed point, which from Proposition 1 means that the IVPFO (1) has a unique solution on C 2 ( [ a , b ] ) . □

4. Ulam–Hyers and Ulam–Hyers–Rassias Stabilities

In this section, we analyse the Ulam–Hyers and the Ulam–Hyers–Rassias stabilities of the above class of problems. In fact, since from Proposition 1 we have a new Equation (3) to describe the IVPFO (1) equivalently, we may choose to discuss the stabilities of (1) or (3). Thus, in here, we choose to exhibit, in detail, conditions for the Ulam–Hyers stability of (1) and the Ulam–Hyers–Rassias stability of (3). To this purpose, let us first point out what are the definitions of such stabilities in each of those cases.
Definition 3.
The IVPFO (1) is Ulam–Hyers stable if there exists a real constant k > 0 such that, for each ϵ > 0 and for each solution y C 2 ( [ a , b ] ) of the inequality problem
D a + α y ( t ) + λ ( D a + α 1 y ) ( t ) f ( t , y ( t ) ) ϵ , t [ a , b ] , y ( a ) = y ( a ) = 0 ,
there exists a solution x C 2 ( [ a , b ] ) of the problem (1) (or, equivalently, of (3)) such that
| y ( t ) x ( t ) | k ϵ , t [ a , b ] .
Remark 3.
If we look at what is inside the modulus function in (15) as a single “new” function h, it directly follows that a function y C 2 ( [ a , b ] ) is a solution of the inequality in (15) if and only if there exists a function h C ( [ a , b ] ) (which depends on y) such that
(i) 
| h ( t ) | ϵ , t [ a , b ] ,
(ii) 
y ( a ) = y ( a ) = 0 ,
(iii) 
D a + α y ( t ) + λ ( D a + α 1 y ) ( t ) f ( t , y ( t ) ) = h ( t ) , t [ a , b ] .
Definition 4.
The fractional integral Equation (3) is Ulam–Hyers–Rassias stable with respect to φ : [ a , b ] R + if there exists a real constant k φ > 0 such that, for each ϵ > 0 and for each solution y of
y ( t ) e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , y ( s ) ) d s d u ϵ φ ( t ) , t [ a , b ] ,
there exists a solution x of the problem (3) with
| y ( t ) x ( t ) | k φ ϵ φ ( t ) , t [ a , b ] .

4.1. Ulam–Hyers Stability

As indicated above, we will start by identifying conditions that guarantee the Ulam–Hyers of the IVPFO (1).
Theorem 6.
Let the continuously differentiable function f satisfy the Lipschitz conditions (13)–(14), for all t [ a , b ] , and assume that
K L 1 + L 2 ( b a ) α 1 Γ ( α ) < 1 .
If y C 2 ( [ a , b ] ) satisfies the inequality and initial conditions (15) (with ϵ > 0 ), for all t [ a , b ] , then there exists a unique solution x C 2 ( [ a , b ] ) of the IVPFO (1) such that
| y ( t ) x ( t ) | k ϵ , t [ a , b ] ,
for
k = ( b a ) α 1 α 1 e L 1 1 e λ ( b a ) ( b a ) α 1 λ Γ ( α )
which, in particular, means that the IVPFO (1) is Ulam–Hyers stable.
Proof. 
According to the hypothesis, there exists a unique solution of the IVPFO (1).
Let y C 2 ( [ a , b ] ) be any solution of the inequality of (15). By Remark 3, following the procedure of Proposition 1, one has that
y ( t ) = e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , y ( s ) ) d s d u + e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u h ( s ) d s d u ,
with | h ( t ) | < ϵ . Thus, we have that
| x ( t ) y ( t ) | = e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u ( f ( s , x ( s ) ) f ( s , y ( s ) ) ) d s d u e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u h ( s ) d s d u e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u | f ( s , x ( s ) ) f ( s , y ( s ) ) | d s d u + e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u | h ( s ) | d s d u L 1 e λ t Γ ( α 1 ) a t e λ u a u ( u s ) α 2 | x ( s ) y ( s ) | d s d u + ϵ e λ t Γ ( α 1 ) a t e λ u a u ( u s ) α 2 d s d u L 1 1 e λ ( t a ) λ Γ ( α 1 ) a t a u ( u s ) α 2 | x ( s ) y ( s ) | d s d u + ϵ ( 1 e λ ( t a ) ) ( b a ) α 1 λ Γ ( α ) 1 e λ ( t a ) λ Γ ( α 1 ) ϵ ( b a ) α 1 α 1 + a t a u L 1 ( u s ) α 2 | x ( s ) y ( s ) | d s d u 1 e λ ( b a ) λ Γ ( α 1 ) ϵ ( b a ) α 1 α 1 + a t a u L 1 ( u s ) α 2 | x ( s ) y ( s ) | d s d u .
Thus, according to Theorem 3, we have that
| x ( t ) y ( t ) | ϵ ( b a ) α 1 α 1 e 1 e λ ( b a ) λ Γ ( α 1 ) L 1 ( t a ) α 1 α 1 ϵ ( b a ) α 1 α 1 e L 1 1 e λ ( b a ) ( b a ) α 1 λ Γ ( α ) ,
and we conclude the above claimed inequality and that the IVPFO (1) is Ulam–Hyers stable. □

4.2. Ulam–Hyers–Rassias Stability

We will now consider the Ulam–Hyers–Rassias stability. For that purpose, we consider the space C ( [ a , b ] ) equipped with the Bielecki type metric
d ( x , y ) = sup t [ a , b ] | x ( t ) y ( t ) | σ ( t ) ,
where σ is a non-decreasing continuous function σ : [ a , b ] R + . It is known that ( C ( [ a , b ] ) , d ) is a complete metric space (cf. [25]).
Theorem 7.
Let f : [ a , b ] × R R be a continuous function satisfying the Lipschitz condition
| f ( t , ρ 1 ) f ( t , ρ 2 ) | L | ρ 1 ρ 2 | , ρ 1 , ρ 2 R , t [ a , b ] ,
with L > 0 . Additionally, let σ : [ a , b ] R + be a nondecreasing function and suppose that exist a constant ξ [ 0 , 1 ) such that
e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u σ ( s ) d s d u ξ σ ( t ) , t [ a , b ] .
If y satisfies
y ( t ) e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , y ( s ) ) d s d u ϵ σ ( t ) , t [ a , b ] ,
and L ξ < 1 , then there exist a solution x of the fractional integral Equation (3) such that
| x ( t ) y ( t ) | ϵ σ ( t ) 1 L ξ , t [ a , b ] ,
i.e., under the present conditions, the fractional integral Equation (3) has the Ulam–Hyers– Rassias stability.
Proof. 
Having in mind the fractional integral Equation (3), we will consider (in the framework of the above presented Bielecki type metric) the operator T : C ( [ a , b ] , d ) C ( [ a , b ] , d ) defined by
( T y ) ( t ) = e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , y ( s ) ) d s d u .
Let us first prove that T is strictly contractive in C ( [ a , b ] , d ) . For any v , w C ( [ a , b ] , d ) , we have
d ( T v , T w ) = sup t [ a , b ] e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u ( f ( s , v ( s ) ) f ( s , w ( s ) ) d s d u σ ( t ) L sup t [ a , b ] e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u σ ( s ) v ( s ) w ( s ) σ ( s ) d s d u σ ( t ) L ξ d ( v , w ) .
Consequently, for L ξ < 1 , we have that T is strictly contractive in the present framework, and we have a unique solution x to the equation T y = y .
Let us now identify ϵ as an upper bound for d ( T y , y ) , and use this knowledge. Indeed, from the hypothesis, we have
| y ( t ) T y ( t ) | = y ( t ) e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , y ( s ) ) d s d u < ϵ σ ( t ) ,
which allows us to conclude that
d ( x , y ) 1 1 L ξ d ( y , T y ) ϵ 1 L ξ ,
and so
| x ( t ) y ( t ) | ϵ 1 L ξ σ ( t ) , t [ a , b ] .
The Ulam–Hyers stability is a particular case of the Ulam–Hyers–Rassias stability in the sense that instead of having a function φ controlling the differences in the last stability, we simply have a constant k in the first one. Thus, attending that
e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u d s d u ( b a ) α 1 λ Γ ( α ) ( 1 e λ ( b a ) ) , t [ a , b ] ,
and proceeding in an identical way to the proof of Theorem 7, we would pass from an upper bound that depends on a function (of the variable t) to an upper bound in the form of a constant, which is here directly concluded (following the proof of Theorem 7) in the next result:
Corollary 1.
Let f : [ a , b ] × R R be a continuous function satisfying the Lipschitz condition
| f ( t , ρ 1 ) f ( t , ρ 2 ) | L | ρ 1 ρ 2 | , ρ 1 , ρ 2 R , t [ a , b ] ,
with L > 0 . Let
η = ( b a ) α 1 λ Γ ( α ) ( 1 e λ ( b a ) ) .
If L η < 1 and y satisfies
y ( t ) e λ t Γ ( α 1 ) a t a u ( u s ) α 2 e λ u f ( s , y ( s ) ) d s d u ϵ , t [ a , b ] ,
then there exist a solution x of the fractional integral Equation (3) such that
| x ( t ) y ( t ) | ϵ 1 L η , t [ a , b ] ,
i.e., under the above conditions, the fractional integral Equation (3) has the Ulam–Hyers stability.
Remark 4.
Please, note that the constants k in (18) and 1 1 L η in (20) cannot be compared for all the values of the parameters. Consider, for example, the following cases. Admit that L = L 1 = L 2 = 1 20 , α = 7 4 and consider two intervals, one of amplitude equal 1 and another one with amplitude 0.8 . With these values, we have that, for λ ] 2 , 0 [ ] 0 , 5 [ condition (17) is verified and also, L η < 1 for η as defined in (19). For the case b a = 1 , it is possible to observe that k > 1 1 L η (cf. Figure 1, where p ( λ ) > q ( λ ) ). For the case b a = 0.8 , we verify that k < 1 1 L η (cf. Figure 2, where p ( λ ) < q ( λ ) ).

4.3. Concrete Examples

Let us now consider some concrete examples to illustrate the above theory.
We start by considering the following IVPFO
( D 3 2 x ) ( t ) + λ ( D 1 2 x ) ( t ) = t 75 ( x ( t ) + sin ( t ) ) , x ( 2 ) = x ( 2 ) = 0 ,
for t [ 2 , 3 ] . Thus, in the previous notation, we have in here α = 3 2 , a = 2 , b = 3 and
f ( t , ρ ) = t 75 ( ρ + sin ( t ) ) ,
being clear that f is a continuously differentiable function.
According to Theorem 4, there exists, at least, one solution of the IVPFO (21). In addition, having in mind that f ( t , x ( t ) ) = 1 75 x ( t ) + t x ( t ) + sin ( t ) + t cos ( t ) , for t [ 2 , 3 ] , one has that
| f ( t , x ( t ) ) f ( t , y ( t ) ) | 1 25 | x ( t ) y ( t ) | , | f ( t , x ( t ) ) f ( t , y ( t ) ) | 1 25 | x ( t ) y ( t ) | + | x ( t ) y ( t ) | .
Following Theorem 5 and its notation, we have in here L 1 = L 2 = 1 25 . Thus, for a = 2 and b = 3 , we obtain that
K L 1 + L 2 ( b a ) α 1 Γ ( α ) < 1 ,
for λ [ 1 , 0 ] [ 0 , 9 ] (cf. (12) and Figure 3). Thus, for these cases of λ , the IVPFO (21) admits a unique solution in C 2 ( [ 2 , 3 ] ) . Moreover, from Theorem 6, we also know that for those λ the IVPFO (21) is Ulam–Hyers stable.
The example we have just analyzed allows us to see that there really are classes of problems, dependent on λ , in which the conditions required in Theorem 5 are met, and there are still other cases (for different parameters λ ) in which this is not the case. In view of this, and keeping in mind that the conditions of Theorem 5 are just sufficient conditions, an open analysis eventually involves obtaining other weaker conditions according to which the uniqueness of solution for those classes of problems can still be guaranteed. The same can be envisaged for Theorem 6 and its sufficient conditions to guarantee the stability of the Ulam–Hyers type.
Let us now investigate the Ulam–Hyers–Rassias stability of
x ( t ) = e λ t Γ ( 3 2 1 ) 2 t 2 u ( u s ) 3 2 2 e λ u t 75 ( y ( s ) + sin ( s ) ) d s d u ,
for t [ 2 , 3 ] and λ = 3 .
Letting σ ( t ) = e t , we have that σ is a non-decreasing function and
e 3 t Γ ( 1 2 ) 2 t 2 u ( u s ) 3 2 2 e 3 u σ ( s ) d s d u 1 5 σ ( t ) , t [ 2 , 3 ]
(cf. Figure 4). Thus, for the notation of Theorem 7, we have L = 1 25 , ξ = 1 5 and so L ξ = 1 125 < 1 .
Take y ( t ) = 1 10 ( t 2 ) 2 . We have that y C 2 ( [ 2 , 3 ] ) and y ( 2 ) = y ( 2 ) = 0 . We have that (cf. Figure 5)
y ( t ) e 3 t Γ ( 3 2 ) 2 t 2 u ( u s ) 1 2 e 3 u s 75 ( y ( s ) + sin ( s ) ) d s d u 1 200 σ ( t ) , t [ 2 , 3 ] .
Thus, according to Theorem 7, the problem (22) is Ulam–Hyers–Rassias stable with respect to σ ( t ) = e t and
| y ( t ) x ( t ) | 5 e t 992 , t [ 2 , 3 ] .
Moreover, we can also observe that
y ( t ) e 3 t Γ ( 3 2 ) 2 t 2 u ( u s ) 1 2 e 3 u s 75 ( y ( s ) + sin ( s ) ) d s d u 1 10 , t [ 2 , 3 ] .
Thus, applying Corollary 1 and the respective notation, we have that ϵ = 1 10 . Additionally, η = 1 e 3 3 Γ 3 2 0.36 and we conclude that
| x ( t ) y ( t ) | 0.1 .
In this last example, it is relevant to emphasize the importance of the function σ in the whole process, with special predominance, from the outset, in the determination of the exhibited upper bounds. In this case, we chose to work with the exponential function, and this had expected consequences given the growth that the function presents. Incidentally, the importance of the choice and the impact that the σ function has is well evidenced by the fact that the same problem can be Ulam–Hyers–Rassias stable for a given σ 1 function and not Ulam–Hyers–Rassias stable for another σ 2 function. Thus, it is precisely for this reason that the Ulam–Hyers–Rassias stability is determined depending on the chosen σ function (and it is also for this reason that this is explicitly mentioned in the name of this type of stability).
Let us now consider the following different IVPFO
( D 6 5 x ) ( t ) + λ ( D 1 5 x ) ( t ) = t 10 x ( t ) e t , t [ 0 , 1 ] , x ( 0 ) = x ( 0 ) = 0 .
Accordingly to the previous notations, we have now α = 6 5 , a = 0 , b = 1 and f ( t , x ( t ) ) = t 10 x ( t ) e t . It is clear that f is a continuously differentiable function in [ 0 , 1 ] × R . Thus there exists, at least, one solution of the IVPFO (23) (cf. Theorem 4). Moreover, one has that
| f ( t , x ( t ) ) f ( t , y ( t ) ) | 1 10 | x ( t ) y ( t ) | , | f ( t , x ( t ) ) f ( t , y ( t ) ) | 1 10 | x ( t ) y ( t ) | + | x ( t ) y ( t ) | .
Following Theorem 5, we have L 1 = L 2 = 1 10 . Since a = 0 and b = 1 , we obtain that for λ ] 1 , 0 [ 0 , 5 2 , the condition
K L 1 + L 2 ( b a ) α 1 Γ ( α ) < 1
is verified (cf. Figure 6), which means that the IVPFO (23) admits a unique solution in C 2 ( [ 0 , 1 ] ) when considering those values of λ (cf. Theorem 5).
Thus, for these cases of λ , the IVPFO (23) admits a unique solution in C 2 ( [ 0 , 1 ] ) . Moreover, from Theorem 6, we also know that for those λ the IVPFO (23) is Ulam–Hyers stable.
Let us now analyse the Ulam–Hyers–Rassias stability of
x ( t ) = e λ t Γ ( 6 5 1 ) 0 t 0 u ( u s ) 6 5 2 e λ u s 10 x ( s ) e s d s d u ,
for t [ 0 , 1 ] , λ = 2 , and with respect to σ ( t ) = t . Let x C 2 ( [ 0 , 1 ] ) be the exact solution of the IVPFO (23), and let us consider y ( t ) = sin ( t ) t . It follows that y C 2 ( [ 0 , 1 ] ) and y ( 0 ) = y ( 0 ) = 0 . We have that σ is a nondecreasing function and
e 2 t Γ ( 1 5 ) 0 t 0 u ( u s ) 4 5 e 2 u σ ( s ) d s d u 1 4 σ ( t ) , t [ 0 , 1 ]
(cf. Figure 7).
For the notation of Theorem 7, we have L = 1 25 and ξ = 1 4 , and so L ξ = 1 100 < 1 . Thus,
y ( t ) e 2 t Γ ( 1 5 ) 0 t 0 u ( u s ) 4 5 e 2 u s 10 y ( s ) e s d s d u 7 50 σ ( t ) , t [ 0 , 1 ]
(cf. Figure 8).
Therefore, according to Theorem 7, the problem (24) is Ulam–Hyers–Rassias stable with respect to σ ( t ) = t and
| y ( t ) x ( t ) | 14 99 σ ( t ) , t [ 0 , 1 ] .
In this last example, we deliberately chose σ ( t ) = t to work with y ( t ) = sin ( t ) t , which can be considered not the most ideal choice (which, by the way, can be easily noticed when we look at Figure 8 and see, on the right, the “greatest” distance between the two functions represented there). Anyway, we consider this example important because it emphasizes that the theoretical conditions obtained earlier are robust enough to guarantee stability in less favorable or obvious choices.
Moreover, according to Corollary 1, we can also conclude the Ulam–Hyers stability. Using the respective notation of the Corollary, we have that ϵ = 7 50 , η = 1 e 2 2 Γ 6 5 0.47 and L = 1 25 . Thus, we conclude that | x ( t ) y ( t ) | 7 50 .

5. Conclusions

We conclude this article by summarizing the results obtained. We analyze a class of nonlinear fractional differential equations, with initial conditions, characterized by having the Riemann–Liouville fractional derivative of order α ( 1 , 2 ) . Having made use of distinct fixed-point arguments, we were able to deduce conditions that guarantee the existence and uniqueness of solutions in a frame of adequate spaces, and we also obtained sufficient conditions to have the Ulam–Hyers and Ulam–Hyers–Rassias stabilities of the problems in the analysis (where the use of a Bielecki-type metric and some additional contractive arguments were of crucial importance). In the last section, some examples were included mainly to illustrate that the conditions obtained in the theoretical part really exist and can be considered in particular cases.

Author Contributions

Conceptualization, L.P.C. and A.S.S.; methodology, L.P.C. and A.S.S.; formal analysis, L.P.C. and A.S.S.; investigation, L.P.C. and A.S.S.; writing—original draft preparation, L.P.C. and A.S.S.; writing—review and editing, L.P.C. and A.S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Center for Research and Development in Mathematics and Applications (CIDMA) through the Portuguese Foundation for Science and Technology (FCT—Fundação para a Ciência e a Tecnologia), reference UIDB/04106/2020. Additionally, A. Silva is also funded by national funds (OE), through FCT, I.P., in the scope of the framework contract foreseen in the numbers 4, 5, and 6 of article 23, of the Decree-Law 57/2016, of 29 August, changed by Law 57/2017, of 19 July.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the Referees for their work in reviewing the present investigation and for providing their observations and suggestions.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Castro, L.P.; Silva, A.S. On the solution and Ulam-Hyers-Rassias stability of a Caputo fractional boundary value problem. Math. Biosci. Eng. 2022, 19, 10809–10825. [Google Scholar] [CrossRef] [PubMed]
  2. Fahad, H.M.; Fernandez, A.; Ur, R.M.; Siddiqi, M. Tempered and Hadamard-type fractional calculus with respect to functions. Mediterr. J. Math. 2021, 18, 143. [Google Scholar] [CrossRef]
  3. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2016. [Google Scholar]
  4. Mali, A.D.; Kucche, K.D.; Fernandez, A.; Fahad, H.M. On tempered fractional calculus with respect to functions and the associated fractional differential equations. Math. Methods Appl. Sci. 2022, 45, 11134–11157. [Google Scholar] [CrossRef]
  5. Miller, K.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley: New York, NY, USA, 1993. [Google Scholar]
  6. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  7. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives. Theory and Applications; Translated from the 1987 Russian Original; Gordon and Breach Science Publishers: Yverdon, Switzerland, 1993. [Google Scholar]
  8. Aoki, T. On the stability of the linear transformation in Banach spaces. J. Math. Soc. Jpn. 1950, 2, 64–66. [Google Scholar] [CrossRef]
  9. Castro, L.P.; Simões, A.M. Hyers-Ulam-Rassias stability of nonlinear integral equations through the Bielecki metric. Math. Methods Appl. Sci. 2018, 41, 7367–7383. [Google Scholar] [CrossRef]
  10. Castro, L.P.; Simões, A.M. Different types of Hyers-Ulam-Rassias stabilities for a class of integro-differential equations. Filomat 2017, 31, 5379–5390. [Google Scholar] [CrossRef] [Green Version]
  11. Hyers, D.H. On the stability of the linear functional equation. Proc. Natl. Acad. Sci. USA 1941, 27, 222–224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Marian, D. Semi-Hyers–Ulam–Rassias stability of the convection partial differential equation via Laplace transform. Mathematics 2021, 9, 2980. [Google Scholar] [CrossRef]
  13. Marian, D.; Ciplea, S.A.; Lungu, N. Hyers–Ulam stability of a system of hyperbolic partial differential equations. Mathematics 2022, 10, 2183. [Google Scholar] [CrossRef]
  14. Rassias, T.M. On the stability of the linear mapping in Banach spaces. Proc. Am. Math. Soc. 1978, 72, 297–300. [Google Scholar] [CrossRef]
  15. Sousa, J.V.d.C.; Oliveira, D.S.; Capelas de Oliveira, E. A note on the mild solutions of Hilfer impulsive fractional differential equations. Chaos Solitons Fractals 2021, 147, 110944. [Google Scholar] [CrossRef]
  16. Ulam, S.M. Problems in Modern Mathematics; John Wiley & Sons: New York, NY, USA, 1940. [Google Scholar]
  17. Vanterler da Costa Sousa, J.; Kucche, K.D.; de Oliveira, E.C. Stability of mild solutions of the fractional nonlinear abstract Cauchy problem. Electron Res. Arch. 2022, 30, 272–288. [Google Scholar]
  18. Chai, G. Existence results for boundary value problems of nonlinear fractional differential equations. Comp. Math. Appl. 2011, 62, 2374–2382. [Google Scholar] [CrossRef] [Green Version]
  19. Xu, L.; Dong, Q.; Li, G. Existence and Hyers–Ulam stability for three-point boundary value problems with Riemann–Liouville fractional derivatives and integrals. Adv. Differ. Equations 2018, 2018, 458. [Google Scholar] [CrossRef] [Green Version]
  20. Ahmad, B.; Matar, M.M.; El-Salmy, O.M. Existence of Solutions and Ulam Stability for Caputo Type Sequential Fractional Differential Equations of Order α ∈ (2,3). Int. J. Anal. Appl. 2017, 15, 86–101. [Google Scholar]
  21. Alvan, M.; Darzi, R.; Mahmoodi, A. Existence results for a new class of boundary value problems of nonlinear fractional differential equations. Mathematics 2016, 4, 13. [Google Scholar] [CrossRef] [Green Version]
  22. Bilgici, S.S.; Şan, M. Existence and uniqueness results for a nonlinear singular fractional differential equation of order σ ∈ (1,2). AIMS Math. 2021, 6, 13041–13056. [Google Scholar] [CrossRef]
  23. Bainov, D.; Simeonov, P. Integral Inequalities and Applications; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1992. [Google Scholar]
  24. Bai, Z.; Sun, W. Existence and multiplicity of positive solutions for singular fractional boundary value problems. Comput. Math. Appl. 2012, 63, 1369–1381. [Google Scholar] [CrossRef]
  25. Rolewicz, S. Functional Analysis and Control Theory: Linear Systems; Mathematics and its Applications (East European Series), 29; D. Reidel Publishing Co.: Dordrecht, The Netherlands; PWN–Polish Scientific Publishers: Warsaw, Poland, 1987. [Google Scholar]
Figure 1. The graphs of p ( λ ) = k ( λ ) and q ( λ ) = 1 1 L η ( λ ) for λ [ 2 , 0 ] [ 0 , 5 ] : case b a = 1 .
Figure 1. The graphs of p ( λ ) = k ( λ ) and q ( λ ) = 1 1 L η ( λ ) for λ [ 2 , 0 ] [ 0 , 5 ] : case b a = 1 .
Mathematics 11 00297 g001
Figure 2. The graphs of p ( λ ) = k ( λ ) and q ( λ ) = 1 1 L η ( λ ) for λ [ 2 , 0 ] [ 0 , 5 ] : case b a = 0.8 .
Figure 2. The graphs of p ( λ ) = k ( λ ) and q ( λ ) = 1 1 L η ( λ ) for λ [ 2 , 0 ] [ 0 , 5 ] : case b a = 0.8 .
Mathematics 11 00297 g002
Figure 3. The graphs of z 1 ( λ ) = K ( λ ) L 1 + L 2 Γ 3 2 and z 2 = 1 .
Figure 3. The graphs of z 1 ( λ ) = K ( λ ) L 1 + L 2 Γ 3 2 and z 2 = 1 .
Mathematics 11 00297 g003
Figure 4. The graphs of p 1 ( t ) = 1 5 σ ( t ) = 1 5 e t (the upper one), and p 2 ( t ) = e 3 t Γ ( 1 2 ) 2 t 2 u ( u s ) 3 2 2 e 3 u σ ( s ) d s d u , t [ 2 , 3 ] .
Figure 4. The graphs of p 1 ( t ) = 1 5 σ ( t ) = 1 5 e t (the upper one), and p 2 ( t ) = e 3 t Γ ( 1 2 ) 2 t 2 u ( u s ) 3 2 2 e 3 u σ ( s ) d s d u , t [ 2 , 3 ] .
Mathematics 11 00297 g004
Figure 5. The graphs of q 1 ( t ) = 1 200 σ ( t ) (the upper one) and q 2 ( t ) = y ( t ) e 3 t Γ ( 3 2 ) 2 t 2 u ( u s ) 1 2 e 3 u s 75 ( y ( s ) + sin ( s ) ) d s d u , t [ 2 , 3 ] .
Figure 5. The graphs of q 1 ( t ) = 1 200 σ ( t ) (the upper one) and q 2 ( t ) = y ( t ) e 3 t Γ ( 3 2 ) 2 t 2 u ( u s ) 1 2 e 3 u s 75 ( y ( s ) + sin ( s ) ) d s d u , t [ 2 , 3 ] .
Mathematics 11 00297 g005
Figure 6. The graphs of h 1 ( λ ) = K ( λ ) L 1 + L 2 Γ 6 5 and h 2 = 1 .
Figure 6. The graphs of h 1 ( λ ) = K ( λ ) L 1 + L 2 Γ 6 5 and h 2 = 1 .
Mathematics 11 00297 g006
Figure 7. The graphs of m 1 ( t ) = e 2 t Γ ( 1 5 ) 0 t 0 u ( u s ) 4 5 e 2 u σ ( s ) d s d u and m 2 ( t ) = 1 4 σ ( t ) = t 4 , t [ 0 , 1 ] .
Figure 7. The graphs of m 1 ( t ) = e 2 t Γ ( 1 5 ) 0 t 0 u ( u s ) 4 5 e 2 u σ ( s ) d s d u and m 2 ( t ) = 1 4 σ ( t ) = t 4 , t [ 0 , 1 ] .
Mathematics 11 00297 g007
Figure 8. The graphs of w 1 ( t ) = y ( t ) e 2 t Γ ( 1 5 ) 0 t 0 u ( u s ) 4 5 e 2 u s 10 y ( s ) e s d s d u (the lower one) and w 2 ( t ) = 7 50 σ ( t ) , t [ 0 , 1 ] .
Figure 8. The graphs of w 1 ( t ) = y ( t ) e 2 t Γ ( 1 5 ) 0 t 0 u ( u s ) 4 5 e 2 u s 10 y ( s ) e s d s d u (the lower one) and w 2 ( t ) = 7 50 σ ( t ) , t [ 0 , 1 ] .
Mathematics 11 00297 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Castro, L.P.; Silva, A.S. On the Existence and Stability of Solutions for a Class of Fractional Riemann–Liouville Initial Value Problems. Mathematics 2023, 11, 297. https://doi.org/10.3390/math11020297

AMA Style

Castro LP, Silva AS. On the Existence and Stability of Solutions for a Class of Fractional Riemann–Liouville Initial Value Problems. Mathematics. 2023; 11(2):297. https://doi.org/10.3390/math11020297

Chicago/Turabian Style

Castro, Luís P., and Anabela S. Silva. 2023. "On the Existence and Stability of Solutions for a Class of Fractional Riemann–Liouville Initial Value Problems" Mathematics 11, no. 2: 297. https://doi.org/10.3390/math11020297

APA Style

Castro, L. P., & Silva, A. S. (2023). On the Existence and Stability of Solutions for a Class of Fractional Riemann–Liouville Initial Value Problems. Mathematics, 11(2), 297. https://doi.org/10.3390/math11020297

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

Article Metrics

Back to TopTop