Next Article in Journal
An Iterative Method Based on the Marginalized Particle Filter for Nonlinear B-Spline Data Approximation and Trajectory Optimization
Next Article in Special Issue
On Separate Fractional Sum-Difference Equations with n-Point Fractional Sum-Difference Boundary Conditions via Arbitrary Different Fractional Orders
Previous Article in Journal
Refined Neutrosophy and Lattices vs. Pair Structures and YinYang Bipolar Fuzzy Set
Previous Article in Special Issue
Hyers–Ulam Stability and Existence of Solutions for Differential Equations with Caputo–Fabrizio Fractional Derivative
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability, Existence and Uniqueness of Boundary Value Problems for a Coupled System of Fractional Differential Equations

Eastern Mediterranean University, Gazimagusa 99628, T.R. North Cyprus, Mersin 10, Turkey
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(4), 354; https://doi.org/10.3390/math7040354
Submission received: 20 February 2019 / Revised: 10 April 2019 / Accepted: 10 April 2019 / Published: 16 April 2019

Abstract

:
The current article studies a coupled system of fractional differential equations with boundary conditions and proves the existence and uniqueness of solutions by applying Leray-Schauder’s alternative and contraction mapping principle. Furthermore, the Hyers-Ulam stability of solutions is discussed and sufficient conditions for the stability are developed. Obtained results are supported by examples and illustrated in the last section.

1. Introduction

