Next Article in Journal
Yule–Walker Equations Using a Gini Covariance Matrix for the High-Dimensional Heavy-Tailed PVAR Model
Next Article in Special Issue
Analysis and Computation of Solutions for a Class of Nonlinear SBVPs Arising in Epitaxial Growth
Previous Article in Journal
Fractal Stochastic Processes on Thin Cantor-Like Sets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlocal Sequential Boundary Value Problems for Hilfer Type Fractional Integro-Differential Equations and Inclusions

by
Nawapol Phuangthong
1,
Sotiris K. Ntouyas
2,3,
Jessada Tariboon
1,* and
Kamsing Nonlaopon
4
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
Department of Mathematics, University of Ioannina, 451 10 Ioannina, Greece
3
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
4
Department of Mathematics, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(6), 615; https://doi.org/10.3390/math9060615
Submission received: 16 February 2021 / Revised: 7 March 2021 / Accepted: 10 March 2021 / Published: 15 March 2021
(This article belongs to the Special Issue New Trends on Boundary Value Problems)

Abstract

:
In the present research, we study boundary value problems for fractional integro-differential equations and inclusions involving the Hilfer fractional derivative. Existence and uniqueness results are obtained by using the classical fixed point theorems of Banach, Krasnosel’skiĭ, and Leray–Schauder in the single-valued case, while Martelli’s fixed point theorem, a nonlinear alternative for multivalued maps, and the Covitz–Nadler fixed point theorem are used in the inclusion case. Examples are presented to illustrate our results.

1. Introduction

In last few decades, fractional differential equations with initial/boundary conditions have been studied by many researchers. This is because fractional differential equations describe many real world processes more accurately compared to classical order differential equations. Therefore, the fractional-order models become more practical and realistic compared to the integer-order models. Fractional differential equations arise in lots of engineering and clinical disciplines, including biology, physics, chemistry, economics, signal and image processing, control theory, and so on; see the monographs [1,2,3,4,5,6,7,8].
There exist several different definitions of fractional integrals and derivatives in the literature, such as the Riemann–Liouville and Caputo fractional derivatives, the Hadamard fractional derivative, the Erdeyl–Kober fractional derivative, and so on. Both Riemann–Liouville and Caputo fractional derivatives were generalized by Hilfer in [9]. This generalization, known as the Hilfer fractional derivative of order α and a type β [ 0 , 1 ] , interpolates between the Riemann–Liouville and Caputo derivative, since it can be reduced to the Riemann–Liouville and Caputo fractional derivatives when β = 0 and β = 1 , respectively. See [9,10,11] and the references cited therein for some properties and applications of the Hilfer derivative.
Several authors have studied initial value problems involving Hilfer fractional derivatives, see, for example, [12,13,14] and the references included therein. Boundary value problems for the Hilfer fractional derivative and nonlocal boundary conditions were initiated in [15].
Motivated by the research going in this direction, in the present paper, we study existence and uniqueness of solutions for the following new class of boundary value problems consisting of fractional-order sequential Hilfer-type differential equations supplemented with nonlocal integro-multipoint boundary conditions of the form
H D α , β + k H D α 1 , β x ( t ) = f ( t , x ( t ) , I δ x ( t ) ) , t [ a , b ] ,
x ( a ) = 0 , x ( b ) = i = 1 m 2 ζ i I ϕ i x ( θ i ) ,
where H D α , β denote the Hilfer fractional derivative operator of order α , 1 < α < 2 and parameter β , 0 β 1 , f : [ a , b ] × R × R R is a continuous function, given constant k R , I ψ as the Riemann–Liouville fractional integral of order ψ > 0 , ψ { ϕ i , δ } , a 0 , a < θ 1 < θ 2 < < θ m 2 < b and ζ i R , i = 1 , 2 , , m 2 .
By using standard tools from functional analysis, we study existence and uniqueness results for the sequential boundary value in Problems (1) and (2). We establish the existence of a unique solution via Banach’s fixed point theorem [16], while nonlinear alternatives of Leray–Schauder-type [17] and Krasnosel’skiĭ’s fixed point theorem [18] are applied to obtain the existence results.
After that, we look at the corresponding multivalued problem by studying existence of solutions for a new class of sequential boundary value problems of Hilfer-type fractional differential inclusions with nonlocal integro-multipoint boundary conditions of the form
H D α , β + k H D α 1 , β x ( t ) F ( t , x ( t ) , I δ x ( t ) ) , t [ a , b ] ,
x ( a ) = 0 , x ( b ) = i = 1 m 2 ζ i I ϕ i x ( θ i ) ,
where F : [ a , b ] × R × R P ( R ) is a multivalued map ( P ( R ) is the family of all nonempty subjects of R ).
Existence results for the sequential boundary value Problems (3) and (4) with convex-valued maps are derived by applying a fixed point theorem according to Martelli’s [19] nonlinear alternative for multivalued maps. For possible nonconvex-valued maps, we obtain an existence result by using Covitz–Nadler fixed point theorem [20] for contractive maps.
The organization of this paper is as follows: Section 2 contains some preliminary concepts related to our problem. We present our main work for Problems (1) and (2) in Section 3, while the main results for the multivalued Problems (3) and (4) are presented in Section 4. Examples are constructed to illustrate the main results.

2. Preliminaries

We present the preliminary results from fractional calculus needed in our proofs as follows [2,5].
Definition 1.
The fractional integral of Riemann–Liouville, of a continuous function of order α > 0 is defined by
I α u ( t ) = 1 Γ ( α ) a t ( t s ) α 1 u ( s ) d s , n 1 < α < n ,
where n = [ α ] + 1 ( [ α ] is the integer part of α R ), provided the right-hand side is pointwise defined on ( a , ) .
Definition 2.
For a continuous function, the Riemann–Liouville fractional derivative of order α > 0 is defined by
R L D α u ( t ) : = D n I n α u ( t ) = 1 Γ ( n α ) d d t n a t ( t s ) n α 1 u ( s ) d s , n 1 < α < n ,
where D n = d n / d t n , provided the right-hand side is pointwise defined on ( a , ) .
Definition 3.
For a continuous function, the Caputo fractional derivative of order α > 0 is defined by
C D α u ( t ) : = I n α D n u ( t ) = 1 Γ ( n α ) a t ( t s ) n α 1 d d s n u ( s ) d s , n 1 < α < n ,
provided the right-hand side is pointwise defined on ( a , ) .
Another new definition of fractional derivative was introduced in [9] (see also [11]).
Definition 4.
The Hilfer fractional derivative of order α and parameter β is defined by
H D α , β u ( t ) = I β ( n α ) D n I ( 1 β ) ( n α ) u ( t ) ,
where n 1 < α < n , 0 β 1 , t > a , D = d d t .
Remark 1.
When β = 0 , the Hilfer fractional derivative gives the Riemann–Liouville fractional derivative
H D α , 0 u ( t ) = D n I n α u ( t ) ,
while when β = 1 the Hilfer fractional derivative gives the Caputo fractional derivative
H D α , 1 u ( t ) = I n α D n u ( t ) .
The following basic lemma is used in the sequel.
Lemma 1
([11]). Let f L ( a , b ) , n 1 < α n , n N , 0 β 1 , I ( n α ) ( 1 β ) f A C k [ a , b ] . Then,
I α H D α , β f ( t ) = f ( t ) k = 0 n 1 ( t a ) k ( n α ) ( 1 β ) Γ ( k ( n α ) ( 1 β ) + 1 ) lim t a + d k d t k I ( 1 β ) ( n α ) f ( t ) .

