Next Article in Journal
SRIFA: Stochastic Ranking with Improved-Firefly-Algorithm for Constrained Optimization Engineering Design Problems
Next Article in Special Issue
A Coupled System of Fractional Difference Equations with Nonlocal Fractional Sum Boundary Conditions on the Discrete Half-Line
Previous Article in Journal
Fractional Calculus of Extended Mittag-Leffler Function and Its Applications to Statistical Distribution
Previous Article in Special Issue
On Ulam Stability and Multiplicity Results to a Nonlinear Coupled System with Integral Boundary Conditions
 
 
Correction published on 10 December 2019, see Mathematics 2019, 7(12), 1220.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Existence and Stability Results for a Fractional Order Differential Equation with Non-Conjugate Riemann-Stieltjes Integro-Multipoint Boundary Conditions

1
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
2
Department of Mathematics, University of Ioannina, Ioannina 45110, Greece
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(3), 249; https://doi.org/10.3390/math7030249
Submission received: 20 February 2019 / Revised: 5 March 2019 / Accepted: 5 March 2019 / Published: 11 March 2019

Abstract

:
We discuss the existence and uniqueness of solutions for a Caputo-type fractional order boundary value problem equipped with non-conjugate Riemann-Stieltjes integro-multipoint boundary conditions on an arbitrary domain. Modern tools of functional analysis are applied to obtain the main results. Examples are constructed for the illustration of the derived results. We also investigate different kinds of Ulam stability, such as Ulam-Hyers stability, generalized Ulam-Hyers stability, and Ulam-Hyers-Rassias stability for the problem at hand.

1. Introduction

Fractional calculus played a pivotal role in improving the mathematical modeling of many real-world problems. The extensive application of fractional order (differential and integral) operators indeed reflects the popularity of this branch of mathematical analysis. In contrast to the integer order operators, such operators are nonlocal in nature and do have the capacity to trace the history of the phenomenon under investigation. A detailed account of the use of fractional calculus tools can be found in several scientific disciplines such as, chaos and fractional dynamics [1], evolution in honeycomb lattice via fractional Schrödinger equation [2], financial economics [3], ecology [4], bio-engineering [5], etc. For theoretical development and further application of the topic, see the texts [6,7,8,9].
During the past two decades, the study of fractional order boundary value problems has been one of the hot topics of scientific research. Several researchers contributed to the development of this class of problems by producing a huge number of articles, special issues, monographs, etc. Now the literature on the topic contains a variety of existence and uniqueness results, and analytic and numerical methods of solutions for these problems. In particular, there has been shown a great interest in the formulation and investigation of fractional order boundary value problems involving non-classical (nonlocal and integral) boundary conditions. The nonlocal boundary conditions are found to be of great utility in modeling the changes happening within the domain of the given scientific phenomena, while the concept of integral boundary conditions is applied to model the physical problems, such as blood flow problems on arbitrary structures and ill-posed backward problems. For some recent works on fractional order differential equations involving Riemann-Liouville, Caputo, and Hadamard type fractional derivatives, equipped with classical, nonlocal, and integral boundary conditions, we refer the reader to a series of papers [10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28] and the references cited therein.
In this paper, we study the existence of solutions for a nonlinear Liouville-Caputo-type fractional differential equation on an arbitrary domain:
c D q x ( t ) = f ( t , x ( t ) ) , 3 < q 4 , t [ a , b ] ,
supplemented with non-conjugate Riemann-Stieltjes integro-multipoint boundary conditions of the form:
x ( a ) = i = 1 n 2 α i x ( η i ) + a b x ( s ) d A ( s ) , x ( a ) = 0 , x ( b ) = 0 , x ( b ) = 0 ,
where c D q denotes the Caputo fractional derivative of order q , a < η 1 < η 2 < < η n 2 < b , f : [ a , b ] × R R is a given continuous function, A is a function of bounded variation, and α i R , i = 1 , 2 , , n 2 .
The main emphasis in the present work is to introduce non-conjugate Riemann-Stieltjes integro-multipoint boundary conditions and develop the existence theory for a Caputo-type fractional order boundary value problem equipped with these conditions on an arbitrary domain. Conjugate conditions on the body/fluid interface provide continuity of the thermal fields by specifying the equalities of temperatures and heat fluxes of a body and a flow at the vicinity of interface. The results obtained in this paper may have potential applications in diffraction-free and self-healing optoelectronic devices. Moreover, propagation properties for fractional Schrödinger equation similar to our results are well known theoretically [29].
The rest of the paper is organized as follows. An auxiliary result related to the linear variant of the problems (1) and (2), which plays a key role in the forthcoming analysis, is presented in Section 2. Some basic ideas of fractional calculus are also given in this section. In Section 3, we obtain some existence results for the given problem, while Section 4 contains a uniqueness result for the problem at hand. Ulam stability of different kinds for the problem (1) and (2) is studied in Section 5.

2. Preliminary Material