Fractional calculus is undoubtedly one of the very fast-growing fields of modern mathematics, due to its broad range of applications in various fields of science and its unique efficiency in modeling complex phenomena [1,2]. In particular, fractional differential equations with boundary conditions are widely employed to build complex mathematical models for numerous real-life problems such as blood flow problem, underground water flow, population dynamics, and bioengineering. As an example, consider the following equation that describes a thermostat model
x = g ( t ) f ( t , x ) , x ( 0 ) = 0 , β x ( 1 ) = x ( η ) ,
where t ( 0 , 1 ) , η ( 0 , 1 ] and β is a positive constant. Note that solutions of the above equation with the specified integral boundary conditions are in fact solutions of the one-dimensional heat equation describing a heated bar with a controller at point 1, which increases or reduces heat based on the temperature picked by a sensor at η. A few of the relevant studies on coupled systems of fractional differential equations with integral boundary conditions are briefly reviewed below and for further information on this topic, refer to References [3,4].
In Reference [5], Ntouyas and Obaid used Leray-Schauder’s alternative and Banach’s fixed-point theorem to prove the existence and uniqueness of solutions for the following coupled fractional differential equations with Riemann-Liouville integral boundary conditions:
{ D   c 0 + α u ( t ) = g ( t , u ( t ) , v ( t ) ) , t [ 0 , 1 ] , D   c 0 + β v ( t ) = g ( t , u ( t ) , v ( t ) ) , t [ 0 , 1 ] , u ( 0 ) = γ I p u ( η ) = γ 0 η   ( η s ) p 1 Γ ( p ) u ( s ) d s , 0 < η < 1 , v ( 0 ) = δ I q v ( ζ ) = δ 0 ζ   ( ζ s ) q 1 Γ ( q ) v ( s ) d s , 0 < ζ < 1 .
Here, D   c 0 + α and D   c 0 + β are Caputo fractional derivatives, 0 < α , β 1 , f , g C ( [ 0 , 1 ] × 2 , ) and p , q , γ , δ .
Similarly, Ahmed and Ntouyas [6] employed Banach fixed-point theorem and Leray-Schauder’s alternative to prove the existence and uniqueness of solutions for the following coupled fractional differential system:
{ D   c q x ( t ) = f ( t , x ( t ) , y ( t ) ) ,      t [ 0 , 1 ] ,     1 < q 2 , D   c p y ( t ) = g ( t , x ( t ) , y ( t ) ) ,      t [ 0 , 1 ] ,     1 < q 2 ,
supplemented with coupled and uncoupled slit-strips-type integral boundary conditions, respectively, given by
{ x ( 0 ) = 0 ,     x ( ζ ) = a 0 η y ( s ) d s + b ξ 1 y ( s ) d s ,     0 < η < ζ < ξ < 1 , y ( 0 ) = 0 ,     y ( ζ ) = a 0 η x ( s ) d s + b ξ 1 x ( s ) d s ,     0 < η < ζ < ξ < 1 ,
and
{ x ( 0 ) = 0 ,     x ( ζ ) = a 0 η x ( s ) d s + b ξ 1 x ( s ) d s ,     0 < η < ζ < ξ < 1 , y ( 0 ) = 0 ,     y ( ζ ) = a 0 η y ( s ) d s + b ξ 1 y ( s ) d s ,     0 < η < ζ < ξ < 1 .
Furthermore, Alsulami et al. [7] investigated the following coupled system of fractional differential equations:
{   c D α x ( t ) = f ( t , x ( t ) , y ( t ) ) , t [ 0 , T ] , 1 < α 2 ,   c D β y ( t ) = g ( t , x ( t ) , y ( t ) ) , t [ 0 , T ] , 1 < β 2 ,
subject to the following non-separated coupled boundary conditions:
{ x ( 0 ) = λ 1 y ( T ) , x ( 0 ) = λ 2 y ( T ) , y ( 0 ) = μ 1 x ( T ) , y ( 0 ) = μ 2 x ( T ) .
Note that   c D α and   c D β denote Caputo fractional derivatives of order α and β. Moreover, λ i , μ i , i = 1 ,   2 ,   are real constants with λ i μ i 1 and f , g : [ 0 , T ] × × are appropriately chosen functions. For further details on this topic, refer to References [8,9,10,11,12,13,14,15,16,17,18,19,20,21].
The current paper studies the following coupled system of nonlinear fractional differential equations:
{ D   c α x ( t ) = f ( t , x ( t ) , y ( t ) ) ,      t [ 0 , T ] ,     1 < α 2 ,     c D β y ( t ) = g ( t , x ( t ) , y ( t ) ) ,      t [ 0 , T ] ,     1 < β 2 ,        
supplemented with boundary conditions of the form:
x ( T ) = η y ( ρ ) ,     y ( T ) = ζ x ( μ ) ,     x ( 0 ) = 0 ,     y ( 0 ) = 0 , ρ , μ [ 0 , T ]    
Here,   c D k denotes Caputo fractional derivative of order k   ( k = α , β ) ; and f , g C ( [ 0 , T ] × 2 , ) are given continuous functions. Note that η , ζ are real constants such that T 2 η ζ 0 .
The rest of this paper is organized in the following manner: In Section 2, we briefly review some of the relevant definitions from fractional calculus and prove an auxiliary lemma that will be used later. Section 3 deals with proving the existence and uniqueness of solutions for the given problem, and Section 4 discusses the Hyers-Ulam stability of solutions and presents sufficient conditions for the stability. The paper concludes with supporting examples and obtained results.

2. Preliminaries

We begin this section by reviewing the definitions of fractional derivative and integral [1,2].
Definition 1.
The Riemann-Liouville fractional integral of order τ for a continuous function h is given by
I τ h ( s ) = 1 Γ ( τ ) 0 s h ( t ) ( s t ) 1 τ d t ,      τ > 0 ,
provided that the right-hand side is point-wise defined on [ 0 , ) .
Definition 2.
The Caputo fractional derivatives of order τ for ( h 1 ) —times absolutely continuous function g : [ 0 , ) is defined as
  c D τ g ( s ) = 1 Γ ( h τ ) 0 s ( s t ) h τ 1 g ( h ) ( t ) d t ,      h 1 < τ < h ,      h = [ τ ] + 1 ,
where [ τ ] is the integer part of real number τ .
Here we prove the following auxiliary lemma that will be used in the next section.
Lemma 1.
Let u , v C ( [ 0 , T ] , ) then the unique solution for the problem
{ D   c α x ( t ) = u ( t ) ,     t [ 0 , T ] ,     1 < α 2 , D   c β y ( t ) = v ( t ) ,     t [ 0 , T ] ,     1 < β 2 ,    x ( T ) = η y ( ρ ) ,    y ( T ) = ζ x ( μ ) ,    x ( 0 ) = 0 ,    y ( 0 ) = 0 , ρ , μ [ 0 , T ]   
is
x ( t ) = t Δ ( η T 0 ρ ( ρ s ) β 2 Γ ( β 1 ) v ( s ) d s T 0 T ( T s ) α 1 Γ ( α ) u ( s ) d s + η ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) u ( s ) d s η 0 T ( T s ) β 1 Γ ( β ) v ( s ) d s ) + 0 t ( t s ) α 1 Γ ( α ) u ( s ) d s ,   
and
y ( t ) = t Δ ( η ζ 0 ρ ( ρ s ) β 2 Γ ( β 1 ) v ( s ) d s ζ 0 T ( T s ) α 1 Γ ( α ) u ( s ) d s + T ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) u ( s ) T 0 T ( T s ) β 1 Γ ( β ) v ( s ) d s ) + 0 t ( t s ) β 1 Γ ( β ) v ( s ) d s
where Δ = T 2 η ζ 0 .
Proof. 
General solutions of the fractional differential equations in (3) are known [6] as
x ( t ) = a t + b + 1 Γ ( α ) 0 t ( t s ) α 1 u ( s ) d s , y ( t ) = c t + d + 1 Γ ( β ) 0 t ( t s ) β 1 v ( s ) d s ,
where a , b , c , and   d are arbitrary constants.
Apply conditions x ( 0 ) = 0   and y ( 0 ) = 0 , and we obtain b = d = 0 .
Here
x ( t ) = a + 1 Γ ( α 1 ) 0 t ( t s ) α 2 u ( s ) d s ,
y ( t ) = c + 1 Γ ( β 1 ) 0 t ( t s ) β 2 v ( s ) d s .
Considering boundary conditions
x ( T ) = η y ( ρ ) ,      y ( T ) = ζ x ( μ )
we get
a T + 0 T ( T s ) α 1 Γ ( α ) u ( s ) d s = η c + η 0 ρ ( ρ s ) β 2 Γ ( β 1 ) v ( s ) d s ,
and
c T + 0 T ( T s ) β 1 Γ ( β ) v ( s ) d s = a ζ + ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) u ( s ) d s ,
so
a = 1 T ( η c + η 0 ρ ( ρ s ) β 2 Γ ( β 1 ) v ( s ) d s 0 T ( T s ) α 1 Γ ( α ) u ( s ) d s ) ,
c = 1 T ( a ζ + ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) u ( s ) d s 0 T ( T s ) β 1 Γ ( β ) v ( s ) d s ) .
Hence, by substituting the value of a into c, we obtain the final result for these constants as
c = 1 T ( ζ T [ η c + η 0 ρ ( ρ s ) β 2 Γ ( β 1 ) v ( s ) d s 0 T ( T s ) α 1 Γ ( α ) u ( s ) d s ] + ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) u ( s ) d s 0 T ( T s ) β 1 Γ ( β ) v ( s ) d s ) ,
c ζ η c T 2 = 1 T ( ζ T [ η 0 ρ ( ρ s ) β 2 Γ ( β 1 ) v ( s ) d s 0 T ( T s ) α 1 Γ ( α ) u ( s ) d s ] + ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) u ( s ) d s 0 T ( T s ) β 1 Γ ( β ) v ( s ) d s ) ,
c ( T 2 ζ η T 2 ) = 1 T ( ζ T [ η 0 ρ ( ρ s ) β 2 Γ ( β 1 ) v ( s ) d s 0 T ( T s ) α 1 Γ ( α ) u ( s ) d s ] + ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) u ( s ) d s 0 T ( T s ) β 1 Γ ( β ) v ( s ) d s ) ,
c = T T 2 ζ η ( ζ T [ η 0 ρ ( ρ s ) β 2 Γ ( β 1 ) v ( s ) d s 0 T ( T s ) α 1 Γ ( α ) u ( s ) d s ] + ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) u ( s ) d s 0 T ( T s ) β 1 Γ ( β ) v ( s ) d s ) ,
c = 1 T 2 ζ η ( η ζ 0 ρ ( ρ s ) β 2 Γ ( β 1 ) v ( s ) d s ζ 0 T ( T s ) α 1 Γ ( α ) u ( s ) d s + T ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) u ( s ) d s T 0 T ( T s ) β 1 Γ ( β ) v ( s ) d s )
c = 1 Δ ( η ζ 0 ρ ( ρ s ) β 2 Γ ( β 1 ) v ( s ) d s ζ 0 T ( T s ) α 1 Γ ( α ) u ( s ) d s + T ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) u ( s ) d s T 0 T ( T s ) β 1 Γ ( β ) v ( s ) d s ) ,
and
a = 1 Δ ( η T 0 ρ ( ρ s ) β 2 Γ ( β 1 ) v ( s ) d s T 0 T ( T s ) α 1 Γ ( α ) u ( s ) d s + η ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) u ( s ) d s η 0 T ( T s ) β 1 Γ ( β ) v ( s ) d s ) ,
Substituting the values of a , b , c , and   d in (6) and (7) we get (4) and (5). The converse follows by direct computation. This completes the proof. □