3. Main Results

The following lemma concerns a linear variant of the sequential boundary value Problems (1) and (2).
Lemma 2.
Let h C ( [ a , b ] , R ) and a 0 , 1 < α < 2 , γ = α + 2 β α β . Assume that
Λ : = ( b a ) γ 1 Γ ( γ ) i = 1 m 2 ζ i ( θ i a ) γ + ϕ i 1 Γ ( γ + ϕ i ) 0 .
Then, the function x is a solution of the sequential boundary value problem
H D α , β + k H D α 1 , β x ( t ) = h ( t ) , t [ a , b ] ,
x ( a ) = 0 , x ( b ) = i = 1 m 2 ζ i I ϕ i x ( θ i ) ,
if and only if
x ( t ) = I α h ( t ) k a t x ( s ) d s + ( t a ) γ 1 Λ i = 1 m 2 ζ i I α + ϕ i h ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k a b x ( s ) d s I α h ( b ) .
Proof. 
Assume that x is a solution of the sequential nonlocal boundary value Problems (6) and (7). Applying the operator I α to both sides of Equation (6) and using Lemma 1, we obtain
x ( t ) = c 0 ( t a ) ( 2 α ) ( 1 β ) Γ ( 1 ( 2 α ) ( 1 β ) ) + c 1 ( t a ) 1 ( 2 α ) ( 1 β ) Γ ( 2 ( 2 α ) ( 1 β ) ) k I 1 x ( t ) + I α h ( t ) = c 0 ( t a ) γ 2 Γ ( γ 1 ) + c 1 ( t a ) γ 1 Γ ( γ ) k I 1 x ( t ) + I α h ( t ) ,
where c 0 , c 1 R and ( 1 β ) ( 2 α ) = 2 γ .
From the boundary condition x ( a ) = 0 , we obtain c 0 = 0 . Then, we get
x ( t ) = c 1 ( t a ) γ 1 Γ ( γ ) k I 1 x ( t ) + I α h ( t ) .
From the boundary condition x ( b ) = i = 1 m 2 ζ i I ϕ i x ( θ i ) , we find
c 1 = Γ ( γ ) Λ i = 1 m 2 ζ i I α + ϕ i h ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k a b x ( s ) d s I α h ( b ) .
Substituting the value of c 1 into (9), we obtain the solution (8). Conversely, it is easily shown that the solution x given by (8) satisfies the sequential nonlocal boundary value Problems (6) and (7). This completes the proof. □
Let C ( [ a , b ] , R ) denote the Banach space endowed with the sup-norm x = sup { | x ( t ) | : t [ a , b ] } . In view of Lemma 2, we define an operator A : C ( [ a , b ] , R ) C ( [ a , b ] , R ) by
( A x ) ( t ) = I α f x δ ( t ) k a t x ( s ) d s + ( t a ) γ 1 Λ i = 1 m 2 ζ i I α + ϕ i f x δ ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k a b x ( s ) d s I α f x δ ( b ) ,
where
f x δ ( t ) = f ( t , x ( t ) , I δ x ( t ) ) .
It is obvious that the sequential nonlocal boundary value problem has a solution if and only if the operator A has fixed points.
To simplify the computations, we use the following notations:
Ω = ( b a ) γ 1 | Λ | i = 1 m 2 | ζ i | ( θ i a ) α + ϕ i Γ ( α + ϕ i + 1 ) + ( b a ) α Γ ( α + 1 ) + ( b a ) α Γ ( α + 1 )
and
Ω 1 = ( b a ) γ 1 | Λ | | k | i = 1 m 2 | ζ i | ( θ i a ) ϕ i + 1 Γ ( ϕ i + 2 ) + | k | ( b a ) + | k | ( b a ) .
By using classical fixed point theorems, we prove in the next subsections, for the sequential boundary value Problems (1) and (2), our main existence and uniqueness results.

3.1. Existence and Uniqueness Result for Problems (1) and (2)