We begin this section with some basic definitions of fractional calculus [6]. Later we prove an auxiliary lemma, which plays a key role in defining a fixed-point problem associated with the given problem.
Definition 1.
Let g be a locally integrable real-valued function on a < t < b + . The Riemann-Liouville fractional integral I a p of order p R ( p > 0 ) for the function g is defined as
I a p g t = g K p ( t ) = 1 Γ p a t t s p 1 g s d s ,
where K p ( t ) = t p 1 Γ ( p ) , Γ denotes the Euler gamma function.
Definition 2.
The Caputo derivative of fractional order p for an ( m 1 ) -times absolutely continuous function g : [ a , ) R is defined as
c D p g ( t ) = 1 Γ ( m p ) a t ( t t ¯ ) m p 1 g ( m ) ( t ¯ ) d t ¯ , m 1 < p m , m = [ p ] + 1 ,
where [ p ] denotes the integer part of the real number p.
Lemma 1.
[6] The general solution of the fractional differential equation c D q x ( t ) = 0 , m 1 < q < m , t [ a , b ] is
x ( t ) = ω 0 + ω 1 ( t a ) + ω 2 ( t a ) 2 + + ω m 1 ( t a ) m 1 ,
where ω i R , i = 0 , 1 , , m 1 . Furthermore,
I q c D q x ( t ) = x ( t ) + i = 0 m 1 ω i ( t a ) i .
Lemma 2.
Let
γ 1 = A 1 ( b a ) 2 3 + 2 A 3 3 ( b a ) A 2 0 .
For f ^ C ( [ a , b ] , R ) , the unique solution of the linear equation
c D q x ( t ) = f ^ ( t ) , 3 < q 4 , t [ a , b ] ,
supplemented with the boundary conditions (2) is given by
x ( t ) = a t ( t s ) q 1 Γ ( q ) f ^ ( s ) d s g 1 ( t ) a b ( b s ) q 1 Γ ( q ) f ^ ( s ) d s g 2 ( t ) a b ( b s ) q 2 Γ ( q 1 ) f ^ ( s ) d s + g 3 ( t ) a b a s ( s u ) q 1 Γ ( q ) f ^ ( u ) d u d A ( s ) ,
where
g 1 ( t ) = λ 1 ( t a ) 2 A 1 γ 1 + ( t a ) 3 λ 4 , g 2 ( t ) = λ 2 + ( t a ) 2 γ 2 γ 1 + ( t a ) 3 λ 5 ,
g 3 ( t ) = λ 3 + ( t a ) 2 γ 1 + ( t a ) 3 λ 6 ,
γ 2 = A 3 + ( b a ) 3 A 1 3 ( b a ) 2 , λ 1 = 1 + ( b a ) 2 A 1 3 γ 1 , λ 2 = ( b a ) 3 ( b a ) 2 γ 2 3 γ 1 , λ 3 = ( b a ) 2 3 γ 1 ,
λ 4 = 2 A 1 3 ( b a ) γ 1 , λ 5 = γ 1 2 ( b a ) γ 2 3 ( b a ) 2 γ 1 , λ 6 = 2 3 ( b a ) γ 1 , A 1 = 1 i = 1 n 2 α i a b d A ( s ) , A 2 = i = 1 n 2 α i ( η i a ) 2 + a b ( s a ) 2 d A ( s ) ,
A 3 = i = 1 n 2 α i ( η i a ) 3 + a b ( s a ) 3 d A ( s ) .
Proof. 
Applying the integral operator I q to both sides of (4) and using Lemma 1, we get
x ( t ) = a t ( t s ) q 1 Γ ( q ) f ^ ( s ) d s + c 0 + c 1 ( t a ) + c 2 ( t a ) 2 + c 3 ( t a ) 3 ,
where c i R , i = 0 , 1 , 2 , 3 are unknown arbitrary constants. Differentiating (10) with respect to t, we have
x ( t ) = a t ( t s ) q 2 Γ ( q 1 ) f ^ ( s ) d s + c 1 + 2 c 2 ( t a ) + 3 c 3 ( t a ) 2 .
Using the boundary conditions (2) in (10) and (11), we obtain c 1 = 0 and
c 0 + ( b a ) 2 c 2 + ( b a ) 3 c 3 = I 1 ,
2 ( b a ) c 2 + 3 ( b a ) 2 c 3 = I 2 ,
A 1 c 0 A 2 c 2 A 3 c 3 = I 3 ,
where A i ( i = 1 , 2 , 3 ) are given by (9) and
I 1 = a b ( b s ) q 1 Γ ( q ) f ^ ( s ) d s , I 2 = a b ( b s ) q 2 Γ ( q 1 ) f ^ ( s ) d s , I 3 = a b a s ( s u ) q 1 Γ ( q ) f ^ ( u ) d u d A ( s ) .
Solving (12) and (13), for c 0 and c 3 in terms of c 2 , we get
c 0 = I 1 ( b a ) 3 I 2 ( b a ) 2 3 c 2 ,
c 3 = 1 3 ( b a ) 2 I 2 2 3 ( b a ) c 2 .
Substituting (16) and (17) in (14) yields
c 2 = γ 2 γ 1 I 2 A 1 γ 1 I 1 + 1 γ 1 I 3 , γ 1 0 ,
where γ 1 and γ 2 are defined by (3) and (7) respectively. Using (18) in (16) and (17), we find that
c 0 = λ 1 I 1 + λ 2 I 2 + λ 3 I 3 ,
c 3 = λ 4 I 1 + λ 5 I 2 + λ 6 I 3 .
Inserting the values of c 0 , c 1 , c 2 and c 3 in (10) together with notations (6), we obtain the solution (5). The converse of the lemma can be proved by direct computation. □

3. Existence Results

