Next Article in Journal / Special Issue
Some New Fractional Estimates of Inequalities for LR-p-Convex Interval-Valued Functions by Means of Pseudo Order Relation
Previous Article in Journal
Correction: Bates et al. Vector Fields and Differential Forms on the Orbit Space of a Proper Action. Axioms 2021, 10, 118
Previous Article in Special Issue
A Novel Numerical Method for Solving Fractional Diffusion-Wave and Nonlinear Fredholm and Volterra Integral Equations with Zero Absolute Error
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sequential Riemann–Liouville and Hadamard–Caputo Fractional Differential Systems with Nonlocal Coupled Fractional Integral Boundary Conditions

by
Chanakarn Kiataramkul
1,†,
Weera Yukunthorn
2,†,
Sotiris K. Ntouyas
3,4,† and
Jessada Tariboon
1,*
1
Intelligent and Nonlinear Dynamic Innovations Research Center, Department of Mathematics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand
2
Faculty of Science and Technology, Kanchanaburi Rajabhat University, Kanchanaburi 71000, Thailand
3
Department of Mathematics, University of Ioannina, 451 10 Ioannina, Greece
4
Nonlinear Analysis and Applied Mathematics (NAAM)-Research Group, Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 30 June 2021 / Revised: 29 July 2021 / Accepted: 30 July 2021 / Published: 31 July 2021
(This article belongs to the Special Issue Fractional Calculus - Theory and Applications)

Abstract

:
In this paper, we initiate the study of existence of solutions for a fractional differential system which contains mixed Riemann–Liouville and Hadamard–Caputo fractional derivatives, complemented with nonlocal coupled fractional integral boundary conditions. We derive necessary conditions for the existence and uniqueness of solutions of the considered system, by using standard fixed point theorems, such as Banach contraction mapping principle and Leray–Schauder alternative. Numerical examples illustrating the obtained results are also presented.

1. Introduction