Our first result is an existence and uniqueness result, based on Banach’s fixed point theorem [16].
Theorem 1.
Assume that
( H 1 )  
there exists a constant L > 0 such that
| f ( t , x 1 , x 2 ) f ( t , y 1 , y 2 ) | L ( | x 1 y 1 | + | x 2 y 2 | )
for each t [ a , b ] and x i , y i R , i = 1 , 2 .
If
L L 1 Ω + Ω 1 < 1 ,
where Ω and Ω 2 are defined by (11) and (12), respectively, and L 1 = 1 + ( b a ) δ Γ ( δ + 1 ) , then the sequential boundary value Problems (1) and (2) have a unique solution on [ a , b ] .
Proof. 
Consider the operator A defined in (10). The sequential boundary value Problems (1) and (2) are then transformed into a fixed point problem x = A x . We shall show that A has a unique fixed point by applying the Banach contraction mapping principle.
We set sup t [ a , b ] | f ( t , 0 , 0 ) | = M < , and choose r > 0 such that
r M Ω 1 L L 1 Ω Ω 1 .
Now, we show that A B r B r , where B r = { x C ( [ a , b ] , R ) : x r } . By using the assumption ( H 1 ) , we have
f x δ ( t ) : = | f ( t , x ( t ) , I δ x ( t ) ) | | f ( t , x ( t ) , I δ x ( t ) ) f ( t , 0 , 0 , 0 ) | + | f ( t , 0 , 0 , 0 ) | L ( | x ( t ) | + | I δ x ( t ) | ) + M L x + ( b a ) δ Γ ( δ + 1 ) x + M = L x 1 + ( b a ) δ Γ ( δ + 1 ) + M = L L 1 x + M .
For any x B r , we have
| ( A x ) ( t ) | sup t [ a , b ] { I α | f x δ ( t ) | + | k | I 1 | x ( t ) | + ( t a ) γ 1 | Λ | ( i = 1 m 2 | ζ i | I α + ϕ i | f x δ ( θ i ) | + | k | i = 1 m 2 | ζ i | I ϕ i + 1 | x ( θ i ) | + | k | I 1 | x ( b ) | + I α | f x δ ( b ) | ) } ( L L 1 x + M ) I α ( 1 ) ( t ) + | k | I 1 | x ( t ) | + ( t a ) γ 1 | Λ | { ( L L 1 x + M ) i = 1 m 2 | ζ i | I α + ϕ i ( 1 ) ( θ i ) + | k | i = 1 m 2 | ζ i | I ϕ i + 1 | x ( θ i ) | + | k | I 1 | x ( b ) | + ( L L 1 x + M ) I α ( 1 ) ( b ) } ( L L 1 x + M ) ( b a ) α Γ ( α + 1 ) + | k | ( b a ) x + ( b a ) γ 1 | Λ | { ( L L 1 x + M ) i = 1 m 2 | ζ i | ( θ i a ) α + ϕ i Γ ( α + ϕ i + 1 ) + | k | x i = 1 m 2 | ζ i | ( θ i a ) ϕ i + 1 Γ ( ϕ i + 2 ) + | k | ( b a ) x + ( L L 1 x + M ) ( b a ) α Γ ( α + 1 ) } = ( b a ) γ 1 | Λ | i = 1 m 2 | ζ i | ( θ i a ) α + ϕ i Γ ( α + ϕ i + 1 ) + ( b a ) α Γ ( α + 1 ) + ( b a ) α Γ ( α + 1 ) ( L L 1 x + M ) + ( b a ) γ 1 | Λ | | k | i = 1 m 2 | ζ i | ( θ i a ) ϕ i + 1 Γ ( ϕ i + 2 ) + | k | ( b a ) + | k | ( b a ) x ( L L 1 r + M ) Ω + Ω 1 r r ,
which implies that A B r B r .
Next, we let x , y C ( [ a , b ] , R ) . Then, for t [ a , b ] , we have
| ( A x ) ( t ) ( A y ) ( t ) | I α | f x δ ( t ) f y δ ( t ) | + | k | I 1 | x ( t ) y ( t ) | + ( t a ) γ 1 | Λ | { i = 1 m 2 | ζ i | I α + ϕ i | f x δ ( θ i ) f y δ ( θ ) | + | k | i = 1 m 2 | ζ i | I ϕ i + 1 | x ( θ i ) y ( θ ) | + | k | I 1 | x ( b ) y ( b ) | + I α | f x δ ( b ) f y δ ( b ) | } ( b a ) γ 1 | Λ | i = 1 m 2 | ζ i | ( θ i a ) α + ϕ i Γ ( α + ϕ i + 1 ) + ( b a ) α Γ ( α + 1 ) + ( b a ) α Γ ( α + 1 ) L L 1 x y + ( b a ) γ 1 | Λ | | k | i = 1 m 2 | ζ i | ( θ i a ) ϕ i + 1 Γ ( ϕ i + 2 ) + | k | ( b a ) + | k | ( b a ) x y = ( L L 1 Ω + Ω 1 ) x y ,
which implies that A x A y ( L L 1 Ω + Ω 1 ) x y . As L L 1 Ω + Ω 1 < 1 , A is a contraction. Thus, by the Banach’s contraction mapping principle, we deduce that A has a unique fixed point, which is the unique solution of the sequential boundary value Problems (1) and (2). The proof is finished. □
Example 1.
Consider the sequential nonlocal boundary value problem with Hilfer fractional differential equation of the form:
H D 5 3 , 3 4 + 1 10 H D 2 3 , 3 4 x ( t ) = f ( t , x ( t ) , I 2 3 x ( t ) ) , t 1 4 , 7 4 , x 1 4 = 0 , x 7 4 = 1 5 I 1 4 x 1 2 + 2 7 I 1 3 x 3 4 + 3 11 I 1 2 x 1 + 4 13 I 3 2 x 5 4 + 5 16 I 5 4 x 3 2 ,
where
f ( t , x ( t ) , I 2 3 x ( t ) ) = 4 4 t + 159 x 2 ( t ) + 2 | x ( t ) | 1 + | x ( t ) | + 4 4 t + 83 sin | I 2 3 x ( t ) | + 1 3 e t .
Here, α = 5 / 3 , β = 3 / 4 , k = 1 / 10 , δ = 2 / 3 , a = 1 / 4 , b = 7 / 4 , m = 7 , ζ 1 = 1 / 5 , ζ 2 = 2 / 7 , ζ 3 = 3 / 11 , ζ 4 = 4 / 13 , ζ 5 = 5 / 16 , ϕ 1 = 1 / 4 , ϕ 2 = 1 / 3 , ϕ 3 = 1 / 2 , ϕ 4 = 3 / 2 , ϕ 5 = 5 / 4 , θ 1 = 1 / 2 , θ 2 = 3 / 4 , θ 3 = 1 , θ 4 = 5 / 4 , θ 5 = 3 / 2 , and γ = ( 5 / 3 ) + ( 2 ( 5 / 3 ) ) ( 3 / 4 ) = 23 / 12 . From these constants, we can verify that Λ 0.8651080717 , Ω 3.928448219 , Ω 1 0.4944426735 , and L 1 2.451539773 .
Since
| f ( t , w 1 , w 2 ) f ( t , z 1 , z 2 ) | 1 20 ( | w 1 z 1 | + | w 2 z 2 | ) , t 1 4 , 7 4 , w i , z i R , i = 1 , 2 ,
we set L = 1 / 20 , which leads to L L 1 Ω + Ω 1 0.9759800261 < 1 . Then, inequality (13) holds. Therefore, by applying Theorem 1, the Hilfer fractional differential equation with nonlocal Conditions (15) and (16) has a unique solution on [ 1 / 4 , 7 / 4 ] .

3.2. Existence Results