3. Existence and Uniqueness of Solutions

Consider the space C ( [ 0 , T ] , )   endowed with norm x = s u p 0 t T | x ( t ) | . Consequently, the product space C ( [ 0 , T ] , ) ×   C ( [ 0 , T ] , )   is a Banach Space (endowed with ( x , y ) = x + y ).
In view of Lemma 1, we define the operator G : C ( [ 0 , T ] , ) ×   C ( [ 0 , T ] , ) C ( [ 0 , T ] , ) ×   C ( [ 0 , T ] , ) as:
G ( x , y ) ( t ) = ( G 1 ( x , y ) ( t ) , G 2 ( x , y ) ( t ) ) ,
where
G 1 ( x , y ) ( t ) = t Δ ( η T 0 ρ ( ρ s ) β 2 Γ ( β 1 ) g ( s , x ( s ) , y ( s ) ) d s T 0 T ( T s ) α 1 Γ ( α ) f ( s , x ( s ) , y ( s ) ) d s + η ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) f ( s , x ( s ) , y ( s ) ) d s η 0 T ( T s ) β 1 Γ ( β ) g ( s , x ( s ) , y ( s ) ) d s ) + 0 t ( t s ) α 1 Γ ( α ) f ( s , x ( s ) , y ( s ) ) d s ,
and
G 2 ( x , y ) ( t ) = t Δ ( η ζ 0 ρ ( ρ s ) β 2 Γ ( β 1 ) g ( s , x ( s ) , y ( s ) ) d s ζ 0 T ( T s ) α 1 Γ ( α ) f ( s , x ( s ) , y ( s ) ) d s + T ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) f ( s , x ( s ) , y ( s ) ) d s T 0 T ( T s ) β 1 Γ ( β ) g ( s , x ( s ) , y ( s ) ) d s ) + 0 t ( t s ) β 1 Γ ( β ) f ( s , x ( s ) , y ( s ) ) d s ,
Here we establish the existence of the solutions for the boundary value problem (1) and (2) by using Banach’s contraction mapping principle.
Theorem 1.
Assume f , g : C ( [ 0 , T ] × 2 are jointly continuous functions and there exist constants ϕ , ψ , such that   x 1 , x 2 , y 1 , y 2 , t [ 0 , T ] , we have
| f ( t , x 1 , x 2 ) f ( t , y 1 , y 2 ) | ϕ ( | x 2 x 1 | + | y 2 y 1 | ) ,
| g ( t , x 1 , x 2 ) f ( t , y 1 , y 2 ) | ψ ( | x 2 x 1 | + | y 2 y 1 | ) ,
where
ϕ ( Q 1 + Q 3 ) + ψ ( Q 2 + Q 4 ) < 1 ,
then the BVP (1) and (2) has a unique solution on [ 0 , T ] . Here
Q 1 = T | Δ | ( T α + 1 Γ ( α + 1 ) + | η ζ | μ α 1 Γ ( α ) ) + T α Γ ( α + 1 ) , Q 2 = T | Δ | ( | η | T ρ β 1 Γ ( β ) + | η | T β Γ ( β + 1 ) ) ,      Q 3 = T | Δ | ( | ζ | T α Γ ( α + 1 ) + T | ζ | μ α 1 Γ ( α ) ) ,      Q 4 = T | Δ | ( | η ζ | ρ β 1 Γ ( β ) + T β + 1 Γ ( β + 1 ) ) + T β Γ ( β + 1 ) .
Proof. 
Define s u p 0 t T | f ( t , 0 , 0 ) | = f 0 < ,   s u p 0 t T | g ( t , 0 , 0 ) | = g 0 < and Ω ε = { ( x , y ) C ( [ 0 , T ] , ) ×   C ( [ 0 , T ] , ) : ( x , y ) ε } , and ε > 0 , such that
ε ( Q 1 + Q 3 ) f 0 + ( Q 2 + Q 4 ) g 0 1 [ ϕ ( Q 1 + Q 3 ) + ψ ( Q 2 + Q 4 ) ] .
Firstly, we show that G Ω ε Ω ε .
By our assumption, for ( x , y ) Ω ε , t [ 0 , T ] , we have
| f ( t , x ( t ) , y ( t ) ) | | f ( t , x ( t ) , y ( t ) ) f ( t , 0 , 0 ) | + | f ( t , 0 , 0 ) | , ϕ ( | x ( t ) | + | y ( t ) | ) + f 0 ϕ ( x + y ) + f 0 , ϕ ε + f 0 ,
and
| g ( t , x ( t ) , y ( t ) ) | ψ ( | x ( t ) | + | y ( t ) | ) + g 0 ψ ( x + y ) + g 0 , ψ ε + g 0 ,
which lead to
| G 1 ( x , y ) ( t ) | T | Δ | ( | η | T 0 ρ ( ρ s ) β 2 Γ ( β 1 ) d s ( ψ ( x + y ) + g 0 ) + T 0 T ( T s ) α 1 Γ ( α ) d s ( ϕ ( x + y ) + f 0 ) + | η ζ | 0 μ ( μ s ) α 2 Γ ( α 1 ) d s ( ϕ ( x + y ) + f 0 ) + | η | 0 T ( T s ) β 1 Γ ( β ) d s ( ψ ( x + y ) + g 0 ) ) + s u p 0 t T 0 t ( t s ) α 1 Γ ( α ) d s ( ϕ ( x + y ) + f 0 ) ( ϕ ( x + y ) + f 0 ) [ T | Δ | ( T α + 1 Γ ( α + 1 ) + | η ζ | μ α 1 Γ ( α ) ) + T α Γ ( α + 1 ) ] + ( ψ ( x + y ) + g 0 ) [ T | Δ | ( | η | T ρ β 1 Γ ( β ) + | η | T β Γ ( β + 1 ) ) ] ( ϕ ( x + y ) + f 0 ) Q 1 + ( ψ ( x + y ) + g 0 ) Q 2 ( ϕ ε + f 0 ) Q 1 + ( ψ ε + g 0 ) Q 2 .   
In a similar manner:
| G 2 ( x , y ) ( t ) | ( ϕ ( x + y ) + f 0 ) Q 3 + ( ψ ( x + y ) + g 0 ) Q 4   ( ϕ ε + f 0 ) Q 3 + ( ψ ε + g 0 ) Q 4 .   
Hence,
G 1 ( x , y ) ( ϕ ε + f 0 ) Q 1 + ( ψ ε + g 0 ) Q 2 ,
and
G 2 ( x , y ) ( ϕ ε + f 0 ) Q 3 + ( ψ ε + g 0 ) Q 4 .
Consequently,
G ( x , y ) ( ϕ ε + f 0 ) ( Q 1 + Q 3 ) + ( ψ ε + g 0 ) ( Q 2 + Q 4 ) ε   .
and we get G ( x , y ) ε that is G Ω ε Ω ε .
Now let ( x 1 , y 1 ) , ( x 2 , y 2 ) C ( [ 0 , T ] , ) ×   C ( [ 0 , T ] , ) , t [ 0 , T ] .
Then we have
| G 1 ( x 1 , y 1 ) ( t ) G 1 ( x 2 , y 2 ) ( t ) | T | Δ | ( | η | T 0 ρ ( ρ s ) β 2 Γ ( β 1 ) d s ψ ( x 2 x 1 + y 2 y 1 ) + T 0 T ( T s ) α 1 Γ ( α ) d s ϕ ( x 2 x 1 + y 2 y 1 ) + | η ζ | 0 μ ( μ s ) α 2 Γ ( α 1 ) d s ϕ ( x 2 x 1 + y 2 y 1 ) + | η | 0 T ( T s ) β 1 Γ ( β ) d s ψ ( x 2 x 1 + y 2 y 1 ) ) + s u p 0 t T 0 t ( t s ) α 1 Γ ( α ) d s ϕ ( x 2 x 1 + y 2 y 1 ) ,
G 1 ( x 1 , y 1 ) G 1 ( x 2 , y 2 ) Q 1 ϕ ( x 2 x 1 + y 2 y 1 ) + Q 2 ψ ( x 2 x 1 + y 2 y 1 ) .
and likewise
G 2 ( x 1 , y 1 ) G 2 ( x 2 , y 2 ) Q 3 ϕ ( x 2 x 1 + y 2 y 1 ) + Q 4 ψ ( x 2 x 1 + y 2 y 1 ) .
From (11) and (12) we have
G ( x 1 , y 1 ) G ( x 2 , y 2 ) ( ϕ ( Q 1 + Q 3 ) + ψ ( Q 2 + Q 4 ) ) ( x 2 x 1 + y 2 y 1 ) .
Since ϕ ( Q 1 + Q 3 ) + ψ ( Q 2 + Q 4 ) < 1 , therefore, the operator G is a contraction operator. Hence, by Banach’s fixed-point theorem, the operator G   has a unique fixed point, which is the unique solution of the BVP (1) and (2). This completes the proof. □
Next we will prove the existence of solutions by applying the Leray-Schauder alternative.
Lemma 2.
“(Leray-Schauder alternative [7], p. 4) Let F : E E be a completely continuous operator (i.e., a map restricted to any bounded set in E is compact). Let E ( F ) = { x E : x = λ F ( x )   f o r   s o m e   0 < λ < 1 } . Then either the set E ( F ) is unbounded or F has at least one fixed point)”.
Theorem 2.
Assume f , g : C ( [ 0 , T ] × 2 are continuous functions and there exist θ 1 , θ 2 , λ 1 , λ 2 0 where θ 1 , θ 2 , λ 1 , λ 2 are real constants and θ 0 , λ 0 > 0 such that x i , y i ,   ( i = 1 , 2 ) , we have
| f ( t , x 1 , x 2 ) | θ 0 + θ 1 | x 1 | + θ 2 | x 2 | ,
| g ( t , x 1 , x 2 ) | λ 0 + λ 1 | x 1 | + λ 2 | x 2 | ,
If
( Q 1 + Q 3 ) θ 1 + ( Q 2 + Q 4 ) λ 1 < 1 ,
and
( Q 1 + Q 3 ) θ 2 + ( Q 2 + Q 4 ) λ 2 < 1 ,
where Q i , i = 1 , 2 , 3 , 4 are defined in (10), then the problem (1) and (2) has at least one solution.
Proof. 
This proof will be presented in two steps.
Step 1: We will show that G : C ( [ 0 , T ] , ) ×   C ( [ 0 , T ] , ) C ( [ 0 , T ] , ) ×   C ( [ 0 , T ] , ) is completely continuous. The continuity of the operator G holds by the continuity of the functions f , g .
Let B C ( [ 0 , T ] , ) ×   C ( [ 0 , T ] , ) be bounded. Then there exists positive constants k 1 , k 2   such that
| f ( t , x ( t ) , y ( t ) ) | k 1 ,      | g ( t , x ( t ) , y ( t ) ) | k 2 ,      t [ 0 , T ] .
Then ( x , y ) B , and we have
| G 1 ( x , y ) ( t ) | Q 1 k 1 + Q 2 k 2 ,
which implies
G 1 ( x , y ) Q 1 k 1 + Q 2 k 2 ,
and similarly
G 2 ( x , y ) Q 3 k 1 + Q 4 k 2 .
Thus, from the above inequalities, it follows that the operator G is uniformly bounded, since
G ( x , y ) ( Q 1 + Q 3 ) k 1 + ( Q 2 + Q 4 ) k 2 .
Next, we will show that operator G is equicontinuous. Let ω 1 , ω 2 [ 0 , T ] with ω 1 < ω 2 . This yields
| G 1 ( x , y ) ( ω 2 ) G 1 ( x , y ) ( ω 1 ) | ω 2 ω 1 | Δ | ( | η | T 0 ρ ( ρ s ) β 2 Γ ( β 1 ) | g ( s , x ( s ) , y ( s ) ) | d s + T 0 T ( T s ) α 1 Γ ( α ) | f ( s , x ( s ) , y ( s ) ) | d s + | η ζ | 0 μ ( μ s ) α 2 Γ ( α 1 ) | f ( s , x ( s ) , y ( s ) ) | d s + | η | 0 T ( T s ) β 1 Γ ( β ) | g ( s , x ( s ) , y ( s ) ) | d s ) + | 0 ω 2 ( ω 2 s ) α 1 Γ ( α ) f ( s , x ( s ) , y ( s ) ) d s 0 ω 1 ( ω 1 s ) α 1 Γ ( α ) f ( s , x ( s ) , y ( s ) ) d s | ω 2 ω 1 | Δ | ( | η | T k 2 0 ρ ( ρ s ) β 2 Γ ( β 1 ) d s + T k 1 0 T ( T s ) α 1 Γ ( α ) d s + | η ζ | k 1 0 μ ( μ s ) α 2 Γ ( α 1 ) d s + | η | k 2 0 T ( T s ) β 1 Γ ( β ) d s ) + | 0 ω 1 ( ( ω 2 s ) α 1 Γ ( α ) ( ω 1 s ) α 1 Γ ( α ) ) f ( s , x ( s ) , y ( s ) ) d s | + | ω 1 ω 2 ( ω 2 s ) α 1 Γ ( α ) f ( s , x ( s ) , y ( s ) ) d s | , ω 2 ω 1 | Δ | ( k 2 | η | T ρ β 1 Γ ( β ) + k 1 T α + 1 Γ ( α + 1 ) + k 1 | η ζ | μ α 1 Γ ( α ) + k 2 | η | T β Γ ( β + 1 ) ) + k 1 Γ ( α ) ( 0 ω 1 ( ( ω 2 s ) α 1 ( ω 1 s ) α 1 ) d s + ω 1 ω 2 ( ω 2 s ) α 1 d s ) .
And we obtain
| G 1 ( x , y ) ( ω 2 ) G 1 ( x , y ) ( ω 1 ) | ω 2 ω 1 | Δ | ( k 2 | η | T ρ β 1 Γ ( β ) + k 1 T α + 1 Γ ( α + 1 ) + k 1 | η ζ | μ α 1 Γ ( α ) + k 2 | η | T β Γ ( β + 1 ) ) + k 1 Γ ( α + 1 ) [ ω 2 α ω 1 α ] .
Hence, we have G 1 ( x , y ) ( ω 2 ) G 1 ( x , y ) ( ω 1 ) 0 independent of x and y as ω 2 ω 1 . Furthermore, we obtain
| G 2 ( x , y ) ( ω 2 ) G 2 ( x , y ) ( ω 1 ) | ω 2 ω 1 | Δ | ( k 2 | η ζ | ρ β 1 Γ ( β ) + k 1 | ζ | T α Γ ( α + 1 ) + k 1 T | ζ | μ α 1 Γ ( α ) + k 2 T β + 1 Γ ( β + 1 ) ) + k 2 Γ ( β + 1 ) [ ω 2 β ω 1 β ] ,
which implies that G 2 ( x , y ) ( ω 2 ) G 2 ( x , y ) ( ω 1 ) 0 independent of x and y as ω 2 ω 1 .
Therefore, operator G ( x , y ) is equicontinuous, and thus G ( x , y ) is completely continuous.
Step 2: (Boundedness of operator)
Finally, we will show that Z = { ( x , y ) C ( [ 0 , T ] , ) ×   C ( [ 0 , T ] , ) : ( x , y ) = h G ( x , y ) , h [ 0 , 1 ] } is bounded. Let ( x , y ) , with ( x , y ) = h G ( x , y ) for any t [ 0 , T ] , we have
x ( t ) = h G 1 ( x , y ) ( t ) ,     y ( t ) = h G 2 ( x , y ) ( t ) .
Then
| x ( t ) | Q 1 ( θ 0 + θ 1 | x ( t ) | + θ 2 | y ( t ) | ) + Q 2 ( λ 0 + λ 1 | x ( t ) | + λ 2 | y ( t ) | ) ,
and
| y ( t ) | Q 3 ( θ 0 + θ 1 | x ( t ) | + θ 2 | y ( t ) | ) + Q 4 ( λ 0 + λ 1 | x ( t ) | + λ 2 | y ( t ) | ) .
Hence,
x Q 1 ( θ 0 + θ 1 x + θ 2 y ) + Q 2 ( λ 0 + λ 1 x + λ 2 y ) ,
and
y Q 3 ( θ 0 + θ 1 x + θ 2 y ) + Q 4 ( λ 0 + λ 1 x + λ 2 y ) ,
which implies
x + y ( Q 1 + Q 3 ) θ 0 + ( Q 2 + Q 4 ) λ 0 + ( ( Q 1 + Q 3 ) θ 1 + ( Q 2 + Q 4 ) λ 1 ) x + ( ( Q 1 + Q 3 ) θ 2 + ( Q 2 + Q 4 ) λ 2 ) y .
Therefore,
( x , y ) ( Q 1 + Q 3 ) θ 0 + ( Q 2 + Q 4 ) λ 0 Q 0 ,
where Q 0 = m i n { 1 ( Q 1 + Q 3 ) θ 1 ( Q 2 + Q 4 ) λ 1 , 1 ( Q 1 + Q 3 ) θ 2 ( Q 2 + Q 4 ) λ 2 } . This proves that Z   is bounded and hence by Leray-Schauder alternative theorem, operator G   has at least one fixed point. Therefore, the BVP (1) and (2) has at least one solution on [ 0 , T ] . This completes the proof. □

4. Hyers-Ulam Stability

In this section, we will discuss the Hyers-Ulam stability of the solutions for the BVP (1) and (2) by means of integral representation of its solution given by
x ( t ) = G 1 ( x , y ) ( t ) , y ( t ) = G 2 ( x , y ) ( t ) ,
where G 1 and G 2 are defined by (8) and (9).
Define the following nonlinear operators N 1 , N 2 C ( [ 0 , T ] , ) × C ( [ 0 , T ] , ) C ( [ 0 , T ] , ) ;
D   c α x ( t ) f ( t , x ( t ) , y ( t ) ) = N 1 ( x , y ) ( t ) ,     t [ 0 , T ] ,     D   c β y ( t ) g ( t , x ( t ) , y ( t ) ) = N 2 ( x , y ) ( t ) ,     t [ 0 , T ] .    
For some ε 1 , ε 2 > 0 , we consider the following inequality:
N 1 ( x , y ) ε 1 ,      N 2 ( x , y ) ε 2 .
Definition 3.
([8,9]).The coupled system (1) and (2) is said to be Hyers-Ulam stable, if there exist M 1 , M 2 > 0 , such that for every solution ( x * , y * ) C ( [ 0 , T ] , ) × C ( [ 0 , T ] , ) of the inequality (13), there exists a unique solution ( x , y ) C ( [ 0 , T ] , ) × C ( [ 0 , T ] , ) of problems (1) and (2) with
( x , y ) ( x * , y * ) M 1 ε 1 + M 2 ε 2 .
Theorem 3.
Let the assumptions of Theorem 1 hold. Then the BVP (1) and (2) is Hyers-Ulam-stable.
Proof. 
Let ( x , y ) C ( [ 0 , T ] , ) × C ( [ 0 , T ] , ) be the solution of the problems (1) and (2) satisfying (8) and (9). Let ( x * , y * ) be any solution satisfying (13):
D   c α x * ( t ) = f ( t , x * ( t ) , y * ( t ) ) + N 1 ( x * , y * ) ( t ) ,     t [ 0 , T ] ,     D   c β y * ( t ) = g ( t , x * ( t ) , y * ( t ) ) + N 2 ( x * , y * ) ( t ) ,     t [ 0 , T ] .    
So
x * ( t ) = G 1 ( x * , y * ) ( t ) + t Δ ( η T 0 ρ ( ρ s ) β 2 Γ ( β 1 ) N 2 ( x * , y * ) ( s ) d s T 0 T ( T s ) α 1 Γ ( α ) N 1 ( x * , y * ) ( s ) d s + η ζ 0 μ ( μ s ) α 2 Γ ( α 1 ) N 1 ( x * , y * ) ( s ) d s η 0 T ( T s ) β 1 Γ ( β ) N 2 ( x * , y * ) ( s ) d s ) + 0 t ( t s ) α 1 Γ ( α ) N 1 ( x * , y * ) ( s ) d s ,
It follows that
| G 1 ( x * , y * ) ( t ) x * ( t ) | T | Δ | ( | η | T 0 ρ ( ρ s ) β 2 Γ ( β 1 ) d s ε 2 + T 0 T ( T s ) α 1 Γ ( α ) d s ε 1 + | η ζ | 0 μ ( μ s ) α 2 Γ ( α 1 ) d s ε 1 + | η | 0 T ( T s ) β 1 Γ ( β ) d s ε 2 ) + 0 T ( T s ) α 1 Γ ( α ) d s ε 1 ,   [ T | Δ | ( T α + 1 Γ ( α + 1 ) + | η ζ | μ α 1 Γ ( α ) ) + T α Γ ( α + 1 ) ]   ε 1 + T | Δ | ( | η | T ρ β 1 Γ ( β ) + | η | T β Γ ( β + 1 ) ) ε 2   , Q 1 ε 1 + Q 2 ε 2 .
Similarly,
| G 1 ( x * , y * ) ( t ) x * ( t ) | T | Δ | ( | ζ | T α Γ ( α + 1 ) + T | ζ | μ α 1 Γ ( α ) ) ε 1 + [ T | Δ | ( | η ζ | ρ β 1 Γ ( β ) + T β + 1 Γ ( β + 1 ) ) + T β Γ ( β + 1 ) ] ,   Q 3 ε 1 + Q 4 ε 2 ,
where Q i , i = 1 ,   2 ,   3 ,   4 are defined in (10).
Therefore, we deduce by the fixed-point property of operator G , that is given by (8) and (9), which
| x ( t ) x * ( t ) | = | x ( t ) G 1 ( x * , y * ) ( t ) + G 1 ( x * , y * ) ( t ) x * ( t ) | | G 1 ( x , y ) ( t ) G 1 ( x * , y * ) ( t ) | + | G 1 ( x * , y * ) ( t ) x * ( t ) | ( Q 1 ϕ + Q 2 ψ ) ( x , y ) ( x * , y * ) + Q 1 ε 1 + Q 2 ε 2 ,
and similarly
| y ( t ) y * ( t ) | = | y ( t ) G 2 ( x * , y * ) ( t ) + G 2 ( x * , y * ) ( t ) y * ( t ) | | G 2 ( x , y ) ( t ) G 2 ( x * , y * ) ( t ) | + | G 2 ( x * , y * ) ( t ) y * ( t ) | ( Q 3 ϕ + Q 4 ψ ) ( x , y ) ( x * , y * ) + Q 3 ε 1 + Q 4 ε 2 ,
From (14) and (15) it follows that
( x , y ) ( x * , y * ) ( Q 1 ϕ + Q 2 ψ + Q 3 ϕ + Q 4 ψ ) ( x , y ) ( x * , y * ) + ( Q 1 + Q 3 ) ε 1 + ( Q 2 + Q 4 ) ε 2 ,
( x , y ) ( x * , y * ) ( Q 1 + Q 3 ) ε 1 + ( Q 2 + Q 4 ) ε 2 1 ( ( Q 1 + Q 3 ) ϕ + ( Q 2 + Q 4 ) ψ ) , M 1 ε 1 + M 2 ε 2 .
with
M 1 = ( Q 1 + Q 3 ) 1 ( ( Q 1 + Q 3 ) ϕ + ( Q 2 + Q 4 ) ψ ) ,
M 2 = ( Q 2 + Q 4 ) 1 ( ( Q 1 + Q 3 ) ϕ + ( Q 2 + Q 4 ) ψ ) .
Thus, sufficient conditions for the Hyers-Ulam stability of the solutions are obtained. □

5. Examples

Example 1.
Consider the following coupled system of fractional differential equations
{   c D 3 2 x ( t ) = 1 6 π 81 + t 2 ( | x ( t ) | 3 + | x ( t ) | + | y ( t ) | 5 + | x ( t ) | ) ,   c D 7 4 y ( t ) = 1 12 π 64 + t 2 ( s i n ( x ( t ) ) + s i n ( y ( t ) ) ) , x ( 1 ) = 2 y ( 1 ) ,     y ( 1 ) = x ( 1 / 2 ) ,     x ( 0 ) = 0 ,     y ( 0 ) = 0 ,
α = 3 2 , β = 7 4 , T = 1 , η = 2 , ζ = 1 , μ = 1 2 , ρ = 1 .
Using the given data, we find that Δ = 3 , Q 1 = 1.269 , Q 2 = 1.1398 , Q 3 = 0.5167 , Q 4 = 1.554 , ϕ = 1 54 π , ψ = 1 48 π .
It is clear that
f ( t , x ( t ) , y ( t ) ) = 1 6 π 81 + t 2 ( | x ( t ) | 3 + | x ( t ) | + | y ( t ) | 5 + | x ( t ) | ) ,
and
g ( t , x ( t ) , y ( t ) ) = 1 12 π 64 + t 2 ( s i n ( x ( t ) ) + s i n ( y ( t ) ) ) ,
are jointly continuous functions and Lipschitz function with ϕ = 1 54 π , ψ = 1 48 π . Moreover,
1 54 π ( 1.269 + 0.5167 ) + 1 48 π ( 1.1398 + 1.554 ) = 0.0283 < 1 .
Thus, all the conditions of Theorem 1 are satisfied, then problem (16) has a unique solution on [ 0 , 1 ] , which is Hyers-Ulam-stable.
Example 2.
Consider the following system of fractional differential equation
{   c D 5 / 3 x ( t ) = 1 80 + t 4 + | x ( t ) | 120 ( 1 + y 2 ( t ) ) + 1 4 2500 + t 2 e 3 t cos ( y ( t ) ) ,   t [ 0 , 1 ]   c D 6 / 5 y ( t ) = 1 16 + t 2 cos t + 1 150 e 3 t s i n ( y ( t ) ) + 1 180 x ( t ) ,   t [ 0 , 1 ]    x ( 1 ) = 3 y ( 1 / 3 ) ,    y ( 1 ) = x ( 1 ) ,    x ( 0 ) = 0 ,    y ( 0 ) = 0 ,    
α = 5 3 , β = 6 5 , T = 1 , η = 3 , ζ = 1 , μ = 1 , ρ = 1 / 3 .
Using the given data, we find that Δ = 3 , Q 1 = 1.269 , Q 2 = 1.1398 , Q 3 = 0.5167 , Q 4 = 1.554 , ϕ = 1 54 π , ψ = 1 48 π .
It is clear that
| f ( t , x , y ) | 1 80 + 1 120 | x | + 1 200 | y | ,
| g ( t , x , y ) | 1 4 + 1 180 | x | + 1 150 | y | .
Thus, θ 0 = 1 80 , θ 1 = 1 120 , θ 2 = 1 200 , λ 0 = 1 4 , λ 1 = 1 180 , λ 2 = 1 150 .
Note that ( Q 1 + Q 3 ) θ 1 + ( Q 2 + Q 4 ) λ 1 = 0.0298 < 1   and ( Q 1 + Q 3 ) θ 2 + ( Q 2 + Q 4 ) λ 2 = 0.0269 < 1 , and hence by Theorem 2, problem (17) has at least one solution on [ 0 , 1 ] .

6. Conclusions

In this paper, the existence, uniqueness and the Hyers-Ulam stability of solutions for a coupled system of nonlinear fractional differential equations with boundary conditions were established and discussed.
Future studies may focus on different concepts of stability and existence results to a neutral time-delay system/inclusion, time-delay system/inclusion with finite delay.

Author Contributions

The authors have made the same contribution. All authors read and approved the final manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors wish to thank the anonymous reviewers for their valuable comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  2. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations. North-Holland Mathematics Studies; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  3. Chalishajar, D.; Raja, D.S.; Karthikeyan, K.; Sundararajan, P. Existence results for nonautonomous impulsive fractional evolution equations. Res. Nonlinear Anal. 2018, 1, 133–147. [Google Scholar]
  4. Chalishajar, D.; Kumar, A. Existence, uniqueness and Ulam’s stability of solutions for a coupled system of fractional differential equations with integral boundary conditions. Mathematics 2018, 6, 96. [Google Scholar] [CrossRef]
  5. Ntouyas, S.K.; Obaid, M. A coupled system of fractional differential equations with nonlocal integral boundary conditions. Adv. Differ. Equ. 2012, 2012, 130–139. [Google Scholar] [CrossRef]
  6. Ahmad, B.; Ntouyas, S.K. A Coupled system of nonlocal fractional differential equations with coupled and uncoupled slit-strips-type integral boundary conditions. J. Math. Sci. 2017, 226, 175–196. [Google Scholar] [CrossRef]
  7. Alsulami, H.H.; Ntouyas, S.K.; Agarwal, R.P.; Ahmad, B.; Alsaedi, A. A study of fractional-order coupled systems with a new concept of coupled non-separated boundary conditions. Bound. Value Probl. 2017, 2017, 68–74. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Bai, Z.; Feng, T. Existence results for a coupled system of nonlinear fractional three-point boundary value problems at resonance. Comput. Math. Appl. 2011, 61, 1032–1047. [Google Scholar] [CrossRef]
  9. Granas, A.; Dugundji, J. Fixed Point Theory; Springer: New York, NY, USA, 2005. [Google Scholar]
  10. Hyers, D.H. On the stability of the linear functional equation. Proc. Nat. Acad. Sci. USA 1941, 27, 222. [Google Scholar] [CrossRef] [PubMed]
  11. Rus, I.A. Ulam stabilities of ordinary differential equations in a Banach space. Carpathian J. Math. 2010, 103–107. [Google Scholar]
  12. Cabada, A.; Wang, G. Positive solutions of nonlinear fractional differential equations with integral boundary value conditions. J. Math. Anal. Appl. 2012, 389, 403–411. [Google Scholar] [CrossRef]
  13. Graef, J.R.; Kong, L.; Wang, M. Existence and uniqueness of solutions for a fractional boundary value problem on a graph. Fract. Calc. Appl. Anal. 2014, 17, 499–510. [Google Scholar] [CrossRef]
  14. Ahmad, B.; Nieto, J.J. Existence results for a coupled system of nonlinear fractional differential equations with three-point boundary conditions. Comput. Math. Appl. 2009, 58, 1838–1843. [Google Scholar] [CrossRef]
  15. Su, X. Boundary value problem for a coupled system of nonlinear fractional differential equations. Appl. Math. Lett. 2009, 22, 64–69. [Google Scholar] [CrossRef]
  16. Wang, J.; Xiang, H.; Liu, Z. Positive solution to nonzero boundary values problem for a coupled system of nonlinear fractional differential equations. Int. J. Differ. Equ. 2010, 10, 12. [Google Scholar] [CrossRef]
  17. Ahmad, B.; Ntouyas, S.K.; Alsaedi, A. On a coupled system of fractional differential equations with coupled nonlocal and integral boundary conditions. Chaos Solitons Fractals 2016, 83, 234–241. [Google Scholar] [CrossRef]
  18. Zhai, C.; Xu, L. Properties of positive solutions to a class of four-point boundary value problem of Caputo fractional differential equations with a parameter. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 2820–2827. [Google Scholar] [CrossRef]
  19. Ahmad, B.; Ntouyas, S.K. Existence results for a coupled system of Caputo type sequential fractional differential equations with nonlocal integral boundary conditions. Appl. Math. Comput. 2015, 266, 615–622. [Google Scholar] [CrossRef]
  20. Tariboon, J.; Ntouyas, S.K.; Sudsutad, W. Coupled systems of Riemann-Liouville fractional differential equations with Hadamard fractional integral boundary conditions. J. Nonlinear Sci. Appl. 2016, 9, 295–308. [Google Scholar] [CrossRef]
  21. Mahmudov, N.I.; Bawaneh, S.; Al-Khateeb, A. On a coupled system of fractional differential equations with four point integral boundary conditions. Mathematics 2019, 7, 279. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Mahmudov, N.I.; Al-Khateeb, A. Stability, Existence and Uniqueness of Boundary Value Problems for a Coupled System of Fractional Differential Equations. Mathematics 2019, 7, 354. https://doi.org/10.3390/math7040354

AMA Style

Mahmudov NI, Al-Khateeb A. Stability, Existence and Uniqueness of Boundary Value Problems for a Coupled System of Fractional Differential Equations. Mathematics. 2019; 7(4):354. https://doi.org/10.3390/math7040354

Chicago/Turabian Style

Mahmudov, Nazim I, and Areen Al-Khateeb. 2019. "Stability, Existence and Uniqueness of Boundary Value Problems for a Coupled System of Fractional Differential Equations" Mathematics 7, no. 4: 354. https://doi.org/10.3390/math7040354

APA Style

Mahmudov, N. I., & Al-Khateeb, A. (2019). Stability, Existence and Uniqueness of Boundary Value Problems for a Coupled System of Fractional Differential Equations. Mathematics, 7(4), 354. https://doi.org/10.3390/math7040354

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