We provide a detailed analysis of the shock formation process for the non-isentropic 2d Euler equations in azimuthal symmetry. We prove that from an open set of smooth and generic initial data, solutions of the Euler equations form a first singularity or gradient blow-up. This first singularity is termed a Hölder C13 pre-shock, and our analysis provides the first detailed description of this cusp solution. The novelty of this work relative to [
Citation: Isaac Neal, Steve Shkoller, Vlad Vicol. A characteristics approach to shock formation in 2D Euler with azimuthal symmetry and entropy[J]. Communications in Analysis and Mechanics, 2025, 17(1): 188-236. doi: 10.3934/cam.2025009
[1] | Thabet Abdeljawad, Pshtiwan Othman Mohammed, Hari Mohan Srivastava, Eman Al-Sarairah, Artion Kashuri, Kamsing Nonlaopon . Some novel existence and uniqueness results for the Hilfer fractional integro-differential equations with non-instantaneous impulsive multi-point boundary conditions and their application. AIMS Mathematics, 2023, 8(2): 3469-3483. doi: 10.3934/math.2023177 |
[2] | Weerawat Sudsutad, Chatthai Thaiprayoon, Sotiris K. Ntouyas . Existence and stability results for $ \psi $-Hilfer fractional integro-differential equation with mixed nonlocal boundary conditions. AIMS Mathematics, 2021, 6(4): 4119-4141. doi: 10.3934/math.2021244 |
[3] | Mohamed Adel, M. Elsaid Ramadan, Hijaz Ahmad, Thongchai Botmart . Sobolev-type nonlinear Hilfer fractional stochastic differential equations with noninstantaneous impulsive. AIMS Mathematics, 2022, 7(11): 20105-20125. doi: 10.3934/math.20221100 |
[4] | Muneerah AL Nuwairan, Ahmed Gamal Ibrahim . The weighted generalized Atangana-Baleanu fractional derivative in banach spaces- definition and applications. AIMS Mathematics, 2024, 9(12): 36293-36335. doi: 10.3934/math.20241722 |
[5] | Xiaoming Wang, Rizwan Rizwan, Jung Rey Lee, Akbar Zada, Syed Omar Shah . Existence, uniqueness and Ulam's stabilities for a class of implicit impulsive Langevin equation with Hilfer fractional derivatives. AIMS Mathematics, 2021, 6(5): 4915-4929. doi: 10.3934/math.2021288 |
[6] | J. Vanterler da C. Sousa, E. Capelas de Oliveira, F. G. Rodrigues . Ulam-Hyers stabilities of fractional functional differential equations. AIMS Mathematics, 2020, 5(2): 1346-1358. doi: 10.3934/math.2020092 |
[7] | Murugesan Manigandan, R. Meganathan, R. Sathiya Shanthi, Mohamed Rhaima . Existence and analysis of Hilfer-Hadamard fractional differential equations in RLC circuit models. AIMS Mathematics, 2024, 9(10): 28741-28764. doi: 10.3934/math.20241394 |
[8] | Thabet Abdeljawad, Sabri T. M. Thabet, Imed Kedim, Miguel Vivas-Cortez . On a new structure of multi-term Hilfer fractional impulsive neutral Levin-Nohel integrodifferential system with variable time delay. AIMS Mathematics, 2024, 9(3): 7372-7395. doi: 10.3934/math.2024357 |
[9] | Muneerah Al Nuwairan, Ahmed Gamal Ibrahim . Nonlocal impulsive differential equations and inclusions involving Atangana-Baleanu fractional derivative in infinite dimensional spaces. AIMS Mathematics, 2023, 8(5): 11752-11780. doi: 10.3934/math.2023595 |
[10] | Ugyen Samdrup Tshering, Ekkarath Thailert, Sotiris K. Ntouyas . Existence and stability results for a coupled system of Hilfer-Hadamard sequential fractional differential equations with multi-point fractional integral boundary conditions. AIMS Mathematics, 2024, 9(9): 25849-25878. doi: 10.3934/math.20241263 |
We provide a detailed analysis of the shock formation process for the non-isentropic 2d Euler equations in azimuthal symmetry. We prove that from an open set of smooth and generic initial data, solutions of the Euler equations form a first singularity or gradient blow-up. This first singularity is termed a Hölder C13 pre-shock, and our analysis provides the first detailed description of this cusp solution. The novelty of this work relative to [
A non-instantaneous impulsive differential equation is due to Hernándaz et. al.[1], and is used to describe impulsive action, which stays active on a finite time interval. Hilfer [2] introduced a fractional derivative, which is a generalization for Riemann-Liouville fractional derivative and Caputo fractional derivative. Many works have been appeared studying various models involving fractional differential with instantaneous and non-instantaneous impulses and providing solutions to those models. For example, Saravanakumar et al. [3] analyzed the existence of mild solution of non instantaneous impulsive for Hilfer fractional stochastic differential equations driven by fractional Brownian motion, . Shu et al.[4] presented a right formula of mild solutions to a fractional semilinear evolution equation generated by a sectorial operator, and its order belongs to the intervals (0,1) and (1,2), Wang et al. [5] studied the global attracting solutions to non-instantaneous impulsive differential inclusions containing Hilfer fractional, and Ngo et al.[6] presented a formula of solution for a non-instantaneous impulsive differential equation containing ψ−Hilfer derivative with lower limit of the fractional derivative at zero. For more works on non-instantaneous impulsive differential equations and inclusions, we refer to [7,8,9,10,11,12,13].
Moreover, Ulam problem [14] has been attracted by many researchers. We highlight some recent works on the existence and Hyers-Ulam stability of solutions for fractional differential equations. Guo et al.[15] investigated the existence and Hyers-Ulam stability of mild solution for an impulsive Riemann-Liouville fractional neutral functional stochastic differential equation with infinite delay of order between one and two, Guo et al.[16] proved the existence and Hyers-Ulam stability of the almost periodic solution to fractional differential equations with impulse involving fractional Brownian, Wang et al. [11] presented the generalized Ulam-Hyers stability for a non-instantaneous impulsive differential inclusions containing the Caputo derivative and Vanterler et al. [17] studied, in finite dimensional Banach spaces, the stability of a Volterra integro-differential equation containing ψ− Hilfer derivative in the sense of Ulam-Hyers. More recently, Vanterler et al.[18] investigated, in finite dimensional Banach spaces, the δ−Ulam-Hyers-Rassias stability for a non-instantaneous impulsive fractional differential equation containing ψ−Hilfer derivative, Benchohra et al. [19] established, in finite dimensional spaces, the existence and stability of solutions for an implicit fractional differential equations with Riemann-Liouville fractional derivative, and Kumar et al. [20] studied the existence and stability of solution for a fractional differential equation with non-instantaneous integral impulses. Very recently, Ben Mahlouf et al.[21] given sufficient conditions to guarantee the existence and stability of solutions for generalized nonlinear fractional differential equations of order α∈(1,2), Elsayed et al.[22] established the existence and stability of boundary value problem for differential equation with Hilfer-Katugampola fractional derivative. For more papers on Ulam-Hyers stability of solutions, we refer to [23,24,25,26,27,28,29,30,31].
It is worth noting that, when the considered problem contains non-instantaneous impulses, there are two approaches in the literature, one by keeping the lower limit of the fractional derivative at zero [6,11,17,18], and the other by switching it at the impulsive points [5,10] Motivated by the above cited work, we prove two existence results of solutions, in infinite dimensional Banach spaces, for a non-instantaneous impulsive fractional differential inclusion involving ψ−Hilfer derivative with delay and we switch the lower limit of the fractional derivative at the impulsive points, and then we study the ψ−generalized Ulam-Hyers stability.
Let E be a real Banach space, J=[0,b], b>0, J∗=(0,b], r>0, 0<ϑ<1, 0≤ν≤1, μ=ϑ+ν−ϑν, Ψ:[−r,0]→E a continuous function except a finite number of discontinuity points s≠0 such that Ψ(0)=0, all values x(s+), and x(s−) are finite, ψ∈C1([0,b],R) be increasing, ψ´(ϱ)>0,;∈J, and Dϑ,ν,ψsi+ be the ψ−Hilfer derivative with lower limit at si of orderϑ and type ν. Moreover, 0=s0< ϱ1<s1<ϱ2<⋯<ϱn<sn<ϱn+1=b, and I1−μ,ψs+ix(s+i)=limϱ→s+i I1−μ,ψs+ix(ϱ) , F:J×E→2E−{ϕ} is a multifunction, and gi:[ϱi,si]×E⟶Ei=1,2,…,n. Finally, for any ϱ∈J, τ(ϱ):H→Θ, x(θ)=x(ϱ+θ),θ∈[−r,0]x∈H, where Θ and H will be introduced in the next section. In this paper, we establish existence results of solutions of the following ψ−Hilfer fractional non-instantaneous impulsive differential inclusions with delay:
{Dϑ,ν,ψs+ix(ϱ)∈F(ϱ,τ(ϱ)x)),a.e.ϱ∈∪i=ni=0(si,ϱi+1]x(ϱ)=Ψ(ϱ),ϱ∈[−r,0]I1−μ,ψ0+x(0+)=Ψ(0),x(ϱ+i)=gi(ϱi,x(ϱ−i)),i=1,......n,x(ϱ)=gi(ϱ,x(ϱ−i)),ϱ∈(ϱisi],i=1,......n,I1−μ,ψs+ix(s+i)=gi(si,x(ϱ−i)),i=1,......n. | (1) |
Then, we investigate the ψ-generalized Ulam-Hyers stability of Problem (1). To achieve our aim, we present a relation between the solutions of this problem and the corresponding fractional integral equation (Lemma5).
To make a comparison between the present paper objectives and other relevant recent papers, we refer to the following:
1-Abbas et al.[31] proved the existence of solutions and studied Ulam-Hyers-Rassias stability of problem (1) in the absence of both delay and impulses effect, E=R, ψ(ϱ)=lnϱ,ϱ∈[1,b], F is a single-valued function.
2-Benchohra [19] investigated the existence and Ulam stability for nonlinear implicit differential equations with Riemann-Liouville derivative, which is including in Hilfer derivative
3-Ngo et al.[6] presented a formula of solution for a non-instantaneous impulsive differential equation containing ψ−Hilfer derivative with lower limit of the fractional derivative at zero and in the absence of delay.
4-Vanterler el al.[17] established the existence and stability of solutions for Problem (1) in the absence of delay. and when E=R, the lower limit of the fractional derivative at zero F is a single-valued function
5-Wang et al. [10] considered Problem (1) in the absence of delay, when ψ(ϱ)=ϱ and without studying the stability of solutions.
6-Wang et al.[11] consider a non-instantaneous impulsive semilinear differential inclusions containing Caputo derivative and in the absence of delay.
To clarify the novelty and contribution of this study, we refer to, in this paper, we present a relation between a solution of Problem (1) and the corresponding fractional integral equation (Lemma5), provide two methods to demonstrate the existence of solutions for Problem (1), then, investigate the ψ−generalized Ulam-Hyers stability of solutions. Because our considered problem contains ψ−Hilfer fractional derivative, non-instantaneous impulses with the lower limit of the fractional derivative switches at the impulsive points, presence of delay, and the right hand side is a multi-valued function, therefore, this study generalize recent results, as it is shown above, such as [6,10,11,17,19,31]. In addition, there isn't work in the literature, on ψ−Hilfer fractional non-instantaneous impulsive differential inclusions, in infinite dimensional spaces, in the presence of delay, and the lower limit of the fractional derivative switches at the impulsive points. Moreover, the technique presented in this paper can be used to study the existence and Ulam-Hyers stability of solutions or mild solutions for the problems considered in [3,4,15,16,20,21,22] to the case when, there are impulses and delay on the system, the right hand side is a multi-valued function and involving ψ−Hilfer fractional derivative.
In section 2, we prove some properties for ψ−fractional integral and ψ−fractional derivative, then we present, in Lemma 3.3, a relation between any solution of problem (1) and the corresponding fractional integral equation. In section 3, we prove an existence result of Problem (1). In section 4, we give another existence result of (1), then we investigate the ψ−generalized Ulam-Hyers stability of solutions. In the last section, examples are given to clarify the possibility of applicability of our assumptions.
Let Pck(E) be the family of non-empty convex and compact subsets of E. Since the given problem containing Hilfer derivative we need to the the spaces:
C1−μ,ψ(J,E):={x∈C(J∗,E):(ψ(.)−ψ(a))1−μx(.)∈C([a,b],E)}, |
and
Cn1−μ,ψ(J,E):={x∈Cn−1(J,E):[1ψ′(ϱ)ddϱ]n x∈C1−μ,ψ(J,E),n∈N. |
Obviously C1−μ,ψ(J,E) and Cn1−μ,ψ(J,E) are Banach spaces with norms
||x||C1−μ,ψ(J,E):=supϱ∈J|| (ψ(ϱ) −ψ(0))1−μx(ϱ)||, |
and
||x||Cn1−μ,ψ(J,E):=k=n−1∑k=1||x(k)||C(J,E)+||[1ψ′(ϱ)ddϱ]n x||C1−μ,ψ(J,E). |
Because Problem (1) involving impulses effect we recall the Banach space:
PC1−μ,ψ(J,E):={x:J∗→E, (ψ(.)−ψ(sk))1−μx(.)∈C(Jk,E),k=0,1,..,n,limϱ→s+k(ψ(ϱ)−ψ(sk))1−μx(ϱ)exists,x∈C(ϱi,E), andlimϱ→ϱ+ix(ϱ)exist,i=1,2,..,n}, |
endowed with the norm
||x||PC1−μ,ψ(J,E):=max{supϱ∈¯Jkk=0,1,,,,n,(ψ(ϱ)−ψ(sk))1−μ||x(ϱ)||E, supϱ∈¯ϱii=1,...,n||x(ϱ)||E}, |
where Jk=(sk,ϱk+1], ¯Jk=[sk,ϱk+1] (k=0,1,..,n),ϱi=(ϱi,si] and ¯ϱi=[ϱi,si] (i=1,2,..,n),
Next, the function χPC1−μ(J,E):Pb(PC1−μ,ψ(J,E))→[0,∞),given by by
χPC1−μ,ψ(J,E)(D):=max{maxk=0,1,..,nχC(¯Jk,E)(D|¯Jk),maxi=1,..,nχC(¯ϱi,E)(D|¯ϱi)} |
is a measure of noncompactness on PC1−μ(J,E), where
D∣¯Jk:={h∗∈C(¯Jk,E):h∗(ϱ)=(ψ(ϱ)−ψ(sk))1−μh(ϱ),ϱ∈Jk,h∗(sk)=limϱ→s+kh∗(ϱ),h∈D}, |
and
D∣¯ϱi:={h∗∈C(¯ϱi,E):h∗(ϱ)=h(ϱ),ϱ∈ϱi,h∗(ϱi)=h(ϱ+i),h∈D}. |
In the sequel, Iq,ψa+ denotes to the ψ−Riemann-Liouville fractional integral operator of order q with the lower limit at a, Dϑ,ψa+f to the ψ−Riemann-Liouville fractional derivative operator of order ϑ with the lower limit at a and cDϑ,ψa+f to the ψ−Caputo fractional derivative of order ϑ with the lower limit at a for f ∈AC1,ψ([a,b],E), where
AC1,ψ([a,b],E):={x:J→E,[1ψ′(ϱ)ddϱ] x∈AC(J,E)}. |
If f ∈C1,ψ([a,b],E):={x∈C(J,E):[1ψ′(ϱ)ddϱ] x∈C(J,E)}, then
cDϑ,ψa+f(ϱ):=I1−ϑ,ψa+[1ψ′(ϱ)ddϱ] f(ϱ),ϱ∈[a,b]. |
If ψ(ϱ)=ϱ, we obtain the Caputo fractional derivative, and if ψ(ϱ)=lnϱ, we obtain the Caputo-Hadamard fractional derivative. The following remark and more information about ψ−fractional integral and derivative can be found in [32,33,34,35]
Remark 2.1. If q=1, then I1,ψa+f(ϱ)=∫ϱaψ′(s)f(s) ds, and hence
1ψ′(ϱ)ddϱ I1,ψa+f(ϱ)=1ψ′(ϱ)ddϱ∫ϱa ψ′(s)f(s) ds=f(ϱ),a.e.for ϱ∈[a,b]. |
In the following lemma we give an important for Iϑ,ψa+, which we need later.
Lemma 2.1. Let 0<η≤ϑ. Then Iϑ,ψa+is bounded from Cη,ψ([a,b],E) into C([a,b],E) and
Iϑ,ψa+f(a)=limϱ→a+Iϑ,a+f(ϱ)=0,∀f∈Cη,ψ([a,b],E). |
Proof. The assumption f∈Cη,ψ([a,b],E) leads to (ψ(ϱ)−ψ(a))ηf(ϱ) is continuous on [a,b], and hence there is M>0 such that ||(ψ(ϱ)−ψ(a))ηf(ϱ)||≤M,∀ϱ∈ [a,b]. As a consequence, ||f(ϱ)||≤(ψ(ϱ)−ψ(a))−ηM,∀ϱ∈ (a,b]. Then
||Iϑ,ψa+f(ϱ)||=||∫ϱa (ψ(ϱ)−ψ(s))ϑ−1ψ′(s)Γ(ϑ)f(s) ds||≤MΓ(ϑ) |Iϑ,ψa+(ψ(ϱ)−ψ(a))−η|. |
In virtu of Lemma2 in[35], Iϑ,ψa+(ψ(ϱ)−ψ(a))−η=Γ(1−η)Γ(1−η+ϑ)(ψ(ϱ)−ψ(a))ϑ−η, and hence
||Iϑ,ψa+f(ϱ)||≤MΓ(ϑ)Γ(1−η)Γ(1−η+ϑ)(ψ(ϱ)−ψ(a))ϑ−η |
Since 0<η≤ϑ, we get limϱ→a+Iϑ,ψa+f(ϱ)=0, which means that Iϑ,ψa+f(ϱ) is continuous. Moreover, Iϑ,ψa+is bounded from Cη,ψ([a,b],E) into C([a,b],E).
Let us recall the definition of ψ−Hilfer fractional derivative.
Definition 1. [2] Let f∈L1([a,b],E) be such that I(1−ν)(1−ϑ),ψa+f∈AC1,ψ([a,b],E). The ψ−Hilfer fractional derivative of order 0<ϑ<1and type 0≤ν≤1 and with lower limit at a for a function f:[a,b]→E is defined by
Dϑ,ν,ψa+f(ϱ)=Iν(1−ϑ),ψa+[1ψ′(ϱ)ddϱ] (I(1−ν)(1−ϑ),ψa+f)(ϱ),ϱ∈[a,b], |
Denote
Cμ1−μ,ψ(J,E):={x∈C1−μ,ψ(J,E),Dμ,ψa+x∈C1−μ,ψ(J,E)}, |
Cϑ,ν1−μ,,ψ(J,E):={x∈C1−μ(J,E),Dϑ,ν,ψa+x∈C1−μ,ψ(J,E)}, |
PCμ1−μ,ψ(J,E):={x∈PC1−μ(J,E),Dμ,ψs+kx|Jk∈C1−μ,ψ(Jk,E),k=0,1,..,n}, |
and
PCϑ,ν1−μ,ψ(J,E):={x∈PC1−μ,ψ(J,E),Dϑ,ν,ψs+kx|Jk∈C1−μ,ψ(Jk,E), k=0,1,..,n}. |
Notice that, the operator Dϑ,ν,ψa+ can be written as:
Let
Dϑ,ν,ψa+f(ϱ)=Iν(1−ϑ),ψa+[1ψ′(ϱ)ddϱ](I1−μa+f(ϱ)=Iν(1−ϑ),ψa+Dμ,ψa+f(ϱ), μ=ϑ+ν−ϑν. |
So, if f∈Cμ1−μ,ψ([a,b],E), then by Lemma 2.1, Dϑ,ν,ψa+f(ϱ) exists ϱ∈[a,b].
Remark 2.2. Since Dϑ,ν,ψa+x=(Iν(1−ϑ),ψa+Dμ,ψa+x)(ϱ), it follows from lemma1. that Cμ1−μ,ψ([a,b],E) ⊆Cϑ,ν1−μ,ψ([a,b],E). Similarly, PCμ1−μ,ψ([a,b],E) ⊆PCϑ,ν1−μ,ψ([a,b],E).
Now, since our considered problem contains a delay we need to present the following spaces:
1- The normed space
Θ={z:[−r,0]→E suchthatz hasafinitenumberofdiscontinuitypointss≠0,allvaluesz(s+),and z(s−) arefinite}, |
endowed with the norm:
||z||Θ=∫0−r||z(s)||ds. |
2-The metric space (the space of solutions)
H:={x:[−r,b]→E,x|[−r,0]=Ψ, x|J∗∈PC1−μ,ψ(J,E)}, |
where the metric function is given by:
dH(x,y):=supϱ∈J||x(ϱ)−y(ϱ)||. |
3- The Banach space
H:={x:[−r,b]→Esuchthat x(ϱ)=0,∀ϱ∈[−r,0],x|J∗∈PC1−μ,ψ(J,E)}, |
endowed with the norm:
||x||H:=||x|J∗||PC1−μ,ψ(J,E)+||x|[−r,0]||Θ=||x|J∗||PC1−μ,ψ(J,E). |
Remark 2.3. (i) If x∈H, then x(0−)=0 and x(0+)=limϱ→0+(ψ(ϱ)−ψ(0))1−μx(ϱ).
(ⅱ) If x∈H, then x(0−)=Ψ(0) and x(0+−)=limϱ→0+(ψ(ϱ)−ψ(0))1−μx(ϱ). So, if Ψ(0)=0, then x will be continuous at zero.
It is easily seen that the function:
χH(B):=χPC1−μ(J,E){x|J:x∈B}, |
define a measure of noncompactness on H, where B is a bounded subset of H.
We need to the following fixed point for multi-valued functions.
Lemma 2.2. ([36], Theorem 3.1) Let W be a closed convex subset of a Banach space X and ϱ:W→Pc(W). Suppose that ϱ is closed, ϱ(D) is relatively compact, whenever D is compact, and that, for some x0∈W, one has
B⊆W, B=conv({x0}∪ϱ(B)), ¯B=¯C withC⊆B countable⟹B isrelativelycompact. |
Then, there is a fixed point for ϱ.
In this section, we demonstrate the existence of solutions of Problem (1). For any x∈H let
S1F(.,τ(.)x)={z∈L1(J,E):z(ϱ)∈F(ϱ,τ(.)x),a.e. for ϱ∈J }. |
and
Iν(1−ϑ),ψ(PC1−μ,ψ(J,E))={f:J→E, thereish∈PC1−μ,ψ(J,E)suchthatf(ϱ)=Iν(1−ϑ),ψsi+h(ϱ),ϱ∈Jk, k=0,1,2,..,n}. |
and for any x∈H let
¯x(ϱ):={Ψ(ϱ),ϱ∈[−r,0],x(ϱ),ϱ∈(0,b]. |
In order to derive the relation between any solution for Problem (1) and the corresponding fractional integral equation, we need to the following essential Lemmas.
Lemma 3.1. Let 0<ϑ<1,η∈[0,1).If f∈Cη,ψ([a,b],E) and I1−ϑ,ψa+f∈C1η,ψ([a,b],E), then
Iϑ,ψa+Dϑ,ψa+f(ϱ)=f(ϱ)−(ψ(ϱ)−ψ(a))ϑ−1I1−ϑ,ψa+f(a)Γ(ϑ),ϱ∈[a,b]. |
Proof. Since f∈Cη,ψ([a,b],E), then I1−ϑ,ψa+f(ϱ) is defined for ϱ∈[a,b]. Moreover, the assumption I1−ϑ,ψa+f∈C1η,ψ([a,b],E) implies to Dϑ,ψa+f(ϱ)∈Cη,ψ([a,b],E), and hence Iϑ,ψa+Dϑ,ψa+f(ϱ) and Iϑ+1,ψa+Dϑ,ψa+f(ϱ) are well defined for ϱ∈[a,b]. Observe that
Iϑ,ψa+Dϑ,ψa+f(ϱ)=1ψ′(ϱ)ddϱ(∫ϱa ψ′(s)Iϑ,ψa+Dϑ,ψa+f(s)ds)=1ψ′(ϱ)ddϱI1,ψa+Iϑ,ψa+Dϑ,ψa+f(ϱ)=1ψ′(ϱ)ddϱIϑ+1,ψa+Dϑ,ψa+f(ϱ),ϱ∈[a,b]. |
Now,
Iϑ+1,ψa+Dϑ,ψa+f(ϱ)=1Γ(ϑ+1)∫ϱa (ψ(ϱ)−ψ(s))ϑψ′(s)(1ψ′(s)dds)(I1−ϑ,ψa+f(s))ds=1Γ(ϑ+1)∫ϱa (ψ(ϱ)−ψ(s))ϑdds(I1−ϑ,ψa+f(s))ds. |
By integration by parts, we get
Iϑ+1,ψa+Dϑ,ψa+f(ϱ)=−(ψ(ϱ)−ψ(a))ϑI1−ϑ,ψa+f(a)Γ(ϑ+1)+1Γ(ϑ)∫ϱa (ψ(ϱ)−ψ(s))ϑ−1ψ′(s)I1−ϑ,ψa+f(s)ds=−(ψ(ϱ)−ψ(a))ϑI1−ϑ,ψa+f(a)Γ(ϑ+1)+Iϑ,ψa+I1−ϑ,ψa+f(ϱ)=I1,ψa+f(ϱ)−(ψ(ϱ)−ψ(a))ϑI1−ϑ,ψa+f(a)Γ(ϑ+1). |
It follows from Remark 2.1, that
Iϑ,ψa+Dϑ,ψa+f(ϱ)=1ψ′(ϱ)ddϱI1,ψa+f(ϱ)−1ψ′(ϱ)ddϱ(ψ(ϱ)−ψ(a))ϑ(I1−ϑ,ψa+f(a))Γ(ϑ+1)]=f(ϱ)−(ψ(ϱ)−ψ(a))ϑ−1I1−ϑ,ψa+f(a)Γ(ϑ). |
Lemma 3.2. Let α∈(0,1), η>0 and f∈Cη,ψ([a,b],E). Then Dα,ψa+Iα,ψa+f(ϱ)=f(ϱ),a.e. for ϱ∈[a,b].
Proof. In view of Remark 2.1, we get for a.e.ϱ∈[a,b]
Dα,ψa+Iα,ψa+f(ϱ)=[1ψ′(ϱ)ddϱ] I1−α,ψa+Iα,ψa+f(ϱ)=[1ψ′(ϱ)ddϱ] I1,ψa+f(ϱ)=f(ϱ). |
Now, we give in the following lemma the relation between any solution for Problem (1) and the corresponding fractional integral equation.
Lemma 3.3. Let 0<ϑ<1, 0≤ν≤1, μ=ϑ+ν−ϑν and x∈PCν(1−ϑ),ψ1−μ,ψ([0,b],E), gi:(ϱi,si]→E,i=1,......n, is continuous. The following hold.
(1) The function y:(0,b]→E given by
y(ϱ)={(ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+x(ϱ),ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,x(ϱ−i))+Iϑ,ψsi+x(ϱ),ϱ∈(siϱi+1],i=1,......n, | (2) |
belongs to PCμ1−μ,ψ([0,b],E), Dϑ,ν,ψsi+y(ϱ) exists for any ϱ∈(si ϱi+1], i=0,1,......n, and verifies the ψ−Hilfer fractional problem:
{Dϑ,ν,ψs+iy(ϱ)=x(ϱ),ϱ∈∪i=ni=0(si,ϱi+1]I1−μ0+y(0+)=Ψ(0),y(ϱ+i)=gi(ϱi,x(ϱ−i)),i=1,......n,y(ϱ)=gi(ϱ,x(ϱ−i)),ϱ∈(ϱisi],i=1,......n,I1−μs+iy(s+i)=gi(si,x(ϱ−i)),i=1,......n, | (3) |
(2) If y∈PCμ1−μ,ψ([0,b],E) is a solution of (3), then y satisfies (2).
Proof. Let ϱ∈(0,ϱ1]. In view of (2)
y(ϱ)=(ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+x(ϱ). | (4) |
According to Lemma 3.1 in[35], Dμ,ψ0+(ψ(ϱ)−ψ(a))μ−1=0. Then, by applying Dμ,ψa+ to both side of (4), it yields
D0+y(ϱ)=(1ψ′(ϱ)ddϱ)I(1−μ),ψ0+Iϑ,ψ0+x(ϱ)=(1ψ′(ϱ)ddϱ)I1−(μ−ϑ),ψ0+x(ϱ)=Dμ−ϑ,ψ0+x(ϱ)=Dν(1−ϑ),ψ0+x(ϱ). | (5) |
Observe that the assumption x∈PCν(1−ϑ),ψ1−μ,ψ([0,b],E) implies to Dν(1−ϑ),ψ0+x|J0∈C1−μ,ψ([0,ϱ1],E). It follows from (5) that Dμ,ψa+y∈C1−μ,ψ([0,ϱ1],E) and hence y|J0∈Cμ1−μ,ψ([0,ϱ1],E). Consequently y|J0∈Cϑ,ν1−μ,ψ([0,ϱ1],E) (see, Remark 2), and this assures that Dϑ,ν,ψ0+y(ϱ) is well defined for ϱ∈J0. Now, since Dν(1−ϑ),ψ0+x|J0∈C1−μ,ψ([0,ϱ1],E), then (1ψ′(ϱ)ddϱ)I1−ν(1−ϑ),ψ0+x|J0∈C1−μ,ψ([0,ϱ1],E). So, I1−ν(1−ϑ),ψ0+x|J0∈C11−μ,ψ([0,ϱ1],E). As a result from Lemma (3) one obtains
Dϑ,ν,ψ0+y(ϱ)=Iν(1−ϑ)0+Dμ,ψ0+y(ϱ)=Iν(1−ϑ)0+Dν(1−ϑ),ψ0+x(ϱ)=x(ϱ)−(ψ(ϱ)−ψ(0))ν(1−ϑ)−1Γ(ν(1−ϑ))(I1−ν(1−ϑ),ψ0+x)(0). |
Since I1−ν(1−ϑ),ψ0+x∈C1−μ,ψ([0,b],E) and 1−μ<1−ν(1−ϑ), then, by Lemma(1) (I1−ν(1−ϑ),ψ0+x)(0)=0. So, Dϑ,ν,ψ0+y(ϱ)=x(ϱ), ϱ∈[0,ϱ1]. It remains to demonstrate that y satisfies I1−μ0+y(0+)=ya. To do this, apply I1−μ,ψa+ to both side of (4)
I1−μ,ψ0+y(ϱ)=Ψ(0)Γ(μ)I1−μ,ψ0+(ψ(ϱ)−ψ(0))μ−1(0)+I1−μ,ψ0+Iϑ,ψa+x(ϱ)=Ψ(0)+I1−μ+ϑ,ψ0+Iϑ,ψa+x(ϱ). |
Because x∈C1−μ,ψ([0,ϱ1],E) and 1−μ<1−μ+ϑ, then, by Lemma(1),I1−μ+ϑ,ψ0+Iϑ,ψ0+x(0)=0. So, I1−μ,ψ0+y(0)=Ψ(0).
Similarly, we can show that for ϱ∈(si ϱi+1],i=1,......n,we have Dϑ,ν,ψs+iy(ϱ)=x(ϱ) and I1−μs+iy(s+i)=gi(si,x(ϱ−i)), and hence y is a solution for(3).
(2) Let y∈PCμ1−μ,ψ(J,E) be a solution of (3). Let i=0. Then y∈C1−μ,ψ(J0,E) and (1ψ′(ϱ)ddϱ)I1−μ,ψ0+y|J0∈C1−μ,ψ(J0,E). Thus, I1−μ,ψsiy|J0∈C11−μ,ψ(J0,E). By applying Lemma (3) it yields
Iμ,ψ0+Dμ,ψ0+y(ϱ)=y(ϱ)−(ψ(ϱ)−ψ(0))μ−1Γ(μ)(I1−μ,ψ0+x)(0)=y(ϱ)−(ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0),ϱ∈(0,ϱ1). | (6) |
Next, applying Iϑ,ψ0+ to both side of the equation Dϑ,ν,ψ0+y(ϱ)=x(ϱ), we get from (6)
Iϑ,ψ0+x(ϱ)=Iϑ,ψ0+Dϑ,ν,ψ0+y(ϱ)=Iϑ,ψ0+Iμ−ϑ,ψ0+Dμ,ψ0+y(ϱ)=Iμ,ψ0+Dμ,ψ0+y(ϱ)=y(ϱ)−(ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0),ϱ∈(0,ϱ1]. |
So,
y(ϱ)=(ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+x(ϱ),ϱ∈(0,ϱ1]. |
Similarly, we can show that for i=1,......n,
y(ϱ)=(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,x(ϱ−i))+Iϑ,ψ0+x(ϱ),ϱ∈(siϱi+1]. |
Now, based on Lemma 3.3, we can give the concept of solutions for problem (1).
Definition 2. A function ¯x∈H is called a mild solution of (1) if there is f∈S1F(,τ(ϱ)¯x) such that ¯x satisfies the integral equation
¯x(ϱ)={Ψ(ϱ),ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+f(ϱ),ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,x(ϱ−i))+Iϑ,ψs+if(ϱ),ϱ∈(siϱi+1],i=1,......n. |
Remark 3.1. The function ¯x is not necessarily continuous at the points si,i=0,1,......n. But if Ψ(0)=0, then it will be continuous at s0, and if gi(si,x(ϱ−i))=0, it will be continuous at si,i=1,......n.
In the following, we present our first existence result of solutions for Problem (1).
Theorem 3.1. Let F:J×Θ→Pck(E) be a multifunction, Ψ∈Θ, and gi:[ϱi,si]×E→E (i=1,2,....n). We assume the following conditions
(F1) For every z∈H, S1F(.,τ(.)z) is not empty subset of Iν(1−ϑ),ψ(PC1−μ,ψ(J,E)) and for almost every ϱ∈J,z⟶F(ϱ,τ(ϱ)z) is upper semicontinuous.
(F2) There is a φ∈L1q(I,R+),(0<q<ϑ) such that for any z∈H
‖F(ϱ,τ(ϱ)¯z‖≤φ(ϱ)(1+||z||H),a.e,ϱ∈J |
(F3) There is a ς∈L1q(I,R+),(0<q<ϑ) with the property that for any bounded subset D⊆H, any k=0,1,2,..,n, and a.e., for ϱ∈Jk
χE(F(ϱ,{τ(ϱ)¯x:x∈D))≤(ψ(ϱ)−ψ(sk))1−μς(ϱ)supθ∈[−r,0]χE{(τ(ϱ))¯x(θ):x∈D}, |
and
2κ1Γ(ϑ) η||ς||L1q(I,R+)<1, | (7) |
where η=δq(ψ(b)−ψ(0))ϑ−q(ϑ−q1−q)1−q, δ=maxs∈Jψ′(s), κ1=(ψ(b)−ψ(0))1−μ and χ is the Hausdorff measure of noncompactness on E.
(H1) for every i=1,2,…,n, gi:[ϱi,si]×E→E is uniformly continuous on bounded sets and for any ϱ∈J, gi (ϱ,.) maps any bounded subset of E to a relatively compact subset and there exists a positive constant hi such that for any x∈E
‖gi(ϱ,x))‖≤hi(ψ(ϱi)−ψ(si−1))1−μ‖x‖, ϱ∈(ϱi,si],x∈E. |
Then Problem (1) has a mild solution provided that
κ1Γ(ϑ)||φ||L1q(J,R+)η+h∗+h∗Γ(μ)<1, | (8) |
where, h∗=∑i=ni=0hi.
Proof. We define a multioperator Φ:H→P(H) as follows: let x∈H, then due to (F1) there is f ∈S1F(.,τ(.)¯x), and hence we can define y∈Φ(x) if and only if
y(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+f(ϱ),ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(0))μ−1Γ(μ)gi(si,x(ϱ−i))+Iϑ,ψs+if(ϱ),ϱ∈(siϱi+1],i=1,......n. | (9) |
Let us clarify that a point x is a fixed point for Φ if and only if ¯x is a solution for (1). Let x be a fixed point to Φ. Then
x(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+f(ϱ),ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(0))μ−1Γ(μ)gi(si,x(ϱ−i))+Iϑ,ψs+if(ϱ),ϱ∈(siϱi+1],i=1,......n. | (10) |
where f ∈S1F(.,τ(.)¯x). Therefore,
¯x(ϱ)={Ψ(ϱ),ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+f(ϱ)ds,ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(0))μ−1Γ(μ)gi(si,x(ϱ−i))+Iϑ,ψs+if(ϱ),ϱ∈(siϱi+1],i=1,......n. |
This means ¯x satisfies (4), and hence it is a solution for (1). Similarly, it is easy to see that if ¯x satisfies (4), then x is a fixed point for Φ. So we prove, by application Lemma 3.3, that Φ has a fixed point. Obviously the values of Φ are convex.
Step1. We demonstrate that there is a n∈N with Φ(Bn)⊆Bn,where Bn={x∈H:‖x‖H≤n}. Suppose that for any natural number n, there are xn,yn∈H with yn∈Φ(xn),‖xn‖H≤n and ‖yn‖H>n. Then, there is fn∈ S1F(.,τ(.)¯xn), n≥1 such that
yn(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+Iϑ,ψ0+fn(ϱ)∈(0,ϱ1]gi(ϱ,xn(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,xn(ϱ−i))+Iϑ,ψs+ifn(ϱ),ϱ∈(siϱi+1],i=1,......n. |
In view of (F2), we get for almost ϱ∈J
||fn(ϱ)||≤φ(ϱ)(1+‖xn‖H)≤φ(ϱ)(1+n). | (11) |
Then, if for almost ϱ∈(0,ϱ1], we get from Holder's inequality and (12)
(ψ(ϱ)−ψ(0))1−μ||yn (ϱ)||≤Ψ(0)Γ(μ)+(1+n)(ψ(b)−ψ(0))1−μΓ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)φ(ϱ)ds |
Now, from Holder's inequality
∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)φ(ϱ)ds≤||φ||L1q(J,R+)[∫ϱ0(ψ(ϱ)−ψ(s))ϑ−11−q(ψ′(s))11−qds]1−q≤||φ||L1q(J,R+)[∫ϱ0(ψ(ϱ)−ψ(s))ϑ−11−qψ′(s)(ψ′(s))q1−qds]1−q≤δq||φ||L1q(J,R+)[∫ϱ0(ψ(ϱ)−ψ(s))ϑ−11−qψ′(s)ds]1−q≤δq(ψ(b)−ψ(0))ϑ−q(ϑ−q1−q)1−q=η. | (12) |
Then,
(ψ(ϱ)−ψ(0))1−μ||yn (ϱ)||≤Ψ(0)Γ(μ)+(1+n)(ψ(b)−ψ(0))1−μΓ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)φ(ϱ)ds≤Ψ(0)Γ(μ)+κ1Γ(ϑ)||φ||L1q(J,R+)(1+n)η. | (13) |
If ϱ∈(ϱi si],=1,2,…,n, then by (H1)
||yn(ϱ)||=supϱ∈(ϱisi]||gi(ϱ,xn(ϱ−i))||≤h∗(ψ(ϱi) −ψ(si−1))1−μ||xn(ϱ−i)||≤h∗||xn||H≤hn. | (14) |
Similar as in (13), we get for almost ϱ∈(si ϱi+1],i=1,2,…,n.
(ψ(ϱ)−ψ(si))1−μ||yn(ϱ)||≤||gi(si,xn(ϱ−i))||Γ(μ)+ κ1Γ(ϑ)||φ||L1q(J,R+)(1+n)η≤h∗nΓ(μ)+ κ1Γ(ϑ)||φ||L1q(J,R+)(1+n)η. |
It follows from this inequality, (13) and (14), that
n<||yn||H≤Ψ(0)Γ(μ)+h∗n+h∗nΓ(μ)+κ1Γ(ϑ)||φ||L1q(J,R+)(1+n)η. |
By dividing both side by n and taking the limit as n→∞, one obtains
1<κ1Γ(ϑ)||φ||L1q(J,R+)η+h∗+h∗Γ(μ), |
which contradicts with (8). Thus, there is a natural number n0 such that Φ(Bn0)⊆Φ(Bn0).
Step2. Φ is closed on Bn0.
Let (xn)n≥1, (yn)n≥1be two sequences in Bn0 with xn→x in Bn0,yn→y in H and yn∈Φ(xn);n≥1. Then, there is fn ∈S1F(.,τ(.)¯xn) such that
yn(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)fn(s)ds,ϱ∈(0,ϱ1]gi(ϱ,xn(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,xn(ϱ−i))+1Γ(ϑ)∫ϱsi(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)fn(s),ϱ∈(siϱi+1],i=1,......n. | (15) |
According to (11), it follows that ||fn(ϱ)||≤φ(ϱ)(1+n),a.e.ϱ∈J, and hence (fn)n≥1 is bounded in L1q(J,E). Since L1q(J,E) is reflexive, then, by using Mazur's Lemma, there is sequence, (zj)j≥1, of convex combinations of (fn)n≥1 converging strongly to f in L1(J,E) as j→∞. Notice that, by (F2) again, for every ϱ∈J,s∈(0,ϱ] and every n≥1
||(ϱ−s)ϑ−1fn(s)||≤ |ϱ−s)|ϑ−1φ(ϱ)(1+n)∈L1((0,ϱ],R+). |
Let
˜yn(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)zn(s)ds,ϱ∈(0,ϱ1]gi(ϱ,xn(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,xn(ϱ−i))+1Γ(ϑ)∫ϱsi(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)zn(s),ϱ∈(siϱi+1],i=1,......n. | (16) |
Clearly ˜yn(ϱ)→y(ϱ), ϱ∈J and zn(ϱ)→f(ϱ), for almost ϱ∈J. Also, τ(ϱ)¯xn→τ(ϱ)¯x;ϱ∈J, and hence, by the upper semicontinuity of F(ϱ,.);a.e. ϱ∈J, it follows that f(ϱ)∈F(ϱ,τ(ϱ)¯x),a.e.[[37]. Theorem1, Sec. 4, Ch.1 ]. Therefore, by, the uniform continuity of gi(si,.) on bounded subsets and by passing to the limit as n→∞ in (15) we obtain, from the Lebesgue dominated convergence theorem y∈Φ(x).
Step3. We show that M∣¯Jk (k=0,1,...,n) and M∣¯ϱi (i=1,2,....,n) are equicontinuous, where
M∣¯Jk={z:¯Jk→E, z(ϱ)=(ψ(ϱ)−ψ(sk))1−μy(ϱ),ϱ∈Jk,z(sk)=limϱ→sk(ψ(ϱ)−ψ(sk))1−μy(ϱ), y∈Φ(x),x∈Bn0}, |
and
M∣¯ϱi={y∗∈C(¯ϱi,E):y∗(ϱ)=y(ϱ),ϱ∈(ϱi,si],y∗(ϱi)=y(ϱ+i),y∈Φ(x),x∈Bn0}. |
Case 1.let z∈M∣¯J0.Then there are x∈Bn0and y∈Φ(x) such that for ϱ∈(0,ϱ1],
y(ϱ)=(ψ(ϱ)−ψ(0))1−μΨ(0)Γ(μ)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s)ds, |
z(ϱ)=(ψ(ϱ)−ψ(0))1−μy(ϱ) and z(0)=limϱ→0+(ψ(ϱ)−ψ(0))1−μy(ϱ), where f ∈S1F(.,τ(.)¯x). It follows
limδ→0+z(0+δ)=limδ→0+z(δ)=limδ→0+(ψ(δ)−ψ(0))1−μy(δ)=z(0). |
Let ϱ1,ϱ2 be two points in (0,ϱ1] be such that ϱ1<ϱ2. Then,
‖z(ϱ2)−z(ϱ1)‖≤||(ψ(ϱ2)−ψ(0))1−μΓ(ϑ)∫ϱ20(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)f(s)ds−(ψ(ϱ1)−ψ(0))1−μΓ(ϑ)∫ϱ10(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)f(s)ds||≤(ψ(ϱ2)−ψ(0))1−μΓ(ϑ)||∫ϱ2ϱ1(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)f(s)ds||+||∫ϱ10(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)||f(s)||ds× |(ψ(ϱ2)−ψ(0))1−μΓ(ϑ)−(ψ(ϱ1)−ψ(0))1−μΓ(ϑ)|+(ψ(ϱ1)−ψ(0))1−μΓ(ϑ)×||∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)−(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)]||f(s)||ds||. |
By the absolute integral of the Lebesgue integral and Holder's inequality, it yields from (12)
limϱ2→ϱ1(ψ(ϱ2)−ψ(0))1−μΓ(ϑ)∫ϱ2ϱ1(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)||f(s)||ds≤κ1Γ(ϑ)(1+n0)limϱ2→ϱ1∫ϱ2ϱ1(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)φ(s)ds≤κ1Γ(ϑ)(1+n0)=0, |
independent of x.
Since ψ is continuous, we get by (12)
||∫ϱ10(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)Γ(ϑ)f(s)ds||× |(ψ(ϱ2)−ψ(0))1−μ−(ψ(ϱ1)−ψ(0))1−μ|≤(1+n0)∫ϱ10(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)Γ(ϑ)φ(s)ds×limϱ2→ϱ1|(ψ(ϱ2)−ψ(0))1−μ−(ψ(ϱ1)−ψ(0))1−μ|=0 |
Moreover,
limϱ2→ϱ1(ψ(ϱ1)−ψ(0))1−μΓ(ϑ)×||∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)−(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)]f(s)ds||≤κ1Γ(ϑ)(1+n0)×limϱ2→ϱ1∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)−(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)]φ(s)ds≤κ1Γ(ϑ)(1+n0)||φ||L1q(I,R+)×limϱ2→ϱ1[∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)−(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)]11−qds]1−q |
Put ¯ω=ϑ−11−q∈(−1,0), then for s<ϱ1<ϱ2, we have (ψ(ϱ1)−ψ(s))¯ω≥(ψ(ϱ2)−ψ(s))¯ω. By applying Lemma 3.1 in [43] and taking into account 1−q∈(0,1), we get
|[(ψ(ϱ1)−ψ(s))¯ω]1−q−[(ψ(ϱ2)−ψ(s))¯ω]1−q| ≤[(ψ(ϱ1)−ψ(s))¯ω−(ψ(ϱ2)−ψ(s))¯ω]1−q. |
Then,
|(ψ(ϱ1)−ψ(s))ϑ−1−(ψ(ϱ2)−ψ(s))ϑ−1|≤ [(ψ(ϱ1)−ψ(s))¯ω−(ψ(ϱ2)−ψ(s))¯ω]1−q. |
This leads to
|(ψ(ϱ1)−ψ(s))ϑ−1−(ψ(ϱ2)−ψ(s))ϑ−1|11−q≤(ψ(ϱ1)−ψ(s))¯ω−(ψ(ϱ2)−ψ(s))¯ω. |
Therefore,
limϱ2→ϱ1[∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)−(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)]11−qds]1−q≤limϱ2→ϱ1[∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1−(ψ(ϱ1)−ψ(s))ϑ−1]11−q(ψ′(s))11−qds]1−q≤limϱ2→ϱ1[∫ϱ10[(ψ(ϱ1)−ψ(s))¯ω−(ψ(ϱ2)−ψ(s))¯ω](ψ′(s))11−qds]1−q≤δqlimϱ2→ϱ1[∫ϱ10[(ψ(ϱ1)−ψ(s))¯ω−(ψ(ϱ2)−ψ(s))¯ω]ψ′(s)ds]1−q=δqlimϱ2→ϱ1[(ψ(ϱ1)−ψ(s))¯ω+1¯ω+1−(ψ(ϱ2)−ψ(s))¯ω+1¯ω+1]1−q=0. |
Thus,
limϱ2→ϱ1(ψ(ϱ1)−ψ(0))1−μΓ(ϑ)×||∫ϱ10[(ψ(ϱ2)−ψ(s))ϑ−1ψ′(s)−(ψ(ϱ1)−ψ(s))ϑ−1ψ′(s)]f(s)ds||=0 |
independent of x.
Case 2. Let i∈{1,2,..,n} be fixed and set Ki={x(ϱ−i):x∈Bn0}.ObviouslyKi is bounded subset of E.
Let y∈M∣ϱi. Then, there is x∈Bn0 such that
y(ϱ)=gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si] |
It follows, from the uniform continuity of gi on the bounded set [ϱi,si]×Ki, that for ϱ1,ϱ2∈(ϱi,si]
limϱ2→ϱ1||y(ϱ2)−y(ϱ1)||=limϱ2→ϱ1||gi(ϱ2,x(ϱ−i))−gi(ϱ1,x(ϱ−i))||=0, |
independent of x. When ϱ=ϱi, let δ>0 be such that ϱi+δ∈(ϱi,si] and λ>0 with ϱi<λ<ϱi+δ≤si. Then, we have
limδ→0+‖y∗(ϱi+δ)−y∗(ϱi)‖=limδ→0+limλ→ϱ+i‖y(ϱi+δ)−y(λ)‖=0. |
Case 3. let k=1,…,n be fixed, z∈M∣¯Jk. Then, there are x∈Bn0and y∈Φ(x) such that for ϱ∈ϱ∈(sk,ϱk+1],
y(ϱ)=(ψ(ϱ)−ψ(sk))μ−1Γ(μ)gk(sk,xn(ϱ−k))+1Γ(ϑ)∫ϱsk(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)zn(s), |
z(ϱ)=(ψ(ϱ)−ψ(sk))1−μy(ϱ) and z(sk)=limϱ→sk+(ψ(ϱ)−ψ(sk))1−μy(ϱ), where f ∈S1F(.,τ(.)¯x). It follows
limδ→0+z(sk+δ)=limδ→0+(ψ(sk+δ)−ψ(sk))1−μy(sk+δ)=limϱ→sk+(ψ(ϱ)−ψ(sk))1−μy(ϱ)=z(sk). |
Next, let ϱ1,ϱ2∈(sk,ϱk+1](ϱ1< ϱ2).By arguing as in case 1, one can show that
limϱ2→ϱ1‖z(ϱ+δ)−z(ϱ)‖=0, |
independent of x.
Step 4. Let K⊆Bn0, K=conv({0}∪Φ(K)), ¯Z=¯C with C⊆K countable. We have to show that K is relatively compact in H Let D={yn:n≥1}⊆Φ(K) with C⊆conv({x0}∪D), xn∈K with yn∈Φ(xn).This means that, there is fn∈S1F(.,τ(.)¯xn) such that the relation (15) holds. Observe that, from (F3) it holds for a.e.s∈J0
χE{fn(s):n≥1}≤(ψ(ϱ)−ψ(0))1−μς(ϱ)χ{F(s,τ(s)¯xn):n≥1}≤ζ(s)(ψ(ϱ)−ψ(0))1−μsupθ∈[−r,0] χ{¯xn(s+θ):n≥1}≤ζ(s)(ψ(ϱ)−ψ(0))1−μsupδ∈[−r,s] χ{¯xn(δ):n≥1}. |
Because supδ∈[−r,0] χ{¯xn(δ):n≥1}=supδ∈[−r,0] χ{Ψ(δ)}=0. Thus,
χE{fn(s):n≥1}≤ζ(s)(ψ(ϱ)−ψ(0))1−μ supδ∈[0,s]χ{xn(δ):n≥1}=ζ(s)supδ∈[0,s]χ{(ψ(ϱ)−ψ(0))1−μxn(δ):n≥1}≤ζ(s)χPC1−μ,ψ(J,E){xn:n≥1}≤ζ(s)χH(Bn−1). | (17) |
According to the definition χH(Z), one obtains
χH(Z)=χH(¯Z)=χH(¯C)=χH(C))≤χH(conv({x0}∪D))=χH(D)=max{maxk=0,1,..,nχC(¯Jk,E)(D|¯Jk),maxi=1,..,nχC(¯ϱi,E)(D|¯ϱi)}. | (18) |
In view of Step 3, D|¯Ji and D|¯ϱi are equicontinuous, consequently
χH(K)≤max{maxi=0,1,..,nmaxϱ∈¯Jkχ{y∗n(ϱ):n≥1},maxi=1,..,nmaxϱ∈¯ϱiχ{y∗n(ϱ):n≥1}}, | (19) |
Let ϱ∈(0,ϱ1]. Then,
y∗n(ϱ)=Ψ(0)Γ(μ)+(ψ(ϱ)−ψ(0))1−μΓ(ϑ)∫ϱ0(ψ(ϱ)−ψ(0))ϑ−1ψ′(s)fn(s)ds, |
which yields with (17)
χ{y∗n(ϱ):n≥1}≤(ψ(ϱ)−ψ(0))1−μΓ(ϑ)∫ϱ0(ψ(ϱ)−ψ(0))ϑ−1ψ′(s)χ{fn(s):n≥1}ds≤(ψ(b)−ψ(0))1−μΓ(ϑ)χH(Bn−1)∫ϱ0(ψ(ϱ)−ψ(0))ϑ−1ψ′(s)ζ(s)ds≤κ1Γ(ϑ) η||ς||L1q(I,R+)χH(Bn−1). | (20) |
Notice that
χ{y∗n(0):n≥1}=χ{limϱ→0+y∗n(ϱ):n≥1}=χ{Ψ(0)Γ(μ)}=0. | (21) |
Moreover, from the fact that ||xn||H≤n0,the set {xn(ϱ−i):n≥1} is bounded for every i=1,2,...,n and hence from the assumption (H1), we get
χ{gi(ϱ,xn(ϱ−i)),n≥1}=0,ϱ∈(ϱi,si],i=1,...,n. | (22) |
and
χ{gi(ϱi,xn(ϱ−i)),n≥1}=χ{gi(si,xn(ϱ−i)):n≥1}=0. | (23) |
Similarly,
χ{y∗n(si):n≥1}=χ{limϱ→s+i(ψ(ϱ)−ψ(si))1−μyn(ϱ),:n≥1}=χ{1Γ(μ)gi(si,xn(ϱ−i)):n≥1}=0,i=1,2,..,n. | (24) |
By arguing as in (20), one can show that
maxϱ∈¯Jkχ{y∗n(ϱ):n≥1}≤κ1Γ(ϑ) η||ς||L1q(I,R+)χH(Bn−1),k=1,2,..n. | (25) |
From the relations (7), (19)–(25), it follows that
χH(K)≤χH(K)2κ1Γ(ϑ) η||ς||L1q(I,R+) <χH(K). |
So, χPC1−μ(J,E)(K)=0, and hence K is relatively compact.
Step5. Φ maps compact sets into relatively compact sets.
Let B be a compact subset of Bn0, {yn,:n≥1}⊆Φ(B) Then, there is xn∈B,n≥1, such that yn ∈Φ(xn). So, there is fn∈S1F(.,¯xn(.)) such that (15) holds. We have to show that the set Z={yn:n≥1} is relatively compact in H. Note that, since B is compact in H, then from (F3) we get for a.e. s∈J0,
χE{fn(s):n≥1}≤(ψ(ϱ)−ψ(0))1−μς(ϱ)χ{F(s,τ(s)¯xn):n≥1}≤ζ(s)(ψ(ϱ)−ψ(0))1−μsupθ∈[−r,0] χ{¯xn(s+θ):n≥1}≤ζ(s)(ψ(ϱ)−ψ(0))1−μsupδ∈[−r,s] χ{¯xn(δ):n≥1}≤ζ(s)(ψ(ϱ)−ψ(0))1−μ supδ∈[0,s]χ{xn(δ):n≥1}=ζ(s)supδ∈[0,s]χ{(ψ(ϱ)−ψ(0))1−μxn(δ):n≥1}≤ζ(s)χPC1−μ,ψ(J,E){xn:n≥1}≤ζ(s)χH(B)=0. |
By the same reasons, one can show that for a.e. s∈Jk, k=1,2,..,n
χE{fn(s):n≥1}=0. |
By arguing as in the previous step one can show that Z is relatively compact, and hence Φ(B) is relatively compact.
Now, by applying Lemma 8, there is x∈H andf ∈S1F(.,τ(.)¯x) such that
x(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s)ds,ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(0))μ−1Γ(μ)gi(si,x(ϱ−i))+1Γ(ϑ)∫ϱsi(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s),ϱ∈(siϱi+1],i=1,......n. |
Next, in view of (F1), there is h∈PC1−μ,ψ(J,E) such that f(ϱ)=Iν(1−ϑ),ψsi+h(ϱ),ϱ∈Jk,k=0,..,n, and hence, from Lemma 3.3,
Dν(1−ϑ),ψsi+f(ϱ)=Dν(1−ϑ),ψsi+Iν(1−ϑ)ψsi+h(ϱ)=h(ϱ),ϱ∈Jk,k=0,..,n. |
This yields that f∈PCν(1−ϑ)1−μ,ψ(J,E). Then, the function
¯x(ϱ)={Ψ(ϱ),ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s)ds,ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,x(ϱ−i))+1Γ(ϑ)∫ϱsi(ψ(ϱ)−ψ(s))ϑ−1ψ′(s),ϱ∈(siϱi+1],i=1,......n. | (26) |
belongs to H and in view of Lemma(6) it is a solution for (1). This completes the proof.
Remark 3.2. Theorem 1 remains true if condition (F2) is replaced by the following assumption:
(F2)∗There is a φ∈L1q(I,R+) such that for any z∈Θ
‖F(ϱ,z)‖≤{φ(ϱ) (1+(ψ(ϱ)−ψ(si))1−μ‖z(0)‖E),a.e.ϱ∈∪i=ni=0(siϱi+1],φ(ϱ) (1+‖z(0)‖E),a.e.ϱ∈∪i=ni=1(ϱi,si]. |
In fact, condition (F2) is used only to prove relation (11). We show this relation by using (F2)∗. Let fn∈ S1F(.,τ(.)¯xn), n≥1 Then, by (F2)∗, for almost ϱ∈∪i=ni=0(si ϱi+1] we get
||fn(ϱ)||≤φ(ϱ)(1+(ψ(ϱ)−ψ(si))1−μ||τ(ϱ)¯xn(0)||)≤φ(ϱ)(1+(ψ(ϱ)−ψ(si))1−μ||¯xn(ϱ)||)≤φ(ϱ)(1+(ψ(ϱ)−ψ(si))1−μ||xn(ϱ)||≤φ(ϱ)(1+‖xn‖H). |
and for almost ϱ∈∪i=ni=1(ϱi si]
||fn(ϱ)||≤φ(ϱ)(1+||τ(ϱ)¯xn(0)||E)≤φ(ϱ)(1+¯xn(ϱ))≤φ(ϱ)(1+||xn(ϱ)||)≤φ(ϱ)(1+‖xn‖H). |
In this section, we give another version for the existence of solutions and investigate the generalized ψ−generalized Ulam-Hyers stability of problem (1). For basic information about multivalued weakly Picard operators we refer to [38].
Definition 3. [39,40] A An increasing function ς:[0,∞)→[0,∞) is called comparison function if limn→∞ςn(s)=0, ∀s∈[0,∞),where ςn(s)=ςn−1(ς(s)).If in addition ∞∑n=1ςn(ϱ)<∞,∀s∈(0,∞), then ς:[0,∞)→[0,∞) is called strictly comparison.
Remark 4.1. [41]
1- If ς:[0,∞)→[0,∞) is comparisons, then ς(s)<s,∀s>0,ς(0)=0 and ς is continuous at 0.
2-the functions ς1(s)=cϱ;c∈[0,1) and ς2(s)=ss+1:s∈[0,∞) are strictly comparison
Let M:={y∈H:y|J∈PCμ1−μ,ψ(J,E), Dϑ,ν,ψsi+y(ϱ) exists for any ϱ∈(si ϱi+1], i=0,1,......n,}
Definition 4. [12] Problem (1) is called ψ−generalized Ulam- Hyers stable if there is a continuous function θ:[0,∞)→[0,∞) and θ(0)=0 such that for each ϵ>0 and each solution y∈M of the inequality
{d(Dϑ,ν,ψs+iy(ϱ),F(ϱ,τ(ϱ)y))≤ϵ,a.e.ϱ∈(si,ϱi+1],i=0,1,....,n,||y(ϱ)−gi(ϱ,y(ϱ−i))||≤ϵ,ϱ∈(ϱisi],i=1,......n,,||y(ϱ+i)−gi(ϱi,y(ϱ−i))||≤ϵ,i=1,......n, | (27) |
there is a solution x∈H for (1) with
supϱ∈J||x(ϱ)−y(ϱ)||≤θ(ϵ). | (28) |
In the following theorem, we establish the existence and generalized ψ− Ulam-Hyers stability of solutions Problem (1).
Theorem 4.1. Let F:J×Θ→Pck(E), Ψ∈Θ andgi:[ϱi,si]×E→E (i=1,2,....n)be such that:
(F4) For every z∈H, S1F(.,τ(.)z) is a non-empty subset of Iν(1−ϑ),ψ(PC1−μ,ψ(J,E)).
(F5) There is a function σ∈L1q(I,R+),0<q<ϑ and a strict comparison function ς:[0,∞)→[0,∞) such that
(ⅰ) For every x1,x2∈H
h(F(ϱ,τ(ϱ)¯x1),F(ϱ,τ(ϱ)¯x2))≤σ(ϱ)ς(‖x1−x2‖)H,∀ϱ∈J. |
(ⅱ)
Sup{ ||y|| : y∈F(ϱ,τ(ϱ)¯x0)}≤σ(ϱ),fora.e.ϱ∈J, |
where ¯x0(s)=Ψ(s);s∈[−r,0] and ¯x0(s)=0;s∈(0,b].
(H2)For anyi=1,2,....n, there ξi >0 such that for any ϱ∈[ϱi,si]
‖gi(ϱ,x)−gi(ϱ,y)‖≤ξi ς(‖x−y‖E),∀x,y∈E. |
Then, there is solution for Problem(1) provided that
ξΓ(μ)+ ηδq(ψ(b)−ψ(0))1−μΓ(ϑ)||σ||L1q(I,R+)<1, | (29) |
where ξ=∑i=ni=1ξi.Moreover, if, in addition there, is a c>1such that ς(cϱ)≤c ς(ϱ);ϱ∈[0,∞) and ϱ=0 is a point of uniform convergence for the series ∞∑n=1ςn(ϱ), then problem (1) is a ψ−generalized Ulam-Hyers stable.
Proof. Condition (F4) allows to define a multifunction Φ:H→P(H) as in (9) Note that by (F5), for every x∈H, and for a.e.ϱ∈J
||F(ϱ,τ(ϱ)¯x)||=h(F(ϱ,τ(ϱ)¯x),{0})≤h(F(ϱ,τ(ϱ)¯x),F(ϱ,τ(ϱ)¯x0))+h(F(ϱ,τ(ϱ)¯x0),{0})≤σ(ϱ)ς(‖x‖H)+σ(ϱ)≤σ(ϱ)(1+ς(‖x‖H)),for a.e.ϱ∈J. |
Then, as the reasons of Theorem 1, we can show that the values of Φ are closed. We show that Φ is ς−contraction. Let x2,x1∈H and y1∈Φ(x1). Then, there exists f∈S1F(.,τ(.)¯x1)such that
y1(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s)ds,ϱ∈(0,ϱ1]gi(ϱ,x(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,x(ϱ−i))+1Γ(ϑ)∫ϱsi(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s),ϱ∈(siϱi+1],i=1,......n. |
Define Π:J→2E as
Π(ϱ)={u∈E:||f (ϱ)−u || ≤σ(ϱ)ς(‖x1−x2‖H)}. |
From (F5)(i), we have
h(F(ϱ,τ(ϱ)¯x1),F(ϱ,τ(ϱ)¯x2))≤σ(ϱ)ς(‖x1−x2‖H),ϱ∈J, |
Since the values of F are compact, there exists uϱ ∈F(ϱ,τ(ϱ)x2) such that
||uϱ−f(ϱ)|| ≤σ(ϱ)ς(‖x1−x2‖H), |
as a consequence Π(ϱ);ϱ∈J is not empty. Furthermore, because f,σ,x1,x2 are measurable and E is separable, it follows from [[42], Theorem Ⅲ.41], the multifunction s→ Π(s)∩F(s,τ(s)¯x2) is measurable and since its values are non-empty and compact, there is h∈S1F(ϱ,τ(ϱ)¯x2) such that
||h(ϱ)−f(ϱ)|| ≤ σ(ϱ)ς(‖x1−x2‖H),a.e.ϱ∈J. |
Set
y2(ϱ)={0,ϱ∈[−r,0](ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)h(s)ds,ϱ∈(0,ϱ1]gi(ϱ,x2(ϱ−i)),ϱ∈(ϱi,si],i=1,......n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,x2(ϱ−i))+1Γ(ϑ)∫ϱsi(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)h(s),ϱ∈(siϱi+1],i=1,......n. |
Notice that, y2∈ Φ(x2) and if ϱ∈J0,, we get from Holder's inequality
(ψ(ϱ)−ψ(0))1−μ||y2(ϱ)−y1(ϱ)|| ≤(ψ(ϱ)−ψ(0))1−μΓ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)σ(s) ||h(s)−f(s)||ds≤(ψ(ϱ)−ψ(0))1−μΓ(ϑ)ς(‖x1−x2‖H)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)σ(s)ds≤(ψ(ϱ)−ψ(0))1−μΓ(ϑ) ς(‖x1−x2‖H)||σ||L1q(I,R+)δq(ψ(b)−ψ(0))ϑ−q(ϑ−q1−q)1−q≤ηδq(ψ(b)−ψ(0))1−μΓ(ϑ)||σ||L1q(I,R+)ς(‖x1−x2‖H). | (30) |
Similarly, if ϱ∈(si ϱi+1],i=1,......,n, we get
(ψ(ϱ)−ψ(si))||y2(ϱ)−y1(ϱ)||≤1Γ(μ)||gi(si,x1(ϱ−i))−gi(si,x2(ϱ−i))||+ηδq(ψ(b)−ψ(0))1−μΓ(ϑ)||σ||L1q(I,R+)ς(‖x1−x2‖H).≤ξΓ(μ)ς(‖x1−x2‖H)+ηδq(ψ(b)−ψ(0))1−μΓ(ϑ)||σ||L1q(I,R+)ς(‖x1−x2‖H).≤ς(‖x1−x2‖H)[ξΓ(μ)+ ηδq(ψ(b)−ψ(0))1−μΓ(ϑ)||σ||L1q(I,R+)], | (31) |
where ξ=∑i=ni=1ξi. Furthermore, If ϱ∈(ϱi,si],i=1,......n, then by (H2)one obtains
||y2(ϱ)−y1(ϱ)||≤||gi(ϱ,x(ϱ−i))−gi(ϱ,x(ϱ−i))||≤ξς(||x1(ϱ−i)−x2(ϱ−i)||) ≤ξ ς(‖x1−x2‖H). | (32) |
It follows from (30)–(32),
hH(Φ(x2),Φ(x1))<ς(‖x1−x2‖H)[ξΓ(μ)+ ηδq(ψ(ϱ)−ψ(0))1−μΓ(ϑ)||σ||L1q(I,R+)], |
where, hH is the Hausdorff distance in H. This inequality and (29) imply that Φ is ς−contraction and thus by Theorem 3.1(i) in [41], Φ has a fixed point and as in Theorem 1, this fixed point is a solution for Problem (1).
Now, in order to demonstrate that Problem(1) is a ψ−generalized Ulam-Hyers stable. Let ϵ>0 and y∈M be a solution of inequality (27). Because the values of F are compact, there f∈S1F(.,τ(ϱ)y) such that
||Dϑ,ν,ψs+iy(ϱ)−f(ϱ)||=d(Dϑ,ν,ψs+iy(ϱ),F(ϱ,τ(ϱ)y)),a.e.ϱ∈(si,ϱi+1],i=0,1,....,n. |
Then, for almost everywhere ϱ∈∪i=ni=0(si,ϱi+1]
Dϑ,ν,ψs+iy(ϱ)=w(ϱ)+f(ϱ), |
where w∈PCν(1−ϑ)1−μ,ψ(J,E) and ||w(ϱ)|≤ϵ,∀ϱ∈J. Furthermore, according to (F4), there is h∈PC1−μ,ψ(J,E) such that f(ϱ)=Iν(1−ϑ),ψsi+h(ϱ),ϱ∈Jk, k=0,1,2,..,n, and hence, from Lemma 3.2,
Dν(1−ϑ),ψsi+f(ϱ)=Dν(1−ϑ),ψsi+Iν(1−ϑ)ψsi+h(ϱ)=h(ϱ),ϱ∈Jk, k=0,1,2,..,n. |
This yields that f∈PCν(1−ϑ)1−μ,ψ(J,E). Therefore, y(ϱ)=Ψ(ϱ),ϱ∈[−r,0] and
{Dϑ,ν,ψs+iy(ϱ)=f(ϱ)+ϑ(ϱ),a.e.ϱ∈(si,ϱi+1],i=0,1,....,n,y(ϱ)=gi(ϱ,y(ϱ−i))+ϵ,ϱ∈(ϱisi],i=1,......n,,y(ϱ+i)=gi(ϱi,y(ϱ−i))+ϵ,i=1,......n, | (33) |
In view of the second assertion of Lemma 3.3, relation (33) one obtains
y(ϱ)={Ψ(ϱ),ϱ∈[−r,0] (ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1f(s)ds+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ϑ(ϱ)ds,ϱ∈(0,ϱ1],gi(ϱ,y(ϱ−i))+ϵ,ϱ∈(ϱisi],i=1,2,…,n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)y(si)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)(f(s))ds+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)ϑ(ϱ)ds,ϱ∈(siϱi+1],i=1,2,…,n. | (34) |
Let
z(ϱ)={0,ϱ∈[−r,0] (ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s)ds,ϱ∈(0,ϱ1],gi(ϱ,y(ϱ−i)),ϱ∈(ϱisi],i=1,2,…,n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,y(ϱ−i))+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)(f(s))ds|+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)ϑ(ϱ)ds,ϱ∈(siϱi+1],i=1,2,…,n. |
Obviously, z∈Φ(y∗), where
y∗(ϱ)={0,ϱ∈[−r,0] (ψ(ϱ)−ψ(0))μ−1Γ(μ)Ψ(0)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)f(s)ds+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ϑ(ϱ)ds,ϱ∈(0,ϱ1],gi(ϱ,y(ϱ−i))+ϵ,ϱ∈(ϱisi],i=1,2,…,n,(ψ(ϱ)−ψ(si))μ−1Γ(μ)y(si)+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)(f(s))ds|+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)ϑ(ϱ)ds,ϱ∈(siϱi+1],i=1,2,…,n. |
Moreover, from (33) and (34) we get for ϱ∈(0,ϱ1]
||z(ϱ)−y∗(ϱ)||≤1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)||ϑ(ϱ)||ds≤ϵ(ψ(b)−ψ(0))ϑΓ(ϑ+1). |
Ifϱ∈(ϱi si], then||¯z(ϱ)−y(ϱ)||≤ϵ.If ϱ∈(si ϱi+1], then
||z(ϱ)−y∗(ϱ)||≤||(ψ(ϱ)−ψ(si))μ−1Γ(μ)gi(si,y(ϱ−i))(ψ(ϱ)−ψ(si))μ−1Γ(μ)y(si)||+1Γ(ϑ)∫ϱ0(ψ(ϱ)−ψ(s))ϑ−1ψ′(s)||ϑ(ϱ)||ds≤(ψ(ϱ)−ψ(si))μ−1Γ(μ)||gi(si,y(ϱ−i))−y(si)||+ϵ(ψ(b)−ψ(0))ϑΓ(ϑ+1)≤ϵ(ψ(b)−ψ(0))μ−1Γ(μ)+ϵ(ψ(b)−ψ(0))ϑΓ(ϑ+1)=ϵ[(ψ(b)−ψ(0))μ−1Γ(μ)+(ψ(b)−ψ(0))ϑΓ(ϑ+1)+1] |
As a consequence,
supϱ∈J||y∗(ϱ)−z(ϱ)||≤ϵ[(ψ(b)−ψ(0))μ−1Γ(μ)+(ψ(b)−ψ(0))ϑΓ(ϑ+1)+1]. | (35) |
On the other hand, from the facts that Φ is ς- contraction, ς is strictly comparison, ς(cϱ)≤c ς(ϱ)(c>1) for everyϱ∈[0,∞) and ϱ=0 is a point of uniform convergence for the series ∞∑n=1ςn(ϱ), it yields from Theorem3.1(ii)[41], that Φ is ρ−multivalued weakly Picard operator, where ρ(ϱ)=ϱ+∞∑n=1 ςn(ϱ). Then, the function Φ∞:Graph(Φ)→Fix(Φ) is well defined and
||y∗−Φ∞(y∗,z)||H≤ρ(||y∗−z||H). | (36) |
Put x=Φ∞(y∗,z). So, x∈Φ(x) and from (35) and (36) we get
dH(y,¯x)≤||y∗−Φ∞(y∗,z)||H≤ρ(||y∗−z||H)=ρ(ϵ[(ψ(b)−ψ(0))μ−1Γ(μ)+(ψ(b)−ψ(0))ϑΓ(ϑ+1)+1])=θ(ϵ), |
where θ(ϱ)=ρ(ϱ[(ψ(b)−ψ(0))μ−1Γ(μ)+(ψ(b)−ψ(0))ϑΓ(ϑ+1)+1]). Consequently, Problem (1) is ψ− generalized Ulam-Hyers stable.
In this section we give examples to clarify the possibility of applicability our assumptions.
EXAMPLE 1 Let E be a Hilbert space, K a non-empty convex compact subset of E,r=1,s0=0,ϱ1=1,s1=2,ϱ2=3, ϑ=12,ν=14and μ=ϑ+ν−ϑν=58. Let F:J×Θ→Pck(E) be defined as follows:
F(ϱ,z)={(ψ(ϱ)−ψ(0))38||z(0)||K, ϱ∈[0,1],||z(0)||K,ϱ∈(1,2],(ψ(ϱ)−ψ(2))38||z(0)||K,ϱ∈(2,3]. | (37) |
Clearly, for almost every ϱ∈J,z⟶F(ϱ,τ(ϱ)z) is upper semicontinuous. Set λ=sup{||x||:x∈K}. Then, for any (ϱ,z)∈J×Θ and any y∈F(ϱ,z) we have
||y(ϱ)||≤{(1+(ψ(ϱ)−ψ(0))38||z(0)||)λ, ϱ∈[0,1],||z(0)||λ,ϱ∈(1,2],(1+(ψ(ϱ)−ψ(2))38||z(0)||)λ,ϱ∈(2,3]. |
Then condition (F2)of Theorem1 is verified with φ(ϱ)=λ;ϱ∈J. Moreover, if f∈S1F(,.z);z∈Θ, then
f(ϱ)={(ψ(ϱ)−ψ(0))38||z(0)||y0, ϱ∈[0,1],y1,ϱ∈(1,2],(ψ(ϱ)−ψ(2))38y2,ϱ∈(2,3]. | (38) |
We define
hx1,x2,x3(ϱ)={Γ(118)Γ(108)(ψ(ϱ)−ψ(0))28x1, ϱ∈[0,1],x2,ϱ∈(1,2],Γ(118)Γ(108)(ψ(ϱ)−ψ(2))28x3,ϱ∈(2,3]. | (39) |
Obviously hx1,x2,x3∈PC1−μ,ψ(J,E) and from Lemma 2.2, we get Iν(1−ϑ),ψs+ihx1,x2,x3(ϱ)=f(ϱ);i=0,1.
Now, let D be a bounded subset of H, z1, z2∈D, ϱ∈(0,1], x∈F(ϱ.τ(ϱ)¯z1) and y∈F(ϱ.τ(ϱ)¯z2). Then x=(ψ(ϱ)−ψ(0))38||z1(ϱ)||x∗ and y=(ψ(ϱ)−ψ(0))38||z2(ϱ)||y∗, where x∗,y∗∈K. It follows that
||x−y||≤(ψ(ϱ)−ψ(0))38||x∗||z1(ϱ)||−||z2(ϱ)||y∗||, |
As a result,
infy∈F(ϱ.τ(ϱ)¯z2)||x−y||≤(ψ(ϱ)−ψ(0))38||x∗||z1(ϱ)||−||z2(ϱ)||x∗||=||x∗||(ψ(ϱ)−ψ(0))38|||z1(ϱ)||−||z2(ϱ)||≤λ((ψ(ϱ)−ψ(0))38||z1(ϱ)−z2(ϱ)||, |
which means that
h(F(ϱ.τ(ϱ)¯z1),F(ϱ.τ(ϱ)¯z2)≤λ((ψ(ϱ)−ψ(0))38||z1(ϱ)−z2(ϱ)||. | (40) |
Then, for ϱ∈(0,ϱ1]
χE(F(ϱ,{τ(ϱ)¯x:x∈D))≤λ(ψ(ϱ)−ψ(0))1−μς(ϱ)χE{x(ϱ):x∈D}≤λ(ψ(ϱ)−ψ(0))1−μς(ϱ)supθ∈[−r,0]χE{(τ(ϱ))¯x(θ):x∈D}, |
Similarly, one can show that if for ϱ∈(s2,ϱ3]
χE(F(ϱ,D))≤λ(ψ(ϱ)−ψ(s2))1−μsupθ∈[−r,0]χE{(τ(ϱ))¯x(θ):x∈D} |
and consequently, by choosing λ small enough such that the relation (7) becomes realized with ς(ϱ)=λ;ϱ∈J. Next, let g1:[1,2]×E→E such that
g(ϱ,x)=ρ(ψ(1)−ψ(0))38ϱ(x), | (41) |
where, ρ is a positive number andϱ is a linear, bounded and compact operator on E. So, condition (H1) is satisfied. As a consequence, from Theorem (1), the problem (1) has a solution where F and g are given by (47) and (50) and Ψ∈Θ provided that
3λκ1Γ(12)η+ρ+ρΓ(58)<1, | (42) |
where, η=δq(ψ(b)−ψ(0))ϑ−q(ϑ−q1−q)1−q, δ=maxs∈Jψ′(s), q=14 and κ1=(ψ(b)−ψ(0))1−μ. By choosing ρ small enough relation (42) becomes realized.
EXAMPLE 2 Let E, K, r, s0, ϱ1, s1, ϱ2, ϑ, ν, μ, F:J×Θ→Pck(E) and g1:[1,2]×E→E be as in EXAMPLE 1. Then assumption (F4) is satisfied. Moreover, in view of (40) for any z1,z1∈H
h(F(ϱ.τ(ϱ)¯z1),F(ϱ.τ(ϱ)¯z2)≤λ((ψ(ϱ)−ψ(0))38||z1(ϱ)−z2(ϱ)||≤λ||z1−z2||H=σ(ϱ)ς(||z1−z2||H), |
where, σ(ϱ)=2λ and ς(ϱ)=ϱ2;ϱ∈J. Observe that ς is strictly comparison and ς(cϱ)=c ς(ϱ), for every c>0, every ϱ∈[0,∞) and ϱ=0 is a point of uniform convergence for the series ∞∑n=1ςn(ϱ), and hence assumption (F5)is satisfied. Now, for any ϱ∈[1,2] and any x,y∈E
||g1(ϱ,x)−g1(ϱ,y)||≤ρ(ψ(1)−ψ(0))38||ϱ(x)−ϱ(y)||≤ρ||ϱ||(ψ(1)−ψ(0))38||x−y||E=ξς(||x−y||E), |
where ξ=2ρ||ϱ||(ψ(1)−ψ(0))38. It follows that (H2) is satisfied. By applying Theorem 2, Problem (1) has a solution and it is ψ−generalized Ulam-Hyers stable, where F and g are given by (37) and (40) and Ψ∈Θ provided that
2ρ||ϱ||(ψ(1)−ψ(0))38Γ(μ)+2ρ||ϱ||(ψ(1)−ψ(0))38Γ(μ)+η||σ||L1q(I,R+)<1. |
By choosing ρ small enough this inequality becomes realized.
A relation between a solution of the considered problem and the corresponding fractional integral equation is given, then two existence results of solutions for an ψ-Hilfer fractional non-instantaneous impulsive differential inclusion in the presence of delay in an infinite dimensional Banah spaces are obtained. Moreover, by using the multivalued weakly Picard operator theory, the stability of solutions for the considered problem in the sense of generalized Ulam-Hyers is studied. This work generalizes many recent results in the literature, for example [6,10,11,17,19,31]. Moreover, our technique can be used to study the existence and Ulam-Hyers stability of solutions or mild solutions for the problems considered in [3,4,15,16,20,21,22] to the case when, there are impulses and delay on the system, the right hand side is a multi-valued function and involving ψ−Hilfer fractional derivative. There are many directions for future work, for example: Generalize the obtained results in [3,4,15,16] when, the considered problems in these works involving ψ-Hilfer fractional derivative.
The authors acknowledge the Deanship of Scientific Research at King Faisal University for the financial support under Nasher Track (Grant No.186042).
The authors declare that they have no conflicts of interest.
All authors contributed equality and read and approved the final version of the manuscript.
[1] |
T. Buckmaster, T. D. Drivas, S. Shkoller, V. Vicol, Simultaneous development of shocks and cusps for 2D Euler with azimuthal symmetry from smooth data, Ann. PDE, 8 (2022). https://doi.org/10.1007/s40818-022-00141-6 doi: 10.1007/s40818-022-00141-6
![]() |
[2] |
T. Buckmaster, S. Shkoller, V. Vicol, Formation of shocks for 2D isentropic compressible Euler, Commun. Pur. Appl. Math., 75 (2022), 2069–2120. https://doi.org/10.1002/cpa.21956 doi: 10.1002/cpa.21956
![]() |
[3] |
P. D. Lax, Development of singularities of solutions of nonlinear hyperbolic partial differential equations, J. Math. Phys., 5 (1964), 611–613. https://doi.org/10.1063/1.1704154 doi: 10.1063/1.1704154
![]() |
[4] |
T. P. Liu, Development of singularities in the nonlinear waves for quasi-linear hyperbolic partial differential equations, J. Differ. Equations, 33 (1979), 92–111. https://doi.org/10.1016/0022-0396(79)90082-2 doi: 10.1016/0022-0396(79)90082-2
![]() |
[5] |
T. C. Sideris, Formation of singularities in three-dimensional compressible fluids, Commun. Math. Phys., 101 (1985), 475–485. https://doi.org/10.1007/BF01210741 doi: 10.1007/BF01210741
![]() |
[6] |
G. Chen, R. Pan, S. Zhu, Singularity formation for the compressible Euler equations, SIAM J. Math. Anal., 49 (2017), 2591–2614. https://doi.org/10.1137/16M1062818 doi: 10.1137/16M1062818
![]() |
[7] |
G. Chen, Optimal density lower bound on nonisentropic gas dynamics, J. Differ. Equations, 268 (2020), 4017–4028. https://doi.org/10.1016/j.jde.2019.10.017 doi: 10.1016/j.jde.2019.10.017
![]() |
[8] | C. Dafermos, Hyperbolic Conservation Laws in Continuum Physics, Springer, 2005. https://doi.org/10.1007/978-3-662-49451-6 |
[9] | G. Q. Chen, D. Wang, The Cauchy problem for the Euler equations for compressible fluids, in Handbook of mathematical fluid dynamics, Vol. I (eds. S. Friedlander, and D. Serre), North-Holland, Amsterdam (2002), 421–543. https://doi.org/10.1016/s1874-5792(02)80012-x |
[10] | T. P. Liu, Shock Waves, American Mathematical Society, Providence, RI, 2021. |
[11] |
J. Rauch, BV estimates fail for most quasilinear hyperbolic systems in dimensions greater than one, Commun. Math. Phys., 106 (1986), 481–484. https://doi.org/10.1007/BF01207258 doi: 10.1007/BF01207258
![]() |
[12] | M. P. Lebaud, Description de la formation d'un choc dans le p-système, Journal de mathématiques pures et appliquées, 73 (1994), 523–566. |
[13] |
S. Chen, L. Dong, Formation and construction of shock for p-system, Sci. China Ser. A-Math., 44 (2001), 1139–1147. https://doi.org/10.1007/BF02877431 doi: 10.1007/BF02877431
![]() |
[14] |
D. X. Kong, Formation and propagation of singularities for 2×2 quasilinear hyperbolic systems, Trans. Amer. Math. Soc., 354 (2002), 3155–3179. https://doi.org/10.1090/S0002-9947-02-02982-3 doi: 10.1090/S0002-9947-02-02982-3
![]() |
[15] | T. Buckmaster, T. D. Drivas, S. Shkoller, V. Vicol, Formation and development of singularities for the compressible Euler equations, in Proceedings of the International Congress of Mathematicians 2022, EMS Press, Berlin, 5 (2023), 3636–3659. https://doi.org/10.4171/ICM2022/87 |
[16] | D. Christodoulou, The formation of shocks in 3-dimensional fluids, European Mathematical Society, 2007. https://doi.org/10.4171/031 |
[17] | D. Christodoulou, The shock development problem, European Mathematical Society, 2019. https://doi.org/10.4171/192 |
[18] |
H. Yin, Formation and construction of a shock wave for 3-d compressible Euler equations with the spherical initial data, Nagoya Math. J., 175 (2004), 125–164. https://doi.org/10.1017/S002776300000893X doi: 10.1017/S002776300000893X
![]() |
[19] |
D. Christodoulou, A. Lisibach, Shock development in spherical symmetry, Ann. PDE, 2 (2016). https://doi.org/10.1007/s40818-016-0009-1 doi: 10.1007/s40818-016-0009-1
![]() |
[20] |
J. Luk, J. Speck, Shock formation in solutions to the 2d compressible Euler equations in the presence of non-zero vorticity, Invent. Mathe., 214 (2018), 1–169. https://doi.org/10.1007/s00222-018-0799-8 doi: 10.1007/s00222-018-0799-8
![]() |
[21] |
J. Luk, J. Speck, The stability of simple plane-symmetric shock formation for 3d compressible Euler flow with vorticity and entropy, Anal PDE, 17 (2024), 831–941. https://doi.org/10.2140/apde.2024.17.831 doi: 10.2140/apde.2024.17.831
![]() |
[22] |
T. Buckmaster, S. Shkoller, V. Vicol, Formation of point shocks for 3D compressible Euler, Comm. Pur. Appl. Math., 76 (2023), 2073–2191. https://doi.org/10.1002/cpa.22068 doi: 10.1002/cpa.22068
![]() |
[23] |
T. Buckmaster, S. Shkoller, V. Vicol, Shock formation and vorticity creation for 3D Euler, Comm. Pur. Appl. Math., 76 (2023), 1965–2072. https://doi.org/10.1002/cpa.22067 doi: 10.1002/cpa.22067
![]() |
[24] |
S. Shkoller, V. Vicol, The geometry of maximal development and shock formation for the Euler equations in multiple space dimensions, Invent. Math., 237 (2024), 871–1252. https://doi.org/10.1007/s00222-024-01269-x doi: 10.1007/s00222-024-01269-x
![]() |
[25] | L. Abbrescia, J. Speck, The emergence of the singular boundary from the crease in 3d compressible Euler flow, preprint, arXiv: 2207.07107. |
[26] | L. D. Landau, E. M. Lifshitz, Fluid Mechanics, Vol. 6, Elsevier, 1987. https://doi.org/10.1016/C2013-0-03799-1 |
[27] | R. J. Walker, Algebraic Curves, Springer-Verlag, 1950. |
[28] | E. Brieskorn, H. Knörrer, Plane algebraic curves, Springer Basel, 1986. https://doi.org/10.1007/978-3-0348-5097-1 |
[29] |
S. Kumagai, An implicit function theorem: comment, J Optim Theory Appl, 31 (1980), 285–288. https://doi.org/10.1007/BF00934117 doi: 10.1007/BF00934117
![]() |
1. | Thabet Abdeljawad, Pshtiwan Othman Mohammed, Hari Mohan Srivastava, Eman Al-Sarairah, Artion Kashuri, Kamsing Nonlaopon, Some novel existence and uniqueness results for the Hilfer fractional integro-differential equations with non-instantaneous impulsive multi-point boundary conditions and their application, 2023, 8, 2473-6988, 3469, 10.3934/math.2023177 | |
2. | Zainab Alsheekhhussain, Ahmed Gamal Ibrahim, Rabie A. Ramadan, Existence of $ S $-asymptotically $ \omega $-periodic solutions for non-instantaneous impulsive semilinear differential equations and inclusions of fractional order $ 1 < \alpha < 2 $, 2023, 8, 2473-6988, 76, 10.3934/math.2023004 | |
3. | Ramasamy Arul, Panjayan Karthikeyan, Kulandhaivel Karthikeyan, Palanisamy Geetha, Ymnah Alruwaily, Lamya Almaghamsi, El-sayed El-hady, On Ψ-Hilfer Fractional Integro-Differential Equations with Non-Instantaneous Impulsive Conditions, 2022, 6, 2504-3110, 732, 10.3390/fractalfract6120732 | |
4. | M. Latha Maheswari, K. S. Keerthana Shri, Karthik Muthusamy, Existence results for coupled sequential ψ-Hilfer fractional impulsive BVPs: topological degree theory approach, 2024, 2024, 1687-2770, 10.1186/s13661-024-01901-y | |
5. | Zainab Alsheekhhussain, Ahmad Gamal Ibrahim, Mohammed Mossa Al-Sawalha, Yousef Jawarneh, The Existence of Solutions for w-Weighted ψ-Hilfer Fractional Differential Inclusions of Order μ ∈ (1, 2) with Non-Instantaneous Impulses in Banach Spaces, 2024, 8, 2504-3110, 144, 10.3390/fractalfract8030144 | |
6. | Feryal Aladsani, Ahmed Gamal Ibrahim, Existence and Stability of Solutions for p-Proportional ω-Weighted κ-Hilfer Fractional Differential Inclusions in the Presence of Non-Instantaneous Impulses in Banach Spaces, 2024, 8, 2504-3110, 475, 10.3390/fractalfract8080475 | |
7. | Zainab Alsheekhhussain, Ahmed Gamal Ibrahim, M. Mossa Al-Sawalha, Khudhayr A. Rashedi, Mild Solutions for w-Weighted, Φ-Hilfer, Non-Instantaneous, Impulsive, w-Weighted, Fractional, Semilinear Differential Inclusions of Order μ ∈ (1, 2) in Banach Spaces, 2024, 8, 2504-3110, 289, 10.3390/fractalfract8050289 | |
8. | F. M. Ismaael, On a New Class of Impulsive η-Hilfer Fractional Volterra-Fredholm Integro-Differential Equations, 2023, 17, 1823-8343, 691, 10.47836/mjms.17.4.10 |