Two existence results are presented in this subsection. The first is based on the well-known Krasnosel’skiĭ’s fixed point theorem [18].
Theorem 2.
Let f : [ a , b ] × R × R R be a continuous function satisfying the following assumption:
( H 2 )  
| f ( t , x , y ) | φ ( t ) , ( t , x , y ) [ a , b ] × R × R , and φ C ( [ a , b ] , R + ) .
Then, if
Ω 1 < 1 ,
the sequential boundary value Problems (1) and (2) have at least one solution on [ a , b ] .
Proof. 
We set sup t [ a , b ] φ ( t ) = φ and choose ρ > 0 such that
ρ φ Ω 1 Ω 1 ,
(where Ω , Ω 1 are defined by (11) and (12), respectively), and consider B ρ = { x C ( [ a , b ] , R ) : x ρ } . On B ρ , we define the operators A 1 , A 2 by
( A 1 x ) ( t ) = I α f x δ ( t ) + ( t a ) γ 1 Λ i = 1 m 2 ζ i I α + ϕ i f x δ ( θ i ) I α f x δ ( b ) , t [ a , b ]
and
( A 2 x ) ( t ) = k I 1 x ( t ) + ( t a ) γ 1 Λ k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) , t [ a , b ] .
For any x , y B ρ , we have
| ( A 1 x ) ( t ) + ( A 2 y ) ( t ) | { I α | f x δ ( t ) | + | k | I 1 | x ( t ) | + ( t a ) γ 1 | Λ | ( i = 1 m 2 | ζ i | I α + ϕ i | f x δ ( θ i ) | + | k | i = 1 m 2 | ζ i | I ϕ i + 1 | x ( θ i ) | + | k | I 1 | x ( b ) | + I α | f x δ ( b ) | } ( b a ) γ 1 | Λ | i = 1 m 2 | ζ i | ( θ i a ) α + ϕ i Γ ( α + ϕ i + 1 ) + ( b a ) α Γ ( α + 1 ) + ( b a ) α Γ ( α + 1 ) φ + ( b a ) γ 1 | Λ | | k | i = 1 m 2 | ζ i | ( θ i a ) ϕ i + 2 Γ ( ϕ i + 1 ) + | k | ( b a ) + | k | ( b a ) x φ Ω + Ω 1 ρ ρ .
This shows that A 1 x + A 2 y B ρ . As we have proved with little difficulty, using (17), operator A 2 is a contraction mapping.
A 1 is continuous, since f is continuous. In addition, A 1 is uniformly bounded on B ρ as
A 1 x ( b a ) γ 1 | Λ | i = 1 m 2 | ζ i | ( θ i a ) α + ϕ i Γ ( α + ϕ i + 1 ) + ( b a ) α Γ ( α + 1 ) + ( b a ) α Γ ( α + 1 ) φ .
Now, we prove the compactness of the operator A 1 . Let t 1 , t 2 [ a , b ] with t 1 < t 2 . Then, we have
| ( A 1 x ) ( t 2 ) ( A 1 x ) ( t 1 ) | = 1 Γ ( α ) a t 1 [ ( t 2 s ) α 1 ( t 1 s ) α 1 ] f ( s , x ( s ) , I δ x ( s ) ) d s + t 1 t 2 ( t 2 s ) α 1 f ( s , x ( s ) , I δ x ( s ) ) d s + ( t 2 a ) γ 1 ( t 1 a ) γ 1 | Λ | i = 1 m 2 ζ i I α + ϕ i | f x δ ( θ i ) | + I α | f x δ ( b ) | ϕ Γ ( α + 1 ) [ 2 ( t 2 t 1 ) α + | ( t 2 a ) α ( t 1 a ) α | ] + ( t 2 a ) γ 1 ( t 1 a ) γ 1 | Λ | i = 1 m 2 ζ i I α + ϕ i | f x δ ( θ i ) | + I α | f x δ ( b ) | ,
which is independent of x and tends to zero as t 2 t 1 0 . Hence, A 1 is equicontinuous. The operator A 1 is compacted on B ρ by the Arzel a ´ –Ascoli theorem. Thus, all the assumptions of Krasnosel’skiĭ’s fixed point theorem are satisfied. So, by its conclusion, the sequential boundary value Problems (1) and (2) have at least one solution on [ a , b ] . The proof is finished. □
Our second existence result is proved by using Leray–Schauder’s Nonlinear Alternative [17].
Theorem 3.
Assume that (17) holds. Moreover, we suppose that
( H 3
there exists a continuous, nondecreasing, subhomogeneous (that is, ψ ( μ x ) μ ψ ( x ) for all μ 1 and x R + ) function ψ : [ 0 , ) ( 0 , ) and a function p C ( [ a , b ] , R + ) such that
| f ( t , u , v ) | p ( t ) ψ ( | u | + | v | ) for each ( t , u , v ) [ a , b ] × R × R ;
( H 4
there exists a constant K > 0 such that
( 1 Ω 1 ) K L 1 ψ ( K ) p Ω > 1 ,
where Ω and Ω 1 are defined by (11) and (12), respectively, and L 1 is defined in Theorem 1.
Then, the sequential boundary value Problems (1) and (2) have at least one solution on [ a , b ] .
Proof. 
Consider the operator A given by (10). In the first step, we shall show that bounded sets (balls) are mapped by A into bounded sets in C ( [ a , b ] , R ) . For r > 0 , let a bounded ball in C ( [ a , b ] , R ) given by B r = { x C ( [ a , b ] , R ) : x r } . Then, for t [ a , b ] , we have
| ( A x ) ( t ) | sup t [ a , b ] I α | f x δ ( t ) | + | k | I 1 | x ( t ) | + ( t a ) γ 1 | Λ | ( i = 1 m 2 | ζ i | I α + ϕ i | f x δ ( θ i ) | + | k | i = 1 m 2 | ζ i | I ϕ i + 1 | x ( θ i ) | + | k | I 1 | x ( b ) | + I α | f x δ ( b ) | ) } ( b a ) γ 1 | Λ | i = 1 m 2 | ζ i | ( θ i a ) α + ϕ i Γ ( α + ϕ i + 1 ) + ( b a ) α Γ ( α + 1 ) + ( b a ) α Γ ( α + 1 ) ψ L 1 x p + ( b a ) γ 1 | Λ | | k | i = 1 m 2 | ζ i | ( θ i a ) ϕ i + 1 Γ ( ϕ i + 2 ) + | k | ( b a ) + | k | ( b a ) x L 1 ψ ( x ) p ( b a ) γ 1 | Λ | i = 1 m 2 | ζ i | ( θ i a ) α + ϕ i Γ ( α + ϕ i + 1 ) + ( b a ) α Γ ( α + 1 ) + ( b a ) α Γ ( α + 1 ) + ( b a ) γ 1 | Λ | | k | i = 1 m 2 | ζ i | ( θ i a ) ϕ i + 1 Γ ( ϕ i + 2 ) + | k | ( b a ) + | k | ( b a ) x ,
and consequently,
A x L 1 ψ ( r ) p Ω + Ω 1 r .
Next, we will show that bounded sets are mapped by A into equicontinuous sets of C ( [ a , b ] , R ) . Let τ 1 , τ 2 [ a , b ] with τ 1 < τ 2 and x B r . Then, we have
| ( A x ) ( τ 2 ) ( A x ) ( τ 1 ) | 1 Γ ( α ) a τ 1 [ ( τ 2 s ) α 1 ( τ 1 s ) α 1 ] f x δ ( s ) d s + τ 1 τ 2 ( τ 2 s ) α 1 f x δ ( s ) d s + r ( t 2 t 1 ) + ( τ 2 a ) γ 1 ( τ 1 a ) γ 1 | Λ | i = 1 m 2 ζ i I α + ϕ i | f x δ ( θ i ) | + I α | f x δ ( b ) | + | k | i = 1 m 2 | ζ i | I ϕ i + 1 | x ( θ i ) | + | k | I 1 | x ( b ) | p L 1 ψ ( r ) Γ ( α + 1 ) [ 2 ( τ 2 τ 1 ) α + | ( τ 2 a ) α ( τ 1 a ) α | ] + r ( τ 2 τ 1 ) + ( τ 2 a ) γ 1 ( τ 1 a ) γ 1 | Λ | i = 1 m 2 ζ i I α + ϕ i | f x δ ( θ i ) | + I α | f x δ ( b ) | + | k | i = 1 m 2 | ζ i | I ϕ i + 1 | x ( θ i ) | + | k | I 1 | x ( b ) | .
In the above inequality, the right-hand side is independent of x B r tends to zero as τ 2 τ 1 0 . By the Arzel a ´ –Ascoli theorem, the operator A : C ( [ a , b ] , R ) C ( [ a , b ] , R ) is completely continuous.
Now, to finish the proof, we must prove that the set of all solutions to equations x = λ A x for λ ( 0 , 1 ) is bounded.
Let x be a solution. By using the computations in the first step, we have for t [ a , b ] ,
| x ( t ) | L 1 ψ ( x ) p Ω + Ω 1 x ,
and consequently,
( 1 Ω 1 ) x L 1 ψ ( x ) p Ω 1 .
By ( H 4 ) , there exists K such that x K . Consider the set
U = { x C ( [ a , b ] , R ) : x < K } .
We have proved that the operator A : U ¯ C ( [ a , b ] , R ) is completely continuous. There is no x U such that x = λ A x for some λ ( 0 , 1 ) , by the choice of U . By the nonlinear alternative of Leray–Schauder type ([17]), we have as result that A has a fixed point x U ¯ , which is a solution of the sequential boundary value Problems (1) and (2). The proof is completed. □
Example 2.
Consider the sequential nonlocal Hilfer fractional boundary value problem of the form:
H D 7 5 , 3 8 + 1 9 H D 2 5 , 3 8 x ( t ) = f ( t , x ( t ) , I 3 2 x ( t ) ) , t 1 5 , 8 5 , x 1 5 = 0 , x 8 5 = 2 11 I 1 3 x 2 5 + 4 13 I 2 3 x 3 5 + 6 17 I 4 3 x 4 5 + 8 19 I 5 3 x 1 + 10 23 I 7 3 x 6 5 + 12 29 I 8 3 x 7 5 .
Here, α = 7 / 5 ,   β = 3 / 8 ,   k = 1 / 9 ,   δ = 3 / 2 ,   a = 1 / 5 ,   b = 8 / 5 ,   m = 8 ,   ζ 1 = 2 / 11 ,   ζ 2 = 4 / 13 ,   ζ 3 = 6 / 17 ,   ζ 4 = 8 / 19 ,   ζ 5 = 10 / 23 ,   z e t a 6 = 12 / 29 ,   ϕ 1 = 1 / 3 ,   ϕ 2 = 2 / 3 ,   ϕ 3 = 4 / 3 ,   ϕ 4 = 5 / 3 ,   ϕ 5 = 7 / 3 ,   ϕ 6 = 8 / 3 ,   θ 1 = 2 / 5 ,   θ 2 = 3 / 5 ,   θ 3 = 4 / 5 ,   θ 4 = 1 ,   θ 5 = 6 / 5 ,   θ 6 = 7 / 5 and γ = ( 7 / 5 ) + ( 2 ( 7 / 5 ) ) ( 3 / 8 ) = 65 / 40 . From the given data, we find that Λ 0.8346353019 , Ω 3.404450305 , Ω 1 0.4283560839 , and L 1 2.246108486 .
(i) Let the nonlinear function f ( t , x , I 3 2 x ) in Problem (19) be defined by
f ( t , x , I 3 2 x ) = t 3 e x 4 + 1 2 t 2 x 12 1 + x 12 sin 6 I 3 2 x + t tan 1 I 3 2 x + 1 4 .
Then, we see that the function f satisfies condition ( H 2 ) in Theorem 2 as
| f ( t , w , z ) | t 3 + 1 2 t 2 + π 2 t + 1 4 : = φ ( t ) , ( t , w , z ) 1 5 , 8 5 × R 2 .
By the benefit of the conclusion in Theorem 2, it implies that the nonlocal Hilfer fractional boundary value Problems (19) and (20), with f given by (20), have at least one solution on [ 1 / 5 , 8 / 5 ] .
(ii) Let the nonlinear function f ( t , x , I 3 2 x ) in Problem (19) be given by
f ( t , x , I 3 2 x ) = 5 5 t + 74 | x | 23 1 + x 22 + 1 + 5 5 t + 79 | I 3 2 x | cos 8 | I 3 2 x | + 5 5 t + 84 1 | x | + I 3 2 x + 1 .
Thus, we can compute that the above function f satisfies inequality
| f ( t , w , z ) | 5 5 t + 74 | w | + | z | + 1 + 1 | w | + | z | + 1 .
Setting p ( t ) = 5 / ( 5 t + 74 ) and ψ ( u ) = u + 1 + ( 1 / ( u + 1 ) ) , u = | w | + | z | , we obtain p = 1 / 15 and ψ ( μ u ) μ ψ ( u ) , μ 1 , which is a subhomogeneous function. Then, there exists a constant K > 9.060160134 satisfying condition ( H 4 ) in Theorem 3. Therefore, all assumptions of Theorem 3 are fulfilled. Then, the boundary value Problems (19)–(21), with f given by (21), have at least one solution on [ 1 / 5 , 8 / 5 ] .

4. Existence Results for Problems (3) and (4)

In this section, we will establish existence results for a new class of sequential boundary value problems of Hilfer-type fractional differential inclusions with nonlocal integro-multipoint boundary Conditions (3) and (4).
In the following, by P p , we denote the set of all nonempty subsets of X that have the property “p”, where “p” will be bounded (b), closed (cl), convex (c), compact (cp), etc. Thus, P c l ( X ) = { Y P ( X ) : Y   is   closed } , P c p ( X ) = { Y P ( X ) : Y   is   compact } , P c p , c ( X ) = { Y P ( X ) : Y   is   compact   and   convex } , and P b , c l , c ( X ) = { Y P ( X ) : Y   is   bounded ,   closed ,   and convex } .
For details on multivalued maps, we refer the interested reader to the books by Castaing and Valadier [21], Deimling [16], Gorniewicz [22], and Hu and Papageorgiou [23]. For a brief summary of the needed results of multivalued analysis for fractional differential equations, we refer to [7,24,25].
For the sequential boundary value Problems (3) and (4), we will give the definition of its solution.
Definition 5.
A function x C 2 ( [ a , b ] , R ) is a solution of the sequential boundary value Problems (3) and (4) if there exists a function v L 1 ( J , R ) with v F ( t , x , y ) a.e. on [ a , b ] such that x satisfies the sequential fractional differential equation D α x ( t ) = v ( t ) on [ a , b ] and the nonlocal integro-multipoint boundary Condition (4).

4.1. The Upper Semicontinuous Case

Let us discuss first the case when the multivalued F has convex values and we give an existence result based on Martelli’s fixed point theorem [19]. For the reader’s convenience, we state the following form of Martelli’s fixed point, which is the multivalued version of Schaefer’s fixed point theorem.
Lemma 3.
(Martelli’s fixed point theorem) [19] Let X be a Banach space and T : X P b , c l , c ( X ) be a completely continuous multivalued map. If the set E = { x X : λ x T ( x ) , λ > 1 } is bounded, then T has a fixed point.
Theorem 4.
Assume that (17) holds. In addition, we assume that
( A 1 )  
F : [ a , b ] × R × R P c p , c ( R ) is L 1 -Carathéodory i.e.,
(i) 
t F ( t , x , y ) is measurable for each y R ;
(ii) 
( x , y ) F ( t , x , y ) is upper semicontinuous for almost all t [ a , b ] ;
(iii) 
for each r > 0 , there exists ϕ r L 1 ( J , R + ) such that
F ( t , x , y ) = sup { | v | : v F ( t , x , y } ϕ r ( t )
for all x , y R with x , y r and for a.e. t [ a , b ] ;
( A 2 )  
there exists a function q C ( [ a , b ] , R ) such that
F ( t , x , y ) q ( t ) , for a . e . t [ a , b ] and each x , y R .
Then, the sequential boundary value Problems (3) and (4) have at least one solution on [ a , b ] .
Proof. 
We consider the multivalued map N : C ( [ a , b ] , R ) P ( C ( [ a , b ] , R ) ) in order to transform Problems (3) and (4) into a fixed point problem:
N ( x ) = h C ( [ a , b ] , R ) : h ( t ) = I α v ( t ) k I 1 x ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α v ( b ) ] , v S F , x , y .
Obviously, the solutions of the sequential boundary value Problems (3) and (4) are fixed points of N . We will show that the operator N satisfies all conditions of Martelli’s fixed point theorem (Lemma 3). The proof is constructed in several steps.
Step 1. N ( x ) is convex for each x C ( [ a , b ] , R ) .
Indeed, if h 1 , h 2 belongs to N ( x ) , then there exist v 1 , v 2 S F , x such that for each t [ a , b ] , we have
h i ( t ) = I α v i ( t ) k I 1 x ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v i ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α v i ( b ) ] , i = 1 , 2 .
Let 0 λ 1 . Then, for each t [ a , b ] , we have
[ λ h 1 + ( 1 λ ) h 2 ] ( t ) = I α [ λ v 1 ( s ) + ( 1 λ ) v 2 ( s ) ] ( t ) k I 1 x ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i [ λ v 1 ( s ) + ( 1 λ ) v 2 ( s ) ] ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α [ λ v 1 ( s ) + ( 1 λ ) v 2 ( s ) ] ( b ) ] , i = 1 , 2 .
Then, we have
λ h 1 + ( 1 λ ) h 2 N ( x ) ,
since F has convex values.
Step 2.Bounded sets are mapped by N into bounded sets in C ( [ a , b ] , R ) .
For r > 0 , let a bounded ball in C ( [ a , b ] , R ) , defined by B r = { x C ( [ a , b ] , R ) : x r } . Then, for each h N ( x ) , x B r , there exists v S F , x such that
h ( t ) = I α v ( t ) k I 1 x ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α v ( b ) ] .
Then, for t [ a , b ] , we have
| h ( t ) | I α | v ( t ) | + | k | I 1 | x ( t ) | + ( t a ) γ 1 | Λ | { i = 1 m 2 | ζ i | I α + ϕ i | v ( θ i ) | + | k | i = 1 m 2 | ζ i | I ϕ i + 1 | x ( θ i | + | k | I 1 | x ( b ) | + I α | v ( b ) | } ( b a ) γ 1 | Λ | i = 1 m 2 | ζ i | ( θ i a ) α + ϕ i Γ ( α + ϕ i + 1 ) + ( b a ) α Γ ( α + 1 ) + ( b a ) α Γ ( α + 1 ) q + ( b a ) γ 1 | Λ | | k | i = 1 m 2 | ζ i | ( θ i a ) ϕ i + 1 Γ ( ϕ i + 2 ) + | k | ( b a ) + | k | ( b a ) x ,
and consequently,
N ( x ) q Ω + Ω 1 r .
Step 3.Bounded sets are mapped by N into equicontinuous sets of C ( [ a , b ] , R ) .
Let x B r and h N ( x ) . For each t [ a , b ] , there exists a function v S F , x such that
h ( t ) = I α v ( t ) k I 1 x ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α v ( b ) ] .
Let t 1 , t 2 [ a , b ] , t 1 < t 2 . Then, we have
h ( t 2 ) h ( t 1 ) 1 Γ ( α ) a t 1 [ ( t 2 s ) α 1 ( t 1 s ) α 1 ] v ( s ) d s + t 1 t 2 ( t 2 s ) α 1 v ( s ) d s + r ( t 2 t 1 ) + ( t 2 a ) γ 1 ( t 1 a ) γ 1 | Λ | i = 1 m 2 ζ i I α + ϕ i | v ( θ i ) | + I α | v ( b ) | + | k | i = 1 m 2 | ζ i | I ϕ i + 1 | x ( θ i ) | + | k | I 1 | x ( b ) | q Γ ( α + 1 ) [ 2 ( t 2 t 1 ) α + | ( t 2 a ) α ( t 1 a ) α | ] + r ( t 2 t 1 ) + ( τ 2 a ) γ 1 ( τ 1 a ) γ 1 | Λ | i = 1 m 2 ζ i I α + ϕ i | v ( θ i ) | + I α | v ( b ) | + | k | i = 1 m 2 | ζ i | I ϕ i + 1 | x ( θ i ) | + | k | I 1 | x ( b ) | .
The right side of the above inequality clearly tends to zero independently of x B r as t 1 t 2 . As a consequence of Steps 1–3 together with Arzel a ´ –Ascoli’s theorem, we conclude that N : C ( [ a , b ] , R ) P ( C ( [ a , b ] , R ) ) is completely continuous.
Now, we show that the operator N is upper semicontinuous. To prove this, [16] [Proposition 1.2] is enough to show that N has a closed graph.
Step 4.N has a closed graph.
Let x n x * , h n N ( x n ) and h n h * . We will show that h * N ( x * ) . Now, h n N ( x n ) implies that there exists v n S F , x n such that for each t [ a , b ] ,
h n ( t ) = I α v n ( t ) k I 1 x ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v n ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α v n ( b ) ] .
We show that there exists v * S F , x * such that for each t [ a , b ] ,
h * ( t ) = I α v * ( t ) k I 1 x ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v * ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α v * ( b ) ] .
Consider the continuous linear operator Θ : L 1 ( [ a , b ] , R ) C ( [ a , b ] ) by
v Θ ( v ) ( t ) = I α v ( t ) k I 1 x ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α v ( b ) ] .
Observe that h n ( t ) h * ( t ) 0 as n , and thus, it follows from a closed graph theorem [26] that Θ S F , x is a closed graph operator. Moreover, we have
h n Θ ( S F , x n ) .
Since x n x * , the closed graph theorem [26] implies that
h * ( t ) = I α v * ( t ) k I 1 x ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v * ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α v * ( b ) ] ,
for some v * S F , x * .
Step 5.We show the boundedness of the set
E = { x C ( [ a , b ] , R ) : λ x N ( x ) , λ > 1 . }
.
Let x E , then, λ x N ( x ) for some λ > 1 and there exists a function v S F , x such that
x ( t ) = 1 λ I α v ( t ) k 1 λ I 1 x ( t ) + 1 λ ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α v ( b ) ] .
For each t [ a , b ] , we have from Step 2 that
x q Ω + Ω 1 x ,
and consequently,
x q Ω 1 Ω 1 ,
which means that the set E is bounded. By using Lemma 3, we conclude that N has at least one fixed point and consequently, the sequential boundary value Problems (3)–(4) have a solution on [ a , b ] . □
Our second existence result is proved via Leray–Schauder nonlinear alternative for multivalued maps [17].
Theorem 5.
Assume that (17) and ( A 1 ) hold. In addition we assume that:
( A 3 )  
there exists a continuous, nondecreasing, subhomogeneous function ψ : [ 0 , ) ( 0 , ) and a function p C ( [ a , b ] , R + ) such that
F ( t , x , y ) P : = sup { | y | : y F ( t , x ) } p ( t ) ψ ( | x | + | y | ) f o r   e a c h ( t , x , y ) [ a , b ] × R × R ;
( A 4 )  
there exists a constant M > 0 such that
( 1 Ω 1 ) M L 1 ψ ( M ) p Ω > 1 .
Then the sequential boundary value Problems (3) and (4) has at least one solution on [ a , b ] .
Proof. 
Consider the operator N defined in the proof of Theorem 4. Let x λ N ( x ) for some λ ( 0 , 1 ) . We show there exists U C ( [ a , b ] , R ) , U is an open set, with x N ( x ) for all x U and for any λ ( 0 , 1 ) . Let λ ( 0 , 1 ) and x λ N ( x ) . Then, there exists v L 1 ( [ a , b ] , R ) with v S F , x such that, for t J , we have
x ( t ) = I α v ( t ) k I 1 x ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α v ( b ) ] .
In view of ( A 3 ) , we have for each t [ a , b ] , as in Theorem 3, that
| x ( t ) | L 1 ψ ( x ) p Ω + Ω 1 x
or
( 1 Ω 1 ) x L 1 ψ ( x ) p Ω 1 .
In view of ( A 4 ) , there exists M such that x M . Consider the set
U = { x C ( J , R ) : x < M } .
Proceeding as in the proof of Theorem 4, we have that the operator N : U ¯ P ( C ( [ a , b ] , R ) ) is a compact, upper semicontinuous, multivalued map with convex closed values. There is no x U such that x λ N ( x ) for some λ ( 0 , 1 ) , by the choice of U. Thus, N has a fixed point x U ¯ , which is a solution of the sequential boundary value Problems (3) and (4), by the nonlinear alternative of Leray–Schauder type ([17]). This completes the proof. □

