Next Article in Journal
Coincidence Point Results for Multivalued Suzuki Type Mappings Using θ-Contraction in b-Metric Spaces
Next Article in Special Issue
A Mollification Regularization Method for the Inverse Source Problem for a Time Fractional Diffusion Equation
Previous Article in Journal
The Impact of Sinusoidal Surface Temperature on the Natural Convective Flow of a Ferrofluid along a Vertical Plate
Previous Article in Special Issue
Positive Solutions for a System of Fractional Integral Boundary Value Problems of Riemann–Liouville Type Involving Semipositone Nonlinearities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On a Generalized Langevin Type Nonlocal Fractional Integral Multivalued Problem

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, 451 10 Ioannina, Greece
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(11), 1015; https://doi.org/10.3390/math7111015
Submission received: 29 August 2019 / Revised: 17 October 2019 / Accepted: 21 October 2019 / Published: 25 October 2019

Abstract

:
We establish sufficient criteria for the existence of solutions for a nonlinear generalized Langevin-type nonlocal fractional-order integral multivalued problem. The convex and non-convex cases for the multivalued map involved in the given problem are considered. Our results rely on Leray–Schauder nonlinear alternative for multivalued maps and Covitz and Nadler’s fixed point theorem. Illustrative examples for the main results are included.

1. Introduction

Fractional calculus is the extension of classical calculus which deals with differential and integral operators of fractional order. It has evolved into a significant and popular branch of mathematical analysis owing to its extensive applications in the mathematical modeling of applied and technical problems. The literature on fractional calculus is now much enriched and covers a wide range of interesting results, for instance [1,2,3,4,5,6]. For a comprehensive treatment of Hadamard-type fractional differential equations and inclusions, we refer the reader to the text [7].
The Langevin equation is found to be an effective tool to describe stochastic problems in fluctuating situations. A modified type of this equation is used in various functional approaches for fractal media. A variety of boundary value problems involving the Langevin equation have been investigated by several authors. In [8], existence and uniqueness results for a nonlinear Langevin equation involving two fractional orders supplemented with three-point boundary conditions were obtained. An impulsive boundary value problem for a nonlinear Langevin equation involving two different fractional derivatives was investigated in [9]. Some existing results for Langevin fractional differential inclusions with two indices were derived in [10]. In [11], the authors proved the existence of and uniqueness results for an anti-periodic boundary value problem of a system of Langevin fractional differential equations. In [12], the authors investigated a nonlinear fractional Langevin equation with anti-periodic boundary conditions by applying coupled fixed point theorems. In a recent work [13], the authors obtained some existence results for a fractional Langevin equation with nonlinearity depending on Riemann–Liouville fractional integral, and complemented with nonlocal multi-point and multi-strip boundary conditions.
In the present paper, we study the existence of solutions for a nonlinear generalized Langevin type nonlocal fractional-order integral multivalued problem given by
c ρ D a + α ( c ρ D a + β + λ ) x ( t ) F ( t , x ( t ) ) , t J : = [ a , T ] , λ R , x ( a ) = 0 , x ( η ) = 0 , x ( T ) = μ ρ I a + γ x ( ξ ) , a < η < ξ < T , μ R ,
where c ρ D a + α , c ρ D a + β denote the Caputo-type generalized fractional differential operators of order 1 < α 2 , 0 < β < 1 , ρ > 0 , respectively, F : [ a , T ] × R P ( R ) is a multi-valued map ( P ( R ) is the family of all nonempty subsets of R ), ρ I a + γ is the generalized fractional integral operator of order γ > 0 and ρ > 0 . Here we emphasize that the single-valued analogue of the problem (1) was discussed in [14].
The rest of the paper is arranged as follows. The background material related to our work is outlined in Section 3. The existence results for the problem (1) are presented in Section 3. The first result for the problem (1), associated with the convex valued mutivalued map, is derived with the aid of Leray–Schauder nonlinear alternative for multivalued maps, while the result for non-convex valued map for the problem (1) is proved by applying a fixed point theorem due to Covitz and Nadler. Section 4 contains the illustrative examples for the main results. We summarize the work established in this paper, and its implications, in the last section.

2. Preliminaries