Fractional differential equations have played a very important role in almost all branches of applied sciences because they are considered a valuable tool to model many real world problems. For details and applications, we refer the reader to monographs [1,2,3,4,5,6,7,8,9,10,11]. The study of coupled systems of fractional differential equations is important as such systems appear in various problems in applied sciences, see [12,13,14,15,16].
On the other hand, multi-term fractional differential equations also gained considerable importance in view of their occurrence in the mathematical models of certain real world problems, such as behavior of real materials [17], continuum and statistical mechanics [18], an inextensible pendulum with fractional damping terms [19], etc.
Fractional differential equations have several kinds of fractional derivatives, such as Riemann–Liouville fractional derivative, Caputo fractional derivative, Hadamard fractional derivative, and so on. In the literature, there are many papers studying existence and uniqueness results for boundary value problems and coupled systems of fractional differential equations and used mixed types of fractional derivatives, see [20,21,22,23,24,25,26,27,28,29]. In [23], the following boundary value problem is considered:
R L D q [ C D r x ( t ) g ( t , x ( t ) ) ] = f ( t , x ( t ) ) , 0 < t < T , x ( η ) = ϕ ( x ) , I p x ( T ) = h ( x ) ,
where R L D q , C D r are Riemann–Liouville and Caputo fractional derivatives of orders q , r ( 0 , 1 ) , respectively, I p is the Riemann–Liouville fractional integral of order p > 0 , f , g : J × R R are given continuous functions and ϕ , h : C ( J , R ) R are two given functionals.
In [24], the authors initiated the study of a coupled system of sequential mixed Caputo and Hadamard fractional differential equations supplemented with coupled separated boundary conditions. To be more precisely, in [24], existence and uniqueness results are established for the following couple system:
C D p 1 H D q 1 x ( t ) = f ( t , x ( t ) , y ( t ) ) , t [ a , b ] , H D q 2 C D p 2 y ( t ) = g ( t , x ( t ) , y ( t ) ) , t [ a , b ] , α 1 x ( a ) + α 2 C D p 2 y ( a ) = 0 , β 1 x ( b ) + β 2 C D p 2 y ( b ) = 0 , α 3 y ( a ) + α 4 H D q 1 x ( a ) = 0 , β 3 y ( b ) + β 4 H D q 1 x ( b ) = 0 ,
where C D p i and H D q i are notations of the Caputo and Hadamard fractional derivatives of orders p i and q i , respectively, 0 < p i , q i 1 , i = 1 , 2 , f , g : [ a , b ] × R × R R are nonlinear continuous functions, a > 0 , α i R { 0 } , β i R , i = 1 , , 4 .
In [25], the existence and uniqueness of solutions for neutral fractional order coupled systems containing mixed Caputo and Riemann–Liouville sequential fractional derivatives were studied, complemented with nonlocal multi-point and Riemann–Stieltjes integral multi-strip conditions of the form:
c D q ( R L D p x ( t ) + f ( t , x ( t ) ) = g ( t , x ( t ) , y ( t ) ) , t ( 0 , 1 ) , c D q 1 ( R L D p 1 y ( t ) + f 1 ( t , y ( t ) ) = g 1 ( t , x ( t ) , y ( t ) ) , t ( 0 , 1 ) , x ( 0 ) = 0 , b x ( 1 ) = a 0 1 y ( s ) d H ( s ) + i = 1 n α i ξ i η i y ( s ) d s , y ( 0 ) = 0 , b 1 y ( 1 ) = a 1 0 1 x ( s ) d H ( s ) + j = 1 m β j θ j ζ j x ( s ) d s ,
where R L D p , R L D p 1 , and c D q , c D q 1 denote the Riemann–Liouville and Caputo fractional derivatives of order p , p 1 and q , q 1 , respectively, 0 < p , p 1 , q , q 1 1 , with 1 < p + q 2 , 1 < p 1 + q 1 2 , f , f 1 and g , g 1 are given continuous functions, 0 < ξ i < η i < 1 , 0 < θ j < ζ j < 1 , α i , β j R , i = 1 , 2 , , n , j = 1 , 2 , , m , a , a 1 , b , b 1 R , and H ( · ) is a function of bounded variation.
To the best of the authors’ knowledge, there are some papers dealing with sequential mixed type fractional derivatives, but we not find in the literature papers dealing with coupled systems with sequential Riemann–Liouville and Hadamard–Caputo fractional differential equations. Motivated by this fact, and to fill this gap, in the present paper, we investigate the existence and uniqueness of solutions for the following coupled system of sequential Riemann–Liouville and Hadamard–Caputo fractional differential equations supplemented with nonlocal coupled fractional integral boundary conditions
R L D p 1 H C D q 1 x ( t ) = f ( t , x ( t ) , y ( t ) ) , t [ 0 , T ] , R L D p 2 H C D q 2 y ( t ) = g ( t , x ( t ) , y ( t ) ) , t [ 0 , T ] , H C D q 1 x ( 0 ) = 0 , x ( T ) = i = 1 m α i R L I β i y ( ξ i ) , H C D q 2 y ( 0 ) = 0 , y ( T ) = j = 1 k λ j R L I δ j x ( η j ) ,
where R L D p r and H C D q r are the Riemann–Liouville and Hadamard–Caputo fractional derivatives of orders p r and q r , respectively, 0 < p r , q r < 1 , r = 1 , 2 , the nonlinear continuous functions f , g : [ 0 , T ] × R 2 R , R L I ϕ is the Riemann–Liouville fractional integral of orders ϕ > 0 , ϕ { β i , δ j } and given constants α i , λ j R , ξ i , η j ( 0 , T ) , i = 1 , , m , j = 1 , , k .
Let us compare the coupled system (4) with the coupled system (2) studied in [24].
(i)
In (2), we studied a coupled system consisting by mixed Caputo and Hadamard fractional derivatives, while, in (4), we consider mixed Riemann–Liouville and Hadamard–Caputo fractional derivatives.
(ii)
In (2), the coupled system was subjected to coupled separated boundary conditions, while, in (4), the coupled system is subjected to nonlocal coupled fractional integral boundary conditions.
(iii)
In both problems (4) and (2), the same method is used to establish the existence and uniqueness results, and based on standard fixed point theorems, but their presentation in the framework of mixed coupled Caputo and Hadamard and Riemann–Liouville and Hadamard–Caputo fractional derivatives is new.
We also notice that the conditions H C D q 1 x ( 0 ) = 0 and H C D q 2 y ( 0 ) = 0 are necessary for the well-posedness of the problem.
By using standard tools from fixed point theory in the present study, we establish existence and uniqueness results for the coupled system (4). The Banach contraction mapping principle is used to obtain the existence and uniqueness result, while an existence result is derived via the Leray–Schauder alternative.
The rest of the paper is organized as follows. In Section 2, some basic definitions and lemmas from fractional calculus are recalled. In addition, an auxiliary lemma, concerning a linear variant of (4), which plays a key role in obtaining the main results, is proved. The main results are presented in Section 3, which also include examples illustrating the basic results. We emphasize that our results are new and significantly enhance the existing literature on the topic, and, as far as we know, they are the first results concerning a coupled system with sequential mixed Riemann–Liouville and Hadamard–Caputo fractional derivatives.

2. Preliminaries

In this section, we introduce some notations and definitions of fractional calculus [2,30] and present preliminary results needed in our proofs later.
Definition 1.
The Riemann–Liouville fractional derivative of order p > 0 of a continuous function f : ( 0 , ) R is defined by
R L D p f ( t ) = 1 Γ ( n p ) d d t n 0 t ( t s ) n p 1 f ( s ) d s , n 1 < p < n ,
where n = [ p ] + 1 , [ p ] denotes the integer part of a real number p and Γ is the Gamma function defined by Γ ( p ) = 0 e s s p 1 d s .
Definition 2.
The Riemann–Liouville fractional integral of order p of a function f : ( 0 , ) R , is defined as
R L I p f ( t ) = 1 Γ ( q ) 0 t ( t s ) p 1 f ( s ) d s , p > 0 ,
provided the right side is pointwise defined on R + .
Definition 3.
For an at least n-times differentiable function g : ( 0 , ) R , the Hadamard–Caputo derivative of fractional order q > 0 is defined as
H C D q g ( t ) = 1 Γ ( n q ) 0 t log t s n q 1 δ n g ( s ) d s s , n 1 < q < n , n = [ q ] + 1 ,
where δ = t d d t and log ( · ) = log e ( · ) .
Definition 4.
The Hadamard fractional integral of order q > 0 is defined as
H I q g ( t ) = 1 Γ ( α ) 0 t log t s q 1 g ( s ) d s s ,
provided the integral exists.
Lemma 1
(see [2]). Let p > 0 . Then, for y C ( 0 , T ) L ( 0 , T ) , it holds that
R L I p R L D p y ( t ) = y ( t ) + c 1 t p 1 + c 2 t p 2 + + c n t p n ,
where c i R , i = 1 , 2 , , n and n 1 < p < n .
Lemma 2
([30]). Let u A C δ n [ 0 , T ] or C δ n [ 0 , T ] and q C , where X δ n [ 0 , T ] = { g : [ 0 , T ] C : δ n 1 g ( t ) X [ 0 , T ] } . Then, we have
H I q ( H C D q ) u ( t ) = u ( t ) + c 0 + c 1 log t + c 2 log t 2 + + c n 1 log t n 1 ,
where c i R , i = 0 , 1 , 2 , , n 1 ( n = [ q ] + 1 ) .
Lemma 3
([2], p. 113). Let q > 0 and β > 0 be given constants. Then, the following formula
H I q t β = β q t β ,
holds.
Next, the integral equations are obtained by transformation of a linear variant of problem (4). For convenience in computation, we set some constants
Ω 1 = i = 1 m α i ξ i β i Γ ( β i + 1 ) , Ω 2 = j = 1 k λ j η j δ j Γ ( δ j + 1 )
and Λ = Ω 1 Ω 2 1 0 .
Lemma 4.
Let f * , g * C ( [ a , b ] , R ) be two given functions. Then, the linear system equivalent to problem (4) of sequential Riemann–Liouville and Hadamard–Caputo fractional differential equations
R L D p 1 H C D q 1 x ( t ) = f * ( t ) , t [ 0 , T ] , R L D p 2 H C D q 2 y ( t ) = g * ( t ) , t [ 0 , T ] , H C D q 1 x ( 0 ) = 0 , x ( T ) = i = 1 m α i R L I β i y ( ξ i ) , H C D q 2 y ( 0 ) = 0 , y ( T ) = j = 1 k λ j R L I δ j x ( η j ) ,
can be written into integral equations as
x ( t ) = 1 Λ i = 1 m α i R L I β i H I q 2 R L I p 2 g * ( ξ i ) + 1 Λ H I q 1 R L I p 1 f * ( T ) + Ω 1 Λ H I q 2 R L I p 2 g * ( T ) Ω 1 Λ j = 1 k λ j R L I δ j H I q 1 R L I p 1 f * ( η j ) + H I q 1 R L I p 1 f * ( t ) ,
and
y ( t ) = Ω 2 Λ i = 1 m α i R L I β i H I q 2 R L I p 2 g * ( ξ i ) + Ω 2 Λ H I q 1 R L I p 1 f * ( T ) + 1 Λ H I q 2 R L I p 2 g * ( T ) 1 Λ j = 1 k λ j R L I δ j H I q 1 R L I p 1 f * ( η j ) + H I q 2 R L I p 2 g * ( t ) .
Proof. 
For t [ 0 , T ] and by taking the Riemann–Liouville fractional integral of order p 1 to the first equation of (5), we obtain
H C D q 1 x ( t ) = c 1 t p 1 1 + R L I p 1 f * ( t ) , c 1 R .
Similarly, for the second equation of (5), we have
H C D q 2 y ( t ) = d 1 t p 2 1 + R L I p 2 g * ( t ) , d 1 R .
Since 0 < p r < 1 , r = 1 , 2 , the conditions H C D q 1 x ( 0 ) = 0 and H C D q 2 y ( 0 ) = 0 imply c 1 = 0 and d 1 = 0 , respectively. Applying the Hadamard fractional integral of orders q 1 and q 2 to (8) and (9), respectively, and substituting the values of c 1 , d 1 , we get
x ( t ) = c 0 + H I q 1 R L I p 1 f * ( t ) ,
and
y ( t ) = d 0 + H I q 2 R L I p 2 g * ( t ) .
Now, we consider the terms
i = 1 m α i R L I β i y ( ξ i ) = d 0 i = 1 m α i ξ i β i Γ ( β i + 1 ) + i = 1 m α i R L I β i H I q 2 R L I p 2 g * ( ξ i )
and
j = 1 k λ j R L I δ j x ( η j ) = c 0 j = 1 k λ j η j δ j Γ ( δ j + 1 ) + j = 1 k λ j R L I δ j H I q 1 R L I p 1 f * ( η j ) .
Consequently, by (10)–(13) and boundary fractional integral conditions in (5), it follows that
c 0 = 1 Λ i = 1 m α i R L I β i H I q 2 R L I p 2 g * ( ξ i ) + 1 Λ H I q 1 R L I p 1 f * ( T ) + Ω 1 Λ H I q 2 R L I p 2 g * ( T ) Ω 1 Λ j = 1 k λ j R L I δ j H I q 1 R L I p 1 f * ( η j ) ,
and
d 0 = Ω 2 Λ i = 1 m α i R L I β i H I q 2 R L I p 2 g * ( ξ i ) + Ω 2 Λ H I q 1 R L I p 1 f * ( T ) + 1 Λ H I q 2 R L I p 2 g * ( T ) 1 Λ j = 1 k λ j R L I δ j H I q 1 R L I p 1 f * ( η j ) ,
Substituting the values of c 0 and d 0 in (10) and (11), we obtain integral equations in (6) and (7), respectively, as desired.
The converse follows by direct computation. This completes the proof. □
Next, we establish formulas for multiple fractional integrals of Riemann–Liouville and Hadamard types.
Lemma 5.
Let a , b , c > 0 be constants. Then, we have
( i )
H I b R L I a ( 1 ) ( t ) = a b t a Γ ( a + 1 ) .
( i i )
R L I c H I b R L I a ( 1 ) ( t ) = a b Γ ( a + c + 1 ) t a + c .
Proof. 
Since R L I a ( 1 ) = t a Γ ( a + 1 ) , we have
H I b R L I a ( 1 ) ( t ) = 1 Γ ( a + 1 ) H I b t a = a b t a Γ ( a + 1 ) ,
by using Lemma 3, and (i) is proved. To prove (ii), taking the Riemann–Liouville fractional integral of order c > 0 in (14), we have
R L I c H I b R L I a ( 1 ) ( t ) = a b Γ ( a + 1 ) R L I c t a = a b Γ ( a + c + 1 ) t a + c ,
from R L I c t a = Γ ( a + 1 ) Γ ( a + c + 1 ) t a + c . The proof is completed. □
Corollary 1.
Let constants p r , q r , r = 1 , 2 , β i , ξ i , δ j , η j be defined in problem (4). Then, from Lemma 5, we have
H I q 1 R L I p 1 1 ( T ) = p 1 q 1 T p 1 Γ ( p 1 + 1 ) , H I q 2 R L I p 2 1 ( T ) = p 2 q 2 T p 2 Γ ( p 2 + 1 ) , R L I β i H I q 2 R L I p 2 1 ( ξ i ) = p 2 q 2 Γ ( p 2 + β i + 1 ) ξ i p 2 + β i , R L I δ j H I q 1 R L I p 1 1 ( η j ) = p 1 q 1 Γ ( p 1 + δ j + 1 ) η j p 1 + δ j ,
which will be used in the next section.

3. Main Results

Let C = C ( [ 0 , T ] , R ) be the Banach space of all continuous functions from [ 0 , T ] to R . Let X = { x ( t ) : x ( t ) C 2 ( [ 0 , T ] , R ) } be the space endowed with the norm x = sup { | x ( t ) | , t [ 0 , T ] } . Obviously, ( X , · ) is a Banach space. Next, we set Y = { y ( t ) : y ( t ) C 2 ( [ 0 , T ] , R ) } with the norm y = sup { | y ( t ) | , t [ 0 , T ] } . The product space ( X × Y , ( x , y ) ) is Banach space with the norm ( x , y ) = x + y .
In the following, for brevity, we use the subscript notation
h x , y ( t ) = h ( t , x ( t ) , y ( t ) ) , h { f , g } ,
in fractional integral as
R L I p h x , y ( ϕ ) = 1 Γ ( p ) a ϕ ( ϕ s ) p 1 h ( s , x ( s ) , y ( s ) ) d s ,
where ϕ { t , T , ξ i , η j } . In addition, we use it in multiple fractional integrations.
In view of Lemma 4, we define the operator P : X × Y X × Y by
P ( x , y ) ( t ) = P 1 ( x , y ) ( t ) P 2 ( x , y ) ( t ) ,
where
P 1 ( x , y ) ( t ) = 1 Λ i = 1 m α i R L I β i H I q 2 R L I p 2 g x , y ( ξ i ) + 1 Λ H I q 1 R L I p 1 f x , y ( T ) + Ω 1 Λ H I q 2 R L I p 2 g x , y ( T ) Ω 1 Λ j = 1 k λ j R L I δ j H I q 1 R L I p 1 f x , y ( η j ) + H I q 1 R L I p 1 f x , y ( t )
and
P 2 ( x , y ) ( t ) = Ω 2 Λ i = 1 m α i R L I β i H I q 2 R L I p 2 g x , y ( ξ i ) + Ω 2 Λ H I q 1 R L I p 1 f x , y ( T ) + 1 Λ H I q 2 R L I p 2 g x , y ( T ) 1 Λ j = 1 k λ j R L I δ j H I q 1 R L I p 1 f x , y ( η j ) + H I q 2 R L I p 2 g x , y ( t ) .
For computational convenience, we set
M 1 = 1 + | Λ | | Λ | p 1 q 1 T p 1 Γ ( p 1 + 1 ) + | Ω 1 | | Λ | p 1 q 1 j = 1 k | λ j | η j p 1 + δ j Γ ( p 1 + δ j + 1 ) , M 2 = | Ω 1 | | Λ | p 2 q 2 T p 2 Γ ( p 2 + 1 ) + 1 | Λ | p 2 q 2 i = 1 m | α i | ξ i p 2 + β i Γ ( p 2 + β i + 1 ) , M 3 = | Ω 2 | | Λ | p 1 q 1 T p 1 Γ ( p 1 + 1 ) + 1 | Λ | p 1 q 1 j = 1 k | λ j | η j p 1 + δ j Γ ( p 1 + δ j + 1 ) , M 4 = 1 + | Λ | | Λ | p 2 q 2 T p 2 Γ ( p 2 + 1 ) + | Ω 2 | | Λ | p 2 q 2 i = 1 m | α i | ξ i p 2 + β i Γ ( p 2 + β i + 1 ) .
In the first result, Banach’s contraction mapping principle is used to prove existence and uniqueness of solutions of system (4).
Theorem 1.
Suppose that f , g : [ 0 , T ] × R 2 R are continuous functions. In addition, we assume that f , g satisfies the Lipchitz condition:
( H 1 )
there exist constants m i , n i , i = 1 , 2
| f ( t , u 1 , v 1 ) f ( t , u 2 , v 2 ) | m 1 | u 1 u 2 | + m 2 | v 1 v 2 |
and
| g ( t , u 1 , v 1 ) g ( t , u 2 , v 2 ) | n 1 | u 1 u 2 | + n 2 | v 1 v 2 | ,
for all t [ 0 , T ] and u i , v i R , i = 1 , 2 . Then, the system (4) has a unique solution on [ 0 , T ] , if
( M 1 + M 3 ) ( m 1 + m 2 ) + ( M 2 + M 4 ) ( n 1 + n 2 ) < 1 .
Proof. 
Let us define sup t [ 0 , T ] f ( t , 0 , 0 ) = N 1 < and sup t [ 0 , T ] g ( t , 0 , 0 ) = N 2 < . Choose a constant r > 0 satisfying
r > ( M 1 + M 3 ) N 1 + ( M 2 + M 4 ) N 2 1 [ ( M 1 + M 3 ) ( m 1 + m 2 ) + ( M 2 + M 4 ) ( n 1 + n 2 ) ] .
At first, we shall show that the set P B r B r , where a ball B r = { ( x , y ) X × Y : ( x , y ) r } . For ( x , y ) B r , and using
| f x , y | | f x , y f 0 , 0 | + | f 0 , 0 | m 1 x + m 2 y + N 1 ,
and
| g x , y | | g x , y g 0 , 0 | + | g 0 , 0 | n 1 x + n 2 y + N 2 ,
we get relations
| P 1 ( x , y ) ( t ) | 1 | Λ | i = 1 m | α i | R L I β i H I q 2 R L I p 2 | g x , y | ( ξ i ) + 1 | Λ | H I q 1 R L I p 1 | f x , y | ( T ) + | Ω 1 | | Λ | H I q 2 R L I p 2 | g x , y | ( T ) + | Ω 1 | | Λ | j = 1 k | λ j | R L I δ j H I q 1 R L I p 1 | f x , y | ( η j ) + H I q 1 R L I p 1 | f x , y | ( T ) 1 | Λ | i = 1 m | α i | R L I β i H I q 2 R L I p 2 1 ( ξ i ) ( n 1 x + n 2 y + N 2 ) + 1 | Λ | H I q 1 R L I p 1 1 ( T ) ( m 1 x + m 2 y + N 1 ) + | Ω 1 | | Λ | H I q 2 R L I p 2 1 ( T ) ( n 1 x + n 2 y + N 2 ) + | Ω 1 | | Λ | j = 1 k | λ j | R L I δ j H I q 1 R L I p 1 1 ( η j ) ( m 1 x + m 2 y + N 1 ) + H I q 1 R L I p 1 1 ( T ) ( m 1 x + m 2 y + N 1 ) = 1 | Λ | p 2 q 2 i = 1 m | α i | ξ i p 2 + β i Γ ( p 2 + β i + 1 ) ( n 1 x + n 2 y + N 2 ) + 1 | Λ | p 1 q 1 T p 1 Γ ( p 1 + 1 ) ( m 1 x + m 2 y + N 1 ) + | Ω 1 | | Λ | p 2 q 2 T p 2 Γ ( p 2 + 1 ) ( n 1 x + n 2 y + N 2 ) + | Ω 1 | | Λ | p 1 q 1 j = 1 k | λ j | η j p 1 + δ j Γ ( p 1 + δ j + 1 ) ( m 1 x + m 2 y + N 1 ) + p 1 q 1 T p 1 Γ ( p 1 + 1 ) ( m 1 x + m 2 y + N 1 ) = M 1 ( m 1 x + m 2 y + N 1 ) + M 2 ( n 1 x + n 2 y + N 2 ) = ( M 1 m 1 + M 2 n 1 ) x + ( M 1 m 2 + M 2 n 2 ) y + M 1 N 1 + M 2 N 2 [ M 1 ( m 1 + m 2 ) + M 2 ( n 1 + n 2 ) ] r + M 1 N 1 + M 2 N 2 .
Therefore, we deduce that
P 1 ( x , y ) [ M 1 ( m 1 + m 2 ) + M 2 ( n 1 + n 2 ) ] r + M 1 N 1 + M 2 N 2 .
In a similar way of computation, we get
| P 2 ( x , y ) ( t ) | | Ω 2 | | Λ | i = 1 m | α i | R L I β i H I q 2 R L I p 2 1 ( ξ i ) ( n 1 x + n 2 y + N 2 ) + | Ω 2 | | Λ | H I q 1 R L I p 1 1 ( T ) ( m 1 x + m 2 y + N 1 ) + 1 | Λ | H I q 2 R L I p 2 1 ( T ) ( n 1 x + n 2 y + N 2 ) + 1 | Λ | j = 1 k | λ j | R L I δ j H I q 1 R L I p 1 1 ( η j ) ( m 1 x + m 2 y + N 1 ) + H I q 2 R L I p 2 1 ( T ) ( n 1 x + n 2 y + N 2 ) = | Ω 2 | | Λ | p 2 q 2 i = 1 m | α i | ξ i p 2 + β i Γ ( p 2 + β i + 1 ) ( n 1 x + n 2 y + N 2 ) + | Ω 2 | | Λ | p 1 q 1 T p 1 Γ ( p 1 + 1 ) ( m 1 x + m 2 y + N 1 ) + 1 | Λ | p 2 q 2 T p 2 Γ ( p 2 + 1 ) ( n 1 x + n 2 y + N 2 ) + 1 | Λ | p 1 q 1 j = 1 k | λ j | η j p 1 + δ j Γ ( p 1 + δ j + 1 ) ( m 1 x + m 2 y + N 1 ) + p 2 q 2 T p 2 Γ ( p 2 + 1 ) ( n 1 x + n 2 y + N 2 ) = M 3 ( m 1 x + m 2 y + N 1 ) + M 4 ( n 1 x + n 2 y + N 2 ) ,
which yields
P 2 ( x , y ) [ M 3 ( m 1 + m 2 ) + M 4 ( n 1 + n 2 ) ] r + M 3 N 1 + M 4 N 2 .
Then, we conclude that
P ( x , y ) [ M 1 ( m 1 + m 2 ) + M 2 ( n 1 + n 2 ) ] r + M 1 N 1 + M 2 N 2 + [ M 3 ( m 1 + m 2 ) + M 4 ( n 1 + n 2 ) ] r + M 3 N 1 + M 4 N 2 r ,
which leads to P B r B r .
In the next step, we will show that the P is a contraction operator. For any ( x 1 , y 1 ) , ( x 2 , y 2 ) X × Y , we have
| P 1 ( x 1 , y 1 ) ( t ) P 1 ( x 2 , y 2 ) ( t ) | 1 | Λ | i = 1 m | α i | R L I β i H I q 2 R L I p 2 | g x 1 , y 1 g x 2 , y 2 | ( ξ i ) + 1 | Λ | H I q 1 R L I p 1 | f x 1 , y 1 f x 2 , y 2 | ( T ) + | Ω 1 | | Λ | H I q 2 R L I p 2 | g x 1 , y 1 g x 2 , y 2 | ( T ) + | Ω 1 | | Λ | j = 1 k | λ j | R L I δ j H I q 1 R L I p 1 | f x 1 , y 1 f x 2 , y 2 | ( η j ) + H I q 1 R L I p 1 | f x 1 , y 1 f x 2 , y 2 | ( T ) 1 | Λ | i = 1 m | α i | R L I β i H I q 2 R L I p 2 1 ( ξ i ) ( n 1 x 1 x 2 + n 2 y 1 y 2 ) + 1 | Λ | H I q 1 R L I p 1 1 ( T ) ( m 1 x 1 x 2 + m 2 y 1 y 2 ) + | Ω 1 | | Λ | H I q 2 R L I p 2 1 ( T ) ( n 1 x 1 x 2 + n 2 y 1 y 2 ) + | Ω 1 | | Λ | j = 1 k | λ j | R L I δ j H I q 1 R L I p 1 1 ( η j ) ( m 1 x 1 x 2 + m 2 y 1 y 2 ) + H I q 1 R L I p 1 1 ( T ) ( m 1 x 1 x 2 + m 2 y 1 y 2 ) = M 1 ( m 1 x 1 x 2 + m 2 y 1 y 2 ) + M 2 ( n 1 x 1 x 2 + n 2 y 1 y 2 ) = ( M 1 m 1 + M 2 n 1 ) x 1 x 2 + ( M 1 m 2 + M 2 n 2 ) y 1 y 2 .
Then, we get the result that
P 1 ( x 1 , y 1 ) P 1 ( x 2 , y 2 ) M 1 ( m 1 + m 2 ) + M 2 ( n 1 + n 2 ) x 1 x 2 + y 1 y 2 .
In addition, we have
| P 2 ( x 1 , y 1 ) ( t ) P 2 ( x 2 , y 2 ) ( t ) | | Ω 2 | | Λ | i = 1 m | α i | R L I β i H I q 2 R L I p 2 1 ( n 1 x 1 x 2 + n 2 y 1 y 2 ) + | Ω 2 | | Λ | H I q 1 R L I p 1 1 ( T ) ( m 1 x 1 x 2 + m 2 y 1 y 2 ) + 1 | Λ | H I q 2 R L I p 2 1 ( T ) ( n 1 x 1 x 2 + n 2 y 1 y 2 ) + 1 | Λ | j = 1 k | λ j | R L I δ j H I q 1 R L I p 1 1 ( η j ) ( m 1 x 1 x 2 + m 2 y 1 y 2 ) + H I q 2 R L I p 2 1 ( T ) ( n 1 x 1 x 2 + n 2 y 1 y 2 ) = M 3 ( m 1 x 1 x 2 + m 2 y 1 y 2 ) + M 4 ( n 1 x 1 x 2 + n 2 y 1 y 2 ) = ( M 3 m 1 + M 4 n 1 ) x 1 x 2 + ( M 3 m 2 + M 4 n 2 ) y 1 y 2 ,
which yields
P 2 ( x 1 , y 1 ) P 2 ( x 2 , y 2 ) M 3 ( m 1 + m 2 ) + M 4 ( n 1 + n 2 ) x 1 x 2 + y 1 y 2 .
The above results in (18) and (19) imply
P ( x 1 , y 1 ) P ( x 2 , y 2 ) [ ( M 1 + M 3 ) ( m 1 + m 2 ) + ( M 2 + M 4 ) ( n 1 + n 2 ) ] × x 1 x 2 + y 1 y 2 .
Since ( M 1 + M 3 ) ( m 1 + m 2 ) + ( M 2 + M 4 ) ( n 1 + n 2 ) < 1 , then the operator P is a contraction. From the benefits of Banach’s fixed point theorem, the operator P has a unique fixed point, which is the unique solution of (4) on [ 0 , T ] . The proof is completed. □
The Leray–Schauder alternative is applied to our second existence result.
Lemma 6.
(Leray–Schauder alternative) [31]. Let Q : U U be a completely continuous operator. Let
μ ( Q ) = { x U : x = θ Q ( x ) f o r s o m e 0 < θ < 1 } .
Then, either the set μ ( Q ) is unbounded, or Q has at least one fixed point.
Theorem 2.
Suppose that there exist constants a r , b r 0 for r = 1 , 2 and a 0 , b 0 > 0 . In addition, for any u , v R , we assume that
| f ( t , u , v ) | a 0 + a 1 | u | + a 2 | v | , | g ( t , u , v ) | b 0 + b 1 | u | + b 2 | v | .
If ( M 1 + M 3 ) a 1 + ( M 2 + M 4 ) b 1 < 1 and ( M 1 + M 3 ) a 2 + ( M 2 + M 4 ) b 2 < 1 , then (4) has at least one solution on [ 0 , T ] .
Proof. 
The first task of the proof is to show that the operator P : X × Y X × Y is completely continuous. The continuity of the functions f , g on [ 0 , T ] × R × R can be used to claim that the operator P is continuous. Now, we let Φ be the bounded subset of X × Y . Then, there exist positive constants G 1 and G 2 such that
| f ( t , x , y ) | G 1 , | g ( t , x , y ) | G 2 , ( x , y ) Φ .
For any ( x , y ) Φ , we have
| P 1 ( x , y ) ( t ) | 1 | Λ | i = 1 m | α i | R L I β i H I q 2 R L I p 2 | g x , y | ( ξ i ) + 1 | Λ | H I q 1 R L I p 1 | f x , y | ( T ) + | Ω 1 | | Λ | H I q 2 R L I p 2 | g x , y | ( T ) + | Ω 1 | | Λ | j = 1 k | λ j | R L I δ j H I q 1 R L I p 1 | f x , y | ( η j ) + H I q 1 R L I p 1 | f x , y | ( T ) 1 | Λ | p 2 q 2 i = 1 m | α i | ξ i p 2 + β i Γ ( p 2 + β i + 1 ) G 2 + 1 | Λ | p 1 q 1 T p 1 Γ ( p 1 + 1 ) G 1 + | Ω 1 | | Λ | p 2 q 2 T p 2 Γ ( p 2 + 1 ) G 2 + | Ω 1 | | Λ | p 1 q 1 j = 1 k | λ j | η j p 1 + δ j Γ ( p 1 + δ j + 1 ) G 1 + p 1 q 1 T p 1 Γ ( p 1 + 1 ) G 1 ,
which leads to
P 1 ( x , y ) G 1 M 1 + G 2 M 2 .
Furthermore, we get
P 2 ( x , y ) | Ω 2 | | Λ | i = 1 m | α i | R L I β i H I q 2 R L I p 2 1 G 2 + | Ω 2 | | Λ | H I q 1 R L I p 1 1 ( T ) G 1 + 1 | Λ | H I q 2 R L I p 2 1 ( T ) G 2 + 1 | Λ | j = 1 k | λ j | R L I δ j H I q 1 R L I p 1 1 ( η j ) G 1 + H I q 2 R L I p 2 1 ( T ) G 2 = G 1 M 3 + G 2 M 4 .
Therefore, from above two results, we deduce that the set P Φ is uniformly bounded. The next is to prove that the set P Φ is equicontinuous. Choosing two points τ 1 , τ 2 [ 0 , T ] such that τ 1 < τ 2 , we have, for any ( x , y ) Φ , that
| P 1 ( x , y ) ( τ 2 ) P 1 ( x , y ) ( τ 1 ) | = H I q 1 R L I p 1 f x , y ( τ 2 ) H I q 1 R L I p 1 f x , y ( τ 1 ) G 1 H I q 1 R L I p 1 1 ( τ 2 ) H I q 1 R L I p 1 1 ( τ 1 ) = G 1 p 1 q 1 Γ ( p 1 + 1 ) τ 2 p 1 τ 1 p 1 ,
which implies
| P 1 ( x , y ) ( τ 2 ) P 1 ( x , y ) ( τ 1 ) | 0 , as τ 1 τ 2 .
In addition, we obtain
| P 2 ( x , y ) ( τ 2 ) P 2 ( x , y ) ( τ 1 ) | = H I q 2 R L I p 2 g x , y ( τ 2 ) H I q 2 R L I p 2 g x , y ( τ 1 ) G 2 H I q 2 R L I p 2 1 ( τ 2 ) H I q 2 R L I p 2 1 ( τ 1 ) = G 2 p 2 q 2 Γ ( p 2 + 1 ) τ 2 p 2 τ 1 p 2 .
Then,
| P 2 ( x , y ) ( τ 2 ) P 2 ( x , y ) ( τ 1 ) | 0 , as τ 1 τ 2 .
Thus, the set P Φ is equicontinuous. By taking into account the Arzelá-Ascoli theorem, the set P Φ is relatively compact. Then, operator P is completely continuous.
Finally, we will claim that the set μ = { ( x , y ) X × Y : ( x , y ) = θ P ( x , y ) , 0 θ 1 } is bounded. For any ( x , y ) μ , then ( x , y ) = θ P ( x , y ) . Hence, for t [ a , b ] , we have
x ( t ) = θ P 1 ( x , y ) ( t ) and y ( t ) = θ P 2 ( x , y ) ( t ) .
Therefore, we obtain
x ( a 0 + a 1 x + a 2 y ) M 1 + ( b 0 + b 1 x + b 2 y ) M 2 , y ( a 0 + a 1 x + a 2 y ) M 3 + ( b 0 + b 1 x + b 2 y ) M 4 ,
which lead to
x + y ( M 1 + M 3 ) a 0 + ( M 2 + M 4 ) b 0 + ( M 1 + M 3 ) a 1 + ( M 2 + M 4 ) b 1 x + ( M 1 + M 3 ) a 2 + ( M 2 + M 4 ) b 2 y .
Thus, the following inequality holds:
( x , y ) ( M 1 + M 3 ) a 0 + ( M 2 + M 4 ) b 0 M * ,
where M * = min { 1 ( M 1 + M 3 ) a 1 ( M 2 + M 4 ) b 1 , 1 ( M 1 + M 3 ) a 2 ( M 2 + M 4 ) b 2 } . Hence, the set μ is a bounded set. Then, by using Lemma 6, the operator P has at least one fixed point. Therefore, we conclude that problem (4) has at least one solution on [ 0 , T ] . The proof is complete. □
If a r , b r = 0 , r = 1 , 2 , in Theorem 2, we have following corollary.
Corollary 2.
Assume that | f ( t , x , y ) | a 0 and | g ( t , x , y ) | b 0 , where a 0 , b 0 > 0 , ( t , x , y ) [ 0 , T ] × R 2 . Then, problem (4) has at least one solution on [ 0 , T ] .
Next, we present examples to illustrate our results.
Example 1.
Consider the following sequential Riemann–Liouville and Hadamard–Caputo fractional differential system with coupled fractional integral boundary conditions of the form
R L D 1 5 H C D 4 5 x ( t ) = f ( t , x ( t ) , y ( t ) ) , t 0 , 7 / 4 , R L D 2 5 H C D 3 5 y ( t ) = g ( t , x ( t ) , y ( t ) ) , t 0 , 7 / 4 , H C D 4 5 x ( 0 ) = 0 , x 7 4 = 1 3 R L I 3 4 y 1 2 + 2 7 R L I 5 4 y 5 4 , H C D 3 5 y ( 0 ) = 0 , y 7 4 = 3 11 R L I 1 2 x 1 4 + 4 17 R L I 7 8 x 3 4 + 5 19 R L I 11 8 x 3 2 .
Here, p 1 = 1 / 5 , p 2 = 2 / 5 , q 1 = 4 / 5 , q 2 = 3 / 5 , T = 7 / 4 , m = 2 , α 1 = 1 / 3 , α 2 = 2 / 7 , β 1 = 3 / 4 , β 2 = 5 / 4 , ξ 1 = 1 / 2 , ξ 2 = 5 / 4 , k = 3 , λ 1 = 3 / 11 , λ 2 = 4 / 17 , λ 3 = 5 / 19 , δ 1 = 1 / 2 , δ 2 = 7 / 8 , δ 3 = 11 / 8 , η 1 = 1 / 4 , η 2 = 3 / 4 , η 3 = 3 / 2 . Form all constants, we find that Ω 1 0.5489581728 , Ω 2 0.7217268652 , | Λ | 0.6038021388 , M 1 13.82028787 , M 2 3.420721316 , M 3 9.093047627 , M 4 7.354860071 .
Let the two nonlinear Lipschitzian functions f, g : [ 0 , 7 / 4 ] × R 2 R be defined by
f ( t , x , y ) = 1 12 ( t + 12 ) x 2 + 2 | x | 1 + | x | + e t sin y 15 ( 3 t + 5 ) + 1 2 ,
g ( t , x , y ) = cos π t 6 ( 2 t + 9 ) tan 1 x + 1 36 ( 4 t + 7 ) 3 y 2 + 4 | y | 1 + | y | + 3 4 .
From (22)–(23), we see that
| f ( t , x 1 , y 1 ) f ( t , x 2 , y 2 ) | 1 72 | x 1 x 2 | + 1 75 | y 1 y 2 |
and
| g ( t , x 1 , y 1 ) g ( t , x 2 , y 2 ) | 1 54 | x 1 x 2 | + 1 63 | y 1 y 2 | ,
for all x r , y r R , r = 1 , 2 , we obtain ( M 1 + M 3 ) ( 1 / 72 + 1 / 75 ) + ( M 2 + M 4 ) ( 1 / 54 + 1 / 63 ) 0.9943406888 < 1 . From the benefits of Theorem 1, the problem of a sequential Riemann–Liouville and Hadamard–Caputo fractional differential system with coupled fractional integral boundary conditions (21) with f and g given by (22)–(23), respectively, has a unique solution on [ 0 , 7 / 4 ] .
Example 2.
Consider the sequential Riemann–Liouville and Hadamard–Caputo fractional differential system with coupled fractional integral boundary conditions of the Example 1, where the nonlinear functions f, g : [ 0 , 7 / 4 ] × R 2 R are defined by
f ( t , x , y ) = 2 e t 13 + 1 2 ( 5 t + 23 ) x 16 1 + | x | 15 + cos π t 3 ( 2 t + 15 ) y sin 2 x ,
g ( t , x , y ) = 4 t 3 + x e y 2 2 ( 4 t + 11 ) + | y | 19 cos 4 x 3 ( 3 t + 8 ) ( 1 + y 18 ) .
It is easy to obtain that | f ( t , x , y ) | ( 2 / 13 ) + ( 1 / 46 ) | x | + ( 1 / 45 ) | y | and | g ( t , x , y ) | ( 7 / 3 ) + ( 1 / 22 ) | x | + ( 1 / 24 ) | y | . By setting a 0 = 2 / 13 , a 1 = 1 / 46 , a 2 = 1 / 45 , b 0 = 7 / 3 , b 1 = 1 / 22 and b 2 = 1 / 24 , we can find that ( M 1 + M 3 ) a 1 + ( M 2 + M 4 ) b 1 0.9879151432 < 1 and ( M 1 + M 3 ) a 2 + ( M 2 + M 4 ) b 2 0.9581677912 < 1 . The conclusion of Theorem 2 can be implied that system (21) with f and g given by (24)–(25), respectively, has at least one solution on [ 0 , 7 / 4 ] .
Example 3.
Consider the sequential Riemann–Liouville and Hadamard–Caputo fractional differential system with coupled fractional integral boundary conditions of the Example 1, where the nonlinear functions f, g : [ 0 , 7 / 4 ] × R 2 R are given by
f ( t , x , y ) = 1 2 ( 1 + cos 2 t ) + | x | e t ( 1 + | x | ) + 2 π tan 1 y ,
g ( t , x , y ) = 1 4 ( 3 + sin 2 π t ) + e x 4 + 3 y 22 1 + y 22 .
We can check that | f ( t , x , y ) | 3 , | g ( t , x , y ) | 5 for all x , y R . Using the Corollary 2, the problem (21) with f and g given by (26) and (27), respectively, has at least one solution on [ 0 , 7 / 4 ] .

4. Conclusions

In this paper, we studied a new system of sequential fractional differential equations which consists of mixed fractional derivatives of Riemann–Liouville and Hadamard–Caputo types, supplemented with nonlocal coupled fractional integral boundary conditions. To the best of our knowledge, this is the first system of this type that appeared in the literature. After proving a basic lemma, helping us to transform the considered system into a fixed point problem, we use the standard tools from functional analysis to establish existence and uniqueness results. We use a Banach contraction mapping principle to derive the uniqueness result and Leray–Schauder alternative to obtain an existence result. The obtained results are well illustrated by numerical examples. The obtained results enrich the existing literature on sequential systems of fractional differential equations. Other cases of fractional systems with other types of mixed fractional derivatives or other types of boundary conditions can be studied using the methodology of this paper.

Author Contributions

Conceptualization, C.K., W.Y., S.K.N., and J.T.; methodology, C.K., W.Y., S.K.N., and J.T.; formal analysis, C.K., W.Y., S.K.N., and J.T.; funding acquisition, J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by King Mongkut’s University of Technology North Bangkok, Contract No. KMUTNB-61-KNOW-034.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Diethelm, K. The Analysis of Fractional Differential Equations; Lecture Notes in Mathematics; Springer: New York, NY, USA, 2010. [Google Scholar]
  2. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of the Fractional Differential Equations; North-Holland Mathematics Studies; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  3. VLakshmikantham; Leela, S.; Devi, J.V. Theory of Fractional Dynamic Systems; Cambridge Scientific Publishers: Cottenham, UK, 2009. [Google Scholar]
  4. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Differential Equations; John Wiley: NewYork, NY, USA, 1993. [Google Scholar]
  5. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  6. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives; Gordon and Breach Science: Yverdon, Switzerland, 1993. [Google Scholar]
  7. Ahmad, B.; Alsaedi, A.; Ntouyas, S.K.; Tariboon, J. Hadamard-Type Fractional Differential Equations, Inclusions and Inequalities; Springer: Cham, Switzerland, 2017. [Google Scholar]
  8. Ahmad, B.; Ntouyas, S.K. Nonlocal Nonlinear Fractional-Order Boundary Value Problems; World Scientific: Singapore, 2021. [Google Scholar]
  9. Zhou, Y. Basic Theory of Fractional Differential Equations; World Scientific: Singapore, 2014. [Google Scholar]
  10. Tarasov, V.E. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields, and Media; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  11. Uchaikin, V.; Sibatov, R. Fractional Kinetics in Space: Anomalous Transport Models; World Scientific: Singapore, 2018. [Google Scholar]
  12. Ding, Y.; Wang, Z.; Ye, H. Optimal control of a fractional-order HIV-immune system with memory. IEEE Trans. Control. Syst. Technol. 2012, 20, 763–769. [Google Scholar] [CrossRef]
  13. Faieghi, M.; Kuntanapreeda, S.; Delavari, H.; Baleanu, D. LMI-based stabilization of a class of fractional-order chaotic systems. Nonlinear Dynam. 2013, 72, 301–309. [Google Scholar] [CrossRef]
  14. Zhang, F.; Li, G.C.C.; Kurths, J. Chaos synchronization in fractional differential systems. Phys. Eng. Sci. 2013, 371, 20120155. [Google Scholar] [CrossRef] [PubMed]
  15. Javidi, M.; Ahmad, B. Dynamic analysis of time fractional order phytoplankton-toxic phytoplankton–zooplankton system. Ecol. Model. 2015, 318, 8–18. [Google Scholar] [CrossRef]
  16. Carvalho, A.; Pinto, C.M.A. A delay fractional order model for the co-infection of malaria and HIV/AIDS. Int. J. Dyn. Control 2017, 5, 168–186. [Google Scholar] [CrossRef]
  17. Torvik, P.J.; Bagley, R.L. On the appearance of the fractional derivative in the behavior of real materials. J. Appl. Mech. 1984, 51, 294–298. [Google Scholar] [CrossRef]
  18. Mainardi, F. Some basic problems in continuum and statistical mechanics. In Fractals and Fractional Calculus in Continuum Mechanics; Carpinteri, A., Mainardi, F., Eds.; Springer: Berlin/Heidelberg, Germany, 1997; pp. 291–348. [Google Scholar]
  19. Yin, C.; Liu, F.; Anh, V. Numerical simulation of the nonlinear fractional dynamical systems with fractional damping for the extensible and inextensible pendulum. J. Algorithm Comput. Technol. 2007, 1, 427–447. [Google Scholar] [CrossRef] [Green Version]
  20. Ahmad, B.; Ntouyas, S.K. Existence and uniqueness of solutions for Caputo-Hadamard sequential fractional order neutral functional differential equations. Electron. J. Differ. Equ. 2017, 36, 1–11. [Google Scholar]
  21. Promsakon, C.; Phuangthong, N.; Ntouyas, S.K.; Tariboon, J. Nonlinear sequential Riemann–Liouville and Caputo fractional differential equations with generalized fractional integral conditions. Adv. Differ. Equ. 2018, 2018, 385. [Google Scholar] [CrossRef] [Green Version]
  22. Tariboon, J.; Cuntavepanit, A.; Ntouyas, S.K.; Nithiarayaphaks, W. Separated boundary value problems of sequential Caputo and Hadamard fractional differential equations. J. Funct. Spaces 2018, 2018, 6974046. [Google Scholar] [CrossRef]
  23. Asawasamrit, S.; Phuangthong, N.; Ntouyas, S.K.; Tariboon, J. Nonlinear sequential Riemann–Liouville and Caputo fractional differential equations with nonlocal and integral boundary conditions. Int. J. Anal. Appl. 2019, 17, 47–63. [Google Scholar]
  24. Asawasamrit, S.; Ntouyas, S.K.; Tariboon, J.; Nithiarayaphaks, W. Coupled systems of sequential Caputo and Hadamard fractional differential equations with coupled separated boundary conditions. Symmetry 2018, 10, 701. [Google Scholar] [CrossRef] [Green Version]
  25. Ahmad, B.; Agarwal, R.P.; Alsaedi, A.; Ntouyas, S.K.; Alruwaily, Y. Fractional order coupled systems for mixed fractional derivatives with nonlocal multi-point and Riemann-Stieltjes integral-multi-strip conditions. Dynam. Syst. Appl. 2020, 29, 71–86. [Google Scholar]
  26. Henderson, J.; Luca, R.; Tudorache, A. On a system of fractional differential equations with coupled integral boundary conditions. Fract. Calc. Appl. Anal. 2015, 18, 361–386. [Google Scholar] [CrossRef]
  27. Wang, J.R.; Zhang, Y. Analysis of fractional order differential coupled systems. Math. Methods Appl. Sci. 2015, 38, 3322–3338. [Google Scholar] [CrossRef]
  28. Ahmad, B.; Luca, R. Existence of solutions for a sequential fractional integro-differential system with coupled integral boundary conditions. Chaos Solitons Fractals 2017, 104, 378–388. [Google Scholar] [CrossRef]
  29. Tariboon, J.; Ntouyas, S.K.; Ahmad, B.; Alsaedi, A. Existence results for sequential Riemann–Liouville and Caputo fractional differential inclusions with generalized fractional integral conditions. Mathematics 2020, 8, 1044. [Google Scholar] [CrossRef]
  30. Jarad, F.; Abdeljawad, T.; Baleanu, D. Caputo-type modification of the Hadamard fractional derivatives. Adv. Differ. Equ. 2012, 2012, 142. [Google Scholar] [CrossRef] [Green Version]
  31. Granas, A.; Dugundji, J. Fixed Point Theory; Springer: New York, NY, USA, 2003. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kiataramkul, C.; Yukunthorn, W.; Ntouyas, S.K.; Tariboon, J. Sequential Riemann–Liouville and Hadamard–Caputo Fractional Differential Systems with Nonlocal Coupled Fractional Integral Boundary Conditions. Axioms 2021, 10, 174. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms10030174

AMA Style

Kiataramkul C, Yukunthorn W, Ntouyas SK, Tariboon J. Sequential Riemann–Liouville and Hadamard–Caputo Fractional Differential Systems with Nonlocal Coupled Fractional Integral Boundary Conditions. Axioms. 2021; 10(3):174. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms10030174

Chicago/Turabian Style

Kiataramkul, Chanakarn, Weera Yukunthorn, Sotiris K. Ntouyas, and Jessada Tariboon. 2021. "Sequential Riemann–Liouville and Hadamard–Caputo Fractional Differential Systems with Nonlocal Coupled Fractional Integral Boundary Conditions" Axioms 10, no. 3: 174. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms10030174

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