4.2. The Lipschitz Case

In this subsection, we establish the existence of solutions for the sequential boundary value Problems (3) and (4) with a possible nonconvex-valued right-hand side by using a fixed point theorem for multivalued maps suggested by Covitz and Nadler [20].
Theorem 6.
Assume that
( A 4 )
F : [ a , b ] × R × R P c p ( R ) is such that F ( · , x , y ) : [ a , b ] P c p ( R ) is measurable for each x , y R .
( A 5 )
H d ( F ( t , x , y ) , F ( t , x ¯ , y ¯ ) ) m ( t ) ( | x x ¯ | + | y y ¯ | ) for almost all t [ a , b ] and x , y , x ¯ , y ¯ R with m C ( [ a , b ] , R + ) and d ( 0 , F ( t , 0 , 0 ) ) m ( t ) for almost all t [ a , b ] .
Then, if
L 1 Ω m + Ω 1 < 1 ,
the sequential boundary value Problems (3) and (4) have at least one solution on [ a , b ] .
Proof. 
We transform the boundary value Problems (3) and (4) into a fixed point problem by using the operator N : C ( [ a , b ] , R ) P ( C ( [ a , b ] , R ) ) defined in Theorem 4. We will show that the operator N satisfies the conditions of the Covitz–Nadler theorem.
Step I.N is nonempty and closed for every v S F , x .
The set-valued map F ( · , x ( · ) ) admits a measurable selection v : [ a , b ] R by the measurable selection theorem ([21] Theorem III.6).
We have, by the assumption ( A 5 ) , that
| v ( t ) | m ( t ) + m ( t ) ( | x ( t ) | + | I δ x ( t ) | ) m ( t ) + L 1 m ( t ) | x ( t ) | ,
i.e., v L 1 ( [ a , b ] , R ) and hence, F is integrably bounded, which means S F , x . Further, N ( x ) P c l ( C ( [ a , b ] , R ) ) for each x C ( [ a , b ] , R ) . Let { u n } n 0 N ( x ) be such that u n u ( n ) in C ( [ a , b ] , R ) . Then, u C ( [ a , b ] , R ) and there exists v n S F , x n such that, for each t [ a , b ] ,
h n ( t ) = I α v n ( t ) k I 1 x ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v n ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α v n ( b ) ] .
v n converges to v in L 1 ( [ a , b ] , R ) , as F has compact values. Thus, v S F , x and for each t [ a , b ] , we have
u n ( t ) v ( t ) = I α v ( t ) k I 1 x ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α v ( b ) ] .
Hence, u N ( x ) .
Step II. Next, we show that there exists 0 < θ < 1 ( θ = L 1 Ω m + Ω 1 ) such that
H d ( N ( x ) , N ( x ¯ ) ) θ x x ¯ f o r e a c h x , x ¯ C ( J , R ) .
Let x , x ¯ C ( [ a , b ] , R ) and h 1 N ( x ) . Then, there exists v 1 ( t ) F ( t , x ( t ) , y ( t ) ) such that, for each t [ a , b ] ,
h 1 ( t ) = I α v 1 ( t ) k I 1 x ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v 1 ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ( θ i ) + k I 1 x ( b ) I α v 1 ( b ) ] .
By ( A 5 ) , we have
H d ( F ( t , x , y ) , F ( t , x ¯ , y ¯ ) ) m ( t ) ( | x ( t ) x ¯ ( t ) | + | y ( t ) y ¯ ( t ) | ) .
So, there exists w ( t ) F ( t , x ¯ ( t ) , y ¯ ( t ) ) such that
| v 1 ( t ) w | m ( t ) ( | x ( t ) x ¯ ( t ) | + | y ( t ) y ¯ ( t ) | ) , t [ a , b ] .
Define U : J P ( R ) by
U ( t ) = { w R : | v 1 ( t ) w | m ( t ) ( | x ( t ) x ¯ ( t ) | + | y ( t ) y ¯ ( t ) | ) } .
There exists a function v 2 ( t ) , which is a measurable selection for U, since the multivalued operator U ( t ) F ( t , x ¯ ( t ) , y ¯ ( t ) ) is measurable ([21] Proposition III.4 ). So, v 2 ( t ) F ( t , x ¯ ( t ) , y ¯ ( t ) ) and for each t [ a , b ] , we have | v 1 ( t ) v 2 ( t ) | m ( t ) ( | x ( t ) x ¯ ( t ) | + | y ( t ) y ¯ ( t ) | ) .
For each t [ a , b ] , let us define
h 2 ( t ) = I α v 2 ( t ) k I 1 x ¯ ( t ) + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i v 2 ( θ i ) k i = 1 m 2 ζ i I ϕ i + 1 x ¯ ( θ i ) + k I 1 x ¯ ( b ) I α v 2 ( b ) ] .
Thus,
| h 1 ( t ) h 2 ( t ) | = I α | v 2 ( t ) v 1 ( t ) | + | k | ( b a ) x x ¯ + ( t a ) γ 1 Λ [ i = 1 m 2 ζ i I α + ϕ i | v 2 ( θ i ) v 1 ( θ i ) | + | k | i = 1 m 2 | ζ i | ( θ i a ) ϕ i + 1 Γ ( ϕ i + 2 ) x x ¯ + | k | ( b a ) x x ¯ + I α | v 2 ( b ) v 1 ( b ) | ] ( L 1 Ω m + Ω 1 ) x x ¯
Hence,
h 1 h 2 ( L 1 Ω m + Ω 1 ) x x ¯ .
Interchanging the roles of x and x ¯ , we obtain
H d ( N ( x ) , N ( x ¯ ) ) ( L 1 Ω m + Ω 1 ) x x ¯ .
Hence, N is a contraction. By the Covitz–Nadler fixed point theorem [20], N has a fixed point x, which is a solution of sequential boundary value Problems (3) and (4). This finishes the proof. □
Example 3.
Consider the nonlocal sequential Hilfer inclusion boundary value problem of the form:
H D 5 4 , 3 4 + 1 8 H D 1 4 , 3 4 x ( t ) F ( t , x ( t ) , I 4 3 x ( t ) ) , t 1 6 , 3 2 , x 1 6 = 0 , x 3 2 = 1 6 I 1 4 x 1 3 + 3 26 I 1 2 x 1 2 + 5 36 I 3 4 x 2 3 + 7 46 I 5 4 x 5 6 + 9 56 I 3 2 x 1 + 13 66 I 7 4 x 7 6 + 15 76 I 9 4 x 4 3 ,
where
F ( t , x , I 4 3 x ) = 0 , 6 6 t + 71 x 6 e I 4 3 x 2 1 + | x | 5 + 6 6 t + 77 sin I 4 3 x + 3 x 2 1 + x 2 + 6 6 t + 83 1 | x | + I 4 3 x + 3 .
Now, we set constants α = 5 / 4 ,   β = 3 / 4 ,   k = 1 / 8 ,   δ = 4 / 3 ,   a = 1 / 6 ,   b = 3 / 2 ,   m = 9 ,   ζ 1 = 1 / 6 ,   ζ 2 = 3 / 26 ,   ζ 3 = 5 / 36 ,   ζ 4 = 7 / 46 ,   ζ 5 = 9 / 56 ,   ζ 6 = 13 / 66 ,   ζ 7 = 15 / 76 ,   ϕ 1 = 1 / 4 ,   ϕ 2 = 1 / 2 ,   ϕ 3 = 3 / 4 ,   ϕ 4 = 5 / 4 ,   ϕ 5 = 3 / 2 ,   ϕ 6 = 7 / 4 ,   ϕ 7 = 9 / 4 ,   θ 1 = 1 / 3 ,   θ 2 = 1 / 2 ,   θ 3 = 2 / 3 ,   θ 4 = 5 / 6 ,   θ 5 = 1 ,   θ 6 = 7 / 6 ,   θ 7 = 4 / 3 ,   γ = ( 5 / 4 ) + ( 2 ( 5 / 4 ) ) ( 3 / 4 ) = 29 / 16 . From the information, we find that Λ 1.024874752 , Ω 2.992145342 , Ω 1 0.4023896011 , and L 1 2.232550581 .
From (23), we obtain
F ( t , w , z ) P 6 6 t + 71 | w | + | z | + 3 + 1 | w | + | z | + 3 ,
which can be set as p ( t ) = 6 / ( 6 t + 71 ) and ψ ( u ) = u + 3 + ( 1 / ( u + 3 ) ) , u = | w | + | z | . Consequently, we get p = 1 / 12 and ψ ( μ u ) μ ψ ( u ) , μ 1 , which is the subhomogeneous function. Therefore, the condition ( A 3 ) of Theorem 5 is fulfilled. Indeed, there exists a constant M > 41.10633915 , which makes condition ( A 4 ) true. Thus, the boundary value problem of Hilfer fractional differential Inclusions (22) and (23) have at least one solution on [ 1 / 6 , 3 / 2 ] .