Let E = C ( [ a , b ] , R ) denote the Banach space of all continuous functions from [ a , b ] R equipped with the sup-norm x = sup { | x ( t ) | , t [ a , b ] } . For computational convenience, we introduce
Λ = ( b a ) q Γ ( q + 1 ) + g 1 ¯ ( b a ) q Γ ( q + 1 ) + g 2 ¯ ( b a ) q 1 Γ ( q ) + g 3 ¯ a b ( s a ) q Γ ( q + 1 ) d A ( s ) ,
where g 1 ¯ = sup t [ a , b ] | g 1 ( t ) | , g 2 ¯ = sup t [ a , b ] | g 2 ( t ) | , g 3 ¯ = sup t [ a , b ] | g 3 ( t ) | . By Lemma 2, we transform the problems (1) and (2) into an equivalent fixed-point problem as
x = J x ,
where J : E E is defined by
( J x ) ( t ) = a t ( t s ) q 1 Γ ( q ) f ( s , x ( s ) ) d s g 1 ( t ) a b ( b s ) q 1 Γ ( q ) f ( s , x ( s ) ) d s g 2 ( t ) a b ( b s ) q 2 Γ ( q 1 ) f ( s , x ( s ) ) d s + g 3 ( t ) a b a s ( s u ) q 1 Γ ( q ) f ( u , x ( u ) ) d u d A ( s ) ,
where g 1 ( t ) , g 2 ( t ) and g 3 ( t ) are given by (6).
Evidently, the existence of fixed points of the operator J will imply the existence of solutions for the problems (1) and (2).
Now, the platform is set to present our main results. The following known fixed-point theorem [30] will be used in the proof of our first result.
Theorem 1.
Let X be a Banach space. Assume that G : X X is a completely continuous operator and the set P = { x X | x = β G x , 0 < β < 1 } is bounded. Then G has a fixed point in X.
Theorem 2.
Suppose that there exists ϱ C ( [ a , b ] , R + ) such that | f ( t , x ( t ) ) | ϱ ( t ) , t [ a , b ] , x E , with sup t [ a , b ] | ϱ ( t ) | = ϱ . Then the problems (1) and (2) has at least one solution on [ a , b ] .
Proof. 
Observe that continuity of the operator J follows from that of f. Let Φ E be bounded. Then, x Φ together with the given assumption | f ( t , x ( t ) ) | ϱ ( t ) , we get
| ( J x ) | sup t [ a , b ] { a t ( t s ) q 1 Γ ( q ) | f ( s , x ( s ) ) | d s + | g 1 ( t ) | a b ( b s ) q 1 Γ ( q ) | f ( s , x ( s ) ) | d s + | g 2 ( t ) | a b ( b s ) q 2 Γ ( q 1 ) | f ( s , x ( s ) ) | d s + | g 3 ( t ) | a b a s ( s u ) q 1 Γ ( q ) | f ( u , x ( u ) ) | d u d A ( s ) } ϱ ( b a ) q Γ ( q + 1 ) + g 1 ¯ ( b a ) q Γ ( q + 1 ) + g 2 ¯ ( b a ) q 1 Γ ( q ) + g 3 ¯ a b ( s a ) q Γ ( q + 1 ) d A ( s ) = ϱ Λ = M 1 ,
which shows that J is bounded. Next, for a < t 1 < t 2 < b , we have
| ( J x ) ( t 2 ) ( J x ) ( t 1 ) | a t 1 | ( t 2 s ) q 1 ( t 1 s ) q 1 | Γ ( q ) | f ( s , x ( s ) ) | d s + t 1 t 2 | ( t 2 s ) q 1 | Γ ( q ) | f ( s , x ( s ) ) | d s + | g 1 ( t 2 ) g 1 ( t 1 ) | a b | ( b s ) q 1 | Γ ( q ) | f ( s , x ( s ) ) | d s + | g 2 ( t 2 ) g 2 ( t 1 ) | a b | ( b s ) q 2 | Γ ( q 1 ) | f ( s , x ( s ) ) | d s + | g 3 ( t 2 ) g 3 ( t 1 ) | a b a s | ( s u ) q 1 | Γ ( q ) | f ( u , x ( u ) ) | d u d A ( s ) ϱ [ | ( t 2 a ) q ( t 1 a ) q | + 2 ( t 2 t 1 ) q Γ ( q + 1 ) + | g 1 ( t 2 ) g 1 ( t 1 ) | | ( b a ) q | Γ ( q + 1 ) + | g 2 ( t 2 ) g 2 ( t 1 ) | | ( b a ) q 1 | Γ ( q ) + | g 3 ( t 2 ) g 3 ( t 1 ) | Γ ( q + 1 ) a b a s ( s a ) q ) d A ( s ) ] ,
which tends to zero as t 2 t 1 independent of x. Thus, J is equicontinuous on Φ . Hence, by Arzelá-Ascoli theorem, J is relatively compact on Φ . Therefore, J ( Φ ) is a relatively compact subset of E .
Now we consider a set P = { x E | x = β J x , 0 < β < 1 } , and show that the set P is bounded. Let x P , then x = β J x , 0 < β < 1 . For any t [ a , b ] , we have
| x ( t ) | = β | ( J x ) ( t ) | sup t [ a , b ] { a t ( t s ) q 1 Γ ( q ) | f ( s , x ( s ) ) | d s + | g 1 ( t ) | a b ( b s ) q 1 Γ ( q ) | f ( s , x ( s ) ) | d s + | g 2 ( t ) | a b ( b s ) q 2 Γ ( q 1 ) | f ( s , x ( s ) ) | d s + | g 3 ( t ) | a b a s ( s u ) q 1 Γ ( q ) | f ( u , x ( u ) ) | d u d A ( s ) } ϱ ( b a ) q Γ ( q + 1 ) + g 1 ¯ ( b a ) q Γ ( q + 1 ) + g 2 ¯ ( b a ) q 1 Γ ( q ) + g 3 ¯ a b ( s a ) q Γ ( q + 1 ) d A ( s ) = ϱ Λ ,
where Λ given by (19). Thus, x ϱ Λ for any t [ a , b ] . Therefore, the set P is bounded. In consequence, the conclusion of Theorem 1 applies and that the operator J has at least one fixed point. Thus, there exists at least one solution for the problems (1) and (2) on [ a , b ] . □
Example 1.
Consider the fractional boundary value problem
c D 11 3 x ( t ) = 6 e x 2 t 4 + 24 + cos t t 2 + 1 | x | 3 1 + | x | 3 + t 3 + 6 , t [ 1 , 2 ] , x ( 1 ) = i = 1 4 α i x ( η i ) + 1 2 x ( s ) d A ( s ) , x ( 1 ) = 0 , x ( 2 ) = 0 , x ( 2 ) = 0 ,
where q = 11 / 3 , a = 1 , b = 2 , α 1 = 1 / 2 , α 2 = 1 / 6 , α 3 = 1 / 6 , α 4 = 2 , η 1 = 6 / 5 , η 2 = 7 / 5 , η 3 = 8 / 5 , η 4 = 9 / 5 and f ( t , x ) = 6 e x 2 t 4 + 24 + cos t t 2 + 1 | x | 3 1 + | x | 3 + t 3 + 6 .
Clearly, | f ( t , x ) | 6 t 4 + 24 + cos t t 2 + 1 + t 3 + 6 = ϱ ( t ) > 0 . Therefore, there exists at least one solution for the problem (22) on [ 1 , 2 ] by the conclusion of Theorem 2.
Our next existence result is based on the following fixed-point theorem [30].
Theorem 3.
Let Ω be an open bounded subset of a Banach space X with 0 Ω and the operator F : Ω ¯ X is completely continuous satisfying F x x , x Ω . Then the operator F has a fixed point in Ω ¯ .
Theorem 4.
Let | f ( t , x ) | ξ | x | for 0 < | x | < τ , where τ and ξ are positive constants. Then the problems (1) and (2) has at least one solution for small values of ξ.
Proof. 
Let us choose ξ such that
Λ ξ < 1 ,
where Λ given by (19). Define B r 2 = { x E ; x r 2 } and take x E such that x = r 2 , that is, x B r 2 . As we argued in Theorem 2, it can be shown that J is completely continuous and
| ( J x ) ( t ) | sup t [ a , b ] { a t ( t s ) q 1 Γ ( q ) | f ( s , x ( s ) ) | d s + | g 1 ( t ) | a b ( b s ) q 1 Γ ( q ) | f ( s , x ( s ) ) | d s + | g 2 ( t ) | a b ( b s ) q 2 Γ ( q 1 ) | f ( s , x ( s ) ) | d s + | g 3 ( t ) | a b a s ( s u ) q 1 Γ ( q ) | f ( u , x ( u ) ) | d u d A ( s ) } ξ x ( b a ) q Γ ( q + 1 ) + g 1 ¯ ( b a ) q Γ ( q + 1 ) + g 2 ¯ ( b a ) q 1 Γ ( q ) + g 3 ¯ a b ( s a ) q Γ ( q + 1 ) d A ( s ) = ξ Λ x .
Using (23) and the norm x = sup t [ a , b ] | x ( t ) | , we get J x x , x B r 2 . Therefore, conclusion of Theorem 3 applies and hence the problems (1) and (2) has a solution on [ a , b ] . □
Example 2.
Consider the fractional boundary value problem.
c D 27 7 x ( t ) = 2 | x | 100 + t 2 1 + | x | 1 + | x | , t [ 0 , 1 ] , x ( 0 ) = i = 1 4 α i x ( η i ) + 0 1 x ( s ) d A ( s ) , x ( 0 ) = 0 , x ( 1 ) = 0 , x ( 1 ) = 0 ,
where q = 27 / 7 , a = 0 , b = 1 , α 1 = 2 , α 2 = 1 / 6 , α 3 = 1 / 6 , α 4 = 15 / 4 , η 1 = 1 / 7 , η 2 = 2 / 7 , η 3 = 3 / 7 , η 4 = 4 / 7 , and f ( t , x ( t ) ) = 2 | x | 100 + t 2 1 + | x | 1 + | x | . Let us take A ( s ) = s 2 2 . Using the given data, we have that A 1 1.25 , A 2 1.45068 , A 3 0.903110 , γ 1 0.43194 , γ 2 0.115630 , λ 1 1.96464 , λ 2 0.422566 , λ 3 0.771712 , λ 4 1.92928 , λ 5 0.154867 , λ 6 1.54342 , g 1 ¯ 1.96464 , g 2 ¯ 0.422566 , g 3 ¯ 0.771712 , Λ 0.243646 , where Λ is given by (19). Clearly the hypothesis of Theorem 4 is satisfied with ξ = 2 5 . Also ξ Λ 0.097458 < 1 . Therefore, the problem (24) has at least one solution on [ 0 , 1 ] .
In the next result, we apply a fixed-point theorem due to Krasnoselskii [31] to establish the existence of solutions for the problems (1) and (2).
Theorem 5.
(Krasnoselskii [31]) Let M be a closed, convex, bounded and nonempty subset of a Banach space X and let F 1 , F 2 be the operators defined from M to X such that: ( i ) F 1 x + F 2 y M wherever x , y M ; ( i i ) F 1 is compact and continuous; ( i i i ) F 2 is a contraction. Then there exists z M such that z = F 1 z + F 2 z .
Theorem 6.
Assume that f : [ a , b ] × R R is a continuous function such that the following conditions hold:
  • ( H 1 ) | f ( t , x ) f ( t , y ) | L | x y | , L > 0 , t [ a , b ] , x , y R ;
  • ( H 2 ) | f ( t , x ) | μ ( t ) , ( t , x ) [ a , b ] × R , μ C ( [ a , b ] , R + ) .