Define by X c p ( a , b ) the space of all complex-valued Lebesgue measurable functions ϕ on ( a , b ) equipped with the norm:
ϕ X c p = a b | x c ϕ ( x ) | p d x x 1 / p < , c R , 1 p .
Let A C δ n [ a , b ] denote the class of all absolutely continuous functions g possessing δ n 1 -derivative ( δ n 1 g A C ( [ a , b ] , R ) ) , endowed with the norm g A C δ n = k = 0 n 1 δ k g C .
Definition 1.
The left-sided and right-sided generalized fractional integrals for g X c p ( a , b ) of order β > 0 and ρ > 0 , denoted by ρ I a + β g and ρ I b β g respectively, are defined by [15]
( ρ I a + β g ) ( t ) = ρ 1 β Γ ( β ) a t s ρ 1 ( t ρ s ρ ) 1 β g ( s ) d s , < a < t < b < ,
( ρ I b β g ) ( t ) = ρ 1 α Γ ( β ) t b s ρ 1 ( s ρ t ρ ) 1 β g ( s ) d s , < a < t < b < .
Definition 2.
Let β > 0 , n = [ β ] + 1 and ρ > 0 . We define the generalized fractional derivatives, associated with the generalized fractional integrals (2) and (3), for 0 a < t < b < , as follows [16]:
( ρ D a + β g ) ( t ) = t 1 ρ d d t n ( ρ I a + n β g ) ( t ) = ρ β n + 1 Γ ( n β ) t 1 ρ d d t n a t s ρ 1 ( t ρ s ρ ) β n + 1 g ( s ) d s ,
( ρ D b β g ) ( t ) = t 1 ρ d d t n ( ρ I b n β g ) ( t ) = ρ β n + 1 Γ ( n β ) t 1 ρ d d t n t b s ρ 1 ( s ρ t ρ ) β n + 1 g ( s ) d s ,
provided the integrals in the above expressions exist.
Definition 3.
Let g A C δ n [ a , b ] and β > 0 , n = [ β ] + 1 . Then the Caputo-type generalized fractional derivatives c ρ D a + β g and c ρ D b β g are respectively defined via (4) and (5) by [17]
c ρ D a + β g ( x ) = ρ D a + β g ( t ) k = 0 n 1 δ k g ( a ) k ! t ρ a ρ ρ k ( x ) , δ = x 1 ρ d d x ,
c ρ D b β g ( x ) = ρ D b β g ( t ) k = 0 n 1 ( 1 ) k δ k g ( b ) k ! b ρ t ρ ρ k ( x ) , δ = x 1 ρ d d x .
Remark 1.
The left and right generalized Caputo derivatives of order β for g A C δ n [ a , b ] , are respectively given by [17]
c ρ D a + β g ( t ) = 1 Γ ( n β ) a t t ρ s ρ ρ n β 1 ( δ n g ) ( s ) d s s 1 ρ ,
c ρ D b β g ( t ) = 1 Γ ( n β ) t b s ρ t ρ ρ n α 1 ( 1 ) n ( δ n g ) ( s ) d s s 1 ρ .
Lemma 1.
Let g A C δ n [ a , b ] or C δ n [ a , b ] . Then, for β R , the following results hold [17]:
ρ I a + β c ρ D a + β g ( x ) = g ( x ) k = 0 n 1 ( δ k g ) ( a ) k ! x ρ a ρ ρ k ,
ρ I b β c ρ D b β g ( x ) = g ( x ) k = 0 n 1 ( 1 ) k ( δ k g ) ( a ) k ! b ρ x ρ ρ k .
In particular, for 0 < β 1 , we have
ρ I a + β c ρ D a + β g ( x ) = g ( x ) g ( a ) , ρ I b β c ρ D b β g ( x ) = g ( x ) g ( b ) .
We need the following known lemma [14] in the sequel.
Lemma 2.
Let h C ( [ a , T ] , R ) and x A C δ 3 ( J ) . Then the unique solution of linear problem:
c ρ D a + α ( c ρ D a + β + λ ) x ( t ) = h ( t ) , t J : = [ a , T ] , x ( a ) = 0 , x ( η ) = 0 , x ( T ) = μ ρ I a + γ x ( ξ ) , a < η < ξ < T ,
is given by:
x ( t ) = ρ I a + α + β h ( t ) λ ρ I a + β x ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β h ( T ) λ ρ I a + β x ( T ) μ ρ I a + α + β + γ h ( ξ ) + μ λ ρ I a + β + γ x ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β h ( η ) λ ρ I a + β x ( η ) ,
where it is assumed that
Ω = ( T ρ a ρ ) β ( T ρ η ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ η ρ ) γ ( η ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) 0 .

3. Main Results

We begin this section with the definition of a solution for the multi-valued problem (1).
Definition 4.
A function x C ( J , R ) is called a solution of the problem (1) if we can find a function v L 1 ( J , R ) with v ( t ) F ( t , x ) a.e. on J such that x ( a ) = 0 , x ( η ) = 0 , x ( T ) = μ ρ I a + γ x ( ξ ) and
x ( t ) = ρ I a + α + β v ( t ) λ ρ I a + β x ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v ( T ) λ ρ I a + β x ( T ) μ ρ I a + α + β + γ v ( ξ ) + μ λ ρ I a + β + γ x ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v ( η ) λ ρ I a + β x ( η ) .
For the sake of computational convenience, we set
Λ 1 = ( T ρ a ρ ) α + β ρ α + β Γ ( α + β + 1 ) 1 + ζ 1 ρ β + 1 Γ ( β + 2 ) | Ω | + | μ | ( ξ ρ a ρ ) α + β + γ ζ 1 ρ α + 2 β + γ + 1 Γ ( α + β + γ + 1 ) Γ ( β + 2 ) | Ω | + ( η ρ a ρ ) α ζ 2 ρ α + β Γ ( α + β + 1 ) | Ω | ,
Λ 2 = | λ | ( T ρ a ρ ) β ρ β Γ ( β + 1 ) 1 + ζ 1 ρ β + 1 Γ ( β + 2 ) | Ω | + | μ | | λ | ( ξ ρ a ρ ) β + γ ζ 1 ρ 2 β + γ + 1 Γ ( β + γ + 1 ) Γ ( β + 2 ) | Ω | + | λ | ζ 2 ρ β Γ ( β + 1 ) | Ω | ,
where
ζ 1 : = max t [ a , T ] | ( t ρ a ρ ) β ( η ρ t ρ ) | ,
ζ 2 : = max t [ a , T ] | ( t ρ a ρ ) β ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) | .
We define the set of selections of F by S F , x : = { y L 1 ( J , R ) : y ( t ) F ( t , x ( t ) ) o n J } for each x C ( J , R ) .

