Ann.
Acad.
Rom.
Sci.
Ser.
Math.
Appl.
ISSN
2066-6594
Vol.
18,
No.
3/2026
A
NOTE
ON
THE
EXISTENCE
OF
SOLUTIONS
OF
A
FRACTIONAL
DIFFERENTIAL
INCLUSION
OF
HILFER-HADAMARD
TYPE
∗
Aurelian
Cernea
†
Dedicated
to
the
memory
of
Professor
Mihail
Megan
DOI
10.56082/annalsarscimath.2026.3.29
Abstract
A
fractional
differential
inclusion
with
multi-point
boundary
con-
ditions
involving
Hilfer-Hadamard
derivatives
is
considered.
An
exis-
tence
result
is
established
when
the
set-valued
map
is
lower
semicon-
tinuous.
Keywords:
differential
inclusion,
fractional
derivative,
boundary
value
problem.
MSC:
34A60,
34A12,
34A08.
1
Introduction
In
this
note
we
consider
the
following
differential
inclusion
D
α,β
HH
x
(
t
)
∈
F
(
t,
x
(
t
))
a.e.
([1
,
e
])
(1)
with
boundary
conditions
of
the
form
x
(1
+
ε
)
=
0
,
D
1
,
1
HH
x
(
e
)
=
n
i
=1
ν
i
D
1
,
1
HH
x
(
η
i
)
,
(2)
∗
Accepted
for
publication
on
February
05,
2026
†
acernea@fmi.unibuc.ro
,
Faculty
of
Mathematics
and
Computer
Science,
University
of
Bucharest,
Academiei
14,
010014
Bucharest
and
Academy
of
Romanian
Scientists,
Ilfov
3,
050044
Bucharest,
Romania
29
A
Hilfer-Hadamard
fractional
integro-differential
inclusion
30
where
D
α,β
HH
is
the
Hilfer-Hadamard
fractional
derivative
of
order
α
∈
(1
,
2]
and
type
β
∈
[0
,
1],
η
i
∈
(1
,
e
),
ν
i
∈
R
,
i
=
1
,
n
,
ε
∈
(0
,
1)
and
F
:
[1
,
e
]
×
R
→
P
(
R
)
is
a
set-valued
map.
We
recall
that
Hilfer
fractional
derivative
introduced
in
[
14
]
is
an
interpo-
lation
between
Riemann-Liouville
and
Caputo
derivatives.
The
modification
of
Hilfer
fractional
derivative
resulted
in
the
concept
of
Hilfer-Hadamard
fractional
derivative.
Similarly,
the
Hilfer-Hadamard
fractional
derivative
covers
the
cases
of
Riemann-Liouville-Hadamard
and
Caputo-Hadamard
fractional
derivatives.
Several
recent
papers
are
devoted
to
the
study
to
boundary
value
prob-
lems
associated
to
fractional
differential
equations
and
inclusions
defined
by
Hilfer-Hadamard
fractional
derivative
[
1
,
9
,
10
,
18
]
etc.
For
background
facts
concerning
fractional
calculus
and
differential
equations
of
fractional
order
we
refer
to
[
2
,
12
,
15
,
17
]
etc.
The
present
study
is
motivated
by
a
recent
paper
[
3
],
where
two
existence
results
for
problem
(1)-(2)
are
provided.
A
result
is
obtained
in
the
case
when
the
values
of
the
set-valued
map
are
convex
and
the
set-valued
map
is
upper
semicontinuous.
The
other
result
is
proved
in
the
situation
when
the
set-valued
map
is
Lipschitz
in
the
state
variable
with
non
convex
values.
Our
aim
is
to
complete
the
study
in
[
3
]
by
taking
into
account
the
case
when
the
set-valued
map
has
non
convex
values,
but
instead
the
Lipschitz
hypothesis
in
[
3
]
we
assume
that
F
(
.,
.
)
is
lower
semicontinuous.
The
proof
of
our
existence
result
relies
on
a
nonlinear
alternative
of
Leray-Schauder
type
and
on
Bressan-Colombo-Fryszkowski
selection
theorem
for
lower
semi-
continuous
set-valued
maps
with
decomposable
values.
We
emphasize
that
such
kind
of
approach
is
usual
in
the
study
of
differ-
ential
inclusions
(e.g.,
[
5
–
8
,
11
]),
however
its
presentation
in
the
framework
of
problem
(1)-(2)
is
new.
The
paper
is
organized
as
follows:
in
Section
2
we
recall
some
preliminary
facts
that
we
need
in
the
sequel,
in
Section
3
we
prove
our
result.
2
Preliminaries
Let
(
X,
d
)
be
a
metric
space
with
the
corresponding
norm
|
.
|
and
let
I
⊂
R
be
a
compact
interval.
Denote
by
L
(
I
)
the
σ
-algebra
of
all
Lebesgue
measurable
subsets
of
I
,
by
P
(
X
)
the
family
of
all
nonempty
subsets
of
X
and
by
B
(
X
)
the
family
of
all
Borel
subsets
of
X
.
If
A
⊂
I
then
χ
A
(
.
)
:
I
→
{
0
,
1
}
denotes
the
characteristic
function
of
A
.
For
any
subset
A
⊂
X
we
denote
by
A
the
closure
of
A
.
B
r
(0)
and
B
r
(0)
will
be
the
open
and
closed
balls
centered
at
A.
Cernea
31
the
origin
and
of
radius
r
We
recall
that
the
Pompeiu-Hausdorff
generalized
metric
of
the
closed
subsets
A,
B
⊂
X
is
defined
by
d
H
(
A,
B
)
=
max
{
d
∗
(
A,
B
)
,
d
∗
(
B,
A
)
}
,
d
∗
(
A,
B
)
=
sup
{
d
(
a,
B
);
a
∈
A
}
,
where
d
(
x,
B
)
=
inf
y
∈
B
d
(
x,
y
).
As
usual,
if
X
is
a
separable
Banach
space,
we
denote
by
C
(
I,
X
)
the
Banach
space
of
all
continuous
functions
x
(
.
)
:
I
→
X
endowed
with
the
norm
|
x
|
C
=
sup
t
∈
I
|
x
(
t
)
|
and
by
L
1
(
I,
X
)
the
Banach
space
of
all
(Bochner)
integrable
functions
x
(
.
)
:
I
→
X
endowed
with
the
norm
|
x
|
1
=
I
|
x
(
t
)
|
d
t
.
A
subset
D
⊂
L
1
(
I,
X
)
is
said
to
be
decomposable
if
for
any
u
(
·
)
,
v
(
·
)
∈
D
and
any
subset
A
∈
L
(
I
)
one
has
uχ
A
+
vχ
B
∈
D
,
where
B
=
I
\
A
.
Consider
H
:
X
→
P
(
X
)
a
set-valued
map.
H
(
.
)
is
said
to
be
totally
compact
if
H
(
X
)
is
a
compact
subset
of
X
.
H
(
.
)
is
said
to
be
upper
semi-
continuous
if
for
any
x
0
∈
X
,
H
(
x
0
)
is
a
nonempty
and
closed
subset
of
X
and
if
for
each
open
set
D
of
X
containing
H
(
x
0
)
there
exists
an
open
neighborhood
V
0
of
x
0
such
that
H
(
V
0
)
⊂
D
.
H
(
.
)
is
said
to
be
lower
semi-
continuous
if
for
any
open
subset
D
⊂
X
,
the
set
{
x
∈
X
;
H
(
x
)
∩
D
=
∅}
is
open.
H
(
.
)
is
called
completely
continuous
if
it
is
upper
semicontinuous
and
transfers
bounded
sets
into
relatively
compact
sets.
The
result
in
the
next
section
is
based
on
the
following
nonlinear
alter-
native
of
Leray-Schauder
type
(e.g,
[
18
]).
Theorem
1.
Let
D
be
an
open
subset
in
a
normed
linear
space
X
such
that
0
∈
D
and
let
H
:
D
→
P
(
X
)
be
a
completely
continuous
set-valued
map
with
compact
and
convex
values.
Then
either
i)
the
inclusion
x
∈
H
(
x
)
has
a
solution,
or
ii)
there
exists
x
∈
∂D
(the
boundary
of
D
)
such
that
λx
∈
H
(
x
)
for
some
λ
>
1
.
In
fact,
we
use
the
following
consequence.
Corollary
1.
Let
X
be
a
normed
linear
space
and
let
h
:
B
r
(0)
→
X
be
a
completely
continuous
single
valued
map.
Then
either
i)
the
equation
x
=
h
(
x
)
has
a
solution,
or
ii)
there
exists
x
∈
X
with
|
x
|
=
r
and
x
=
λh
(
x
)
for
some
λ
<
1
.
If
F
(
.,
.
)
:
I
×
X
→
P
(
X
)
is
a
set-valued
map
with
compact
values
we
define
S
F
:
C
(
I,
X
)
→
P
(
L
1
(
I,
X
))
by
S
F
(
x
)
:=
{
f
∈
L
1
(
I,
X
);
f
(
t
)
∈
F
(
t,
x
(
t
))
a.e.
(
I
)
}
.
We
say
that
F
(
.,
.
)
is
of
lower
semicontinuous
type
if
S
F
(
.
)
is
lower
semicon-
tinuous
with
nonempty
closed
and
decomposable
values.
We
need,
also,
the
next
result
(
[
4
]).
A
Hilfer-Hadamard
fractional
integro-differential
inclusion
32
Theorem
2.
Let
S
be
a
separable
metric
space
and
H
(
.
)
:
S
→
P
(
L
1
(
I,
X
))
be
a
lower
semicontinuous
set-valued
map
with
closed
and
decomposable
val-
ues.
Then
H
(
.
)
has
a
continuous
selection
(i.e.,
there
exists
a
continuous
mapping
h
(
.
)
:
S
→
L
1
(
I,
X
)
such
that
h
(
s
)
∈
H
(
s
)
∀
s
∈
S
).
Also,
we
recall
some
definitions
from
fractional
calculus.
In
what
follows
I
=
[1
,
e
].
Definition
1.
a)
The
Hadamard
fractional
integral
of
order
q
>
0
of
a
Lebesgue
integrable
function
f
:
[1
,
∞
)
→
R
is
defined
by
I
q
H
f
(
t
)
=
1
Γ(
q
)
t
1
ln
t
s
q
−
1
f
(
s
)
s
ds
provided
the
integral
exists
and
Γ(
.
)
is
the
(Euler’s)
Gamma
function
defined
by
Γ(
q
)
=
∞
0
t
q
−
1
e
−
t
dt
.
b)
The
Hadamard
fractional
derivative
of
order
q
>
0
of
a
function
f
:
[1
,
∞
)
→
R
is
defined
by
D
q
H
f
(
t
)
=
t
d
n
dt
n
(
I
n
−
q
H
f
)(
t
)
,
provided
the
integral
exists
and
n
=
[
q
]
+
1
,
[
q
]
is
the
integer
part
of
q
.
c)
The
Hilfer-Hadamard
fractional
derivative
of
order
α
∈
(
n
−
1
,
n
)
,
n
≥
2
and
type
β
∈
[0
,
1]
of
a
function
f
∈
L
1
(
I,
R
)
is
defined
by
D
α,β
HH
f
(
t
)
=
(
I
β
(
n
−
α
)
H
D
γ
H
f
)(
t
)
,
γ
=
α
+
nβ
−
αβ.
When
β
=
0
the
Hilfer-Hadamard
fractional
derivative
gives
Hadamard
fractional
derivative
and
when
β
=
1
the
Hilfer-Hadamard
fractional
deriva-
tive
gives
Caputo-Hadamard
fractional
derivative.
In
our
case
n
=
[
α
]
+
1
=
2
and
γ
=
α
+
(2
−
α
)
β
.
We
make,
also,
the
notations
∆
=
(
γ
−
2)
µ
2
[ln(1
+
ε
)]
γ
−
1
−
(
γ
−
1)
µ
1
[ln(1
+
ε
)]
γ
−
2
,
µ
1
=
1
−
n
i
=1
ν
i
(ln
η
i
)
γ
−
2
,
µ
2
=
1
−
n
i
=1
ν
i
(ln
η
i
)
γ
−
3
,
Ω
=
(
γ
−
2)
µ
2
ln(1
+
ε
)
−
(
γ
−
1)
µ
1
.
The
next
result
is
proved
in
[
3
].
Lemma
1.
Assume
Ω
=
0
and
let
h
(
.
)
∈
C
(
I,
R
)
.
A.
Cernea
33
Then,
the
solution
of
problem
D
α,β
HH
x
(
t
)
=
h
(
t
)
with
boundary
conditions
(2)
is
given
by
x
(
t
)
=
I
α
H
h
(
t
)
+
(ln
t
)
γ
−
2
((
γ
−
1)
µ
1
−
(
γ
−
2)
µ
2
ln
t
)
∆
I
α
H
h
(1
+
ε
)+
[ln(1
+
ε
)]
γ
−
2
ln(1
+
ε
)(ln
t
)
γ
−
2
−
(ln
t
)
γ
−
1
∆
[
n
i
=1
ν
i
I
α
H
h
(
η
i
)
−
I
α
−
1
H
h
(
e
)]
.
Definition
2.
x
(
.
)
∈
C
(
I,
R
)
is
called
a
solution
of
a
problem
(1)-(2)
if
there
exists
h
(
.
)
∈
L
1
(
I,
R
)
such
that
h
(
t
)
∈
F
(
t,
x
(
t
))
a.e.
(
I
)
and
x
(
.
)
is
given
as
in
Lemma
1.
Remark
1.
If
we
denote
G
1
(
t,
τ
)
=
(ln
t
−
ln
τ
)
α
−
1
Γ(
α
)
τ
χ
[1
,t
]
(
τ
)
,
G
2
(
t,
τ
)
=
(ln
t
)
γ
−
2
[(
γ
−
1)
µ
1
−
(
γ
−
2)
µ
2
ln
t
]
∆
(ln
t
−
ln
τ
)
α
−
1
Γ(
α
)
τ
χ
[1
,
1+
ε
]
(
τ
)
,
G
3
(
t,
τ
)
=
n
i
=1
ν
i
χ
[1
,η
i
]
(
τ
)
(ln
t
−
ln
τ
)
α
−
2
Γ(
α
−
1)
τ
[ln(1
+
ε
)]
γ
−
2
ln(1+
ε
)(ln
t
)
γ
−
2
−
(ln
t
)
γ
−
1
∆
,
G
4
(
t,
τ
)
=
−
(ln
t
−
ln
τ
)
α
−
2
Γ(
α
−
1)
τ
[ln(1
+
ε
)]
γ
−
2
ln(1+
ε
)(ln
t
)
γ
−
2
−
(ln
t
)
γ
−
1
∆
,
G
(
t,
τ
)
=
G
1
(
t,
τ
)
+
G
2
(
t,
τ
)
+
G
3
(
t,
τ
)
+
G
4
(
t,
τ
)
then
the
solution
x
(
.
)
in
Lemma
1
may
be
put
as
x
(
t
)
=
e
1
G
(
t,
s
)
h
(
s
)
ds
.
Moreover,
|G
(
t,
s
)
|
≤
(ln
t
)
α
−
1
Γ(
α
)
+
|
(
γ
−
1)
µ
1
−
(
γ
−
2)
µ
2
ln
t
|
|
∆
|
Γ(
α
)
(ln
t
)
α
+
γ
−
3
+
(1
+
n
i
=1
|
ν
i
|
)
·
[ln(1+
ε
)]
γ
−
2
|
∆
|
Γ(
α
−
1)
|
ln(1
+
ε
)(ln
t
)
α
+
γ
−
2
−
(ln
t
)
α
+
γ
−
3
|
=:
M
(
t
)
∀
t,
τ
∈
I.
Denote
M
:=
max
t
∈
I
M
(
t
).
3
The
result
We
are
working
under
the
following
hypothesis.
Hypothesis
H
i)
F
(
.,
.
)
:
I
×
R
→
P
(
R
)
has
compact
values,
F
(
.,
.
)
is
L
(
I
)
⊗B
(
R
)
measurable
and
x
→
F
(
t,
x
)
is
lower
semicontinuous
for
almost
all
t
∈
I
.
ii)
There
exists
a
(
.
)
∈
L
1
(
I,
R
)
with
a
(
t
)
>
0
a.e.
(
I
)
and
there
exists
a
nondecreasing
function
b
:
[0
,
∞
)
→
(0
,
∞
)
such
that
sup
{|
v
|
;
v
∈
F
(
t,
x
)
}
≤
a
(
t
)
b
(
|
x
|
)
a.e.
(
I
)
,
∀
x
∈
R
.
A
Hilfer-Hadamard
fractional
integro-differential
inclusion
34
Theorem
3.
Assume
that
Ω
=
0
,
Hypothesis
H
is
satisfied
and
there
exists
c
>
0
such
that
c
>
M
|
a
|
1
b
(
c
)
.
(3)
Then
problem
(1)-(2)
has
at
least
one
solution
on
I
.
Proof.
From
Hypothesis
H
it
follows
that
F
(
.,
.
)
is
of
lower
semicontinu-
ous
type
(e.g.,
[
13
]).
Therefore,
by
Theorem
2
applied
with
S
=
C
(
I,
R
)
and
G
(
.
)
=
S
F
(
.
),
we
infer
that
there
exists
a
continuous
mapping
g
(
.
)
:
C
(
I,
R
)
→
L
1
(
I,
R
)
such
that
g
(
x
)
∈
S
F
(
x
)
for
all
x
∈
C
(
I,
R
).
We
consider
the
integral
equation
x
(
t
)
=
e
1
G
(
t,
s
)
g
(
x
(
s
))
ds,
t
∈
I
(4)
in
the
space
X
=
C
(
I,
R
).
If
x
(
.
)
∈
C
(
I,
R
)
is
a
solution
of
the
problem
(4)
then
x
(
.
)
is
a
solution
to
problem
(1)-(2).
Let
c
>
0
that
satisfies
condition
(3)
and
define
the
map
h
:
B
r
(0)
→
C
(
I,
R
)
by
(
h
(
x
))(
t
)
:=
e
1
G
(
t,
s
)
g
(
x
(
s
))
ds.
The
integral
equation
(4)
becomes
the
operatorial
equation
x
(
t
)
=
(
h
(
x
))(
t
)
,
t
∈
I.
It
remains
to
prove
that
h
(
.
)
satisfies
the
hypotheses
of
Corollary
1.
First,
we
show
that
h
(
.
)
is
continuous
on
B
r
(0).
Taking
into
account
Hypotheses
H
ii)
one
has
|
g
(
x
(
t
))
|
≤
a
(
t
)
b
(
|
x
(
t
)
|
)
a.e.
(
I
)
for
all
x
(
.
)
∈
C
(
I,
R
).
Consider
y
n
,
y
∈
B
r
(0)
with
y
n
→
y
as
n
→
∞
.
Then
|
g
(
y
n
(
t
))
|
≤
a
(
t
)
b
(
c
)
a.e.
(
I
)
.
Lebesgue’s
dominated
convergence
theorem
and
the
continuity
of
g
(
.
)
allow
to
write,
for
all
t
∈
I
lim
n
→∞
(
h
(
y
n
))(
t
)
=
e
1
G
(
t,
s
)
g
(
y
n
(
s
))
ds
=
e
1
G
(
t,
s
)
g
(
y
(
s
))
ds
=
(
h
(
y
))(
t
)
which
means
that
h
(
.
)
is
continuous
on
B
r
(0).
Secondly,
we
prove
that
h
(
B
r
(0))
is
compact.
Take
g
∈
S
F
(
x
)
and
set
w
(
t
)
=
e
1
G
(
t,
s
)
g
(
s
)
ds
.
For
any
t
∈
I
,
we
have
|
w
(
t
)
|
≤
e
1
|G
(
t,
s
)
|
.
|
g
(
s
)
|
ds
≤
e
1
|G
(
t,
s
)
|
a
(
s
)
b
(
|
x
(
t
)
|
)
ds
A.
Cernea
35
which
means
that
|
w
|
C
≤
M
|
a
|
1
b
(
r
)
,
i.e.,
h
(
B
r
(0))
is
bounded.
Next,
we
show
that
h
(
.
)
maps
bounded
sets
into
equi-continuous
sets.
Therefore,
let
B
⊂
B
r
(0)
be
a
bounded
set.
As
before,
take
g
∈
S
F
(
x
)
such
that
w
(
t
)
=
e
1
G
(
t,
s
)
g
(
s
)
ds
.
For
every
τ
1
,
τ
2
∈
I
we
have
|
w
(
τ
1
)
−
w
(
τ
2
)
|
≤
|
e
1
G
(
τ
1
,
s
)
g
(
s
)
ds
−
e
1
G
(
τ
2
,
s
)
g
(
s
)
ds
|
≤
e
1
|G
(
τ
1
,
s
)
−
G
(
τ
2
,
s
)
|
.
|
g
(
s
)
|
ds
≤
e
1
|G
(
τ
1
,
s
)
−
G
(
τ
2
,
s
)
|
a
(
s
)
b
(
r
)
ds,
which
means
that
|
w
(
τ
1
)
−
w
(
τ
2
)
|
→
0
as
τ
1
→
τ
2
.
Therefore,
h
(
B
r
(0))
is
an
equi-continuous
set
in
C
(
I,
R
).
Arzela-Ascoli’s
theorem
allows
us
to
deduce
that
h
(
.
)
is
completely
continuous
on
C
(
I,
R
).
So,
one
may
apply
Corollary
1
in
order
to
deduce
that
either
i)
the
equation
x
=
h
(
x
)
has
a
solution
in
B
r
(0),
or
ii)
there
exists
x
∈
X
with
|
x
|
C
=
c
and
x
=
λh
(
x
)
for
some
λ
<
1.
Assume
that
ii)
is
true.
Following
the
same
arguments
as
in
the
second
step
of
our
proof,
we
get
c
=
|
x
|
C
≤
M
|
a
|
1
b
(
c
).
This
is
a
contradiction
with
(3).
Thus,
only
i)
is
valid,
which
means
that
problem
(1)
has
a
solution
with
|
x
|
C
<
c
.
References
[1]
B.
Ahmad
and
S.K.
Ntouyas,
Hilfer-Hadamard
fractional
boundary
value
problems
with
nonlocal
mixed
boundary
conditions,
Fractal
Fract.
5
(2012),
195.
[2]
D.
B˘
aleanu,
K.
Diethelm,
E.
Scalas
and
J.J.
Trujillo,
Fractional
Calcu-
lus
Models
and
Numerical
Methods
,
World
Scientific,
Singapore,
2012.
[3]
F.Z.
Berrabah,
B.
Hedia
and
G.
Mophou,
Fractional
differential
inclu-
sions
with
multipoint
boundary
conditions
involving
Hilfer-Hadamard
derivative,
Dyn.
Contin.
Discrete
Impulsive
Syst.
Ser.
A
31
(2024),
321-348.
[4]
A.
Bressan
and
G.
Colombo,
Extensions
and
selections
of
maps
with
decomposable
values,
Studia
Math.
90
(1988),
69-86.
[5]
A.
Cernea,
On
a
Sturm-Liouville
type
differential
inclusion
of
fractional
order,
Fract.
Differ.
Calc.
7
(2017),
385-393.
[6]
A.
Cernea,
Existence
of
solutions
for
a
certain
boundary
value
problem
associated
to
a
fourth
order
differential
inclusion,
Int.
J.
Anal.
Appl.
14
(2017),
27-33.
A
Hilfer-Hadamard
fractional
integro-differential
inclusion
36
[7]
A.
Cernea,
A
bilocal
problem
associated
to
a
fractional
differential
in-
clusion
of
Caputo-Fabrizio
type,
Univ.
J.
Math.
Appl.
3
(2020),
133-137.
[8]
A.
Cernea,
On
the
solutions
of
a
certain
boundary
value
problem
as-
sociated
to
a
fractional
differential
inclusion,
Discuss.
Math.,
Differ.
Incl.,
Control
Optim.
41
(2021),
5-17.
[9]
A.
Cernea,
A
note
on
a
Hilfer-Hadamard
fractional
integro-differential
inclusion,
Proc.
Int.
E-Conf.
Math.
Stat.
Sc.
A
Sel¸
cuk
Meeting
,
Ed.
T.
Acar,
2022,
42-47.
[10]
A.
Cernea,
On
the
existence
of
solutions
for
a
class
of
Hilfer-Hadamard
fractional
integro-differential
inclusions,
Ann.
Acad.
Rom.
Sci.,
Ser.
Math.
Appl.
17
(2025),
29-38.
[11]
A.
Cernea,
A
note
on
solutions
of
a
fractional
differential
inclusion
with
certain
nonlocal
boundary
conditions,
Fixed
Point
Theory
,
27
(2026),
197-208.
[12]
K.
Diethelm,
The
Analysis
of
Fractional
Differential
Equations
,
Springer,
Berlin,
2016.
[13]
M.
Frignon,
A.
Granas,
Th´
eor`
emes
d’existence
pour
les
inclusions
diff´
erentielles
sans
convexit´
e,
C.
R.
Acad.
Sci.
Paris
310
(1990),
819-
822.
[14]
R.
Hilfer,
Applications
of
Fractional
Calculus
in
Physics
,
World
Scien-
tific,
Singapore,
2000.
[15]
A.
Kilbas,
H.M.
Srivastavaand
J.J.
Trujillo,
Theory
and
Applications
of
Fractional
Differential
Equations
,
Elsevier,
Amsterdam,
2006.
[16]
D.
O’
Regan,
Fixed
point
theory
for
closed
multifunctions,
Arch.
Math.
(Brno)
34
(1998),
191-197.
[17]
I.
Podlubny,
Fractional
Differential
Equations
,
Academic
Press,
San
Diego,
1999.
[18]
C.
Promsakon,
S.K.
Ntouyas
and
J.
Tariboon,
Hilfer-Hadamard
non-
local
integro-multipoint
fractional
boundary
value
problems,
J.
Funct.
Spaces
2021
(2021),
8031524.