Processing math: 62%
Research article Special Issues

Hybrid ideals in near-subtraction semigroups

  • Received: 05 February 2022 Revised: 11 April 2022 Accepted: 21 April 2022 Published: 20 May 2022
  • MSC : 06D72, 20M12

  • The fuzzy set is highly beneficial for expressing people's hesitations in their everyday lives, and it is a great tool for dealing with uncertainty, which can be described precisely and perfectly from the decision-maker's point of view. Soft set theory has been developed in recent years to address real-world issues. Jun et al. merged fuzzy and soft sets to produce hybrid structures. Hybrid structures are soft set and fuzzy set speculations. The concept of hybrid ideals in near-subtraction semigroups is introduced in this paper, and their equivalent results are obtained. Additionally, we demonstrate the concept of hybrid intersection. Moreover, we define the concept of homomorphism of a hybrid structure in a near-subtraction semigroup.

    Citation: S. Meenakshi, G. Muhiuddin, B. Elavarasan, D. Al-Kadi. Hybrid ideals in near-subtraction semigroups[J]. AIMS Mathematics, 2022, 7(7): 13493-13507. doi: 10.3934/math.2022746

    Related Papers:

    [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 ψ-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
  • The fuzzy set is highly beneficial for expressing people's hesitations in their everyday lives, and it is a great tool for dealing with uncertainty, which can be described precisely and perfectly from the decision-maker's point of view. Soft set theory has been developed in recent years to address real-world issues. Jun et al. merged fuzzy and soft sets to produce hybrid structures. Hybrid structures are soft set and fuzzy set speculations. The concept of hybrid ideals in near-subtraction semigroups is introduced in this paper, and their equivalent results are obtained. Additionally, we demonstrate the concept of hybrid intersection. Moreover, we define the concept of homomorphism of a hybrid structure in a near-subtraction semigroup.



    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 s0 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×E2E{ϕ} is a multifunction, and gi:[ϱi,si]×EEi=1,2,,n. Finally, for any ϱJ,  τ(ϱ):HΘ,  x(θ)=x(ϱ+θ),θ[r,0]xH, 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):={xC(J,E):(ψ(.)ψ(a))1μx(.)C([a,b],E)},

    and

    Cn1μ,ψ(J,E):={xCn1(J,E):[1ψ(ϱ)ddϱ]n xC1μ,ψ(J,E),nN.

    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=n1k=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:JE, (ψ(.)ψ(sk))1μx(.)C(Jk,E),k=0,1,..,n,limϱs+k(ψ(ϱ)ψ(sk))1μx(ϱ)exists,xC(ϱ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:={hC(¯Jk,E):h(ϱ)=(ψ(ϱ)ψ(sk))1μh(ϱ),ϱJk,h(sk)=limϱs+kh(ϱ),hD},

    and

    D¯ϱi:={hC(¯ϱi,E):h(ϱ)=h(ϱ),ϱϱi,h(ϱi)=h(ϱ+i),hD}.

    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:JE,[1ψ(ϱ)ddϱ] xAC(J,E)}.

    If f C1,ψ([a,b],E):={xC(J,E):[1ψ(ϱ)ddϱ] xC(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,fCη,ψ([a,b],E).

    Proof. The assumption fCη,ψ([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 fL1([a,b],E) be such that I(1ν)(1ϑ),ψa+fAC1,ψ([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):={xC1μ,ψ(J,E),Dμ,ψa+xC1μ,ψ(J,E)},
    Cϑ,ν1μ,,ψ(J,E):={xC1μ(J,E),Dϑ,ν,ψa+xC1μ,ψ(J,E)},
    PCμ1μ,ψ(J,E):={xPC1μ(J,E),Dμ,ψs+kx|JkC1μ,ψ(Jk,E),k=0,1,..,n},

    and

    PCϑ,ν1μ,ψ(J,E):={xPC1μ,ψ(J,E),Dϑ,ν,ψs+kx|JkC1μ,ψ(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 fCμ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 hasafinitenumberofdiscontinuitypointss0,allvaluesz(s+),and z(s) arefinite},

    endowed with the norm:

    ||z||Θ=0r||z(s)||ds.

    2-The metric space (the space of solutions)

    H:={x:[r,b]E,x|[r,0]=Ψ, x|JPC1μ,ψ(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|JPC1μ,ψ(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 xH, then x(0)=0 and x(0+)=limϱ0+(ψ(ϱ)ψ(0))1μx(ϱ).

    (ⅱ) If xH, 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:xB},

    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 ϱ:WPc(W). Suppose that ϱ is closed, ϱ(D) is relatively compact, whenever D is compact, and that, for some x0W, one has

    BW, B=conv({x0}ϱ(B)), ¯B=¯C withCB countableB  isrelativelycompact.

    Then, there is a fixed point for ϱ.

    In this section, we demonstrate the existence of solutions of Problem (1). For any xH let

    S1F(.,τ(.)x)={zL1(J,E):z(ϱ)F(ϱ,τ(.)x),a.e. for ϱJ }.

    and

    Iν(1ϑ),ψ(PC1μ,ψ(J,E))={f:JE, thereishPC1μ,ψ(J,E)suchthatf(ϱ)=Iν(1ϑ),ψsi+h(ϱ),ϱJk, k=0,1,2,..,n}.

    and for any xH 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 fCη,ψ([a,b],E) and I1ϑ,ψa+fC1η,ψ([a,b],E), then

    Iϑ,ψa+Dϑ,ψa+f(ϱ)=f(ϱ)(ψ(ϱ)ψ(a))ϑ1I1ϑ,ψa+f(a)Γ(ϑ),ϱ[a,b].

    Proof. Since fCη,ψ([a,b],E), then I1ϑ,ψa+f(ϱ) is defined for ϱ[a,b]. Moreover, the assumption I1ϑ,ψa+fC1η,ψ([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 fCη,ψ([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 xPCν(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 yPCμ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 xPCν(1ϑ),ψ1μ,ψ([0,b],E) implies to Dν(1ϑ),ψ0+x|J0C1μ,ψ([0,ϱ1],E). It follows from (5) that Dμ,ψa+yC1μ,ψ([0,ϱ1],E) and hence y|J0Cμ1μ,ψ([0,ϱ1],E). Consequently y|J0Cϑ,ν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|J0C1μ,ψ([0,ϱ1],E), then (1ψ(ϱ)ddϱ)I1ν(1ϑ),ψ0+x|J0C1μ,ψ([0,ϱ1],E). So, I1ν(1ϑ),ψ0+x|J0C11μ,ψ([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+xC1μ,ψ([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 xC1μ,ψ([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 yPCμ1μ,ψ(J,E) be a solution of (3). Let i=0. Then yC1μ,ψ(J0,E) and (1ψ(ϱ)ddϱ)I1μ,ψ0+y|J0C1μ,ψ(J0,E). Thus, I1μ,ψsiy|J0C11μ,ψ(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 ¯xH is called a mild solution of (1) if there is fS1F(,τ(ϱ)¯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]×EE (i=1,2,....n). We assume the following conditions

    (F1) For every zH, S1F(.,τ(.)z) is not empty subset of Iν(1ϑ),ψ(PC1μ,ψ(J,E)) and for almost every ϱJ,zF(ϱ,τ(ϱ)z) is upper semicontinuous.

    (F2) There is a φL1q(I,R+),(0<q<ϑ) such that for any zH

    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 DH,  any k=0,1,2,..,n, and a.e., for ϱJk

     χE(F(ϱ,{τ(ϱ)¯x:xD))(ψ(ϱ)ψ(sk))1μς(ϱ)supθ[r,0]χE{(τ(ϱ))¯x(θ):xD},

    and

    2κ1Γ(ϑ) η||ς||L1q(I,R+)<1, (7)

    where η=δq(ψ(b)ψ(0))ϑq(ϑq1q)1q, δ=maxsJψ(s), κ1=(ψ(b)ψ(0))1μ and χ is the Hausdorff measure of noncompactness on E.   

    (H1) for every i=1,2,,n,   gi:[ϱi,si]×EE 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 xE

    gi(ϱ,x))hi(ψ(ϱi)ψ(si1))1μx, ϱ(ϱi,si],xE.

    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 Φ:HP(H) as follows: let xH, 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 nN with Φ(Bn)Bn,where Bn={xH:xHn}. Suppose that for any natural number n, there are xn,ynH with ynΦ(xn),xnHn and ynH>n. Then, there is fn S1F(.,τ(.)¯xn), n1 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+xnH)φ(ϱ)(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))ϑ11q(ψ(s))11qds]1q||φ||L1q(J,R+)[ϱ0(ψ(ϱ)ψ(s))ϑ11qψ(s)(ψ(s))q1qds]1qδq||φ||L1q(J,R+)[ϱ0(ψ(ϱ)ψ(s))ϑ11qψ(s)ds]1qδq(ψ(b)ψ(0))ϑq(ϑq1q)1q=η. (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) ψ(si1))1μ||xn(ϱi)||h||xn||Hhn. (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)ηhnΓ(μ)+  κ1Γ(ϑ)||φ||L1q(J,R+)(1+n)η.

    It follows from this inequality, (13) and (14), that

    n<||yn||HΨ(0)Γ(μ)+hn+hnΓ(μ)+κ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)n1, (yn)n1be two sequences in Bn0 with xnx in Bn0,yny in H and ynΦ(xn);n1. 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)n1 is bounded in L1q(J,E). Since L1q(J,E) is reflexive, then, by using Mazur's Lemma, there is sequence, (zj)j1, of convex combinations of (fn)n1  converging strongly to f in L1(J,E) as j.  Notice that, by (F2) again, for every ϱJ,s(0,ϱ] and every n1

    ||(ϱ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:¯JkE, z(ϱ)=(ψ(ϱ)ψ(sk))1μy(ϱ),ϱJk,z(sk)=limϱsk(ψ(ϱ)ψ(sk))1μy(ϱ), yΦ(x),xBn0},

    and

    M¯ϱi={yC(¯ϱi,E):y(ϱ)=y(ϱ),ϱ(ϱi,si],y(ϱi)=y(ϱ+i),yΦ(x),xBn0}.

    Case 1.let zM¯J0.Then there are xBn0and 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)]11qds]1q

    Put ¯ω=ϑ11q(1,0), then for s<ϱ1<ϱ2, we have (ψ(ϱ1)ψ(s))¯ω(ψ(ϱ2)ψ(s))¯ω. By applying Lemma 3.1 in [43] and taking into account 1q(0,1), we get

    |[(ψ(ϱ1)ψ(s))¯ω]1q[(ψ(ϱ2)ψ(s))¯ω]1q| [(ψ(ϱ1)ψ(s))¯ω(ψ(ϱ2)ψ(s))¯ω]1q.

    Then,

    |(ψ(ϱ1)ψ(s))ϑ1(ψ(ϱ2)ψ(s))ϑ1| [(ψ(ϱ1)ψ(s))¯ω(ψ(ϱ2)ψ(s))¯ω]1q.

    This leads to

    |(ψ(ϱ1)ψ(s))ϑ1(ψ(ϱ2)ψ(s))ϑ1|11q(ψ(ϱ1)ψ(s))¯ω(ψ(ϱ2)ψ(s))¯ω.

    Therefore,

    limϱ2ϱ1[ϱ10[(ψ(ϱ2)ψ(s))ϑ1ψ(s)(ψ(ϱ1)ψ(s))ϑ1ψ(s)]11qds]1qlimϱ2ϱ1[ϱ10[(ψ(ϱ2)ψ(s))ϑ1(ψ(ϱ1)ψ(s))ϑ1]11q(ψ(s))11qds]1qlimϱ2ϱ1[ϱ10[(ψ(ϱ1)ψ(s))¯ω(ψ(ϱ2)ψ(s))¯ω](ψ(s))11qds]1qδqlimϱ2ϱ1[ϱ10[(ψ(ϱ1)ψ(s))¯ω(ψ(ϱ2)ψ(s))¯ω]ψ(s)ds]1q=δqlimϱ2ϱ1[(ψ(ϱ1)ψ(s))¯ω+1¯ω+1(ψ(ϱ2)ψ(s))¯ω+1¯ω+1]1q=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):xBn0}.ObviouslyKi is bounded subset of E.

    Let yMϱi. Then, there is xBn0 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λϱ+iy(ϱi+δ)y(λ)=0.

    Case 3. let k=1,,n be fixed, zM¯Jk. Then, there are xBn0and 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ϱ1z(ϱ+δ)z(ϱ)=0,

    independent of x.

    Step 4. Let  KBn0, K=conv({0}Φ(K)), ¯Z=¯C with CK countable. We have to show that K is relatively compact in H Let D={yn:n1}Φ(K) with Cconv({x0}D), xnK with ynΦ(xn).This means that, there is fnS1F(.,τ(.)¯xn) such that the relation (15) holds. Observe that, from (F3) it holds for a.e.sJ0

    χE{fn(s):n1}(ψ(ϱ)ψ(0))1μς(ϱ)χ{F(s,τ(s)¯xn):n1}ζ(s)(ψ(ϱ)ψ(0))1μsupθ[r,0] χ{¯xn(s+θ):n1}ζ(s)(ψ(ϱ)ψ(0))1μsupδ[r,s] χ{¯xn(δ):n1}.

    Because supδ[r,0] χ{¯xn(δ):n1}=supδ[r,0] χ{Ψ(δ)}=0. Thus,

    χE{fn(s):n1}ζ(s)(ψ(ϱ)ψ(0))1μ supδ[0,s]χ{xn(δ):n1}=ζ(s)supδ[0,s]χ{(ψ(ϱ)ψ(0))1μxn(δ):n1}ζ(s)χPC1μ,ψ(J,E){xn:n1}ζ(s)χH(Bn1). (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χ{yn(ϱ):n1},maxi=1,..,nmaxϱ¯ϱiχ{yn(ϱ):n1}}, (19)

    Let ϱ(0,ϱ1]. Then,

    yn(ϱ)=Ψ(0)Γ(μ)+(ψ(ϱ)ψ(0))1μΓ(ϑ)ϱ0(ψ(ϱ)ψ(0))ϑ1ψ(s)fn(s)ds,

    which yields with (17)

    χ{yn(ϱ):n1}(ψ(ϱ)ψ(0))1μΓ(ϑ)ϱ0(ψ(ϱ)ψ(0))ϑ1ψ(s)χ{fn(s):n1}ds(ψ(b)ψ(0))1μΓ(ϑ)χH(Bn1)ϱ0(ψ(ϱ)ψ(0))ϑ1ψ(s)ζ(s)dsκ1Γ(ϑ) η||ς||L1q(I,R+)χH(Bn1). (20)

    Notice that

    χ{yn(0):n1}=χ{limϱ0+yn(ϱ):n1}=χ{Ψ(0)Γ(μ)}=0. (21)

    Moreover, from the fact that ||xn||Hn0,the set {xn(ϱi):n1} is bounded for every i=1,2,...,n and hence from the assumption (H1), we get

    χ{gi(ϱ,xn(ϱi)),n1}=0,ϱ(ϱi,si],i=1,...,n. (22)

    and

    χ{gi(ϱi,xn(ϱi)),n1}=χ{gi(si,xn(ϱi)):n1}=0. (23)

    Similarly,

    χ{yn(si):n1}=χ{limϱs+i(ψ(ϱ)ψ(si))1μyn(ϱ),:n1}=χ{1Γ(μ)gi(si,xn(ϱi)):n1}=0,i=1,2,..,n. (24)

    By arguing as in (20), one can show that

    maxϱ¯Jkχ{yn(ϱ):n1}κ1Γ(ϑ) η||ς||L1q(I,R+)χH(Bn1),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. \Phi \ maps compact sets into relatively compact sets.

    Let B\ be a compact subset of B_{n_{0}} , \{y_{n}, :n\geq 1\}\subseteq \Phi (B)\ Then, there is x_{n}\in B, n\geq 1 , \ such that y_{n} \in \Phi (x_{n}). So, there is f_{n}\in S_{F(., \overline{x_{n}}(.))}^{1}\ such that (15) holds. We have to show that the set Z = \{y_{n}:n\geq 1\}\ is relatively compact in \mathcal{H}.\ Note that, since B\ is compact in \mathcal{H}, \ then from (F_{3})\ we get for a.e. s\in J_{0} ,

    \begin{eqnarray*} \chi _{E}\{f_{n}(s) &:&n\geq 1\}\leq (\psi (\mathcal{\varrho })-\psi (0))^{1-\mu }\varsigma (\mathcal{\varrho })\chi \{F(s,\tau (s)\overline{x} _{n}):n\geq 1\} \\ &\leq &\zeta (s)(\psi (\mathcal{\varrho })-\psi (0))^{1-\mu }\sup\limits_{\theta \in \lbrack -r,0]}\ \chi \{\overline{x}_{n}(s+\theta ):n\geq 1\} \\ &\leq &\zeta (s)(\psi (\mathcal{\varrho })-\psi (0))^{1-\mu }\sup\limits_{\delta \in \lbrack -r,s]}\ \chi \{\overline{x}_{n}(\delta ):n\geq 1\} \\ &\leq &\zeta (s)(\psi (\mathcal{\varrho })-\psi (0))^{1-\mu }\ \sup\limits_{\delta \in \lbrack 0,s]}\chi \{x_{n}(\delta ):n\geq 1\} \\ & = &\zeta (s)\sup\limits_{\delta \in \lbrack 0,s]}\chi \{(\psi (\mathcal{\varrho } )-\psi (0))^{1-\mu }x_{n}(\delta ):n\geq 1\} \\ &\leq &\zeta (s)\chi _{PC_{1-\mu ,\psi }(J,E)}\{x_{n}:n\geq 1\} \\ &\leq &\zeta (s)\chi _{\mathcal{H}}(B) = 0. \end{eqnarray*}

    By the same reasons, one can show that for a.e. s\in J_{k} , k = 1, 2, .., n\

    \begin{equation*} \chi _{E}\{f_{n}(s):n\geq 1\} = 0. \end{equation*}

    By arguing as in the previous step one can show that Z is relatively compact, and hence \Phi (B)\ is relatively compact.

    Now, by applying Lemma 8, there is x\in \mathcal{H} and f \in S_{F(., \tau (.)\overline{x})}^{1}\ such that

    \begin{equation*} x(\mathcal{\varrho }) = \left\{ \begin{array}{l} 0,\mathcal{\varrho }\in \lbrack -r,0] \\ \frac{(\psi (\mathcal{\varrho })-\psi (0))^{\mu -1}}{\Gamma (\mu )}\Psi (0) \\ +\frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{ \varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)f(s)ds,\mathcal{\varrho } \in (0,\mathcal{\varrho }_{1}] \\ g_{i}(\mathcal{\varrho },x(\mathcal{\varrho }_{i}^{-})),\mathcal{\varrho } \in (\mathcal{\varrho }_{i},s_{i}],i = 1,......n, \\ \frac{(\psi (\mathcal{\varrho })-\psi (0))^{\mu -1}}{\Gamma (\mu )} g_{i}(s_{i},x(\mathcal{\varrho }_{i}^{-})) \\ +\frac{1}{\Gamma (\vartheta )}\int_{s_{i}}^{\mathcal{\varrho }}(\psi ( \mathcal{\varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)f(s),\mathcal{ \varrho }\in (s_{i}{\rm{ }}\mathcal{\varrho }_{i+1}],i = 1,......n. \end{array} \right. \end{equation*}

    Next, in view of (F_{1}) , there is h\in PC_{1-\mu, \psi }(J, E) such that f(\mathcal{\varrho }) = I_{s_{i}+}^{\nu (1-\vartheta), \psi }h(\mathcal{ \varrho }), \mathcal{\varrho }\in J_{k}, k = 0, .., n, and hence, from Lemma 3.3,

    \begin{equation*} D_{s_{i}+}^{\nu (1-\vartheta ),\psi }f(\mathcal{\varrho }) = D_{s_{i}+}^{\nu (1-\vartheta ),\psi }I_{s_{i}+}^{\nu (1-\vartheta )\psi }h(\mathcal{\varrho } ) = h(\mathcal{\varrho }),\mathcal{\varrho }\in J_{k},k = 0,..,n. \end{equation*}

    This yields that f\in PC_{1-\mu, \psi }^{\nu (1-\vartheta)}(J, E). Then, the function

    \begin{equation} \overline{x}(\mathcal{\varrho }) = \left\{ \begin{array}{l} \Psi (\mathcal{\varrho }),\mathcal{\varrho }\in \lbrack -r,0] \\ \frac{(\psi (\mathcal{\varrho })-\psi (0))^{\mu -1}}{\Gamma (\mu )}\Psi (0)+ \frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{ \varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)f(s)ds,\mathcal{\varrho } \in (0,\mathcal{\varrho }_{1}] \\ g_{i}(\mathcal{\varrho },x(\mathcal{\varrho }_{i}^{-})),\mathcal{\varrho } \in (\mathcal{\varrho }_{i},s_{i}],i = 1,......n, \\ \frac{(\psi (\mathcal{\varrho })-\psi (s_{i}))^{\mu -1}}{\Gamma (\mu )} g_{i}(s_{i},x(\mathcal{\varrho }_{i}^{-})) \\ +\frac{1}{\Gamma (\vartheta )}\int_{s_{i}}^{\mathcal{\varrho }}(\psi ( \mathcal{\varrho })-\psi (s))^{\vartheta -1}\psi \prime (s),\mathcal{\varrho }\in (s_{i}{\rm{ }}\mathcal{\varrho }_{i+1}],i = 1,......n. \end{array} \right. \end{equation} (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 (F_{2})\ is replaced by the following assumption:

    (F_{2})^{\ast } There is a \varphi \in L^{\frac{1}{q}}(I, \mathbb{R} ^{+}) such that for any z\in \Theta

    \begin{equation*} \left\Vert F(\mathcal{\varrho },z)\right\Vert \leq \left\{ \begin{array}{l} \varphi (\mathcal{\varrho })\ (1+(\psi (\mathcal{\varrho })-\psi (s_{i}))^{1-\mu }\left\Vert z(0)\right\Vert _{E}),a.e.\mathcal{\varrho }\in \cup _{i = 0}^{i = n}(s_{i}{\rm{ }}\mathcal{\varrho }_{i+1}], \\ \varphi (\mathcal{\varrho })\ (1+\left\Vert z(0)\right\Vert _{E}),a.e. \mathcal{\varrho }\in \cup _{i = 1}^{i = n}(\mathcal{\varrho }_{i},s_{i}]. \end{array} \right. \end{equation*}

    In fact, condition (F_{2})\ is used only to prove relation (11). We show this relation by using (F_{2})^{\ast } . Let f_{n}\in S_{F(., \tau (.) \overline{x}_{n})}^{1} , n\geq 1 Then, by (F_{2})^{\ast }, \ for almost \mathcal{\varrho }\in \cup _{i = 0}^{i = n}(s_{i} \mathcal{\varrho }_{i+1}] we get \

    \begin{eqnarray*} ||f_{n}(\mathcal{\varrho })|| &\leq &\varphi (\mathcal{\varrho })(1+(\psi ( \mathcal{\varrho })-\psi (s_{i}))^{1-\mu }||\tau (\mathcal{\varrho }) \overline{x}_{n}(0)||) \\ &\leq &\varphi (\mathcal{\varrho })(1+(\psi (\mathcal{\varrho })-\psi (s_{i}))^{1-\mu }||\overline{x}_{n}(\mathcal{\varrho })||) \\ &\leq &\varphi (\mathcal{\varrho })(1+(\psi (\mathcal{\varrho })-\psi (s_{i}))^{1-\mu }||x_{n}(\mathcal{\varrho })|| \\ &\leq &\varphi (\mathcal{\varrho })(1+\left\Vert x_{n}\right\Vert _{\mathcal{ H}}). \end{eqnarray*}

    and for almost \mathcal{\varrho }\in \cup _{i = 1}^{i = n}(\mathcal{\varrho } _{i} s_{i}]

    \begin{eqnarray*} ||f_{n}(\mathcal{\varrho })|| &\leq &\varphi (\mathcal{\varrho })(1+||\tau ( \mathcal{\varrho })\overline{x}_{n}(0)||_{E}) \\ &\leq &\varphi (\mathcal{\varrho })(1+\overline{x}_{n}(\mathcal{\varrho })) \\ &\leq &\varphi (\mathcal{\varrho })(1+||x_{n}(\mathcal{\varrho })||) \\ &\leq &\varphi (\mathcal{\varrho })(1+\left\Vert x_{n}\right\Vert _{\mathcal{ H}}). \end{eqnarray*}

    In this section, we give another version for the existence of solutions and investigate the generalized \ \psi- 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 \varsigma :[0, \infty)\rightarrow \lbrack 0, \infty)\ is called comparison function if \lim_{n\rightarrow \infty }\varsigma ^{n}(s) = 0 , \ \forall s\in \lbrack 0, \infty), where \varsigma ^{n}(s) = \varsigma ^{n-1}(\varsigma (s)). If in addition \sum\limits_{n = 1}^{\infty }\varsigma ^{n}(\mathcal{\varrho }) < \infty, \forall s\in (0, \infty) , then \varsigma :[0, \infty)\rightarrow \lbrack 0, \infty)\ is called strictly comparison.

    Remark 4.1. [41]

    1- If \varsigma :[0, \infty)\rightarrow \lbrack 0, \infty)\ is comparisons, then \varsigma (s) < s, \forall s > 0, \varsigma (0) = 0 and \varsigma \ is continuous at 0.

    2-the functions \varsigma _{1}(s) = c\mathcal{\varrho }; c\in \lbrack 0, 1)\ and \varsigma _{2}(s) = \frac{s}{s+1}:s\in \lbrack 0, \infty)\ are strictly comparison

    Let M: = \{y\in H:y_{|J}\in PC_{1-\mu, \psi }^{\mu }(J, E) , D_{s_{i}+}^{\vartheta, \nu, \psi }y(\mathcal{\varrho }) exists for any \mathcal{\varrho }\in (s_{i} \mathcal{\varrho }_{i+1}] , i = 0, 1, ......n, \}

    Definition 4. [12] Problem (1) is called \psi - generalized Ulam- Hyers stable if there is a continuous function \theta :[0, \infty)\rightarrow \lbrack 0, \infty) and \theta (0) = 0 such that for each \epsilon > 0 and each \ solution y\in M of the inequality

    \begin{equation} \left\{ \begin{array}{l} d(D_{s_{i}^{+}}^{\vartheta ,\nu ,\psi }y(\mathcal{\varrho }),F(\mathcal{ \varrho },\tau (\mathcal{\varrho })y))\leq \epsilon ,{\rm{ }}a.e.{\rm{ }} \mathcal{\varrho }\in (s_{i},{\rm{ }}\mathcal{\varrho }_{i+1}],i = 0,1,....,n, \\ ||y(\mathcal{\varrho })-g_{i}(\mathcal{\varrho },y(\mathcal{\varrho } _{i}^{-}))||\leq \epsilon ,\mathcal{\varrho }\in (\mathcal{\varrho }_{i} {\rm{ }}s_{i}],i = 1,......n,, \\ ||y(\mathcal{\varrho }_{i}^{+})-g_{i}(\mathcal{\varrho }_{i},y(\mathcal{ \varrho }_{i}^{-}))||\leq \epsilon ,i = 1,......n, \end{array} \right. \end{equation} (27)

    there is a solution x\in H for (1) with

    \begin{equation} \sup\limits_{\mathcal{\varrho }\in J}||x(\mathcal{\varrho })-y(\mathcal{\varrho } )||\leq \theta (\epsilon ). \end{equation} (28)

    In the following theorem, we establish the existence and generalized \psi- Ulam-Hyers stability of solutions Problem (1).

    Theorem 4.1. Let F:J\times \Theta \rightarrow P_{ck}(E) , \Psi \in \Theta and g_{i}:[\mathcal{\varrho }_{i}, s_{i}]\times E\rightarrow E\ (i = 1, 2, ....n) be such that:

    (F_{4}) For every z\in H , S_{F(., \tau (.)z)}^{1} is a non-empty subset of I^{\nu (1-\vartheta), \psi }(PC_{1-\mu, \psi }(J, E)).

    (F_{5}) There is a function \ \sigma \in L^{\frac{1}{q}}(I, \mathbb{R} ^{+}), 0 < q < \vartheta and a strict comparison function \varsigma :[0, \infty)\rightarrow \lbrack 0, \infty)\ such that

    (ⅰ) For every x_{1}, x_{2}\in \mathcal{H}

    \begin{equation*} h(F(\mathcal{\varrho },\tau (\mathcal{\varrho })\overline{x}_{1}),F(\mathcal{ \varrho },\tau (\mathcal{\varrho })\overline{x}_{2}))\leq \sigma (\mathcal{ \varrho })\varsigma (\left\Vert x_{1}-x_{2}\right\Vert )_{\mathcal{H}},{\rm{ }}\forall \mathcal{\varrho }\in J. \end{equation*}

    (ⅱ)

    \begin{equation*} Sup\{\ ||y||\ :\ y\in F(\mathcal{\varrho },\tau (\mathcal{\varrho }) \overline{x}_{0})\}\leq \sigma (\mathcal{\varrho }),\rm{ for }a.e.{\rm{ }} \mathcal{\varrho }\in J, \end{equation*}

    where \ \overline{x}_{0}(s) = \Psi (s); s\in \lbrack -r, 0] and \overline{x} _{0}(s) = 0;s\in (0, b].

    (H_{2}) For any i = 1, 2, ....n, there \ \xi _{i}\ > 0\ such that for any \mathcal{\varrho }\in \lbrack \mathcal{\varrho }_{i}, s_{i}]

    \begin{equation*} \left\Vert g_{i}(\mathcal{\varrho },x)-g_{i}(\mathcal{\varrho } ,y)\right\Vert \leq \xi _{i}\ \varsigma (\left\Vert x-y\right\Vert _{E}),\forall x,y\in E. \end{equation*}

    Then, there is solution for Problem(1) provided that

    \begin{equation} \frac{\xi }{\Gamma (\mu )}+\ \eta \delta ^{q}\frac{(\psi (b)-\psi (0))^{1-\mu }}{\Gamma (\vartheta )}||\sigma ||_{L^{\frac{1}{q}}(I, \mathbb{R} ^{+})} < 1, \end{equation} (29)

    where \xi = \sum_{i = 1}^{i = n}\xi _{i}. Moreover, if, in addition there, is a c > 1 such that \varsigma (c\mathcal{\varrho })\leq c\ \varsigma (\mathcal{ \varrho }); \mathcal{\varrho }\in \lbrack 0, \infty) and \mathcal{\varrho } = 0 is a point of uniform convergence for the series \sum\limits_{n = 1}^{ \infty }\varsigma ^{n}(\mathcal{\varrho }) , then problem (1) is a \ \psi- generalized Ulam-Hyers stable.

    Proof. Condition (F_{4})\ allows to define a multifunction \Phi :\mathcal{H} \rightarrow P(\mathcal{H})\ as in (9) Note that by (F_{5}), for every x\in \mathcal{H}, and for a.e. \mathcal{\varrho }\in J

    \begin{eqnarray*} ||F(\mathcal{\varrho },\tau (\mathcal{\varrho })\overline{x})|| & = &h(F( \mathcal{\varrho },\tau (\mathcal{\varrho })\overline{x}),\{0\}) \\ &\leq &h(F(\mathcal{\varrho },\tau (\mathcal{\varrho })\overline{x}),F( \mathcal{\varrho },\tau (\mathcal{\varrho })\overline{x}_{0}))+h(F(\mathcal{ \varrho },\tau (\mathcal{\varrho })\overline{x}_{0}),\{0\}) \\ &\leq &\sigma (\mathcal{\varrho })\varsigma (\left\Vert x\right\Vert _{ \mathcal{H}})+\sigma (\mathcal{\varrho }) \\ &\leq &\sigma (\mathcal{\varrho })(1+\varsigma (\left\Vert x\right\Vert _{ \mathcal{H}})),\rm{for\ }a.e.\mathcal{\varrho }\in J. \end{eqnarray*}

    Then, as the reasons of Theorem 1, we can show that the values of \Phi are closed. We show that \Phi is \ \varsigma - contraction. Let x_{2}, x_{1}\in \mathcal{H}\ and y_{1}\in \Phi (x_{1}). Then, there exists f\in S_{F(., \tau (.)\overline{x}_{1})}^{1} such that

    \begin{equation*} y_{1}(\mathcal{\varrho }) = \left\{ \begin{array}{l} 0,\mathcal{\varrho }\in \lbrack -r,0] \\ \frac{(\psi (\mathcal{\varrho })-\psi (0))^{\mu -1}}{\Gamma (\mu )}\Psi (0)+ \frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{ \varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)f(s)ds,\mathcal{\varrho } \in (0,\mathcal{\varrho }_{1}] \\ g_{i}(\mathcal{\varrho },x(\mathcal{\varrho }_{i}^{-})),\mathcal{\varrho } \in (\mathcal{\varrho }_{i},s_{i}],i = 1,......n, \\ \frac{(\psi (\mathcal{\varrho })-\psi (s_{i}))^{\mu -1}}{\Gamma (\mu )} g_{i}(s_{i},x(\mathcal{\varrho }_{i}^{-})) \\ +\frac{1}{\Gamma (\vartheta )}\int_{s_{i}}^{\mathcal{\varrho }}(\psi ( \mathcal{\varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)f(s),\mathcal{ \varrho }\in (s_{i}{\rm{ }}\mathcal{\varrho }_{i+1}],i = 1,......n. \end{array} \right. \end{equation*}

    Define \ \Pi :J\rightarrow 2^{E} as

    \begin{equation*} \Pi (\mathcal{\varrho }) = \{u\in E:||f\ (\mathcal{\varrho })-u\ ||\ \leq \sigma (\mathcal{\varrho })\varsigma (\left\Vert x_{1}-x_{2}\right\Vert _{ \mathcal{H}})\}. \end{equation*}

    From (F_{5})(i), we have

    \begin{eqnarray*} &&h(F(\mathcal{\varrho },\tau (\mathcal{\varrho })\overline{x}_{1}),F( \mathcal{\varrho },\tau (\mathcal{\varrho })\overline{x}_{2})) \\ &\leq &\sigma (\mathcal{\varrho })\varsigma (\left\Vert x_{1}-x_{2}\right\Vert _{\mathcal{H}}),\mathcal{\varrho }\in J, \end{eqnarray*}

    Since the values of \ F\ are compact, there exists u_{\mathcal{\varrho } } \in F(\mathcal{\varrho }, \tau (\mathcal{\varrho })x_{2}) such that

    \begin{equation*} ||u_{\mathcal{\varrho }}-f(\mathcal{\varrho })||\ \leq \sigma (\mathcal{ \varrho })\varsigma (\left\Vert x_{1}-x_{2}\right\Vert _{\mathcal{H}}), \end{equation*}

    as a consequence \Pi (\mathcal{\varrho }); \mathcal{\varrho }\in J is not empty. Furthermore, because f, \sigma, x_{1}, x_{2}\ are measurable and E\ is separable, it follows from [[42], Theorem Ⅲ.41], the multifunction \ s\rightarrow \Pi (s)\cap F(s, \tau (s)\overline{x}_{2}) is measurable and since its values are non-empty and compact, there is h\in S_{F(\mathcal{\varrho }, \tau (\mathcal{\varrho })\overline{x}_{2})}^{1}\ such that

    \begin{equation*} ||h(\mathcal{\varrho })-f(\mathcal{\varrho })||\ \leq \ \sigma (\mathcal{ \varrho })\varsigma (\left\Vert x_{1}-x_{2}\right\Vert _{\mathcal{H}}),a.e. \mathcal{\varrho }\in J. \end{equation*}

    Set

    \begin{equation*} y_{2}(\mathcal{\varrho }) = \left\{ \begin{array}{l} 0,\mathcal{\varrho }\in \lbrack -r,0] \\ \frac{(\psi (\mathcal{\varrho })-\psi (0))^{\mu -1}}{\Gamma (\mu )}\Psi (0)+ \frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{ \varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)h(s)ds,\mathcal{\varrho } \in (0,\mathcal{\varrho }_{1}] \\ g_{i}(\mathcal{\varrho },x_{2}(\mathcal{\varrho }_{i}^{-})),\mathcal{\varrho }\in (\mathcal{\varrho }_{i},s_{i}],i = 1,......n, \\ \frac{(\psi (\mathcal{\varrho })-\psi (s_{i}))^{\mu -1}}{\Gamma (\mu )} g_{i}(s_{i},x_{2}(\mathcal{\varrho }_{i}^{-})) \\ +\frac{1}{\Gamma (\vartheta )}\int_{s_{i}}^{\mathcal{\varrho }}(\psi ( \mathcal{\varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)h(s),\mathcal{ \varrho }\in (s_{i}{\rm{ }}\mathcal{\varrho }_{i+1}],i = 1,......n. \end{array} \right. \end{equation*}

    Notice that, y_{2}\in \ \Phi (x_{2}) and if \mathcal{\varrho }\in J_{0, }, we get from Holder's inequality

    \begin{eqnarray} &&(\psi (\mathcal{\varrho })-\psi (0))^{1-\mu }||y_{2}(\mathcal{\varrho } )-y_{1}(\mathcal{\varrho })||\ \\ &\leq &\frac{(\psi (\mathcal{\varrho })-\psi (0))^{1-\mu }}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{\varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)\sigma (s)\ ||h(s)-f(s)||ds \\ &\leq &\frac{(\psi (\mathcal{\varrho })-\psi (0))^{1-\mu }}{\Gamma (\vartheta )}\varsigma (\left\Vert x_{1}-x_{2}\right\Vert _{\mathcal{H} })\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{\varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)\sigma (s)ds \\ &\leq &\frac{(\psi (\mathcal{\varrho })-\psi (0))^{1-\mu }}{\Gamma (\vartheta )}\ \varsigma (\left\Vert x_{1}-x_{2}\right\Vert _{\mathcal{H} })||\sigma ||_{L^{\frac{1}{q}}(I, \mathbb{R} ^{+})}\delta ^{q}\frac{(\psi (b)-\psi (0))^{\vartheta -q}}{(\frac{\vartheta -q}{1-q})^{1-q}} \\ &\leq &\eta \delta ^{q}\frac{(\psi (b)-\psi (0))^{1-\mu }}{\Gamma (\vartheta )}||\sigma ||_{L^{\frac{1}{q}}(I, \mathbb{R} ^{+})}\varsigma (\left\Vert x_{1}-x_{2}\right\Vert _{\mathcal{H}}). \end{eqnarray} (30)

    Similarly, if \mathcal{\varrho }\in (s_{i} \mathcal{\varrho } _{i+1}], i = 1, ......, n, we get

    \begin{eqnarray} &&(\psi (\mathcal{\varrho })-\psi (s_{i}))||y_{2}(\mathcal{\varrho })-y_{1}( \mathcal{\varrho })|| \\ &\leq &\frac{1}{\Gamma (\mu )}||g_{i}(s_{i},x_{1}(\mathcal{\varrho } _{i}^{-}))-g_{i}(s_{i},x_{2}(\mathcal{\varrho }_{i}^{-}))|| \\ &&+\eta \delta ^{q}\frac{(\psi (b)-\psi (0))^{1-\mu }}{\Gamma (\vartheta )} ||\sigma ||_{L^{\frac{1}{q}}(I, \mathbb{R} ^{+})}\varsigma (\left\Vert x_{1}-x_{2}\right\Vert _{\mathcal{H}}). \\ &\leq &\frac{\xi }{\Gamma (\mu )}\varsigma (\left\Vert x_{1}-x_{2}\right\Vert _{\mathcal{H}})+\eta \delta ^{q}\frac{(\psi (b)-\psi (0))^{1-\mu }}{\Gamma (\vartheta )}||\sigma ||_{L^{\frac{1}{q}}(I, \mathbb{R} ^{+})}\varsigma (\left\Vert x_{1}-x_{2}\right\Vert _{\mathcal{H}}). \\ &\leq &\varsigma (\left\Vert x_{1}-x_{2}\right\Vert _{\mathcal{H}})[\frac{ \xi }{\Gamma (\mu )}+\ \eta \delta ^{q}\frac{(\psi (b)-\psi (0))^{1-\mu }}{ \Gamma (\vartheta )}||\sigma ||_{L^{\frac{1}{q}}(I, \mathbb{R} ^{+})}], \end{eqnarray} (31)

    where \xi = \sum_{i = 1}^{i = n}\xi _{i}. Furthermore, If \mathcal{\varrho } \in (\mathcal{\varrho }_{i}, s_{i}], i = 1, ......n, \ then by (H_{2})\; one obtains

    \begin{eqnarray} ||y_{2}(\mathcal{\varrho })-y_{1}(\mathcal{\varrho })|| &\leq &||g_{i}( \mathcal{\varrho },x(\mathcal{\varrho }_{i}^{-}))-g_{i}(\mathcal{\varrho },x( \mathcal{\varrho }_{i}^{-}))|| \\ &\leq &\xi \varsigma (||x_{1}(\mathcal{\varrho }_{i}^{-})-x_{2}(\mathcal{ \varrho }_{i}^{-})||)_{\ } \\ &\leq &\xi \ \varsigma (\left\Vert x_{1}-x_{2}\right\Vert _{\mathcal{H}}). \end{eqnarray} (32)

    It follows from (30)–(32),

    \begin{eqnarray*} &&h_{\mathcal{H}}(\Phi (x_{2}),\Phi (x_{1})) \\ & < &\varsigma (\left\Vert x_{1}-x_{2}\right\Vert _{\mathcal{H}})[\frac{\xi }{ \Gamma (\mu )}+\ \eta \delta ^{q}\frac{(\psi (\mathcal{\varrho })-\psi (0))^{1-\mu }}{\Gamma (\vartheta )}||\sigma ||_{L^{\frac{1}{q}}(I, \mathbb{R} ^{+})}], \end{eqnarray*}

    where, h_{\mathcal{H}} is the Hausdorff distance in \mathcal{H} . This inequality and (29) imply that \Phi is \varsigma - contraction and thus by Theorem 3.1(i) in [41], \Phi 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 \ \psi- generalized Ulam-Hyers stable. Let \epsilon > 0 and y\in M\ be a solution of inequality (27). Because the values of F\ are compact, there f\in S_{F(., \tau (\mathcal{\varrho })y)}^{1} such that

    \begin{equation*} ||D_{s_{i}^{+}}^{\vartheta ,\nu ,\psi }y(\mathcal{\varrho })-f(\mathcal{ \varrho })|| = d(D_{s_{i}^{+}}^{\vartheta ,\nu ,\psi }y(\mathcal{\varrho }),F( \mathcal{\varrho },\tau (\mathcal{\varrho })y)),a.e.{\rm{ }}\mathcal{\varrho }\in (s_{i},{\rm{ }}\mathcal{\varrho }_{i+1}],i = 0,1,....,n. \end{equation*}

    Then, for almost everywhere \mathcal{\varrho }\in \cup _{i = 0}^{i = n}(s_{i}, \mathcal{\varrho }_{i+1}]

    \begin{equation*} D_{s_{i}^{+}}^{\vartheta ,\nu ,\psi }y(\mathcal{\varrho }) = w(\mathcal{ \varrho })+f(\mathcal{\varrho }), \end{equation*}

    where \ w\in PC_{1-\mu, \psi }^{\nu (1-\vartheta)}(J, E)\ and ||w(\mathcal{ \varrho })|\leq \epsilon, \forall \mathcal{\varrho }\in J . Furthermore, according to (F_{4}), \ there is h\in PC_{1-\mu, \psi }(J, E) such that f(\mathcal{\varrho }) = I_{s_{i}+}^{\nu (1-\vartheta), \psi }h(\mathcal{\varrho }), \mathcal{\varrho }\in J_{k}, \ k = 0, 1, 2, .., n , and hence, from Lemma 3.2,

    \begin{equation*} D_{s_{i}+}^{\nu (1-\vartheta ),\psi }f(\mathcal{\varrho }) = D_{s_{i}+}^{\nu (1-\vartheta ),\psi }I_{s_{i}+}^{\nu (1-\vartheta )\psi }h(\mathcal{\varrho } ) = h(\mathcal{\varrho }),\mathcal{\varrho }\in J_{k},\ k = 0,1,2,..,n. \end{equation*}

    This yields that f\in PC_{1-\mu, \psi }^{\nu (1-\vartheta)}(J, E). Therefore, y(\mathcal{\varrho }) = \Psi (\mathcal{\varrho }), \mathcal{\varrho }\in \lbrack -r, 0]\ and

    \begin{equation} \left\{ \begin{array}{l} D_{s_{i}^{+}}^{\vartheta ,\nu ,\psi }y(\mathcal{\varrho }) = f(\mathcal{ \varrho })+\vartheta (\mathcal{\varrho }),{\rm{ }}a.e.{\rm{ }}\mathcal{ \varrho }\in (s_{i},{\rm{ }}\mathcal{\varrho }_{i+1}],i = 0,1,....,n, \\ y(\mathcal{\varrho }) = g_{i}(\mathcal{\varrho },y(\mathcal{\varrho } _{i}^{-}))+\epsilon ,\mathcal{\varrho }\in (\mathcal{\varrho }_{i}{\rm{ }} s_{i}],i = 1,......n,, \\ y(\mathcal{\varrho }_{i}^{+}) = g_{i}(\mathcal{\varrho }_{i},y(\mathcal{ \varrho }_{i}^{-}))+\epsilon ,i = 1,......n, \end{array} \right. \end{equation} (33)

    In view of the second assertion of Lemma 3.3, relation (33) one obtains

    \begin{equation} y(\mathcal{\varrho }) = \left\{ \begin{array}{l} \Psi (\mathcal{\varrho }),\mathcal{\varrho }\in \lbrack -r,0]\ \\ \frac{(\psi (\mathcal{\varrho })-\psi (0))^{\mu -1}}{\Gamma (\mu )}\Psi (0)+ \frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{ \varrho })-\psi (s))^{\vartheta -1}f(s)ds \\ +\frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{ \varrho })-\psi (s))^{\vartheta -1}\vartheta (\mathcal{\varrho })ds,\mathcal{ \varrho }\in (0,\mathcal{\varrho }_{1}], \\ g_{i}(\mathcal{\varrho },y(\mathcal{\varrho }_{i}^{-}))+\epsilon ,\mathcal{ \varrho }\in (\mathcal{\varrho }_{i}{\rm{ }}s_{i}],i = 1,2,\ldots ,n, \\ \frac{(\psi (\mathcal{\varrho })-\psi (s_{i}))^{\mu -1}}{\Gamma (\mu )} y(s_{i})+\frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi ( \mathcal{\varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)(f(s))ds \\ +\frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{ \varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)\vartheta (\mathcal{ \varrho })ds,\mathcal{\varrho }\in (s_{i}{\rm{ }}\mathcal{\varrho } _{i+1}],i = 1,2,\ldots ,n. \end{array} \right. \end{equation} (34)

    Let

    \begin{equation*} z(\mathcal{\varrho }) = \left\{ \begin{array}{l} 0,\mathcal{\varrho }\in \lbrack -r,0]\ \\ \frac{(\psi (\mathcal{\varrho })-\psi (0))^{\mu -1}}{\Gamma (\mu )}\Psi (0)+ \frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{ \varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)f(s)ds,\mathcal{\varrho } \in (0,\mathcal{\varrho }_{1}], \\ g_{i}(\mathcal{\varrho },y(\mathcal{\varrho }_{i}^{-})),\mathcal{\varrho } \in (\mathcal{\varrho }_{i}{\rm{ }}s_{i}],i = 1,2,\ldots ,n, \\ \frac{(\psi (\mathcal{\varrho })-\psi (s_{i}))^{\mu -1}}{\Gamma (\mu )} g_{i}(s_{i},y(\mathcal{\varrho }_{i}^{-}))+\frac{1}{\Gamma (\vartheta )} \int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{\varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)(f(s))ds| \\ +\frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{ \varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)\vartheta (\mathcal{ \varrho })ds,\mathcal{\varrho }\in (s_{i}{\rm{ }}\mathcal{\varrho } _{i+1}],i = 1,2,\ldots ,n. \end{array} \right. \end{equation*}

    Obviously, z\in \Phi (y^{\ast }) , where

    \begin{equation*} y^{\ast }(\mathcal{\varrho }) = \left\{ \begin{array}{l} 0,\mathcal{\varrho }\in \lbrack -r,0]\ \\ \frac{(\psi (\mathcal{\varrho })-\psi (0))^{\mu -1}}{\Gamma (\mu )}\Psi (0)+ \frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{ \varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)f(s)ds \\ +\frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{ \varrho })-\psi (s))^{\vartheta -1}\vartheta (\mathcal{\varrho })ds,\mathcal{ \varrho }\in (0,\mathcal{\varrho }_{1}], \\ g_{i}(\mathcal{\varrho },y(\mathcal{\varrho }_{i}^{-}))+\epsilon ,\mathcal{ \varrho }\in (\mathcal{\varrho }_{i}{\rm{ }}s_{i}],i = 1,2,\ldots ,n, \\ \frac{(\psi (\mathcal{\varrho })-\psi (s_{i}))^{\mu -1}}{\Gamma (\mu )} y(s_{i})+\frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi ( \mathcal{\varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)(f(s))ds| \\ +\frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{ \varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)\vartheta (\mathcal{ \varrho })ds,\mathcal{\varrho }\in (s_{i}{\rm{ }}\mathcal{\varrho } _{i+1}],i = 1,2,\ldots ,n. \end{array} \right. \end{equation*}

    Moreover, from (33) and (34) we get for \mathcal{\varrho }\in (0, \mathcal{\varrho }_{1}]

    \begin{eqnarray*} ||z(\mathcal{\varrho })-y^{\ast }(\mathcal{\varrho })|| &\leq &\frac{1}{ \Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{\varrho } )-\psi (s))^{\vartheta -1}\psi \prime (s)||\vartheta (\mathcal{\varrho })||ds \\ &\leq &\frac{\epsilon (\psi (b)-\psi (0))^{\vartheta }}{\Gamma (\vartheta +1) }. \end{eqnarray*}

    If \mathcal{\varrho }\in (\mathcal{\varrho }_{i} s_{i}], then ||\overline{ z}(\mathcal{\varrho })-y(\mathcal{\varrho })||\leq \epsilon. If \ \mathcal{ \varrho }\in (s_{i} \mathcal{\varrho }_{i+1}], then

    \begin{eqnarray*} ||z(\mathcal{\varrho })-y^{\ast }(\mathcal{\varrho })|| &\leq &||\frac{(\psi (\mathcal{\varrho })-\psi (s_{i}))^{\mu -1}}{\Gamma (\mu )}g_{i}(s_{i},y( \mathcal{\varrho }_{i}^{-}))\frac{(\psi (\mathcal{\varrho })-\psi (s_{i}))^{\mu -1}}{\Gamma (\mu )}y(s_{i})|| \\ &&+\frac{1}{\Gamma (\vartheta )}\int_{0}^{\mathcal{\varrho }}(\psi (\mathcal{ \varrho })-\psi (s))^{\vartheta -1}\psi \prime (s)||\vartheta (\mathcal{ \varrho })||ds \\ &\leq &\frac{(\psi (\mathcal{\varrho })-\psi (s_{i}))^{\mu -1}}{\Gamma (\mu ) }||g_{i}(s_{i},y(\mathcal{\varrho }_{i}^{-}))-y(s_{i})|| \\ &&+\frac{\epsilon (\psi (b)-\psi (0))^{\vartheta }}{\Gamma (\vartheta +1)} \\ &\leq &\epsilon \frac{(\psi (b)-\psi (0))^{\mu -1}}{\Gamma (\mu )}+\frac{ \epsilon (\psi (b)-\psi (0))^{\vartheta }}{\Gamma (\vartheta +1)} \\ & = &\epsilon \lbrack \frac{(\psi (b)-\psi (0))^{\mu -1}}{\Gamma (\mu )}+\frac{ (\psi (b)-\psi (0))^{\vartheta }}{\Gamma (\vartheta +1)}+1] \end{eqnarray*}

    As a consequence,

    \begin{equation} \sup\limits_{\mathcal{\varrho }\in J}||y^{\ast }(\mathcal{\varrho })-z(\mathcal{ \varrho })||\leq \epsilon \lbrack \frac{(\psi (b)-\psi (0))^{\mu -1}}{\Gamma (\mu )}+\frac{(\psi (b)-\psi (0))^{\vartheta }}{\Gamma (\vartheta +1)}+1]. \end{equation} (35)

    On the other hand, from the facts that \Phi \ is \varsigma - contraction, \varsigma is strictly comparison, \varsigma (c\mathcal{ \varrho })\leq c\ \varsigma (\mathcal{\varrho })(c > 1)\ for every \mathcal{ \varrho }\in \lbrack 0, \infty) and \mathcal{\varrho } = 0 is a point of uniform convergence for the series \sum\limits_{n = 1}^{\infty }\varsigma ^{n}(\mathcal{\varrho }), it yields from Theorem3.1(ii)[41], \ that \Phi \ is \rho - multivalued weakly Picard operator, where \rho (\mathcal{\varrho }) = \mathcal{\varrho }+\sum\limits_{n = 1}^{\infty }\ \varsigma ^{n}(\mathcal{\varrho }). Then, the function \Phi ^{\infty }:Graph(\Phi)\rightarrow Fix(\Phi)\ is well defined and

    \begin{equation} ||y^{\ast }-\Phi ^{\infty }(y^{\ast },z)||_{\mathcal{H}}\leq \rho (||y^{\ast }-z||_{\mathcal{H}}). \end{equation} (36)

    Put x = \Phi ^{\infty }(y^{\ast }, z).\ So, x\in \Phi (x)\ and from (35) and (36) we get

    \begin{eqnarray*} d_{H}(y,\overline{x}) &\leq &||y^{\ast }-\Phi ^{\infty }(y^{\ast },z)||_{\mathcal{H}}\leq \rho (||y^{\ast }-z||_{\mathcal{H}}) \\ & = &\rho (\epsilon \lbrack \frac{(\psi (b)-\psi (0))^{\mu -1}}{\Gamma (\mu )}+ \frac{(\psi (b)-\psi (0))^{\vartheta }}{\Gamma (\vartheta +1)}+1]) = \theta (\epsilon ), \end{eqnarray*}

    where \theta (\mathcal{\varrho }) = \rho (\mathcal{\varrho }[\frac{(\psi (b)-\psi (0))^{\mu -1}}{\Gamma (\mu)}+\frac{(\psi (b)-\psi (0))^{\vartheta } }{\Gamma (\vartheta +1)}+1]). Consequently, Problem (1) is \psi - 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, s_{0} = 0, \mathcal{\varrho }_{1} = 1, s_{1} = 2, \mathcal{\varrho }_{2} = 3, \vartheta = \frac{1}{2}, \nu = \frac{1}{4} and \mu = \vartheta +\nu -\vartheta \nu = \frac{5}{8} . Let F:J\times \Theta \rightarrow P_{ck}(E) be defined as follows:

    \begin{equation} F(\mathcal{\varrho },z) = \left\{ \begin{array}{l} (\psi (\mathcal{\varrho })-\psi (0))^{\frac{3}{8}}||z(0)||K,\ \mathcal{ \varrho }\in \lbrack 0,1], \\ ||z(0)||K,\mathcal{\varrho }\in (1,{\rm{ }}2], \\ (\psi (\mathcal{\varrho })-\psi (2))^{\frac{3}{8}}||z(0)||K,\mathcal{\varrho }\in (2,3]. \end{array} \right. \end{equation} (37)

    Clearly, for almost every \mathcal{\varrho }\in J, z\longrightarrow F(\mathcal{\varrho }, \tau (\mathcal{\varrho })z)\ is upper semicontinuous. Set \lambda = \sup \{||x||:x\in K\} . Then, for any (\mathcal{\varrho }, z)\in J\times \Theta and any y\in F(\mathcal{\varrho }, z) we have

    \begin{equation*} ||y(\mathcal{\varrho })||\leq \left\{ \begin{array}{l} (1+(\psi (\mathcal{\varrho })-\psi (0))^{\frac{3}{8}}||z(0)||)\lambda ,\ \mathcal{\varrho }\in \lbrack 0,1], \\ ||z(0)||\lambda ,\mathcal{\varrho }\in (1,{\rm{ }}2], \\ (1+(\psi (\mathcal{\varrho })-\psi (2))^{\frac{3}{8}}||z(0)||)\lambda , \mathcal{\varrho }\in (2,3]. \end{array} \right. \end{equation*}

    Then condition (F_{2}) of Theorem1 is verified with \varphi (\mathcal{ \varrho }) = \lambda; \mathcal{\varrho }\in J . Moreover, if f\in S_{F(, .z)}^{1};z\in \Theta , then

    \begin{equation} f(\mathcal{\varrho }) = \left\{ \begin{array}{l} (\psi (\mathcal{\varrho })-\psi (0))^{\frac{3}{8}}||z(0)||y_{0},\ \mathcal{ \varrho }\in \lbrack 0,1], \\ y_{1},\mathcal{\varrho }\in (1,{\rm{ }}2], \\ (\psi (\mathcal{\varrho })-\psi (2))^{\frac{3}{8}}y_{2},\mathcal{\varrho } \in (2,3]. \end{array} \right. \end{equation} (38)

    We define

    \begin{equation} h_{x_{1},x_{2},x_{3}}(\mathcal{\varrho }) = \left\{ \begin{array}{l} \frac{\Gamma (\frac{11}{8})}{\Gamma (\frac{10}{8})}(\psi (\mathcal{\varrho } )-\psi (0))^{\frac{2}{8}}x_{1},\ \mathcal{\varrho }\in \lbrack 0,1], \\ x_{2},\mathcal{\varrho }\in (1,{\rm{ }}2], \\ \frac{\Gamma (\frac{11}{8})}{\Gamma (\frac{10}{8})}(\psi (\mathcal{\varrho } )-\psi (2))^{\frac{2}{8}}x_{3},\mathcal{\varrho }\in (2,3]. \end{array} \right. \end{equation} (39)

    Obviously h_{x_{1}, x_{2}, x_{3}}\in PC_{1-\mu, \psi }(J, E) and from Lemma 2.2, we get I_{s_{i}^{+}}^{\nu (1-\vartheta), \psi }h_{x_{1}, x_{2}, x_{3}}(\mathcal{\varrho }) = f(\mathcal{\varrho }); i = 0, 1 .

    Now, let D\ be a bounded subset of \mathcal{H} , z_{1} , z_{2}\in D , \mathcal{\varrho }\in (0, 1] , \ x\in F(\mathcal{\varrho }.\tau (\mathcal{ \varrho })\overline{z_{1}}) and y\in F(\mathcal{\varrho }.\tau (\mathcal{ \varrho })\overline{z_{2}}) . Then x = (\psi (\mathcal{\varrho })-\psi (0))^{ \frac{3}{8}}||z_{1}(\mathcal{\varrho })||x^{\ast } and y = (\psi (\mathcal{ \varrho })-\psi (0))^{\frac{3}{8}}||z_{2}(\mathcal{\varrho })||y^{\ast } , where x^{\ast }, y^{\ast }\in K . It follows that

    \begin{equation*} ||x-y||\leq (\psi (\mathcal{\varrho })-\psi (0))^{\frac{3}{8}}||x^{\ast }||z_{1}(\mathcal{\varrho })||-||z_{2}(\mathcal{\varrho })||y^{\ast }||\rm{ ,} \end{equation*}

    As a result,

    \begin{eqnarray*} &&\inf\limits_{y\in F(\mathcal{\varrho }.\tau (\mathcal{\varrho })\overline{z_{2}} )}||x-y|| \\ &\leq &(\psi (\mathcal{\varrho })-\psi (0))^{\frac{3}{8}}||x^{\ast }||z_{1}( \mathcal{\varrho })||-||z_{2}(\mathcal{\varrho })||x^{\ast }|| \\ & = &||x^{\ast }||(\psi (\mathcal{\varrho })-\psi (0))^{\frac{3}{8}}|||z_{1}( \mathcal{\varrho })||-||z_{2}(\mathcal{\varrho })|| \\ &\leq &\lambda ((\psi (\mathcal{\varrho })-\psi (0))^{\frac{3}{8}}||z_{1}( \mathcal{\varrho })-z_{2}(\mathcal{\varrho })||, \end{eqnarray*}

    which means that

    \begin{equation} h(F(\mathcal{\varrho }.\tau (\mathcal{\varrho })\overline{z_{1}}),F(\mathcal{ \varrho }.\tau (\mathcal{\varrho })\overline{z_{2}})\leq \lambda ((\psi ( \mathcal{\varrho })-\psi (0))^{\frac{3}{8}}||z_{1}(\mathcal{\varrho })-z_{2}( \mathcal{\varrho })||. \end{equation} (40)

    Then, for \mathcal{\varrho }\in (0, \mathcal{\varrho }_{1}]

    \begin{eqnarray*} \ \chi _{E}(F(\mathcal{\varrho },\{\tau (\mathcal{\varrho })\overline{x} &:&x\in D)) \\ &\leq &\lambda (\psi (\mathcal{\varrho })-\psi (0))^{1-\mu }\varsigma ( \mathcal{\varrho })\chi _{E}\{x(\mathcal{\varrho }):x\in D\} \\ &\leq &\lambda (\psi (\mathcal{\varrho })-\psi (0))^{1-\mu }\varsigma ( \mathcal{\varrho })\sup\limits_{\theta \in \lbrack -r,0]}\chi _{E}\{(\tau (\mathcal{ \varrho }))\overline{x}(\theta ):x\in D\}, \end{eqnarray*}

    Similarly, one can show that if for \mathcal{\varrho }\in (s_{2}, \mathcal{ \varrho }_{3}]

    \begin{equation*} \chi _{E}(F(\mathcal{\varrho },D))\leq \lambda (\psi (\mathcal{\varrho } )-\psi (s_{2}))^{1-\mu }\sup\limits_{\theta \in \lbrack -r,0]}\chi _{E}\{(\tau ( \mathcal{\varrho }))\overline{x}(\theta ):x\in D\} \end{equation*}

    and consequently, by choosing \lambda small enough such that the relation (7) becomes realized with \varsigma (\mathcal{\varrho }) = \lambda; \mathcal{ \varrho }\in J. Next, let g_{1}:[1, 2]\times E\rightarrow E such that

    \begin{equation} g(\mathcal{\varrho },x) = \rho (\psi (1)-\psi (0))^{\frac{3}{8}}\mathcal{ \varrho }(x), \end{equation} (41)

    where, \rho is a positive number and \mathcal{\varrho }\ is a linear, bounded and compact operator on E. So, condition (H_{1})\ 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 \Psi \in \Theta provided that

    \begin{equation} \ \frac{3\lambda \kappa _{1}}{\Gamma (\frac{1}{2})}\eta +\rho +\frac{\rho }{ \Gamma (\frac{5}{8})} < 1, \end{equation} (42)

    where, \eta = \delta ^{q}\frac{(\psi (b)-\psi (0))^{\vartheta -q}}{(\frac{ \vartheta -q}{1-q})^{1-q}} , \delta = \max_{s\in J}\psi \prime (s) , q = \frac{1}{4} and \kappa _{1} = (\psi (b)-\psi (0))^{1-\mu } . By choosing \rho small enough relation (42) becomes realized.

    EXAMPLE 2 Let E , K , \ r, \ s_{0}, \ \mathcal{\varrho }_{1}, \ s_{1}, \ \mathcal{ \varrho }_{2}, \vartheta, \ \nu, \ \mu , F:J\times \Theta \rightarrow P_{ck}(E)\ and g_{1}:[1, 2]\times E\rightarrow E be as in EXAMPLE 1. Then assumption (F_{4}) is satisfied. Moreover, in view of (40) for any \ z_{1}, z_{1}\in \mathcal{H}

    \begin{eqnarray*} h(F(\mathcal{\varrho }.\tau (\mathcal{\varrho })\overline{z_{1}}),F(\mathcal{ \varrho }.\tau (\mathcal{\varrho })\overline{z_{2}}) &\leq &\lambda ((\psi ( \mathcal{\varrho })-\psi (0))^{\frac{3}{8}}||z_{1}(\mathcal{\varrho })-z_{2}( \mathcal{\varrho })|| \\ &\leq &\lambda ||z_{1}-z_{2}||_{\mathcal{H}} \\ & = &\sigma (\mathcal{\varrho })\varsigma (||z_{1}-z_{2}||_{\mathcal{H}}), \end{eqnarray*}

    where, \sigma (\mathcal{\varrho }) = 2\lambda and \varsigma (\mathcal{ \varrho }) = \frac{\mathcal{\varrho }}{2};\mathcal{\varrho }\in J . Observe that \varsigma \ is strictly comparison and \varsigma (c\mathcal{\varrho }) = c\ \varsigma (\mathcal{\varrho }), for \ every c > 0, every \ \mathcal{ \varrho }\in \lbrack 0, \infty) and \mathcal{\varrho } = 0 is a point of uniform convergence for the series \sum\limits_{n = 1}^{\infty }\varsigma ^{n}(\mathcal{\varrho }), and hence assumption (F_{5}) is satisfied. Now, for any \mathcal{\varrho }\in \lbrack 1, 2] and any x, y\in E

    \begin{eqnarray*} ||g_{1}(\mathcal{\varrho },x)-g_{1}(\mathcal{\varrho },y)|| &\leq &\rho (\psi (1)-\psi (0))^{\frac{3}{8}}||\mathcal{\varrho }(x)-\mathcal{\varrho } (y)|| \\ &\leq &\rho ||\mathcal{\varrho }||(\psi (1)-\psi (0))^{\frac{3}{8} }||x-y||_{E} \\ & = &\xi \varsigma (||x-y||_{E}), \end{eqnarray*}

    where \xi = 2\rho ||\mathcal{\varrho }||(\psi (1)-\psi (0))^{\frac{3}{8}} . It follows that (H_{2}) is satisfied. By applying Theorem 2, Problem (1) has a solution and it is \psi - generalized Ulam-Hyers stable, where F and g are given by (37) and (40) and \Psi \in \Theta provided that

    \begin{equation*} \frac{2\rho ||\mathcal{\varrho }||(\psi (1)-\psi (0))^{\frac{3}{8}}}{\Gamma (\mu )}+\frac{2\rho ||\mathcal{\varrho }||(\psi (1)-\psi (0))^{\frac{3}{8}}}{ \Gamma (\mu )}+\eta ||\sigma ||_{L^{\frac{1}{q}}(I, \mathbb{R} ^{+})} < 1. \end{equation*}

    By choosing \rho 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 \psi -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 \psi - 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 \psi -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] S. Anis, M. Khan, Y. B. Jun, Hybrid ideals in semigroups, Cogent Math., 4 (2017), 1–12. https://doi.org/10.1080/23311835.2017.1352117 doi: 10.1080/23311835.2017.1352117
    [2] J. R. Clay, The near-rings on groups of low order, Math. Z., 104 (1968), 364–371. https://doi.org/10.1007/BF01110428 doi: 10.1007/BF01110428
    [3] P. Dheena, G. S. Kumar, On strongly regular near-subtraction semigroups, Commun. Korean Math. Soc., 22 (2007), 323–330. https://doi.org/10.4134/CKMS.2007.22.3.323 doi: 10.4134/CKMS.2007.22.3.323
    [4] B. Elavarasan, Y. B. Jun, Regularity of semigroups in terms of hybrid ideals and hybrid bi-ideals, Kragujev. J. Math., 46 (2022), 857–864.
    [5] B. Elavarasan, G. Muhiuddin, K. Porselvi, Y. B. Jun, Hybrid structures applied to ideals in near-rings, Complex Intell. Syst., 7 (2021), 1489–1498. https://doi.org/10.1007/s40747-021-00271-7 doi: 10.1007/s40747-021-00271-7
    [6] B. Elavarasan, K. Porselvi, Y. B. Jun, Hybrid generalized bi-ideals in semigroups, Int. J. Math. Comput. Sci., 14 (2019), 601–612.
    [7] Y. B. Jun, H. S. Kim, E. H. Roh, Ideal theory of subtraction algebras, Sci. Math. Jpn., 61 (2005), 459–464.
    [8] Y. B. Jun, H. S. Kim, On ideals in subtraction algebras, Sci. Math. Jpn., 65 (2007), 129–134. https://doi.org/10.32219/isms.65.1_129 doi: 10.32219/isms.65.1_129
    [9] Y. B. Jun, M. Sapanci, M. A. Öztürk, Fuzzy ideals in Gamma near-rings, Tr. J. Math., 22 (1998), 449–459.
    [10] Y. B. Jun, S. Z. Song, G. Muhiuddin, Hybrid structures and applications, Ann. Commun. Math., 1 (2018), 11–25.
    [11] K. J. Lee, C. H. Park, Some questions on fuzzifications of ideals in subtraction algebras, Commun. Korean Math. Soc., 22 (2007), 359–363. https://doi.org/10.4134/CKMS.2007.22.3.359 doi: 10.4134/CKMS.2007.22.3.359
    [12] G. Mason, Strongly regular near-rings, Proc. Edinb. Math. Soc., 23 (1980), 27–35. https://doi.org/10.1017/S0013091500003564 doi: 10.1017/S0013091500003564
    [13] P. K. Maji, A. R. Roy, R. Biswas, An application of soft sets in a decision making problem, Comput. Math. Appl., 44 (2002), 1077–1083. https://doi.org/10.1016/S0898-1221(02)00216-X doi: 10.1016/S0898-1221(02)00216-X
    [14] J. D. P. Meldrum, Varieties and d. g. near-rings, Proc. Edinb. Math. Soc., 17 (1971), 271–274. https://doi.org/10.1017/S0013091500027000 doi: 10.1017/S0013091500027000
    [15] D. Molodtsov, Soft set theory–first results, Comput. Math. Appl., 37 (1999), 19–31. https://doi.org/10.1016/S0898-1221(99)00056-5 doi: 10.1016/S0898-1221(99)00056-5
    [16] G. Muhiuddin, J. C. G. John, B. Elavarasan, Y. B. Jun, K. Porselvi, Hybrid structures applied to modules over semirings, J. Intell. Fuzzy Syst., 42 (2022), 2521–2531. https://doi.org/10.3233/JIFS-211751 doi: 10.3233/JIFS-211751
    [17] G. Muhiuddin, J. C. G. John, B. Elavarasan, K. Porselvi, D. Al-Kadi, Properties of k-hybrid ideals in ternary semiring, J. Intell. Fuzzy Syst., 42 (2022), 5799–5807. https://doi.org/10.3233/JIFS-212311 doi: 10.3233/JIFS-212311
    [18] G. Muhiuddin, D. Al-Kadi, W. A. Khan, C. Jana, Hybrid structures applied to subalgebras of BCH-algebras, Secur. Commun. Netw., 2021 (2021), 1–8. https://doi.org/10.1155/2021/8960437 doi: 10.1155/2021/8960437
    [19] G. Muhiuddin, D. Al-Kadi, A. Mahboob, Hybrid structures applied to ideals in BCI-algebras, J. Math., 2020 (2020), 1–7. https://doi.org/10.1155/2020/2365078 doi: 10.1155/2020/2365078
    [20] K. Porselvi, B. Elavarasan, On hybrid interior ideals in semigroups, Probl. Anal. Issues Anal., 8 (2019), 137–146. https://doi.org/10.15393/j3.art.2019.6150 doi: 10.15393/j3.art.2019.6150
    [21] K. Porselvi, B. Elavarasan, Y. B. Jun, Hybrid interior ideals in ordered semigroups, New Math. Nat. Comput., 18 (2022), 1–8. https://doi.org/10.1142/S1793005722500016 doi: 10.1142/S1793005722500016
    [22] G. Muhiuddin, H. Harizavi, Y. B. Jun, Bipolar-valued fuzzy soft hyper BCK ideals in hyper BCK algebras, Discrete Math. Algorithms Appl., 12 (2020), 2050018. https://doi.org/10.1142/S1793830920500184 doi: 10.1142/S1793830920500184
    [23] G. Muhiuddin, Bipolar fuzzy KU-subalgebras/ideals of KU-algebras, Ann. Fuzzy Math. Inform., 8 (2014), 409–418.
    [24] G. Muhiuddin, D. Al-Kadi, A. Mahboob, A. Albjedi, Interval-valued m-polar fuzzy positive implicative ideals in BCK-algebras, Math. Probl. Eng., 2021 (2021), 1–9. https://doi.org/10.1155/2021/1042091 doi: 10.1155/2021/1042091
    [25] G. Muhiuddin, D. Al-Kadi, Interval valued m-polar fuzzy BCK/BCI-algebras, Int. J. Comput. Intell. Syst., 14 (2021), 1014–1021. https://doi.org/10.2991/ijcis.d.210223.003 doi: 10.2991/ijcis.d.210223.003
    [26] G. Muhiuddin, D. Al-Kadi, A. Mahboob, A. Aljohani, Generalized fuzzy ideals of BCI-algebras based on interval valued m-polar fuzzy structures, Int. J. Comput. Intell. Syst., 14 (2021), 1–9. https://doi.org/10.1007/s44196-021-00006-z doi: 10.1007/s44196-021-00006-z
    [27] G. Muhiuddin, K. P. Shum, New types of (\alpha, \beta)-fuzzy subalgebras of BCK/BCI-algebras, Int. J. Math. Comput. Sci., 14 (2019), 449–464.
    [28] G. Muhiuddin, A. M. Al-Roqi, Subalgebras of BCK/BCI-algebras based on (\alpha, \beta)-type fuzzy sets, J. Comput. Anal. Appl., 18 (2015), 1057–1064.
    [29] D. Al-Kadi, G. Muhiuddin, Bipolar fuzzy BCI-implicative ideals of BCI-algebras, Ann. Commun. Math., 3 (2020), 88–96.
    [30] A. Rosenfeld, Fuzzy groups, J. Math. Anal. Appl., 35 (1971), 512–517.
    [31] B. M. Schein, Difference semigroups, Commun. Algebra, 20 (1992), 2153–2169. https://doi.org/10.1080/00927879208824453 doi: 10.1080/00927879208824453
    [32] D. R. P. Williams, Fuzzy ideals in near-subtraction semigroups, Int. Scholarly Sci. Res. Innov., 2 (2008), 625–632.
    [33] L. A. Zadeh, Fuzzy sets, Inform. Control, 8 (1965), 338–353.
    [34] B. Zelinka, Subtraction semigroups, Math. Bohem., 120 (1995), 445–447. https://doi.org/10.21136/MB.1995.126093 doi: 10.21136/MB.1995.126093
  • This article has been cited by:

    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
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2213) PDF downloads(79) Cited by(9)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog