# Extremal Solutions and Relaxation Problems for Fractional Differential Inclusions

Authors: Juan J. Nieto; Abdelghani Ouahab; P. Prakash
Journal: Abstract and Applied Analysis
 (2013)

---

## Abstract

We present the existence of extremal solution and relaxation
problem for fractional differential inclusion with initial conditions.


---

## Body

## 1. Introduction

Differential equations with fractional order have recently proved to be valuable tools in the modeling of many physical phenomena [1–9]. There has also been a significant theoretical development in fractional differential equations in recent years; see the monographs of Kilbas et al. [10], Miller and Ross [11], Podlubny [12], and Samko et al. [13] and the papers of Kilbas and Trujillo [14], Nahušev [15], Podlubny et al. [16], and Yu and Gao [17].Recently, some basic theory for initial value problems for fractional differential equations and inclusions involving the Riemann-Liouville differential operator was discussed, for example, by Lakshmikantham [18] and Chalco-Cano et al. [19].Applied problems requiring definitions of fractional derivatives are those that are physically interpretable for initial conditions containingy(0), y′(0), and so forth. The same requirements are true for boundary conditions. Caputo’s fractional derivative satisfies these demands. For more details on the geometric and physical interpretation for fractional derivatives of both Riemann-Liouville and Caputo types, see Podlubny [12].Fractional calculus has a long history. We refer the reader to [20].Recently fractional functional differential equations and inclusions and impulsive fractional differential equations and inclusions with standard Riemann-Liouville and Caputo derivatives with differences conditions were studied by Abbas et al. [21, 22], Benchohra et al. [23], Henderson and Ouahab [24, 25], Jiao and Zhou [26], and Ouahab [27–29] and in the references therein.In this paper, we will be concerned with the existence of solutions, Filippov’s theorem, and the relaxation theorem of abstract fractional differential inclusions. More precisely, we will consider the following problem:(1)cDαy(t)∈F(t,y(t)),a.e.t∈J≔[0,b],y(0)=y0,y′(0)=y1,(2)cDαy(t)∈extF(t,y(t)),a.e.t∈J≔[0,b],y(0)=y0,y′(0)=y1,
where cDα is the Caputo fractional derivatives, α∈(1,2], F:J×ℝN→𝒫(ℝN) is a multifunction, and extF(t,y) represents the set of extreme points of F(t,y). (𝒫(ℝN) is the family of all nonempty subsets of ℝN.During the last couple of years, the existence of extremal solutions and relaxation problem for ordinary differential inclusions was studied by many authors, for example, see [30–34] and the references therein.The paper is organized as follows. We first collect some background material and basic results from multivalued analysis and give some results on fractional calculus in Sections2 and 3, respectively. Then, we will be concerned with the existence of solution for extremal problem. This is the aim of Section 4. In Section 5, we prove the relaxation problem.


## 2. Preliminaries

The reader is assumed to be familiar with the theory of multivalued analysis and differential inclusions in Banach spaces, as presented in Aubin et al. [35, 36], Hu and Papageorgiou [37], Kisielewicz [38], and Tolstonogov [32].Let(X,∥·∥) be a real Banach space, [0,b] an interval in R, and C([0,b],X) the Banach space of all continuous functions from J into X with the norm
(3)∥y∥∞=sup{∥y(t)∥:0≤t≤b}.A measurable functiony:[0,b]→X is Bochner integrable if ∥y∥ is Lebesgue integrable. In what follows, L1([0,b],X) denotes the Banach space of functions y:[0,b]→X, which are Bochner integrable with norm
(4)∥y∥1=∫0b‍∥y(t)∥dt.
Denote by Lw1([0,b],X) the space of equivalence classes of Bochner integrable function y:[0,b]→X with the norm
(5)∥y∥w=supt∈[0,t]∥∫0t‍y(s)ds∥.
The norm ∥·∥w is weaker than the usual norm ∥·∥1, and for a broad class of subsets of L1([0,b],X), the topology defined by the weak norm coincides with the usual weak topology (see [37, Proposition 4.14, page 195]). Denote by(6)𝒫(X)={Y⊂X:Y≠∅},𝒫cl(X)={Y∈𝒫(X):Yclosed},𝒫b(X)={Y∈𝒫(X):Ybounded},𝒫cv(X)={Y∈𝒫(X):Yconvex},𝒫cp(X)={Y∈𝒫(X):Ycompact}.
A multivalued map G:X→𝒫(X) has convex (closed) values if G(x) is convex (closed) for all x∈X. We say that G is bounded on bounded sets if G(B) is bounded in X for each bounded set B of X (i.e., supx∈B{sup{∥y∥:y∈G(x)}}<∞).Definition 1.
A multifunctionF:X→𝒫(Y) is said to be upper semicontinuous at the point x0∈X, if, for every open W⊆Y such that F(x0)⊂W, there exists a neighborhood V(x0) of x0 such that F(V(x0))⊂W.A multifunction is calledupper semicontinuous (u.s.c. for short) on X if for each x∈X it is u.s.c. at x.Definition 2.
A multifunctionF:X→𝒫(Y) is said to be lower continuous at the point x0∈X, if, for every open W⊆Y such that F(x0)∩W≠∅, there exists a neighborhood V(x0) of x0 with property that F(x)∩W≠∅ for all x∈V(x0).A multifunction is calledlower semicontinuous (l.s.c. for short) provided that it is lower semicontinuous at every point x∈X.Lemma 3 (see [39, Lemma 3.2]).
LetF:[0,b]→𝒫(Y) be a measurable multivalued map and u:[a,b]→Y a measurable function. Then for any measurable v:[a,b]→(0,+∞), there exists a measurable selection fv of F such that for a.e. t∈[a,b],
(7)∥u(t)-fv(t)∥≤d(u(t),F(t))+v(t).First, consider the Hausdorff pseudometric(8)Hd:𝒫(E)×𝒫(E)⟶ℝ+∪{∞},
defined by
(9)Hd(A,B)=max{supa∈Ad(a,B),supb∈Bd(A,b)},
where d(A,b)=infa∈Ad(a,b) and d(a,B)=infb∈Bd(a,b). (𝒫b,cl(E),Hd) is a metric space and (𝒫cl(X),Hd) is a generalized metric space.Definition 4.
A multifunctionF:Y→𝒫(X) is called Hausdorff lower semicontinuous at the point y0∈Y, if for any ϵ>0 there exists a neighbourhood U(y0) of the point y0 such that
(10)F(y0)⊂F(y)+ϵB(0,1),foreveryy∈U(y0),
where B(0,1) is the unite ball in X.Definition 5.
A multifunctionF:Y→𝒫(X) is called Hausdorff upper semicontinuous at the point y0∈Y, if for any ϵ>0 there exists a neighbourhood U(y0) of the point y0 such that
(11)F(y)⊂F(y0)+ϵB(0,1),foreveryy∈U(y0).F is called continuous, if it is Hausdorff lower and upper semicontinuous.Definition 6.
LetX be a Banach space; a subset A⊂L1([0,b],X) is decomposable if, for all u,v∈A and for every Lebesgue measurable set I⊂J, one has
(12)uχI+vχ[0,b]∖I∈A,
where χA stands for the characteristic function of the set A. We denote  by  Dco(L1([0,b],X)) the family of decomposable sets.LetF:[0,b]×X→𝒫(X) be a multivalued map with nonempty closed values. Assign to F the multivalued operator ℱ:C([0,b],X)→𝒫(L1([0,b],X)) defined by
(13)ℱ(y)={v∈L1([0,b],X):v(t)∈F(t,y(t)),vvva.e.t∈[0,b]∈L1([0,b],X):v(t)∈F(t,y(t))}.
The operator ℱ is called the Nemyts’kiĭ operator associated to F.Definition 7.
LetF:[0,b]×X→𝒫(X) be a multivalued map with nonempty compact values. We say that F is of lower semicontinuous type (l.s.c. type) if its associated Nemyts’kiĭ operator ℱ is lower semicontinuous and has nonempty closed and decomposable values.Next, we state a classical selection theorem due to Bressan and Colombo.Lemma 8 (see [40]).
LetX be a separable metric space and let E be a Banach space. Then every l.s.c. multivalued operator N:X→𝒫cl(L1([0,b],E)) with closed decomposable values has a continuous selection; that is, there exists a continuous single-valued function f:X→L1([0,b],E) such that f(x)∈N(x) for every x∈X.Let us introduce the following hypothesis.(ℋ¯1)
F:[0,b]×X→𝒫(X) is a nonempty compact valued multivalued map such that
(a)
the mapping(t,y)↦F(t,y) is ℒ⊗ℬ measurable;(b)
the mappingy↦F(t,y) is lower semicontinuous for a.e. t∈[0,b].Lemma 9 (see, e.g., [41]).
LetF:J×X→𝒫cp(E) be an integrably bounded multivalued map satisfying (ℋ¯1). Then F is of lower semicontinuous type.Define(14)F(K)={f∈L1([0,b],X):f(t)∈Ka.e.t∈[0,b]},vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvK⊂X,
where X is a Banach space.Lemma 10 (see [37]).
LetK⊂X be a weakly compact subset of X. Then F(K) is relatively weakly compact subset of L1([0,b],X). Moreover if K is convex, then F(K) is weakly compact in L1([0,b],X).Definition 11.
A multifunctionF:[0,b]×Y→𝒫wcpcv(X) possesses the Scorza-Dragoni property (S-D property) if for each ϵ>0, there exists a closed set Jϵ⊂[0,b] whose Lebesgue measure μ(Jϵ)≤ϵ and such that F:[0,b]∖Jϵ×Y→X is continuous with respect to the metric dX(·,·).Remark 12.
It is well known that if the mapF:[0,b]×Y→𝒫wcpcv(X) is continuous with respect to y for almost every t∈[0,b] and is measurable with respect to t for every y∈Y, then it possesses the S-D property.In what follows, we present some definitions and properties of extreme points.Definition 13.
LetA be a nonempty subset of a real or complex linear vector space. An extreme point of a convex set A is a point x∈A with the property that if x=λy+(1-λ)z with y,z∈A and λ∈[0,1], then y=x and/or z=x.  ext(A) will denote the set of extreme points of A.In other words, an extreme point is a point that is not an interior point of any line segment lying entirely inA.Lemma 14 (see [42]).
A nonempty compact set in a locally convex linear topological space has extremal points.Let{xn′}n∈ℕ be a denumerable, dense (in σ(X′,X) topology) subset of the set {x∈X:∥x∥≤1}. For any A∈𝒫cpcv(X) and xn′ define the function
(15)dn(A,u)=max{〈y-z,xn′〉:y,z∈A,u=y+z2}.Lemma 15 (see [33]).
u
∈
ext
(
A
) if and only if dn(A,u)=0 for all n≥1.In accordance with Krein-Milman and Trojansky theorem [43], the set ext(SF) is nonempty and co¯(ext(SF))=SF.Lemma 16 (see [33]).
LetF:[0,b]→𝒫wcpcv(X) be a measurable, integrably bounded map. Then
(16)ext¯(SF)⊆SF,
where ext¯(SF) is the closure of set ext(SF) in the topology of the space L1([0,b],X).Theorem 17 (see [33]).
LetF:[0,b]×Y→𝒫wcpcv(X) be a multivalued map that has the S-D property and let it be integrable bounded on compacts from Y. Consider a compact subset K⊂C([0,b],X) and define the multivalued map G:K→L1([0,b],X), by
(17)G(y(·))={f∈L1([0,b],X):f(t)∈F(t,y(t))a.e.on[0,b]},vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvy∈K.
Then for every K compact in C([0,b],X),  ϵ>0 and any continuous selection f:K→L1([0,b],X), there exists a continuous selector g:K→L1([0,b],X) of the map ext(G):K→L1([0,b],X) such that for all y∈C([0,b],X) one has
(18)supt∈[0,b]∥∫0t‍((fy)(s)-(gy)(s))ds∥≤ϵ.For a background of extreme point ofF(t,y(t)) see Dunford-Schwartz [42, Chapter 5, Section 8] and Florenzano and Le Van [44, Chapter 3].


## 3. Fractional Calculus

According to the Riemann-Liouville approach to fractional calculus, the notation of fractional integral of orderα (α>0) is a natural consequence of the well known formula (usually attributed to Cauchy) that reduces the calculation of the n-fold primitive of a  function   f(t) to a single integral of convolution type. In our notation the Cauchy formula reads
(19)Inf(t):=1(n-1)!∫0t‍(t-s)n-1f(s)ds,t>0,n∈ℕ.Definition 18 (see [13, 45]).
The fractional integral of orderα>0 of a function f∈L1([a,b],ℝ) is defined by
(20)Ia+αf(t)=∫at‍(t-s)α-1Γ(α)f(s)ds,
where Γ is the gamma function. When a=0, we write Iαf(t)=f(t)*ϕα(t), where ϕα(t)=t(α-1)/Γ(α) for t>0, and we write ϕα(t)=0 for t≤0 and ϕα→δ(t) as α→0, where δ is the delta function and Γ is the Euler gamma function defined by
(21)Γ(α)=∫0∞‍tα-1e-tdt,α>0.
For consistency, I0=Id (identity operator), that is, I0f(t)=f(t). Furthermore, by Iαf(0+) we mean the limit (if it exists) of Iαf(t) for t→0+; this limit may be infinite.After the notion of fractional integral, that of fractional derivative of orderα(α>0) becomes a natural requirement and one is attempted to substitute α with -α in the above formulas. However, this generalization needs some care in order to guarantee the convergence of the integral and preserve the well known properties of the ordinary derivative of integer order. Denoting by Dn, with n∈ℕ, the operator of the derivative of order n, we first note that
(22)DnIn=Id,InDn≠Id,n∈ℕ,
that is, Dn is the left inverse (and not the right inverse) to the corresponding integral operator Jn. We can easily prove that
(23)InDnf(t)=f(t)-∑k=0n-1‍f(k)(a+)(t-a)kk!,t>0.
As a consequence, we expect that Dα is defined as the left inverse to Iα. For this purpose, introducing the positive integer n such that n-1<α≤n, one defines the fractional derivative of order α>0.Definition 19.
For a functionf given on interval [a,b], the αth Riemann-Liouville fractional-order derivative of f is defined by
(24)Dαf(t)=1Γ(n-α)(ddt)n∫at‍(t-s)-α+n-1f(s)ds,
where n=[α]+1 and [α] is the integer part of α.Defining for consistency,D0=I0=Id, then we easily recognize that
(25)DαIα=Id,α≥0,(26)α>0,Dαtγ=Γ(γ+1)Γ(γ+1-α)tγ-α,α>0,γ∈(-1,0)∪(0,+∞),t>0.
Of course, properties (25) and (26) are a natural generalization of those known when the order is a positive integer.Note the remarkable fact that the fractional derivativeDαf is not zero for the constant function f(t)=1, if α∉ℕ. In fact, (26) with γ=0 illustrates that
(27)Dα1=(t-a)-αΓ(1-α),α>0,t>0.
It is clear that Dα1=0, for α∈ℕ, due to the poles of the gamma function at the points 0,-1,-2,….We now observe an alternative definition of fractional derivative, originally introduced by Caputo [46, 47] in the late sixties and adopted by Caputo and Mainardi [48] in the framework of the theory of Linear Viscoelasticity (see a review in [4]).Definition 20.
Letf∈ACn([a,b]). The Caputo fractional-order derivative of f is defined by
(28)(Dcαf)(t):=1Γ(n-α)∫at‍(t-s)n-α-1fn(s)ds.This definition is of course more restrictive than Riemann-Liouville definition, in that it requires the absolute integrability of the derivative of orderm. Whenever we use the operator D*α we (tacitly) assume that this condition is met. We easily recognize that in general
(29)Dαf(t):=DmIm-αf(t)≠Jm-αDmf(t):=D*αf(t),
unless the function f(t), along with its first n-1 derivatives, vanishes at t=a+. In fact, assuming that the passage of the m-derivative under the integral is legitimate, we recognize that, for m-1<α<m and t>0,
(30)Dαf(t)=Dcαf(t)+∑k=0m-1‍(t-a)k-αΓ(k-α+1)f(k)(a+),
and therefore, recalling the fractional derivative of the power function (26), one has
(31)Dα(f(t)-∑k=0m-1‍(t-a)k-αΓ(k-α+1)f(k)(a+))=D*αf(t).
The alternative definition, that is, Definition 20, for the fractional derivative thus incorporates the initial values of the function and of lower order. The subtraction of the Taylor polynomial of degree n-1 at t=a+ from f(t) means a sort of regularization of the fractional derivative. In particular, according to this definition, the relevant property for which the fractional derivative of a constant is still zero:(32)cDα1=0,α>0.
We now explore the most relevant differences between the two fractional derivatives given in Definitions 19 and 20. From Riemann-Liouville fractional derivatives, we have
(33)Dα(t-a)α-j=0,forj=1,2,…,[α]+1.
From (32) and (33) we thus recognize the following statements about functions which, for t>0, admit the same fractional derivative of order α, with n-1<α≤n,  n∈ℕ:
(34)Dαf(t)=Dαg(t)⟺f(t)=g(t)+∑j=1m‍cj(t-a)α-j,cDαf(t)=Dcαg(t)⟺f(t)=g(t)+∑j=1m‍cj(t-a)n-j.
In these formulas, the coefficients cj are arbitrary constants. For proving all main results we present the following auxiliary lemmas.Lemma 21 (see [10]).
Letα>0 and let y∈L∞(a,b) or C([a,b]). Then
(35)(DcαIαy)(t)=y(t).Lemma 22 (see [10]).
Letα>0 and n=[α]+1. If y∈ACn[a,b] or y∈Cn[a,b], then
(36)(IαcDαy)(t)=y(t)-∑k=0n-1‍y(k)(a)k!(t-a)k.For further readings and details on fractional calculus, we refer to the books and papers by Kilbas [10], Podlubny [12], Samko [13], and Caputo [46–48].


## 4. Existence Result

Definition 23.
A functiony∈C([0,b],ℝN) is called mild solution of problem (1) if there exist f∈L1(J,ℝN) such that
(37)y(t)=y0+ty1+1Γ(α)∫0t‍(t-s)1-αf(s)ds,t∈[0,b],
where f∈SF,y={v∈L1([0,b],ℝN):f(t)∈F(t,y(t))a.e.on[0,b]}.We will impose the following conditions onF.(ℋ1)
The functionF:J×ℝN→𝒫cpcv(ℝN) such that
(a)
for allx∈ℝN,     the  map     t↦F(t,x) is measurable,(b)
for everyt∈[0,b], the multivalued map x→F(t,x) is Hd continuous(ℋ2)
There existp∈L1(J,ℝ+) and a continuous nondecreasing function ψ:[0,∞)→(0,∞) such that
(38)∥F(t,x)∥𝒫=sup{∥v∥:v∈F(t,x)}≤p(t)ψ(∥x∥),vvvvfora.e.t∈[0,b]andeachx∈ℝN,
with
(39)∫0b‍p(s)ds<∫∥y0∥+b∥y1∥∞‍duψ(u).Theorem 24.
Assume that the conditions(ℋ1)-(ℋ2) and then the problem (2) have at least one solution.Proof.
From(ℋ2) there exists M>0 such that ∥y∥∞≤M   for  each   y∈Sc.
Let(40)F1(t,y)={F(t,y)if∥y∥≤M,F(t,My∥y∥)if∥y∥≥M.
We consider
(41)Dcαy(t)∈F1(t,y(t)),a.e.t∈[0,b],y(0)=y0,y′(0)=y1.
It is clear that all the solutions of (41) are solutions of (2).
Set(42)V={f∈L1([0,b],ℝN):∥f(t)∥≤ψ*(t)},vvvvvvvvvvvvvvvvvvvψ*(t)=p(t)ψ(M).
It is clear that V is weakly compact in L1([0,b],ℝN). Remark that for every f∈V, there exists a unique solution L(f) of the following problem:
(43)Dcαy(t)=f(t),a.e.t∈[0,b],y(t)=y0,y′(0)=y1;
this solution is defined by
(44)L(f)(t)=y0+ty1+1Γ(α)∫0t‍(t-s)α-1f(s)ds,vvvvvvvvvvvvvvvvvvvva.e.t∈[0,b].
We claim that L is continuous. Indeed, let fn→f converge in L1([0,b],ℝN), as n→∞, set yn=L(fn), n∈ℕ. It is clear that {yn:n∈ℕ} is relatively compact in C([0,b],ℝN)   and yn converge to y∈C([0,b],ℝN). Let
(45)z(t)=y0+y1t+1Γ(α)∫0t‍(t-s)α-1f(s)ds,t∈[0,b].
Then
(46)∥yn-z∥∞≤bαΓ(α)∫0b‍∥fn(s)-f(s)∥ds⟶0,vvvvvvvvvvvvvvvvvvasn⟶∞.
Hence K=L(V) is compact and convex subset of C([0,b],ℝN). Let SF:K→𝒫clcv(L1([0,b],ℝN)) be the multivalued Nemitsky operator defined by
(47)SF1(y)={f∈L1([0,b],ℝN):f(t)∈F1(t,y(t)),vvva.e.t∈[0,b]([0,b],(t,y(t))ℝN)}:=SF1,y.
It is clear that F1(·,·) is Hd continuous and F1(·,·)∈𝒫wkcpcv(ℝN) and is integrably bounded, then by Theorem 17 (see also Theorem 6.5 in [32] or Theorem 1.1 in [34]), we can find a continuous function g:K→Lw1([0,b],ℝN) such that
(48)g(x)∈extSF1(y)∀y∈K.
From Benamara [49] we know that
(49)extSF1(y)=SextF1(·,y(·))∀y∈K.
Setting N=L∘g and letting y∈K, then
(50)g(y)∈F1(·,y(·))⟹g(y)∈V⟹N(y)=L(g(y))∈K.
Now, we prove that N is continuous. Indeed, let yn∈K converge to y in C([0,b],ℝN).
Then(51)g(yn)convergeweaklytog(y)asn⟶∞.
Since N(yn)=L(g(yn))∈K and g(yn)(·)∈F(t,yn(t)), then
(52)g(yn)(·)∈F(·,B¯M)∈𝒫cp(ℝN).
From Lemma 10, g(yn) converge weakly to y in L1([0,b],ℝN) as n→∞. By the definition of N, we have
(53)N(yn)=y0+y1t+1Γ(α)∫0t‍(t-s)α-1g(yn)(s)ds,vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvt∈[0,b],N(y)=y0+y1t+1Γ(α)∫0t‍(t-s)α-1g(y)(s)ds,vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvt∈[0,b].
Since {N(yn):n∈ℕ}⊂K, then there exists subsequence of N(yn) converge in C([0,b],ℝN). Then
(54)N(yn)(t)⟶N(y)(t),∀t∈[0,b],asn⟶∞.
This proves that N is continuous. Hence by Schauder’s fixed point there exists y∈K such that y=N(y).


## 5. The Relaxed Problem

In this section, we examine whether the solutions of the extremal problem are dense in those of the convexified one. Such a result is important in optimal control theory whether the relaxed optimal state can be approximated by original states; the relaxed problems are generally much simpler to build. For the problem for first-order differential inclusions, we refer, for example, to [35, Theorem 2, page 124] or [36, Theorem 10.4.4, page 402]. For the relaxation of extremal problems we see the following recent references [30, 50].Now we present our main result of this section.Theorem 25.
LetF:[0,b]×ℝN→𝒫(ℝN) be a multifunction satisfying the following hypotheses.(ℋ3)
The functionF:[0,b]×ℝN→𝒫cpcv(ℝN) such that, for all x∈ℝN, the map
(55)t⟼F(t,x)
is measurable.(ℋ4)
There existsp∈L1(J,ℝ+) such that
(56)Hd(F(t,x),F(t,y))≤p(t)∥x-y∥,fora.e.t∈[0,b]andeachx,y∈ℝN,Hd(F(t,0),0)≤p(t)fora.e.t∈[0,b].Then Se¯=Sc.Proof.
By Coviz and Nadlar fixed point theorem, we can easily prove thatSc≠∅, and since F has compact and convex valued, then Sc is compact in C([0,b],ℝN). For more information we see [25, 27–29, 51, 52].
Lety∈Sc; then there exists f∈SF,y such that
(57)y(t)=y0+y1t+1Γ(α)∫0t‍(t-s)α-1f(s)ds,vvvvvvvvvvvvvvvvvvva.e.t∈[0,b].
Let K be a compact and convex set in C([0,b],ℝN) such that Sc⊂K. Given that y*∈K and ϵ>0, we define the following multifunction Uϵ:[0,b]→𝒫(ℝN) by
(58)Uϵ(t)={u∈ℝN:∥f(t)-u∥<d(f(t),F(t,y(t)))+ϵ,===u∈F(t,y*(t))u∈ℝN:∥f(t)-u∥<d(f(t),F(t,y(t)))+ϵ,}.
The multivalued map t→F(t,·) is measurable and x→F(·,x) is Hd continuous. In addition, if F(·,·)  has   compact values, then F(·,·) is graph measurable, and the mapping t→F(t,y(t)) is a measurable multivalued map for fixed y∈C([0,b],ℝN). Then by Lemma 3, there exists a measurable selection v1(t)∈F(t,y(t))a.e.t∈[0,b] such that
(59)∥f(t)-v1(t)∥<d(f(t),F(t,y(t)))+ϵ;
this implies that Uϵ(·)≠∅. We consider Gϵ:K→𝒫(L1(J,ℝN) defined by
(60)Gϵ(y)={f*∈ℱ(y):∥f(t)-f*(t)∥vvv<ϵ+d(f*(t),F(t,y(t)))}.
Since the measurable multifunction F is integrable bounded, Lemma 9 implies that the Nemyts’kiĭ operator ℱ has decomposable values. Hence y→Gϵ(y)¯ is l.s.c. with decomposable values. By Lemma 8, there exists a continuous selection fϵ:C([0,b],ℝN)→L1(J,ℝN) such that
(61)fϵ(y)∈Gϵ(y)¯∀y∈C([0,b],ℝN).
From Theorem 17, there exists function gϵ:K→Lw([0,b],ℝN) such that
(62)gϵ(y)∈extSF(y)=SextF(·,y(·))∀y∈K,∥gϵ(y)-fϵ(y)∥w≤ϵ,∀y∈K.
From (ℋ3) we can prove that there exists M>0 such that
(63)∥y∥∞≤M∀y∈Sc.
Consider the sequence ϵn→0,  as  n→∞, and set gn=gϵn, fn=fϵn. Set
(64)V={f∈L1([0,b],ℝN):∥f(t)∥≤ψ(t)a.e.t∈[0,b]},vvvvvvvvvvvvvvvvvvvvvvvvvvψ(t)=(1+M)p(t).
Let L:V→C([0,b],ℝN) be the map such that each f∈V assigns the unique solution of the problem
(65)Dcαy(t)=f(t),a.e.t∈[0,b],y(0)=y0,y′(0)=y1.
As in Theorem 24, we can prove that L(V) is compact in C([0,b],ℝN) and the operator Nn=L∘gn:K→K is compact; then by Schauder’s fixed point there exists y~n∈K such that y~n∈Se and
(66)y~n(t)=y0+ty1+1Γ(α)∫0t‍(t-s)α-1gn(yn)(s)ds,vvvvvvvvvvvvvvvvvva.e.t∈[0,b],n∈ℕ.
Hence
(67)∥y(t)-y~n(t)∥≤1Γ(α)∥∫0t‍(t-s)α-1[gn(y~n)(s)-f(s)]ds∥≤1Γ(α)∥∫0t‍(t-s)α-1[gn(y~n)(s)-fn(y~n)(s)]ds∥+bαΓ(α)∫0t‍∥fn(y~n)(s)-f(s)∥ds≤1Γ(α)∥∫0t‍(t-s)α-1[gn(y~n)(s)-fn(y~n)(s)]ds∥+bαΓ(α)∫0t‍(ϵn+d(f(s),fn(y~n)(s)))ds≤1Γ(α)∥∫0t‍(t-s)α-1[gn(y~n)(s)-fn(y~n)(s)]ds∥+bαΓ(α)∫0t‍(ϵn+Hd(F(s,y(s)),F(s,y~n(s))))ds≤bα+1Γ(α+1)ϵn+bα+1Γ(α)ϵn+∫0t‍p(s)∥y(s)-y~n(s)∥.
Let y~(·) be a limit point of the sequence y~n(·). Then, it follows that from the above inequality, one has
(68)∥y(t)-y~(t)∥≤∫0t‍p(s)∥y(s)-y~(s)∥ds,
which implies y(·)=y~(·). Consequently, y∈Sc is a unique limit point of y~n(·)∈Se.Example 26.
LetF:J×ℝN→𝒫cpcv(ℝN) with
(69)F(t,y)=B-(f1(t,y),f2(t,y)),
where f1,f2:J×ℝN→ℝN are Carathéodory functions and bounded.
Then (2) is solvable.Example 27.
If, in addition to the conditions onF of Example 26, f1 and f2 are Lipschitz functions, then Se¯=Sc.



---

*Source: 292643-2013-09-17.xml*