Then the problems (1) and (2) has at least one solution on [ a , b ] if
Λ ( b a ) q Γ ( q + 1 ) L < 1 ,
where Λ is defined by (19).
Proof. 
Consider a closed ball B r = { x E : x r } with r Λ μ , sup t [ a , b ] | μ ( t ) | = μ , and Λ is given by (19). Define operators J 1 and J 2 on B r as
( J 1 x ) ( t ) = a t ( t s ) q 1 Γ ( q ) f ( s , x ( s ) ) d s , ( J 2 x ) ( t ) = g 1 ( t ) a b ( t s ) q 1 Γ ( q ) f ( s , x ( s ) ) d s g 2 ( t ) a b ( t s ) q 2 Γ ( q 1 ) f ( s , x ( s ) ) d s + g 3 ( t ) a b a s ( s u ) q 1 Γ ( q ) f ( u , x ( u ) ) d u d A ( s ) .
Please note that J = J 1 + J 2 . For x , y B r , we have
J 1 x + J 2 y sup t [ a , b ] { a t ( t s ) q 1 Γ ( q ) | f ( s , x ( s ) ) | d s + | g 1 ( t ) | a b ( t s ) q 1 Γ ( q ) | f ( s , y ( s ) ) | d s + | ( g 2 ( t ) | a b ( t s ) q 2 Γ ( q 1 ) | f ( s , y ( s ) ) | d s + | g 3 ( t ) | a b a s ( s u ) q 1 Γ ( q ) | f ( u , y ( u ) | ) d u d A ( s ) } μ sup t [ a , b ] { a t ( t s ) q 1 Γ ( q ) d s + | g 1 ( t ) | a b ( t s ) q 1 Γ ( q ) d s + | g 2 ( t ) | a b ( t s ) q 2 Γ ( q 1 ) d s + | g 3 ( t ) | a b a s ( s u ) q 1 Γ ( q ) d u d A ( s ) } μ Λ r ,
where we have used (19). Thus, J 1 x + J 2 y B r . Next we show that J 2 is a contraction. For x , y B r , we have
J 2 x J 2 y = sup t [ a , b ] | ( J 2 x ) ( t ) ( J 2 y ) ( t ) | sup t [ a , b ] { | g 1 ( t ) | a b ( b s ) q 1 Γ ( q ) | f ( s , x ( s ) ) f ( s , y ( s ) ) | d s + | g 2 ( t ) | a b ( b s ) q 2 Γ ( q 1 ) | f ( s , x ( s ) ) f ( s , y ( s ) ) | d s + | g 3 ( t ) | a b a s ( s u ) q 1 Γ ( q ) | f ( u , x ( u ) ) f ( u , y ( u ) ) | d u d A ( s ) } L x y g ¯ 1 ( b a ) q Γ ( q + 1 ) + g ¯ 2 ( b a ) q 1 Γ ( q ) + g ¯ 3 a b ( s a ) q Γ ( q + 1 ) d A ( s ) = L Λ ( b a ) q Γ ( q + 1 ) x y ,
which shows that J 2 is a contraction by the condition (25). Continuity of f implies that the operator J 1 is continuous. Also, J 1 is uniformly bounded on B r as
J 1 x ( b a ) q Γ ( q + 1 ) μ .
Next, we establish that the operator J 1 is compact. Setting S = [ a , b ] × B r , we define sup ( t , x ) S | f ( t , x ) | = M r . For a < t 2 < t 1 < b , we get
| ( J 1 x ) ( t 1 ) ( J 1 x ) ( t 2 ) | = | a t 2 ( t 1 s ) q 1 ( t 2 s ) q 1 Γ ( q ) f ( s , x ( s ) ) d s + t 2 t 1 ( t 1 s ) q 1 Γ ( q ) f ( s , x ( s ) ) d s | M r Γ ( q + 1 ) | ( t 1 a ) q ( t 2 a ) q | + 2 ( t 1 t 2 ) q 0 when t 1 t 2 0 ,
independent of x. Thus, J 1 is equicontinuous on B r . Hence, by Arzelá-Ascoli theorem, J 1 is compact on B r . Therefore, the conclusion of Theorem 5 applies to the problems (1) and (2). □
Remark 1.
By interchanging the role of the operators J 1 and J 2 in Theorem 6, the condition (25) becomes:
L ( b a ) q Γ ( q + 1 ) < 1 .
Example 3.
Consider the fractional differential equation:
c D 11 3 x ( t ) = δ t 2 + 8 t a n 1 x + e 2 t , t [ 1 , 2 ] .
subject to the boundary conditions of Example 1.
Let us take A ( s ) = ( s 1 ) 2 2 + 1 2 . Using the given data (from Example 1), it is found that A 1 = 1 , A 2 1.543333 , A 3 1.245333 , γ 1 0.379778 , γ 2 0.081778 , λ 1 1.877706 , λ 2 0.261556 , λ 3 0.877706 , λ 4 1.755413 , λ 5 0.476887 , λ 6 1.755413 , g 1 ¯ 1.877706 , g 2 ¯ 0.264808 , g 3 ¯ 0.877706 , Λ 0.272140 (Λ is given by (19) and Λ ( b a ) 2 Γ ( q + 1 ) 0.204166 .
Clearly the hypotheses of Theorem 6 are satisfied with L = δ / 3 and μ ( t ) = δ π 2 t 2 + 8 + e 2 t . Also L Λ ( b a ) 2 Γ ( q + 1 ) < 1 for δ < 14.693960 . Therefore, there exists at least one solution for the problem (26) on [ 1 , 2 ] .

4. Uniqueness of Solution

Here, we prove the uniqueness of solutions for the problems (1) and (2).
Theorem 7.
Assume that f : [ a , b ] × R R is a continuous function satisfying the condition ( H 1 ) . Then the problems (1) and (2) has a unique solution on [ a , b ] if
L Λ < 1 ,
where Λ is given by (19).
Proof. 
Setting sup t [ a , b ] | f ( t , 0 ) | = N < , and selecting
r 1 N Λ ( 1 L Λ ) 1 ,
we define B r 1 = { x E : x r 1 } , and show that J B r 1 B r 1 , where the operator J is defined by (21). For x B r 1 ,
| f ( t , x ( t ) ) | = | f ( t , x ( t ) ) f ( t , 0 ) + f ( t , 0 ) | | f ( t , x ( t ) ) f ( t , 0 ) | + | f ( t , 0 ) | L | x ( t ) | + N L x + N L r 1 + N .
Then,
J x sup t [ a , b ] { a t ( t s ) q 1 Γ ( q ) | f ( s , x ( s ) ) | d s + | g 1 ( t ) | a b ( b s ) q 1 Γ ( q ) | f ( s , x ( s ) ) | d s + | g 2 ( t ) | a b ( b s ) q 2 Γ ( q 1 ) | f ( s , x ( s ) ) | d s + | g 3 ( t ) | a b a s ( s u ) q 1 Γ ( q ) | f ( u , x ( u ) ) | d u d A ( s ) } ( L x + N ) ( b a ) q Γ ( q + 1 ) + g 1 ¯ ( b a ) q Γ ( q + 1 ) + g 2 ¯ ( b a ) q 1 Γ ( q ) + g 3 ¯ a b ( s a ) q Γ ( q + 1 ) d A ( s ) ( L r 1 + N ) Λ r 1 .
This shows that J x B r 1 for any x B r 1 . Therefore, J B r 1 B r 1 . Now, we show that J is a contraction. For x , y E and t [ a , b ] , we obtain
( J x ) ( J y ) = sup t [ a , b ] | ( J x ) ( t ) ( J y ) ( t ) | sup t [ a , b ] { a t ( t s ) q 1 Γ ( q ) | f ( s , x ( s ) ) f ( s , y ( s ) ) | d s + | g 1 ( t ) | a b ( b s ) q 1 Γ ( q ) | f ( s , x ( s ) ) f ( s , y ( s ) ) | d s + | g 2 ( t ) | a b ( b s ) q 2 Γ ( q 1 ) | f ( s , x ( s ) ) f ( s , y ( s ) ) | d s + | g 3 ( t ) | a b a s ( s u ) q 1 Γ ( q ) | f ( u , x ( u ) ) f ( u , y ( u ) ) | d u d A ( s ) } L x y ( b a ) q Γ ( q + 1 ) + g 1 ¯ ( b a ) q Γ ( q + 1 ) + g 2 ¯ ( b a ) q 1 Γ ( q ) + g 3 ¯ a b ( s a ) q Γ ( q + 1 ) d A ( s ) = L Λ x y .
By the condition (27), we deduce from the above inequality that J is a contraction. Thus, by the conclusion of Banach fixed-point theorem, the problems (1) and (2) has a unique solution on [ a , b ] . □
Example 4.
Let us take the problem considered in Example 1, and note that L Λ < 1 for δ < 11.023738 . Clearly the hypothesis of Theorem 7 is satisfied. Hence it follows by the conclusion of Theorem 7 that the problem (22) has a unique solution on [ 1 , 2 ] .

5. Ulam Stability

In this section, we discuss the Ulam stability for the problems (1) and (2) by means of integral representation of its solution given by
y ( t ) = a t ( t s ) q 1 Γ ( q ) f ( s , y ( s ) ) d s g 1 ( t ) a b ( b s ) q 1 Γ ( q ) f ( s , y ( s ) ) d s g 2 ( t ) a b ( b s ) q 2 Γ ( q 1 ) f ( s , y ( s ) ) d s + g 3 ( t ) a b a s ( s u ) q 1 Γ ( q ) f ( u , y ( u ) ) d u d A ( s ) .
Here y C ( [ a , b ] , R ) possesses a fractional derivative of order 3 < q 4 and f : [ a , b ] × R R is a continuous function. Then the nonlinear operator Q : C ( [ a , b ] , R ) C ( [ a , b ] , R ) defined by
Q y ( t ) = c D q y ( t ) f ( t , y ( t ) )
is continuous.
Definition 3.
For each ϵ > 0 and for each solution y of (1) and (2) such that
Q y ϵ ,
the problems (1) and (2) is said to be Ulam-Hyers stable if we can find a positive real number ν and a solution x C ( [ a , b ] , R ) of (1) and (2) satisfying the inequality:
x y ν ϵ ,
where ϵ is a positive real number depending on ϵ .
Definition 4.
Let there exists κ C ( R + , R + ) such that for each solution y of (1) and (2), we can find a solution x C ( [ a , b ] , R ) of (1) and (2) such that
| x ( t ) y ( t ) | κ ( ϵ ) , t [ a , b ] .
Then the problems (1) and (2) is said to be generalized Ulam-Hyers stable
Definition 5.
For each ϵ > 0 and for each solution y of (1) and (2), the problems (1) and (2) is called Ulam-Hyers-Rassias stable with respect to σ C ( [ a , b ] , R + ) if
| Q y ( t ) | ϵ σ ( t ) , t [ a , b ] ,
and there exist a real number ν > 0 and a solution x C ( [ a , b ] , R ) of (1) and (2) such that
| x ( t ) y ( t ) | ν ϵ σ ( t ) , t [ a , b ] ,
where ϵ is a positive real number depending on ϵ .
Theorem 8.
Let the assumption ( H 1 ) hold with L Λ < 1 , where Λ is defined by (19). Then the problems (1) and (2) is both Ulam-Hyers and generalized Ulam-Hyers stable.
Proof. 
Let x C ( [ a , b ] , R ) be a solution of (1) and (2) satisfying (21) by Theorem 7. Let y be any solution satisfying (29). Then by Lemma 2, y satisfies the integral equation (28). Furthermore, the equivalence in Lemma 2 implies the equivalence between the operators Q and J I (where I is identity operator) for every solution y C ( [ a , b ] , R ) of (1) and (2) satisfying (27). Therefore, we deduce by the fixed-point property of the operator J (given by (21)) and (29) that
| y ( t ) x ( t ) | = | y ( t ) J y ( t ) + J y ( t ) J x ( t ) | | J x ( t ) J y ( t ) | + | J y ( t ) y ( t ) | L Λ x y + ϵ ,
where ϵ > 0 and L Λ < 1 . In consequence, we get
x y ϵ 1 L Λ .
Fixing ϵ = ϵ 1 L Λ and ν = 1 , we obtain the Ulam-Hyers stability condition. In addition, the generalized Ulam-Hyers stability follows by taking κ ( ϵ ) = ϵ 1 L Λ .  □
Theorem 9.
Assume that ( H 1 ) holds with L < Λ 1 (where Λ is defined by (19)), and there exists a function σ C ( [ a , b ] , R + ) satisfying the condition (30). Then the problems (1) and (2) is Ulam-Hyers-Rassias stable with respect to σ.
Proof. 
Following the arguments employed in the proof of Theorem 8, we have
x y ϵ σ ( t ) ,
where ϵ = ϵ 1 L Λ . This completes the proof. □
Example 5.
Consider the following fractional differential equation
c D 27 7 x ( t ) = 3 t 2 + 16 cos x + x 1 + x , t [ 0 , 1 ] ,
subject to the same data and the boundary conditions given in Example 1 with f ( t , x ( t ) ) = 3 t 2 + 16 cos x + x 1 + x . Obviously | f ( t , x ) f ( t , y ) | 3 8 x y , so, L = 3 / 8 and L Λ 0.102053 < 1 . Then the problem (31) is Ulam-Hyers stable, and generalized Ulam-Hyers stable. In addition, if there exists a continuous and positive function σ = e 1 + t 2 + 5 satisfying the condition (30), then the problem (31) is Ulam-Hyers-Rassias stable with the given value of f ( t , x ) .

6. Conclusions

We have obtained several existence results for a new class of Caputo-type fractional differential equations of order q ( 3 , 4 ] , supplemented with non-conjugate Riemann-Stieltjes integro-multipoint boundary conditions on an arbitrary domain by imposing different kinds of conditions on the nonlinear function involved in the problem. The existence results, relying on different fixed-point theorems, are presented in Section 3. The uniqueness of solution for the given problem is studied in Section 4 with the aid of Banach fixed-point theorem. Section 5 is concerned with different kinds of Ulam stability for the problem at hand. Some new results follow as special cases of the ones presented in this paper. For example, taking A ( s ) = s , our results correspond to the ones for non-conjugate integro-multipoint boundary conditions of the form: x ( a ) = i = 1 n 2 α i x ( η i ) + a b x ( s ) d s , x ( a ) = 0 , x ( b ) = 0 , x ( b ) = 0 and the value of Λ given by (19) takes the following form in this situation:
Λ = ( b a ) q Γ ( q + 1 ) + g 1 ¯ ( b a ) q Γ ( q + 1 ) + g 2 ¯ ( b a ) q 1 Γ ( q ) + g 3 ¯ ( b a ) q + 1 Γ ( q + 2 ) .
Letting A ( s ) = 0 in our results, we get the ones for non-conjugate multipoint boundary conditions of the form: x ( a ) = i = 1 n 2 α i x ( η i ) , x ( a ) = 0 , x ( b ) = 0 , x ( b ) = 0 . In case we take α i = 0 for all i = 1 , , n 2 , our results reduce to the ones with non-conjugate Riemann-Stieltjes boundary conditions. By fixing A ( s ) = 0 and α i = 0 for all i = 1 , , n 2 , our results correspond to a boundary value problem of fractional order q ( 3 , 4 ] with conjugate boundary conditions.

Author Contributions

Formal Analysis, B.A., Y.A., S.K.N. and A.A.

Funding

This project was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah, Saudi Arabia, under grant no. (RG-25-130-38).

Acknowledgments

The authors acknowledge with thanks DSR technical and financial support. The authors also thank the reviewers for their constructive remarks on our work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zaslavsky, G.M. Hamiltonian Chaos and Fractional Dynamics; Oxford University Press: Oxford, UK, 2005. [Google Scholar]
  2. Zhang, D.; Zhang, Y.; Zhang, Z.; Ahmed, N.; Zhang, Y.; Li, F.; Belic, M.R.; Xiao, M. Unveiling the link between fractional Schrödinger equation and light propagation in honeycomb lattice. Ann. Phys. 2017, 529, 1700149. [Google Scholar] [CrossRef]
  3. Fallahgoul, H.A.; Focardi, S.M.; Fabozzi, F.J. Fractional Calculus and Fractional Processes with Applications to Financial Economics. Theory and Application; Elsevier/Academic Press: London, UK, 2017. [Google Scholar]
  4. Javidi, M.; Ahmad, B. Dynamic analysis of time fractional order phytoplankton-toxic phytoplankton-zooplankton system. Ecol. Model. 2015, 318, 8–18. [Google Scholar] [CrossRef]
  5. Magin, R.L. Fractional Calculus in Bioengineering; Begell House Publishers: Danbury, CT, USA, 2006. [Google Scholar]
  6. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematics Studies, 204; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  7. Lakshimikantham, V.; Leela, S.; Devi, J.V. Theory of Fractional Dynamic Systems; Cambridge Academic Publishers: Cambridge, UK, 2009. [Google Scholar]
  8. Diethelm, K. The Analysis of Fractional Differential Equations. An Application—Oriented Exposition Using Differential Operators of Caputo Type; Lecture Notes in Mathematics 2004; Springer: Berlin, Germany, 2010. [Google Scholar]
  9. Ahmad, B.; Alsaedi, A.; Ntouyas, S.K.; Tariboon, J. Hadamard-Type Fractional Differential Equations, Inclusions and Inequalities; Springer: Cham, Switzerland, 2017. [Google Scholar]
  10. Ahmad, B.; Nieto, J.J. Riemann-Liouville fractional integro-differential equations with fractional nonlocal integral boundary conditions. Bound. Value Probl. 2011, 2011, 36. [Google Scholar] [CrossRef] [Green Version]
  11. Liang, S.; Zhang, J. Existence of multiple positive solutions for m-point fractional boundary value problems on an infinite interval. Math. Comput. Model. 2011, 54, 1334–1346. [Google Scholar] [CrossRef]
  12. Bai, Z.B.; Sun, W. Existence and multiplicity of positive solutions for singular fractional boundary value problems. Comput. Math. Appl. 2012, 63, 1369–1381. [Google Scholar] [CrossRef] [Green Version]
  13. Agarwal, R.P.; O’Regan, D.; Stanek, S. Positive solutions for mixed problems of singular fractional differential equations. Math. Nachr. 2012, 285, 27–41. [Google Scholar] [CrossRef]
  14. Ahmad, B.; Ntouyas, S.K.; Alsaedi, A. A study of nonlinear fractional differential equations of arbitrary order with Riemann-Liouville type multistrip boundary conditions. Math. Probl. Eng. 2013, 2013, 320415. [Google Scholar] [CrossRef]
  15. Ahmad, B.; Ntouyas, S.K. Existence results for higher order fractional differential inclusions with multi-strip fractional integral boundary conditions. Electron. J. Qual. Theory Differ. Equ. 2013, 2013, 20. [Google Scholar] [CrossRef]
  16. O’Regan, D.; Stanek, S. Fractional boundary value problems with singularities in space variables. Nonlinear Dyn. 2013, 71, 641–652. [Google Scholar] [CrossRef]
  17. 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]
  18. 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] [Green Version]
  19. Wang, G.; Liu, S.; Zhang, L. Eigenvalue problem for nonlinear fractional differential equations with integral boundary conditions. Abstr. Appl. Anal. 2014, 2014, 916260. [Google Scholar] [CrossRef]
  20. Henderson, J.; Kosmatov, N. Eigenvalue comparison for fractional boundary value problems with the Caputo derivative. Fract. Calc. Appl. Anal. 2014, 17, 872–880. [Google Scholar] [CrossRef] [Green Version]
  21. Zhang, L.; Ahmad, B.; Wang, G. Successive iterations for positive extremal solutions of nonlinear fractional differential equations on a half line. Bull. Aust. Math. Soc. 2015, 91, 116–128. [Google Scholar] [CrossRef]
  22. Henderson, J.; Luca, R. Nonexistence of positive solutions for a system of coupled fractional boundary value problems. Bound. Value Probl. 2015, 2015, 138. [Google Scholar] [CrossRef] [Green Version]
  23. Ntouyas, S.K.; Etemad, S. On the existence of solutions for fractional differential inclusions with sum and integral boundary conditions. Appl. Math. Comput. 2015, 266, 235–243. [Google Scholar] [CrossRef]
  24. Mei, Z.D.; Peng, J.G.; Gao, J.H. Existence and uniqueness of solutions for nonlinear general fractional differential equations in Banach spaces. Indag. Math. 2015, 26, 669–678. [Google Scholar] [CrossRef]
  25. Ahmad, B.; Ntouyas, S.K. Existence results for fractional differential inclusions with Erdelyi-Kober fractional integral conditions. Analele Universitatii Ovidius Constanta-Seria Matematica 2017, 25, 5–24. [Google Scholar] [CrossRef]
  26. Srivastava, H.M. Remarks on some families of fractional-order differential equations. Integr. Transf. Spec. Funct. 2017, 28, 560–564. [Google Scholar] [CrossRef]
  27. Wang, G.; Pei, K.; Agarwal, R.P.; Zhang, L.; Ahmad, B. Nonlocal Hadamard fractional boundary value problem with Hadamard integral and discrete boundary conditions on a half-line. J. Comput. Appl. Math. 2018, 343, 230–239. [Google Scholar] [CrossRef]
  28. Ahmad, B.; Luca, R. Existence of solutions for sequential fractional integro-differential equations and inclusions with nonlocal boundary conditions. Appl. Math. Comput. 2018, 339, 516–534. [Google Scholar] [CrossRef]
  29. Zhang, Y.; Wang, R.; Zhong, H.; Zhang, J.; Belic, M.R.; Zhang, Y. Resonant mode conversions and Rabi oscillations in a fractional Schrödinger equation. Opt. Express 2017, 25, 32401. [Google Scholar] [CrossRef]
  30. Sun, J.X. Nonlinear Functional Analysis and Its Application; Science Press: Bejing, China, 2018. [Google Scholar]
  31. Krasnoselskii, M.A. Two remarks on the method of successive approximations. Uspekhi Matematicheskikh Nauk 1955, 10, 123–127. [Google Scholar]