3.1. The Upper Semicontinuous Case

In the following result, we assume that the multivalued map F is convex-valued and apply Leray–Schauder nonlinear alternative for multivalued maps [18] to prove the existence of solutions for the problem at hand.
Theorem 1.
Assume that:
( A 1 )
F : J × R P c p , c ( R ) is L 1 -Carathéodory, where P c p , c ( R ) = { Y P ( R ) : Y is compact and convex};
( A 2 )
there exist a function P C ( J , R + ) and a continuous nondecreasing function Q : [ 0 , ) ( 0 , ) such that F ( t , x ) P : = sup { | y | : y F ( t , x ) } P ( t ) Q ( | x | ) f o r e a c h ( t , x ) J × R ;
( A 3 )
there exists a constant M > 0 such that
( 1 Λ 2 ) M Λ 1 P Q ( M ) > 1 , Λ 2 < 1 ,
where Λ 1 and Λ 2 are respectively given by (14) and (15).
Then the problem (1) has at least one solution on J.
Proof. 
Let us first convert the problem (1) into a fixed point problem by introducing a multivalued map: N : C ( J , R ) P ( C ( J , R ) ) as
N ( x ) = h C ( J , R ) : h ( t ) = ρ I a + α + β v ( t ) λ ρ I a + β x ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v ( T ) λ ρ I a + β x ( T ) μ ρ I a + α + β + γ v ( ξ ) + μ λ ρ I a + β + γ x ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v ( η ) λ ρ I a + β x ( η ) ,
for v S F , x .
It is clear that fixed points of N are solutions of problem (1). So we need to verify that the operator N satisfies all the conditions of Leray–Schauder nonlinear alternative [18]. This will be done in several steps.
Step 1. N ( x ) is convex for each x C ( J , 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 J , we have
h i ( t ) = ρ I a + α + β v i ( t ) λ ρ I a + β x ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v i ( T ) λ ρ I a + β x ( T ) μ ρ I a + α + β + γ v i ( ξ ) + μ λ ρ I a + β + γ x ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v i ( η ) λ ρ I a + β x ( η ) , i = 1 , 2 .
Let t J and θ ( 0 , 1 ) . Then
[ θ h 1 + ( 1 θ ) h 2 ] ( t ) = ρ I a + α + β [ θ v 1 ( s ) + ( 1 θ ) v 2 ( s ) ] ( t ) λ ρ I a + β x ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β [ θ v 1 ( s ) + ( 1 θ ) v 2 ( s ) ] ( T ) λ ρ I a + β x ( T ) μ ρ I a + α + β + γ [ θ v 1 ( s ) + ( 1 θ ) v 2 ( s ) ] ( ξ ) + μ λ ρ I a + β + γ x ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β [ θ v 1 ( s ) + ( 1 θ ) v 2 ( s ) ] ( η ) λ ρ I a + β x ( η ) .
Since F has convex values ( S F , x is convex), therefore, θ h 1 + ( 1 θ ) h 2 N ( x ) .
Step 2. N ( x ) maps bounded sets (balls) into bounded sets in C ( J , R ) .
Let B r = { x C ( J , R ) : x r } be a bounded ball in C ( J , R ) , where r is a positive number. Then, for each h N ( x ) , x B r , there exists v S F , x such that
h ( t ) = ρ I a + α + β v ( t ) λ ρ I a + β x ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v ( T ) λ ρ I a + β x ( T ) μ ρ I a + α + β + γ v ( ξ ) + μ λ ρ I a + β + γ x ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v ( η ) λ ρ I a + β x ( η ) .
In view of ( H 2 ) , for each t J , we find that
| h ( t ) | ρ I a + α + β | v ( t ) | + | λ | ρ I a + β | x ( t ) | + | ( t ρ a ρ ) β ( η ρ t ρ ) | ρ β + 1 Γ ( β + 2 ) | Ω | { ρ I a + α + β | v ( T ) | + λ ρ I a + β | x ( T ) | + | μ | ρ I a + α + β + γ | v ( ξ ) | + | μ λ | ρ I a + β + γ | x ( ξ ) | } + | ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) + μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) | ρ I a + α + β | v ( η ) | + | λ | ρ I a + β | x ( η ) | P Q ( x ) ( ( T ρ a ρ ) α + β ρ α + β Γ ( α + β + 1 ) 1 + ζ 1 ρ β + 1 Γ ( β + 2 ) | Ω | + | μ | ( ξ ρ a ρ ) α + β + γ ζ 1 ρ α + 2 β + γ + 1 Γ ( α + β + γ + 1 ) Γ ( β + 2 ) | Ω | + ( η ρ a ρ ) α ζ 2 ρ α + β Γ ( α + β + 1 ) | Ω | ) + x ( | λ | ( T ρ a ρ ) β ρ β Γ ( β + 1 ) 1 + ζ 1 ρ β + 1 Γ ( β + 2 ) | Ω | + | μ | | λ | ( ξ ρ a ρ ) β + γ ζ 1 ρ 2 β + γ + 1 Γ ( β + γ + 1 ) Γ ( β + 2 ) | Ω | + | λ | ζ 2 ρ β Γ ( β + 1 ) | Ω | ) = Λ 1 P Q ( x ) + Λ 2 x ,
which leads to h Λ 1 P Q ( r ) + Λ 2 r .
Step 3. N ( x ) maps bounded sets into equicontinuous sets of C ( J , R ) .
Let x be any element in B r and h N ( x ) . Then there exists a function v S F , x such that, for each t J we have
h ( t ) = ρ I a + α + β v ( t ) λ ρ I a + β x ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v ( T ) λ ρ I a + β x ( T ) μ ρ I a + α + β + γ v ( ξ ) + μ λ ρ I a + β + γ x ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v ( η ) λ ρ I a + β x ( η ) .
Let τ 1 , τ 2 J , τ 1 < τ 2 . Then
| h ( t 2 ) h ( t 1 ) | | ρ 1 ( α + β ) Γ ( α + β ) 0 t 1 s ρ 1 [ ( t 2 ρ s ρ ) α + β 1 ( t 1 ρ s ρ ) α + β 1 ] v ( s ) d s + t 1 t 2 s ρ 1 ( t 2 ρ s ρ ) α + β 1 v ( s ) d s + ( t 2 ρ a ρ ) β ( η ρ t 2 ρ ) ρ β + 1 Γ ( β + 2 ) Ω ( t 1 ρ a ρ ) β ( η ρ t 1 ρ ) ρ β + 1 Γ ( β + 2 ) Ω ρ I a + α + β v ( T ) μ ρ I a + α + β + γ v ( ξ ) [ ( t 2 ρ a ρ ) β Ω ( η ρ a ρ ) β ( T ρ a ρ ) β ( T ρ t 2 ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t 2 ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ( t 1 ρ a ρ ) β Ω ( η ρ a ρ ) β ( T ρ a ρ ) β ( T ρ t 1 ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t 1 ρ ) γ ( t 1 ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ] × × ρ I a + α + β v ( η ) | + | λ ρ 1 β Γ ( β ) 0 t 1 s ρ 1 [ ( t 2 ρ s ρ ) β 1 ( t 1 ρ s ρ ) β 1 ] x ( s ) d s + t 1 t 2 s ρ 1 ( t 2 ρ s ρ ) β 1 x ( s ) d s + ( t 2 ρ a ρ ) β ( η ρ t 2 ρ ) ρ β + 1 Γ ( β + 2 ) Ω ( t 1 ρ a ρ ) β ( η ρ t 1 ρ ) ρ β + 1 Γ ( β + 2 ) Ω λ ρ I a + β x ( T ) + μ λ ρ I a + β + γ x ( ξ ) + [ λ ( t 2 ρ a ρ ) β Ω ( η ρ a ρ ) β ( T ρ a ρ ) β ( T ρ t 2 ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t 2 ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) λ ( t 1 ρ a ρ ) β Ω ( η ρ a ρ ) β ( T ρ a ρ ) β ( T ρ t 1 ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t 1 ρ ) γ ( t 1 ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ] × × ρ I a + β x ( η ) | P Q ( r ) ρ α + β Γ ( α + β + 1 ) | t 2 ρ ( α + β ) t 1 ρ ( α + β ) | + 2 ( t 2 ρ t 1 ρ ) α + β + | ( t 2 ρ a ρ ) β ( η ρ t 2 ρ ) ρ β + 1 Γ ( β + 2 ) Ω ( t 1 ρ a ρ ) β ( η ρ t 1 ρ ) ρ β + 1 Γ ( β + 2 ) Ω | × P Q ( r ) ( T ρ a ρ ) α + β ρ α + β Γ ( α + β + 1 ) + | μ | ( ξ ρ a ρ ) α + β + γ ρ α + β + γ Γ ( α β + γ + 1 ) + | ( t 2 ρ a ρ ) β Ω ( η ρ a ρ ) β ( T ρ a ρ ) β ( T ρ t 2 ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t 2 ρ ) γ ( t 2 ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ( t 1 ρ a ρ ) β Ω ( η ρ a ρ ) β ( T ρ a ρ ) β ( T ρ t 1 ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t 1 ρ ) γ ( t 1 ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) | × × P Q ( r ) ( η ρ a ρ ) α + β ρ α + β Γ ( α + β + 1 ) + r ρ β Γ ( β + 1 ) | t 2 ρ β t 1 ρ β | + 2 ( t 2 ρ t 1 ρ ) β + | ( t 2 ρ a ρ ) β ( η ρ t 2 ρ ) ρ β + 1 Γ ( β + 2 ) Ω ( t 1 ρ a ρ ) β ( η ρ t 1 ρ ) ρ β + 1 Γ ( β + 2 ) Ω | | λ | r ( T ρ a ρ ) β ρ β Γ ( β + 1 ) + | μ | ( ξ ρ a ρ ) β + γ ρ β + γ Γ ( β + γ + 1 ) + | λ ( t 2 ρ a ρ ) β Ω ( η ρ a ρ ) β ( T ρ a ρ ) β ( T ρ t 2 ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t 2 ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) λ ( t 1 ρ a ρ ) β Ω ( η ρ a ρ ) β ( T ρ a ρ ) β ( T ρ t 1 ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t 1 ρ ) γ ( t 1 ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) | × × ( η ρ a ρ ) β ρ β Γ ( β + 1 ) 0 when t 1 t 2 , independently of x B r .
Combining the outcome of Steps 1–3 with Arzelá-Ascoli theorem leads to the conclusion that N : C ( J , R ) P ( C ( J , R ) ) is completely continuous.
Next, we show that N has a closed graph. Then it will follow by Proposition 1.2 in [19] that the operator N is u.s.c.
Step 4.N has a closed graph.
Suppose that there exists x n x , h n N ( x n ) and h n h . Then we have to establish that h N ( x ) . Since h n N ( x n ) , there exists v n S F , x n . In consequence, for each t J , we get
h n ( t ) = ρ I a + α + β v n ( t ) λ ρ I a + β x n ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v n ( T ) λ ρ I a + β x n ( T ) μ ρ I a + α + β + γ v n ( ξ ) + μ λ ρ I a + β + γ x n ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v n ( η ) λ ρ I a + β x n ( η ) .
Next we show that there exists v S F , x such that, for each t J ,
h ( t ) = ρ I a + α + β v ( t ) λ ρ I a + β x ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v ( T ) λ ρ I a + β x ( T ) μ ρ I a + α + β + γ v ( ξ ) + μ λ ρ I a + β + γ x ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v ( η ) λ ρ I a + β x ( η ) .
Consider the continuous linear operator Θ : L 1 ( J , R ) C ( J , R ) given by
v Θ ( v ) ( t ) = ρ I a + α + β v ( t ) λ ρ I a + β x ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v ( T ) λ ρ I a + β x ( T ) μ ρ I a + α + β + γ v ( ξ ) + μ λ ρ I a + β + γ x ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v ( η ) λ ρ I a + β x ( η ) .
Notice that h n ( t ) h ( t ) 0 as n . So we deduce by a closed graph result obtained in [20] that Θ S F , x is a closed graph operator. Furthermore, h n Θ ( S F , x n ) . Since x n x , therefore we have
h ( t ) = ρ I a + α + β v ( t ) λ ρ I a + β x ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v ( T ) λ ρ I a + β x ( T ) μ ρ I a + α + β + γ v ( ξ ) + μ λ ρ I a + β + γ x ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v ( η ) λ ρ I a + β x ( η ) ,
for some v S F , x
Step 5.There exists an open set V C ( J , R ) with x θ N ( x ) for any θ ( 0 , 1 ) and all x V .
Take θ ( 0 , 1 ) , x θ N ( x ) and t J . Then we show that there exists v L 1 ( J , R ) with v S F , x such that
x ( t ) = θ ρ I a + α + β v ( t ) θ λ ρ I a + β x ( t ) + θ ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v ( T ) λ ρ I a + β x ( T ) μ ρ I a + α + β + γ v ( ξ ) + μ λ ρ I a + β + γ x ( ξ ) } θ ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v ( η ) λ ρ I a + β x ( η ) .
Using the computations done in Step 2, for each t J , we get
| x ( t ) | Λ 1 P Q ( x ) + Λ 2 x ,
which yields
( 1 Λ 2 ) x Λ 1 P Q ( x ) 1 .
By ( A 3 ) , there exists M such that x M . Define a set
V = { x C ( J , R ) : x < M } .
Observe that the operator N : V ¯ P ( C ( J , R ) ) is a compact multivalued map, u.s.c. with convex closed values. With the given choice of V , it is not possible to find x V satisfying x θ N ( x ) for some θ ( 0 , 1 ) . Consequently, by the nonlinear alternative of Leray–Schauder type [18], the operator N has a fixed point x V ¯ , which corresponds to a solution of the problem (1). This finishes the proof. □

3.2. The Lipschitz Case

Let ( X , d ) denote a metric space induced from the normed space ( X ; · ) . Let H d : P ( X ) × P ( X ) R { } be defined by H d ( A 1 , A 2 ) = max { sup a 1 A 1 d ( a 1 , A 2 ) , sup a 2 A 2 d ( A 1 , a 2 ) } , where d ( A 1 , a 2 ) = inf a 1 A 1 d ( a 1 ; a 2 ) and d ( a 1 , A 2 ) = inf a 2 A 2 d ( a 1 ; a 2 ) . Then ( P b , c l ( X ) , H d ) is a metric space (see [21]), where P b , c l ( X ) = { Y P ( X ) : Y is bounded and closedg},
The following result deals with the non-convex valued case of the problem (1) and is based on Covitz and Nadler’s fixed point theorem [22]: “If N : X P c l ( X ) is a contraction, then F i x N , where P c l ( X ) = { Y P ( X ) : Y is closed } .
Theorem 2.
Assume that
( A 4 )
F : J × R P c p ( R ) is such that F ( · , x ) : J P c p ( R ) is measurable for each x R , where P c p ( R ) = { Y P ( R ) : Y is compact } ;
( A 5 )
H d ( F ( t , x ) , F ( t , x ^ ) ) ϖ ( t ) | x x ^ | for almost all t J and x , x ^ R with ϖ C ( J , R + ) and d ( 0 , F ( t , 0 ) ) ϖ ( t ) for almost all t J .
Then the problem (1) has at least one solution on J if
ϖ Λ 1 + Λ 2 < 1 ,
where Λ 1 and Λ 2 are respectively given by (14) and (15).
Proof. 
Let us verify that the operator N : C ( J , R ) P ( C ( J , R ) ) , defined in the proof of the last theorem, satisfies the hypothesis of Covitz and Nadler fixed point theorem [22]. We establish it in two steps.
Step I. N ( x ) is nonempty and closed for every v S F , x .
Since the set-valued map F ( · , x ( · ) ) is measurable, it admits a measurable selection v : J R by the measurable selection theorem ([23], Theorem III.6). By ( A 5 ) , we have
| v ( t ) | ϖ ( t ) ( 1 + | x ( t ) | ) ,
that is, v L 1 ( J , R ) . So F is integrably bounded. Therefore, S F , x .
Now we establish that N ( x ) is closed for each x C ( J , R ) . Let { u n } n 0 N ( x ) be such that u n u as n in C ( J , R ) . Then u C ( J , R ) and we can find v n S F , x n such that, for each t J ,
u n ( t ) = ρ I a + α + β v n ( t ) λ ρ I a + β x n ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v ( T ) λ ρ I a + β x n ( T ) μ ρ I a + α + β + γ v n ( ξ ) + μ λ ρ I a + β + γ x n ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v n ( η ) λ ρ I a + β x n ( η ) .
As F has compact values, we can pass onto a subsequence (if necessary) to obtain that v n converges to v in L 1 ( J , R ) . So v S F , x . Then, for each t J , we get
u n ( t ) v ( t ) = ρ I a + α + β v ( t ) λ ρ I a + β x ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v ( T ) λ ρ I a + β x ( T ) μ ρ I a + α + β + γ v ( ξ ) + μ λ ρ I a + β + γ x ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v ( η ) λ ρ I a + β x ( η ) ,
which implies that u N ( x ) .
Step II. We establish that there exists 0 < θ ^ < 1 ( θ ^ = Λ 1 ϖ + Λ 2 ) satisfying
H d ( N ( x ) , N ( x ^ ) ) θ ^ x x ^ for each x , x ^ C ( J , R ) .
Let us take x , x ^ C ( J , R ) and h 1 N ( x ) . Then there exists v 1 ( t ) F ( t , x ( t ) ) such that, for each t J ,
h 1 ( t ) = ρ I a + α + β v 1 ( t ) λ ρ I a + β x ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v 1 ( T ) λ ρ I a + β x ( T ) μ ρ I a + α + β + γ v 1 ( ξ ) + μ λ ρ I a + β + γ x ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v 1 ( η ) λ ρ I a + β x ( η ) .
By ( A 5 ) , we have that H d ( F ( t , x ) , F ( t , x ^ ) ) ϖ ( t ) | x ( t ) x ^ ( t ) | . So, there exists w ( t ) F ( t , x ^ ( t ) ) satisfying | v 1 ( t ) w | ϖ ( t ) | x ( t ) x ^ ( t ) | , t J .
Define W : J P ( R ) by
W ( t ) = { w R : | v 1 ( t ) w | ϖ ( t ) | x ( t ) x ^ ( t ) | } .
As the multivalued operator W ( t ) F ( t , x ^ ( t ) ) is measurable by Proposition III.4 [23], we can find a function v 2 ( t ) which is a measurable selection for W . So v 2 ( t ) F ( t , x ^ ( t ) ) and for each t J , we have | v 1 ( t ) v 2 ( t ) | ϖ ( t ) | x ( t ) x ^ ( t ) | . For each t J , we define
h 2 ( t ) = ρ I a + α + β v 2 ( t ) λ ρ I a + β x ^ ( t ) + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β v 2 ( T ) λ ρ I a + β x ^ ( T ) μ ρ I a + α + β + γ v 2 ( ξ ) + μ λ ρ I a + β + γ x ^ ( ξ ) } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β v 2 ( η ) λ ρ I a + β x ^ ( η ) .
As a result, we get
h 1 ( t ) h 2 ( t ) = | ρ I a + α + β [ v 2 ( s ) v 1 ( s ) ] ( t ) λ ρ I a + β [ x ( t ) x ^ ( t ) ] + ( t ρ a ρ ) β ( η ρ t ρ ) ρ β + 1 Γ ( β + 2 ) Ω { ρ I a + α + β [ v 2 ( s ) v 1 ( s ) ] ( T ) λ ρ I a + β [ x ( T ) x ^ ( T ) ] μ ρ I a + α + β + γ [ v 2 ( s ) v 1 ( s ) ] ( ξ ) + μ λ ρ I a + β + γ [ x ( ξ ) x ^ ( ξ ) ] } ( t ρ a ρ ) β Ω ( η ρ a ρ ) β ( ( T ρ a ρ ) β ( T ρ t ρ ) ρ β + 1 Γ ( β + 2 ) μ ( ξ ρ a ρ ) β + γ [ ( β + 1 ) ( ξ ρ t ρ ) γ ( t ρ a ρ ) ] ρ β + γ + 1 Γ ( β + γ + 2 ) ( β + 1 ) ) ρ I a + α + β [ v 2 ( s ) v 1 ( s ) ] ( η ) λ ρ I a + β [ x ( η ) x ^ ( η ) ] | ϖ x x ^ ( ( T ρ a ρ ) α + β ρ α + β Γ ( α + β + 1 ) 1 + ζ 1 ρ β + 1 Γ ( β + 2 ) | Ω | + | μ | ( ξ ρ a ρ ) α + β + γ ζ 1 ρ α + 2 β + γ + 1 Γ ( α + β + γ + 1 ) Γ ( β + 2 ) | Ω | + ( η ρ a ρ ) α ζ 2 ρ α + β Γ ( α + β + 1 ) | Ω | ) + x x ^ ( | λ | ( T ρ a ρ ) β ρ β Γ ( β + 1 ) 1 + ζ 1 ρ β + 1 Γ ( β + 2 ) | Ω | + | μ | | λ | ( ξ ρ a ρ ) β + γ ζ 1 ρ 2 β + γ + 1 Γ ( β + γ + 1 ) Γ ( β + 2 ) | Ω | + | λ | ζ 2 ρ β Γ ( β + 1 ) | Ω | ) = ( ϖ Λ 1 + Λ 2 ) x x ^ .
Hence
h 1 h 2 ( ϖ Λ 1 + Λ 2 ) x x ^ .
Analogously, we can interchange the roles of x and x ^ to get
H d ( N ( x ) , N ( x ^ ) ) ( ϖ Λ 1 + Λ 2 ) x x ^ ,
which implies that N is a contraction by the condition (18). Hence, by the conclusion of Covitz and Nadler fixed point theorem [22], N has a fixed point x, which corresponds to a solution of (1). This finishes the proof. □

4. Examples

We illustrate our main results by presenting a numerical example.
Example 1.
Consider the following problem
c 1 / 3 D 5 / 4 c 1 / 3 D 1 / 4 + 1 / 5 x ( t ) F ( t , x ( t ) ) , t J : = [ 1 , 2 ] , x ( 1 ) = 0 , x ( 3 / 2 ) = 0 , x ( 2 ) = 2 / 7 1 / 3 I 3 / 4 x ( 7 / 4 ) .
Here ρ = 1 / 3 , α = 5 / 4 , β = 1 / 4 , λ = 1 / 5 , a = 1 , T = 2 η = 3 / 2 , μ = 2 / 7 , γ = 3 / 4 , ξ = 7 / 4 . Using the given data, we find that ζ 1 0.082260 , ζ 2 0.232036 , | Ω | 0.293634 , Λ 1 1.336009 and Λ 2 0.673563 , where ζ 1 , ζ 2 , Λ 1 and Λ 2 are given by (16), (17), (14) and (15) respectively.
(i) Let us consider the function
F ( t , x ( t ) ) = 1 t 2 + 63 | x ( t ) | 3 | x ( t ) | | x ( t ) | + 1 + 2 + 1 , e t 9 t + 8 sin x ( t ) + 1 80 .
We note that | F ( t , x ( t ) ) | P ( t ) Q ( x ) , where P ( t ) = 1 t 2 + 63 , Q ( x ) = x + 1 . So the assumption ( A 2 ) holds. Moreover, there exists M > 1.047447394 satisfying ( A 3 ) . Thus the hypothesis of Theorem 1 holds true and hence there exists at least one solution for the problem (19) with F ( t , x ) given by (20) on [ 1 , 2 ] .
(ii) To illustrate Theorem 2 we consider the function
F ( t , x ) = e t 19 + t , 1 ( t + 4 ) 2 x + tan 1 ( x ) + 1 15 .
Clearly H d ( F ( t , x ) , F ( t , x ¯ ) ) ϖ ( t ) | x x ¯ | , where ϖ ( t ) = 2 ( t + 4 ) 2 . Also d ( 0 , F ( t , 0 ) ) ϖ ( t ) for almost all t [ 0 , 1 ] and Λ 1 ϖ + Λ 2 0.7804433180 < 1 . As the hypothesis of Theorem 2 is satisfied, therefore we conclude that the multivalued problem (19) with F ( t , x ) given by (21) has at least one solution on [ 1 , 2 ] .

5. Conclusions

We have introduced a new class of multivalued (inclusions) boundary value problems on an arbitrary domain containing Caputo-type generalized fractional differential operators of different orders and a generalized integral operator. We have considered convex as well as non-convex valued cases for the multi-valued map involved in the given problem. Leray–Schauder nonlinear alternative for multivalued maps plays a central role in proving the existence of solutions for convex valued case of the given problem, while the existence result for the non-convex valued case is based on Covitz and Nadler fixed point theorem. The work presented in this paper is not only new in the given configuration, but will also lead to some new results as special cases. For example, fixing μ = 0 in the obtained results, we obtain the ones for nonlocal three-point boundary conditions: x ( a ) = 0 , x ( η ) = 0 , x ( T ) = 0 , 0 < η < T . For ρ = 1 , our results specialize to the ones for Liouville–Caputo type fractional differential inclusions complemented with nonlocal generalized integral boundary conditions on an arbitrary domain.

Author Contributions

Conceptualization, B.A. and M.A.; methodology, A.A. and S.K.N.; validation, A.A., B.A., M.A. and S.K.N.; formal analysis, A.A., B.A., M.A. and S.K.N.; writing—original draft preparation, M.A.; writing—review and editing, A.A., B.A. and S.K.N.; project administration, B.A.; funding acquisition, A.A.

Funding

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

Acknowledgments

This project was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under grant no. (KEP-PhD-70-130-38). The authors, therefore, acknowledge with thanks DSR technical and financial support.

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. Magin, R.L. Fractional Calculus in Bioengineering; Begell House Publishers: Danbury, CT, USA, 2006. [Google Scholar]
  3. 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]
  4. 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]
  5. Javidi, M.; Ahmad, B. Dynamic analysis of time fractional order phytoplankton-toxic phytoplankton– zooplankton system. Ecol. Model. 2015, 318, 8–18. [Google Scholar] [CrossRef]
  6. 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]
  7. Ahmad, B.; Alsaedi, A.; Ntouyas, S.K.; Tariboon, J. Hadamard-Type Fractional Differential Equations, Inclusions and Inequalities; Springer: Cham, Switzerland, 2017. [Google Scholar]
  8. Ahmad, B.; Nieto, J.J.; Alsaedi, A.; El-Shahed, M. A study of nonlinear Langevin equation involving two fractional orders in different intervals. Nonlinear Anal. Real World Appl. 2012, 13, 599–606. [Google Scholar] [CrossRef]
  9. Wang, G.; Zhang, L.; Song, G. Boundary value problem of a nonlinear Langevin equation with two different fractional orders and impulses. Fixed Point Theory Appl. 2012, 2012, 200. [Google Scholar] [CrossRef] [Green Version]
  10. Ahmad, B.; Ntouyas, S.K. New existence results for differential inclusions involving Langevin equation with two indices. J. Nonlinear Convex Anal. 2013, 14, 437–450. [Google Scholar]
  11. Muensawat, T.; Ntouyas, S.K.; Tariboon, J. Systems of generalized Sturm-Liouville and Langevin fractional differential equations. Adv. Differ. Equ. 2017, 2017, 63. [Google Scholar] [CrossRef] [Green Version]
  12. Fazli, H.; Nieto, J.J. Fractional Langevin equation with anti-periodic boundary conditions. Chaos Solitons Fractals 2018, 114, 332–337. [Google Scholar] [CrossRef]
  13. Ahmad, B.; Alsaedi, A.; Salem, S. On a nonlocal integral boundary value problem of nonlinear Langevin equation with different fractional orders. Adv. Differ. Equ. 2019, 2019, 57. [Google Scholar] [CrossRef] [Green Version]
  14. Ahmad, B.; Alghanmi, M.; Alsaedi, A.; Srivastava, H.; Ntouyas, K. The Langevin equation in terms of generalized Liouville-Caputo derivatives with nonlocal boundary conditions involving a generalized fractional integral. Mathematics 2019, 7, 533. [Google Scholar] [CrossRef]
  15. Katugampola, U.N. New Approach to a generalized fractional integral. Appl. Math. Comput. 2015, 218, 860–865. [Google Scholar] [CrossRef]
  16. Katugampola, U.N. A new approach to generalized fractional derivatives. Bull. Math. Anal. Appl. 2014, 6, 1–15. [Google Scholar]
  17. Jarad, F.; Abdeljawad, T.; Baleanu, D. On the generalized fractional derivatives and their caputo modification. J. Nonlinear Sci. Appl. 2017, 10, 2607–2619. [Google Scholar] [CrossRef]
  18. Granas, A.; Dugundji, J. Fixed Point Theory; Springer: New York, NY, USA, 2005. [Google Scholar]
  19. Deimling, K. Multivalued Differential Equations; De Gruyter: Berlin, Germany, 1992. [Google Scholar]
  20. 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. Astronom. Phys. 1965, 13, 781–786. [Google Scholar]
  21. Kisielewicz, M. Differential Inclusions and Optimal Control; Kluwer: Dordrecht, The Netherlands, 1991. [Google Scholar]
  22. Covitz, H.; Nadler, S.B., Jr. Multivalued contraction mappings in generalized metric spaces. ISRAEL J. Math. 1970, 8, 5–11. [Google Scholar] [CrossRef]
  23. 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]

Share and Cite

MDPI and ACS Style

Alsaedi, A.; Ahmad, B.; Alghanmi, M.; Ntouyas, S.K. On a Generalized Langevin Type Nonlocal Fractional Integral Multivalued Problem. Mathematics 2019, 7, 1015. https://doi.org/10.3390/math7111015

AMA Style

Alsaedi A, Ahmad B, Alghanmi M, Ntouyas SK. On a Generalized Langevin Type Nonlocal Fractional Integral Multivalued Problem. Mathematics. 2019; 7(11):1015. https://doi.org/10.3390/math7111015

Chicago/Turabian Style

Alsaedi, Ahmed, Bashir Ahmad, Madeaha Alghanmi, and Sotiris K. Ntouyas. 2019. "On a Generalized Langevin Type Nonlocal Fractional Integral Multivalued Problem" Mathematics 7, no. 11: 1015. https://doi.org/10.3390/math7111015

APA Style

Alsaedi, A., Ahmad, B., Alghanmi, M., & Ntouyas, S. K. (2019). On a Generalized Langevin Type Nonlocal Fractional Integral Multivalued Problem. Mathematics, 7(11), 1015. https://doi.org/10.3390/math7111015

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