5. Conclusions

In the present paper, we discussed a new class of boundary value problems for sequential fractional differential equations and inclusions involving Hilfer fractional derivatives, supplemented with nonlocal integro-multipoint boundary conditions. In the single-valued case, existence and uniqueness results are established by using the classical fixed point theorems of Banach and Krasnosel’skiĭ and the nonlinear alternative of Leray–Schauder. In the multivalued case and convex-valued multivalued maps, we proved two existence results by applying Martelli’s fixed point theorem and the nonlinear alternative for Kakutani maps. In the case of possible nonconvex-valued maps, we obtained an existence result via the Covitz–Nadler fixed point theorem for contractive maps. The obtained results are well illustrated by numerical examples. The results obtained in the present paper are new and significantly contribute to the existing literature on the topic.

Author Contributions

Conceptualization, S.K.N. and J.T.; methodology, N.P., S.K.N., J.T., and K.N.; formal analysis, N.P., S.K.N., J.T., and K.N.; 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-GRAD-59-03.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

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 Science Inc.: New York, NY, USA, 2006. [Google Scholar]
  3. Lakshmikantham, V.; Leela, S.; Devi, J.V. Theory of Fractional Dynamic Systems; Cambridge Scientific Publishers: Cambridge, 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. Zhou, Y. Basic Theory of Fractional Differential Equations; World Scientific: Singapore, 2014. [Google Scholar]
  9. Hilfer, R. (Ed.) Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  10. Hilfer, R. Experimental evidence for fractional time evolution in glass forming materials. J. Chem. Phys. 2002, 284, 399–408. [Google Scholar] [CrossRef]
  11. Hilfer, R.; Luchko, Y.; Tomovski, Z. Operational method for the solution of fractional differential equations with generalized Riemann-Liouvill fractional derivatives. Frac. Calc. Appl. Anal. 2009, 12, 299–318. [Google Scholar]
  12. Furati, K.M.; Kassim, N.D.; Tatar, N.E. Existence and uniqueness for a problem involving Hilfer fractional derivative. Comput. Math. Appl. 2012, 64, 1616–1626. [Google Scholar] [CrossRef] [Green Version]
  13. Gu, H.; Trujillo, J.J. Existence of mild solution for evolution equation with Hilfer fractional derivative. Appl. Math. Comput. 2015, 257, 344–354. [Google Scholar] [CrossRef]
  14. Wang, J.; Zhang, Y. Nonlocal initial value problems for differential equations with Hilfer fractional derivative. Appl. Math. Comput. 2015, 266, 850–859. [Google Scholar] [CrossRef]
  15. Asawasamrit, S.; Kijjathanakorn, A.; Ntouyas, S.K.; Tariboon, J. Nonlocal boundary value problems for Hilfer fractional differential equations. Bull. Korean Math. Soc. 2018, 55, 1639–1657. [Google Scholar]
  16. Deimling, K. Multivalued Differential Equations; De Gruyter: Berlin, Germany, 1992. [Google Scholar]
  17. Granas, A.; Dugundji, J. Fixed Point Theory; Springer: New York, NY, USA, 2003. [Google Scholar]
  18. Krasnosel’skiĭ, M.A. Two remarks on the method of successive approximations. Uspekhi Mat. Nauk 1955, 10, 123–127. [Google Scholar]
  19. Martelli, M. A Rothe’s theorem for non compact acyclic-valued maps. Boll. Un. Mat. Ital. 1975, 4, 70–76. [Google Scholar]
  20. Covitz, H.; Nadler, S.B., Jr. Multivalued contraction mappings in generalized metric spaces. Israel J. Math. 1970, 8, 5–11. [Google Scholar] [CrossRef]
  21. Castaing, C.; Valadier, M. Convex Analysis and Measurable Multifunctions; Lecture Notes in Mathematics 580; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 1977. [Google Scholar]
  22. Górniewicz, L. Topological Fixed Point Theory of Multivalued Mappings; Mathematics and its Applications, 495; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1999. [Google Scholar]
  23. Hu, S.; Papageorgiou, N. Handbook of Multivalued Analysis, Volume I: Theory; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1997. [Google Scholar]
  24. Wongcharoen, A.; Ntouyas, S.K.; Tariboon, J. Nonlocal boundary value problems for Hilfer type pantograph fractional differential equations and inclusions. Adv. Differ. Equ. 2020, 2020, 279. [Google Scholar] [CrossRef]
  25. Ntouyas, S.K.; Ahmad, B.; Alsaedi, A. Fractional order single-valued and multi-valued problems with integro-multistrip-multipoint boundary conditions. Fractal Fract 2019, 4, 31. [Google Scholar] [CrossRef]
  26. Lasota, A.; Opial, Z. An application of the Kakutani-Ky Fan theorem in the theory of ordinary differential equations. Bull. Acad. Pol. Sci. Ser. Sci. Math. Astron. Phys. 1965, 13, 781–786. [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

Phuangthong, N.; Ntouyas, S.K.; Tariboon, J.; Nonlaopon, K. Nonlocal Sequential Boundary Value Problems for Hilfer Type Fractional Integro-Differential Equations and Inclusions. Mathematics 2021, 9, 615. https://doi.org/10.3390/math9060615

AMA Style

Phuangthong N, Ntouyas SK, Tariboon J, Nonlaopon K. Nonlocal Sequential Boundary Value Problems for Hilfer Type Fractional Integro-Differential Equations and Inclusions. Mathematics. 2021; 9(6):615. https://doi.org/10.3390/math9060615

Chicago/Turabian Style

Phuangthong, Nawapol, Sotiris K. Ntouyas, Jessada Tariboon, and Kamsing Nonlaopon. 2021. "Nonlocal Sequential Boundary Value Problems for Hilfer Type Fractional Integro-Differential Equations and Inclusions" Mathematics 9, no. 6: 615. https://doi.org/10.3390/math9060615

APA Style

Phuangthong, N., Ntouyas, S. K., Tariboon, J., & Nonlaopon, K. (2021). Nonlocal Sequential Boundary Value Problems for Hilfer Type Fractional Integro-Differential Equations and Inclusions. Mathematics, 9(6), 615. https://doi.org/10.3390/math9060615

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