Share and Cite

MDPI and ACS Style

Ahmad, B.; Alruwaily, Y.; Alsaedi, A.; K. Ntouyas, S. Existence and Stability Results for a Fractional Order Differential Equation with Non-Conjugate Riemann-Stieltjes Integro-Multipoint Boundary Conditions. Mathematics 2019, 7, 249. https://doi.org/10.3390/math7030249

AMA Style

Ahmad B, Alruwaily Y, Alsaedi A, K. Ntouyas S. Existence and Stability Results for a Fractional Order Differential Equation with Non-Conjugate Riemann-Stieltjes Integro-Multipoint Boundary Conditions. Mathematics. 2019; 7(3):249. https://doi.org/10.3390/math7030249

Chicago/Turabian Style

Ahmad, Bashir, Ymnah Alruwaily, Ahmed Alsaedi, and Sotiris K. Ntouyas. 2019. "Existence and Stability Results for a Fractional Order Differential Equation with Non-Conjugate Riemann-Stieltjes Integro-Multipoint Boundary Conditions" Mathematics 7, no. 3: 249. https://doi.org/10.3390/math7030249

APA Style

Ahmad, B., Alruwaily, Y., Alsaedi, A., & K. Ntouyas, S. (2019). Existence and Stability Results for a Fractional Order Differential Equation with Non-Conjugate Riemann-Stieltjes Integro-Multipoint Boundary Conditions. Mathematics, 7(3), 249. https://doi.org/10.3390/math7030249

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