Processing math: 79%
Research article Special Issues

Conditional Ulam stability and its application to von Bertalanffy growth model

  • The purpose of this paper is to apply conditional Ulam stability, developed by Popa, Rașa, and Viorel in 2018, to the von Bertalanffy growth model dwdt=aw23bw, where w denotes mass and a>0 and b>0 are the coefficients of anabolism and catabolism, respectively. This study finds an Ulam constant and suggests that the constant is biologically meaningful. To explain the results, numerical simulations are performed.

    Citation: Masakazu Onitsuka. Conditional Ulam stability and its application to von Bertalanffy growth model[J]. Mathematical Biosciences and Engineering, 2022, 19(3): 2819-2834. doi: 10.3934/mbe.2022129

    Related Papers:

    [1] Weerawat Sudsutad, Sotiris K. Ntouyas, Chatthai Thaiprayoon . Nonlocal coupled system for ψ-Hilfer fractional order Langevin equations. AIMS Mathematics, 2021, 6(9): 9731-9756. doi: 10.3934/math.2021566
    [2] Abdelatif Boutiara, Mohammed S. Abdo, Manar A. Alqudah, Thabet Abdeljawad . On a class of Langevin equations in the frame of Caputo function-dependent-kernel fractional derivatives with antiperiodic boundary conditions. AIMS Mathematics, 2021, 6(6): 5518-5534. doi: 10.3934/math.2021327
    [3] Mohamed A. Barakat, Abd-Allah Hyder, Doaa Rizk . New fractional results for Langevin equations through extensive fractional operators. AIMS Mathematics, 2023, 8(3): 6119-6135. doi: 10.3934/math.2023309
    [4] Arjumand Seemab, Mujeeb ur Rehman, Jehad Alzabut, Yassine Adjabi, Mohammed S. Abdo . Langevin equation with nonlocal boundary conditions involving a ψ-Caputo fractional operators of different orders. AIMS Mathematics, 2021, 6(7): 6749-6780. doi: 10.3934/math.2021397
    [5] Abdelkrim Salim, Sabri T. M. Thabet, Imed Kedim, Miguel Vivas-Cortez . On the nonlocal hybrid (k,φ)-Hilfer inverse problem with delay and anticipation. AIMS Mathematics, 2024, 9(8): 22859-22882. doi: 10.3934/math.20241112
    [6] Zheng Kou, Saeed Kosari . On a generalization of fractional Langevin equation with boundary conditions. AIMS Mathematics, 2022, 7(1): 1333-1345. doi: 10.3934/math.2022079
    [7] Bashir Ahmad, Manal Alnahdi, Sotiris K. Ntouyas, Ahmed Alsaedi . On a mixed nonlinear boundary value problem with the right Caputo fractional derivative and multipoint closed boundary conditions. AIMS Mathematics, 2023, 8(5): 11709-11726. doi: 10.3934/math.2023593
    [8] Lakhlifa Sadek, Tania A Lazǎr . On Hilfer cotangent fractional derivative and a particular class of fractional problems. AIMS Mathematics, 2023, 8(12): 28334-28352. doi: 10.3934/math.20231450
    [9] Zaid Laadjal, Fahd Jarad . On a Langevin equation involving Caputo fractional proportional derivatives with respect to another function. AIMS Mathematics, 2022, 7(1): 1273-1292. doi: 10.3934/math.2022075
    [10] Kishor D. Kucche, Sagar T. Sutar, Kottakkaran Sooppy Nisar . Analysis of nonlinear implicit fractional differential equations with the Atangana-Baleanu derivative via measure of non-compactness. AIMS Mathematics, 2024, 9(10): 27058-27079. doi: 10.3934/math.20241316
  • The purpose of this paper is to apply conditional Ulam stability, developed by Popa, Rașa, and Viorel in 2018, to the von Bertalanffy growth model dwdt=aw23bw, where w denotes mass and a>0 and b>0 are the coefficients of anabolism and catabolism, respectively. This study finds an Ulam constant and suggests that the constant is biologically meaningful. To explain the results, numerical simulations are performed.



    Fractional differential equations constitute an interesting area of research in view of their extensive applications in a variety of fields such as economics, financial issues, disease models, physical and chemical processes, etc. For theoretical and applications details of the subject, see the monographs [6,12,16,19,23,24,26] and the papers [1,4,7,13,25,29,30]. An important factor accounting for the popularity of the topic is the nonlocal nature of fractional order differential and integral operators, which provides insight into the past history of the phenomena. On the other hand, this aspect was missing in the models based on the tools of classical (integer order) calculus.

    The Langevin equation (proposed by Langevin in 1908) is an important equation of mathematical physics, which can describe the evolution of physical phenomena in fluctuating environments [9]. For some new developments on the fractional Langevin equation, see, for example, [3,22,28] and references cited therein.

    In the literature there are many papers dealing with boundary value problems for fractional Langevin equation of the form:

    cDβ(cDα+λ)z(t)=ϕ(t,z(t)),t[a,b],

    equipped with different kinds of boundary conditions, such as Dirichlet boundary conditions [2], nonlocal conditions [27], etc. Moreover, the fractional derivatives involved in Langevin equation are both Caputo type or Riemann-Liouville type or Hadamard type.

    In the present paper we introduce a new nonlocal-terminal value problem for Langevin equation containing both Riemann-Liouville and Caputo type fractional derivatives and variable coefficient, supplemented with nonlocal-terminal fractional integro-differential conditions. In precise terms, we consider the following problem:

    {RLDq(CDr+λ(t))x(t)=f(t,x(t)),t[0,T],x(ξ)=αCDνx(η),x(T)=βIpx(ζ),0<ξ,η,ζ<T, (1.1)

    where RLDq denotes the Riemann-Liouville fractional derivative of order q(0,1), CDr, CDν denote the Caputo fractional derivatives of orders r and ν respectively, 0<r<1,0<ν<r, Ip is the Riemann-Liouville fractional integral of order p>0, λC(R+,R), f:[0,T]×R R and α,βR.

    Concerning the motivation for the equation considered in (1.1), we refer the reader to the applications of such equations in physical phenomena exhibiting anomalous diffusion [15] and the height loss over time of the granular material contained in a silo [21]. The nonlocal conditions involved in the problem (1.1) are flux-integral type boundary conditions, which appear in several applications of diffusion processes and computational fluid dynamics (CFD) studies of blood flow problems, for instance, see [5].

    The existence and uniqueness results for the problem (1.1), based on modern methods of functional analysis (fixed point theorems due to Banach, Krasnoselskii and nonlinear alternative of Leray-Schauder type), are obtained in Section 3. In Section 4 we extend our study to the multivalued analogue of the problem (1.1), that is, we investigate the following multivalued problem:

    {RLDq(CDr+λ(t))x(t)F(t,x(t)),t[0,T],x(ξ)=αCDνx(η),x(T)=βIpx(ζ),0<ξ,η,ζ<T, (1.2)

    where F:[0,T]×RP(R) is a multi-valued map, (P(R) is the family of all nonempty subsets of R) and all other constants are as in problem (1.1).

    We derive the existence results for the problem (1.2) with the aid of standard fixed point theorems for multivalued maps. In case of convex valued right-hand side of the inclusion, we use Leray-Schauder nonlinear alternative for multi-valued maps. In case of non-convex valued right hand side of the inclusion, we apply a fixed point theorem for multivalued contractions due to Covitz and Nadler. Examples illustrating the obtained results are presented in Section 5. Some basic concepts of fractional calculus, multivalued analysis and fixed point theory are outlined in Section 2. We also prove a basic result associated with the linear variant of the problem (1.1) in this section. Section 6 contains some interesting observations and a short discussion for the case when Caputo and Riemann-Liouville fractional derivatives are interchanged in the fractional Langevin equation in (1.1).

    Let us recall some basic concepts of fractional calculus, multivalued analysis and state fixed point results related to our work.

    In this subsection, we recall some basic concepts of fractional calculus [16,24] and present known results needed in our forthcoming analysis.

    Definition 2.1. The Riemann-Liouville fractional derivative of order q for a function f:(0,)R is defined by

    RLDqf(t)=1Γ(nq)(ddt)nt0+(ts)nq1f(s)ds,q>0,n=[q]+1,

    where [q] denotes the integer part of the real number q, provided the right-hand side is pointwise defined on (0,).

    Definition 2.2. The Riemann-Liouville fractional integral of order q for a function f:(0,)R is defined by

    Iqf(t)=1Γ(q)t0+(ts)q1f(s)ds,q>0,

    provided the right-hand side is pointwise defined on (0,).

    Definition 2.3. The Caputo derivative of fractional order q for a n-times derivative function f:(0,)R is defined as

    CDqf(t)=1Γ(nq)t0+(ts)nq1(dds)nf(s)ds,q>0,n=[q]+1.

    Lemma 2.1. If α+β>1, then the equation (IαIβu)(t)=(Iα+βu)(t),t[a,b] is satisfied for uL1([a,b],R),0a<b<.

    Lemma 2.2. Let β>α. Then the equation (DαIβu)(t)=(Iβαu)(t),t[a,b] is satisfied for uC([a,b],R).

    Lemma 2.3. Let n=[α]+1] if αN and n=α if αN. Then the following relations hold:

    (i) for k{0,1,2,,n1},Dαtk=0;

    (ii) if β>n then Dαtβ1=Γ(β)Γ(βα)tβα1;

    (iii) Iαtβ1=Γ(β)Γ(β+α)tβ+α1.

    Lemma 2.4. (see [16]) Let q>0. Then for yC(0,T)L(0,T) holds

    RLIq(RLDqy)(t)=y(t)+c1tq1+c2tq2++cntqn,

    where ciR, i=1,2,,n and n1<q<n.

    Lemma 2.5. (see [16]) Let q>0. Then for yC(0,T)L(0,T) holds

    RLIq(CDqy)(t)=y(t)+c0+c1t+c2t2++cn1tn1,

    where ciR, i=0,1,2,,n1 and n=[q]+1.

    Lemma 2.6. Let Λ:=AΓB0 and yC([0,T],R). Then the linear problem

    {RLDq(CDr+λ(t))x(t)=y(t),0<t<T,x(ξ)=αCDνx(η),x(T)=βIpx(ζ),ξ,η,ζ(0,T), (2.1)

    is equivalent to the integral equation

    x(t)=Iq+ry(t)Ir(λ(t)x(t))+λ1(t)(αIq+rνy(η)αIrν(λ(η)x(η))Iq+ry(ξ)+Ir(λ(ξ)x(ξ)))+λ2(t)(βIq+r+py(ζ)βIp+r(λ(ζ)x(ζ))Iq+ry(T)+Ir(λ(T)x(T))),

    where

    A=Γ(q)Γ(q+r)ξq+r1αΓ(q)Γ(q+rν)ηq+rν1,B=Γ(q)Γ(q+r)Tq+r1βΓ(q)Γ(q+r+p)ζq+r+p1,Γ=1β1Γ(1+p)ζp,λ1(t)=ΓΓ(q)ΛΓ(q+r)tq+r1BΛ,λ2(t)=Γ(q)ΛΓ(q+r)tq+r1+AΛ. (2.2)

    Proof. Firstly, we apply the Riemann-Liouville fractional integral of order q to both sides of equation in (2.1), and then use Lemma 2.5 to obtain

    CDrx(t)+λ(t)x(t)=Iqy(t)+c1tq1, (2.3)

    where c1R. Applying Riemann-Liouville fractional integral of order r to both sides of (2.3), we get

    x(t)=Iq+ry(t)Ir(λ(t)x(t))+c1Γ(q)Γ(q+r)tq+r1+c2, (2.4)

    where c2R.

    From (2.4), we have

    CDνx(t)=Iq+rνy(t)Irν(λ(t)x(t))+c1Γ(q)Γ(q+rν)tq+rν1,
    Ipx(t)=Iq+r+py(t)Ip+r(λ(t)x(t))+c1Γ(q)Γ(q+r+p)tq+r+p1+c21Γ(1+p)tp.

    Using the above expressions in the fractional nonlocal-terminal conditions of the problem (2.1), we find that

    c1=1Λ[Γ(αIq+rνy(η)αIrν(λ(η)x(η))Iq+ry(ξ)+Ir(λ(ξ)x(ξ)))(βIq+r+py(ζ)βIp+r(λ(ζ)x(ζ))Iq+ry(T)+Ir(λ(T)x(T)))]c2=1Λ[A(βIq+r+py(ζ)βIp+r(λ(ζ)x(ζ))Iq+ry(T)+Ir(λ(T)x(T)))B(αIq+rνy(η)αIrν(λ(η)x(η))Iq+ry(ξ)+Ir(λ(ξ)x(ξ)))].

    Substituting the values of c1 and c2 into equation (2.4) we have

    x(t)=Iq+ry(t)Ir(λ(t)x(t))+Γ(q)ΛΓ(q+r)tq+r1[Γ(αIq+rνy(η)αIrν(λ(η)x(η))Iq+ry(ξ)+Ir(λ(ξ)x(ξ)))(βIq+r+py(ζ)βIp+r(λ(ζ)x(ζ))Iq+ry(T)+Ir(λ(T)x(T)))]+1Λ[A(βIq+r+py(ζ)βIp+r(λ(ζ)x(ζ))Iq+ry(T)+Ir(λ(T)x(T)))B(αIq+rνy(η)αIrν(λ(η)x(η))Iq+ry(ξ)+Ir(λ(ξ)x(ξ)))]=Iq+ry(t)Ir(λ(t)x(t))+[ΓΛΓ(q)Γ(q+r)tq+r1BΛ]×(αIq+rνy(η)αIrν(λ(η)x(η))Iq+ry(ξ)+Ir(λ(ξ)x(ξ)))+[Γ(q)ΛΓ(q+r)tq+r1+AΛ]×(βIq+r+py(ζ)βIp+r(λ(ζ)x(ζ))Iq+ry(T)+Ir(λ(T)x(T))),

    which yields the required solution. By direct computation we can prove the converse. This ends the proof.

    Remark 2.1. In Lemma 2.6, note that the conditions Λ0 corresponds to non-resonance case.

    Corollary 2.1. (Special case: λ(t)=λ = constant) Let Λ0 and yC([0,T],R). Then the unique solution of the linear problem

    {RLDq(CDr+λ)x(t)=y(t),t[0,T],x(ξ)=αCDνx(η),x(T)=βIpx(ζ),ξ,η,ζ(0,T), (2.5)

    is given by

    x(t)=Iq+ry(t)λIrx(t)+λ1(t)(αIq+rνy(η)αλIrνx(η)Iq+ry(ξ)+λIrx(ξ))+λ2(t)(βIq+r+py(ζ)βλIp+rx(ζ)Iq+ry(T)+λIrx(T)).

    In order to transform the problem (1.1) into a fixed point problem, we introduce an operator A:CC by Lemma 2.6 as follows:

    (Ax)(t)=Iq+rf(s,x(s))(t)Ir(λ(t)x(t))+λ1(t)(αIq+rνf(s,x(s))(η)αIrν(λ(η)x(η))Iq+rf(s,x(s))(ξ)+Ir(λ(ξ)x(ξ)))+λ2(t)(βIq+r+pf(s,x(s))(ζ)βIr+p(λ(ζ)x(ζ))Iq+rf(s,x(s))(T)+Ir(λ(T)x(T))), (3.1)

    where C is the Banach space of continuous functions from [0,T] to R equipped with the norm x=supt[0,T]|x(t)|. Evidently the existence of solutions for the problem (1.1) is related to the existence of fixed points for the operator A.

    For computational convenience, we use the notations:

    Φ=Tq+rΓ(q+r+1)+ˉλ1(|α|ηq+rνΓ(q+rν+1)+ξq+rΓ(q+r+1))+ˉλ2(|β|ζq+r+pΓ(q+r+p+1)+Tq+rΓ(q+r+1)) (3.2)
    Φλ=Ir|λ(T)|+ˉλ1(|α|Irν|λ(η)|+Ir|λ(ξ)|)+ˉλ2(|β|Ip+r|λ(ζ)|+Ir|λ(T)|), (3.3)

    where

    ˉλ1=maxt[0,T]|λ1(t)|=1|Λ|[|Γ|Γ(q)Γ(q+r)Tq+r1+|B|],ˉλ2=maxt[0,T]|λ2(t)|=1|Λ|[Γ(q)Γ(q+r)Tq+r1+|A|]. (3.4)

    Remark 3.1. In the special case of constant function λ(t)=λ, Φλ becomes

    ¯Φλ=|λ|TrΓ(r+1)+ˉλ1|λ|(|α|ηrνΓ(rν+1)+ξrΓ(r+1))+ˉλ2|λ|(|β|ζr+pΓ(r+p+1)+TrΓ(r+1)).

    We are now in a position to give our first existence and uniqueness result, which relies on the contraction mapping principle due to Banach.

    Theorem 3.1. Let f:[0,T]×RR be a continuous function. Assume that:

    (H1) there exists a positive constant L such that

    |f(t,x)f(t,y)|L|xy|,t[0,T],x,yR.

    If

    LΦ+Φλ<1, (3.5)

    where Φ and Φλ are respectively given by (3.2) and (3.3), then the nonlocal-terminal value problem (1.1) has a unique solution on [0,T].

    Proof. The hypothesis of contraction mapping principle will be verified in two steps. In the first step, we show that ABrBr, where A is the operator defined by (3.1), Br={xC:xr}, rMΦ(1LΦΦλ)1, supt[0,T]|f(t,0)|=M<. For any xBr, and taking into account assumption (H1), we obtain

    Axsupt[0,T]{Iq+r|f(s,x(s))|(t)+|λ1(t)|(|α|Iq+rν|f(s,x(s))|(η)+Iq+r|f(s,x(s))|(ξ))+|λ2(t)|[|β|Iq+r+p|f(s,x(s))|(ζ)+Iq+r|f(s,x(s))|(T)]+Ir|λ(t)x(t)|+|λ1(t)|[|α|Irν|λ(η)x(η)|+Ir|λ(ξ)x(ξ)|]+|λ2(t)|[|β|Ir+p|λ(ζ)x(ζ)|+Ir|λ(T)x(T)|]}supt[0,T]{Iq+r(|f(s,x(s))f(s,0)|+|f(s,0)|)(t)+|λ1(t)|[|α|Iq+rν(|f(s,x(s))f(s,0)|+|f(s,0)|)(η)+Iq+r(|f(s,x(s))f(s,0)|+|f(s,0)|)(ξ)]+|λ2(t)|[|β|Iq+r+p(|f(s,x(s))f(s,0)|+|f(s,0)|)(ζ)+Iq+r(|f(s,x(s))f(s,0)|+|f(s,0)|)(T)]+Ir|λ(t)x(t)|+|λ1(t)|[|α|Irν|λ(η)x(η)|+Ir|λ(ξ)x(ξ)|]+|λ2(t)|[|β|Ir+p|λ(ζ)x(ζ)|+Ir|λ(T)x(T)|]}(Lr+M){Tq+rΓ(q+r+1)+ˉλ1[|α|ηq+rνΓ(q+rν+1)+ξq+rΓ(q+r+1)]+ˉλ2[|β|ζq+r+pΓ(q+r+p+1)+Tq+rΓ(q+r+1)]}+x{Ir|λ(T)|+ˉλ1[|α|Irν|λ(η)|+Ir|λ(ξ)|]+ˉλ2[|β|Ir+p|λ(ζ)|+Ir|λ(T)|]}=(Lr+M)Φ+rΦλr,

    which implies that ABrBr. In the next step, it will be shown that the operator A given by (3.1) is a contraction. For t[0,T] and for x,yC, we have

    |Ax(t)Ay(t)|Iq+r|f(s,x(s))f(s,y(s))|(t)+|λ1(t)|[|α|Iq+rν|f(s,x(s))|(η)+Iq+r|f(s,x(s))f(s,y(s))|(ξ)]+|λ2(t)|[|β|Iq+r+p|f(s,x(s))f(s,y(s))|(ζ)+Iq+r|f(s,x(s))f(s,y(s))|(T)]+Ir|λ(t)(x(t)y(t))|+|λ1(t)|[|α|Irν|λ(η)(x(η)y(η))|+Ir|λ(ξ)(x(ξ)y(ξ)|]+|λ2(t)|[|β|Ir+p|λ(ζ)(x(ζ)y(ζ))|+Ir|λ(T)(x(T)y(T))|]Lxy{Tq+rΓ(q+r+1)+ˉλ1[|α|ηq+rνΓ(q+rν+1)+ξq+rΓ(q+r+1)]+ˉλ2[|β|ζq+r+pΓ(q+r+p+1)+Tq+rΓ(q+r+1)]}+xy{Ir|λ(T)|+ˉλ1[|α|Irν|λ(η)|+Ir|λ(ξ)|]+ˉλ2[|β|Ir+p|λ(ζ)|+Ir|λ(T)|]}=(LΦ+Φλ)xy,

    which leads to AxAy(LΦ+Φλ)xy. As LΦ+Φλ<1, therefore A is a contraction. Since the hypothesis of Banach contraction mapping principle is satisfied, therefore we deduce by its conclusion that the operator A has a unique fixed point, which corresponds to a unique solution of the problem (1.1). The proof is completed.

    Next, we prove an existence result for the problem (1.1) by using Krasnoselskii fixed point theorem [18].

    Theorem 3.2. Assume that f:[0,T]×RR is a continuous function satisfying the assumptions:

    (H2) |f(t,x)|δ(t), (t,x)[0,T]×R, δC([0,T],R+) with δ=supt[0,T]|δ(t)|.

    (¯H2) Φλ<1, where Φλ is given by (3.4).

    Then the nonlocal-terminal value problem (1.1) has at least one solution on [0,T].

    Proof. Let us select a positive number ¯r such that ¯rδΦ(1Φλ)1 and define operators A1 and A2 on B¯r={xC:x¯r} as

    (A1x)(t)=Iq+rf(s,x(s))(t)+λ1(t)[αIq+rνf(s,x(s))(η)Iq+rf(s,x(s))(ξ)]+λ2(t)[βIq+r+pf(s,x(s))(ζ)Iq+rf(s,x(s))(T)],t[0,T],(A2x)(t)=Ir(λ(t)x(t))+λ1(t)[αIrν(λ(η)x(η))+Ir(λ(ξ)x(ξ))]+λ2(t)[βIr+p(λ(ζ)x(ζ))+Ir(λ(T)x(T))],t[0,T].

    Observe that Ax=A1x+A2x. For x,yB¯r, we have

    A1x+A2yδ{Tq+rΓ(q+r+1)+ˉλ1[|α|ηq+rνΓ(q+rν+1)+ξq+rΓ(q+r+1)]+ˉλ2[|β|ζq+r+pΓ(q+r+p+1)+Tq+rΓ(q+r+1)]}+y{Ir|λ(T)|+ˉλ1[|α|Irν|λ(η)|+Ir|λ(ξ)|]+ˉλ2[|β|Ir+p|λ(ζ)|+Ir|λ(T)|]}=δΦ+xΦλ¯r.

    This shows that A1x+A2yB¯r. With the aid of the assumption (¯H2), it is easy to show that A2 is a contraction. Further the operator A1 is continuous in view of continuity of f. Further A1 is uniformly bounded on B¯r as

    A1xδΦ.

    Next, we prove the compactness of the operator A1. Let us set sup(t,x)[0,T]×B¯r|f(t,x)|=¯f< and take t1,t2[0,T] with t1<t2. Then we have

    |(A1x)(t2)(A1x)(t1)|¯fΓ(q+r)|t10[(t2s)q+r1(t1s)q+r1ds|+¯fΓ(q+r)|t2t1(t2s)q+r1ds|+|Γ||tq+r12tq+r11||Λ|Γ(q)Γ(q+r)[|α|ηq+rνΓ(q+rν+1)+ξq+rΓ(q+r+1)]+|tq+r12tq+r11||Λ|Γ(q)Γ(q+r)[|β|ζq+r+pΓ(q+r+p+1)+Tq+rΓ(q+r+1)]¯fΓ(q+r+1)[|tq+r2tq+r1|+2(t2t1)q+r]+|Γ||tq+r12tq+r11||Λ|Γ(q)Γ(q+r)[|α|ηq+rνΓ(q+rν+1)+ξq+rΓ(q+r+1)]+|tq+r12tq+r11||Λ|Γ(q)Γ(q+r)[|β|ζq+r+pΓ(q+r+p+1)+Tq+rΓ(q+r+1)],

    which tends to zero as t2t1 independent of x. Thus, A1 is equicontinuous. So A1 is relatively compact on B¯r. Hence, by the Arzelá-Ascoli theorem, A1 is compact on B¯r. Thus the operators A1 and A2 satisfy the hypothesis of Krasnoselskii fixed point theorem [18]. Hence it follows by the conclusion of Krasnoselskii fixed point theorem [18] that the operator A(=A1+A2) has a fixed point, which corresponds to a solution of the problem (1.1) on [0,T]. The proof is completed.

    In the following result, we prove the existence of solutions for the problem (1.1) by means of Leray-Schauder nonlinear alternative [14].

    Theorem 3.3. Let f:[0,T]×RR be a continuous function satisfying the conditions:

    (H3) there exist a continuous nondecreasing functions ψ:[0,)(0,) and a function ϕC([0,T],R+) such that |f(t,x)|ϕ(t)ψ(|x|) for each (t,x)[0,T]×R;

    (H4) there exists a constant N>0 such that

    (1Φλ)NΦϕψ(N)>1,Φλ<1,

    where Φ and Φλ are respectively given by (3.2) and (3.3).

    Then there exists at least one solution for the nonlocal-terminal value problem (1.1) on [0,T].

    Proof. We verify the hypothesis of Leray-Schauder nonlinear alternative [14] in several steps. Let us first show that the operator A, defined by (3.1), maps bounded sets (balls) into bounded sets in C. For a positive number R, let BR={xC:xR} be a bounded ball in C. Then, for t[0,T], we have

    |Ax(t)|Iq+r|f(s,x(s))|(t)+ˉλ1[|α|Iq+rν|f(s,x(s))|(η)+Iq+r|f(s,x(s))|(ξ)]+ˉλ2[|β|Iq+r+p|f(s,x(s))|(ζ)+Iq+r|f(s,x(s))|(T)]+Ir|λ(t)x(t)|+ˉλ1[|α|Irν|λ(η)x(η)|+Ir|λ(ξ)x(ξ)|]+ˉλ2[|β|Ir+p|λ(ζ)x(ζ)|+Ir|λ(T)x(T)|]ϕ(t)ψ(x){Tq+rΓ(q+r+1)+ˉλ1[|α|ηq+rνΓ(q+rν+1)+ξq+rΓ(q+r+1)]+ˉλ2[|β|ζq+r+pΓ(q+r+p+1)+Tq+rΓ(q+r+1)]}+x{Ir|λ(T)|+ˉλ1[|α|Irν|λ(η)|+Ir|λ(ξ)|]+ˉλ2[|β|Ir+p|λ(ζ)|+Ir|λ(T)|]}.

    In view of (3.2) and (3.3), the above inequality takes the form:

    AxΦϕψ(R)+RΦλ.

    Secondly, we show that A maps bounded sets into equicontinuous sets of C. Let ν1,ν2[0,T] with ν1<ν2 and xBR. Then we have

    |(Ax)(ν2)(Ax)(ν1)|Iq+r(|f(s,x(s))(ν2)f(s,x(s))(ν1)|)+|Γ||νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|α|ηq+rνΓ(q+rν+1)+ξq+rΓ(q+r+1)]+|νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|β|ζq+r+pΓ(q+r+p+1)+Tq+rΓ(q+r+1)]+|ν10[(ν2s)r1(ν1s)r1]Γ(r)λ(s)x(s)ds+ν2ν1(ν2s)r1Γ(r)λ(s)x(s)ds|+|Γ||νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|α|Irν|λ(η)|+Ir|λ(ξ)|]+|νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|β|Ir+p|λ(ζ)|+Ir|λ(T)|]ϕψ(R)Γ(q+r+1)[|tq+r2tq+r2|+2(t2t1)q+r]+|Γ||νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|α|ηq+rνΓ(q+rν+1)+ξq+rΓ(q+r+1)]+|νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|β|ζq+r+pΓ(q+r+p+1)+Tq+rΓ(q+r+1)]+|ν10[(ν2s)r1(ν1s)r1]Γ(r)λ(s)x(s)ds+ν2ν1(ν2s)r1Γ(r)λ(s)x(s)ds|+|Γ||νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|α|Irν|λ(η)|+Ir|λ(ξ)|]+|νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|β|Ir+p|λ(ζ)|+Ir|λ(T)|].

    Obviously the right hand side of the above inequality tends to zero independently of xBR as ν2ν1. Therefore it follows by the Arzelá-Ascoli theorem that F:CC is completely continuous.

    In order to complete the hypothesis of the Leray-Schauder nonlinear alternative [14], it will be shown that the set of all solutions to the equation x=θAx is bounded for θ[0,1]. For that, let x be a solution of x=θAx for θ[0,1]. Then, for t[0,T], we apply the strategy used in the first step to obtain

    xΦϕψ(x)+xΦλ,

    where Φ and Φλ are respectively given by (3.3) and (3.4). Consequently, we have

    (1Φλ)xΦϕψ(x)1.

    By the condition (H4), we can find a positive number N such that xN. Introduce a set

    U={xC:x<N}, (3.7)

    and observe that the operator A:¯UC is continuous and completely continuous. With this choice of U, we cannot find xU satisfying the relation x=θAx for some θ(0,1). Therefore, it follows by nonlinear alternative of Leray-Schauder type [14] that the operator A has a fixed point in ¯U. Thus there exists a solution of the problem (1.1) on [0,T]. The proof is complete.

    We begin this section with the definition of a solution for the multi-valued problem (1.2).

    Definition 4.1. A function xC([0,T],R) is said to be a solution of the problem (1.2) if there exists a function vL1([0,T],R) with v(t)F(t,x) a.e. on [0,T] such that x(ξ)=λCDνx(η),x(T)=μIpx(ζ) and

    x(t)=Iq+rv(s)(t)Ir(λ(t)x(t))+λ1(t)[αIq+rνv(s)(η)αIrν(λ(η)x(η))Iq+r1v(s)(ξ)+Ir(λ(ξ)x(ξ))]+λ2(t)[βIq+r+pv(s)(ζ)βIr+p(λ(ζ)x(ζ))Iq+rv(s)(T)+Ir(λ(T)x(T))].

    Our first result, dealing with the convex-valued F, is based on Leray-Schauder nonlinear alternative for multi-valued maps.

    Definition 4.2. A multivalued map F:[0,T]×RP(R) is said to be Carathéodory if (i) tF(t,x) is measurable for each xR and (ii) xF(t,x) is upper semicontinuous for almost all t[0,T]. Further a Carathéodory function F is called L1Carathéodory if (iii) for each ρ>0, there exists φρL1([0,T],R+) such that F(t,x)=sup{|v|:vF(t,x)}φρ(t) for all xR with xρ and for a.e. t[0,T].

    Theorem 4.1. Assume that:

    (A1) F:[0,T]×RPcp,c(R) is L1-Carathéodory, where Pcp,c(R)={YP(R): Yiscompactandconvex};

    (A2) there exist a continuous nondecreasing function Q:[0,)(0,) and a function PC([0,T],R+) such that F(t,x)P:=sup{|y|:yF(t,x)}P(t)Q(|x|)for each(t,x) [0,T]×R;

    (A3) there exists a constant M>0 such that

    (1Φλ)MΦPQ(M)>1,Φλ<1,

    where Φ and Φλ are respectively given by (3.3) and (3.4).

    Then the nonlocal-terminal value problem (1.2) has at least one solution on [0,T].

    Proof. Firstly, we transform the problem (1.2) into a fixed point problem by defining a multi-valued map: N:C([0,T],R)P(C([0,T],R)) as

    N(x)={hC([0,T],R):h(t)={Iq+rv(s)(t)Ir(λ(t)x(t))+λ1(t)[αIq+rνv(s)(η)αIrν(λ(η)x(η))Iq+r1v(s)(ξ)+Ir(λ(ξ)x(ξ))]+λ2(t)[βIq+r+pv(s)(ζ)βIr+p(λ(ζ)x(ζ))Iq+rv(s)(T)+Ir(λ(T)x(T))],}

    for vSF,x.

    It is clear that fixed points of N are solutions of problem (1.2). So we need to verify that the operator N satisfies all the conditions of Leray-Schauder nonlinear alternative [14]. This will be done in several steps.

    Step 1. N(x) is convex for each xC([0,T],R).

    Indeed, if h1,h2 belongs to N(x), then there exist v1,v2SF,x such that, for each t[0,T], we have

    hi(t)=Iq+rvi(s)(t)Ir(λ(t)x(t))+λ1(t)[αIq+rνvi(s)(η)αIrν(λ(η)x(η))Iq+r1vi(s)(ξ)+Ir(λ(ξ)x(ξ))]+λ2(t)[βIq+r+pvi(s)(ζ)βIr+p(λ(ζ)x(ζ))Iq+rvi(s)(T)+Ir(λ(T)x(T))],

    i=1,2. Let 0θ1. Then, for each t[0,T], we have

    [θh1+(1θ)h2](t)=Iq+r[θv1(s)+(1θ)v2(s)](t)Ir(λ(t)x(t))+λ1(t)[αIq+rν[θv1(s)+(1θ)v2(s)](η)αIrν(λ(η)x(η))Iq+r1[θv1(s)+(1θ)v2(s)](ξ)+Ir(λ(ξ)x(ξ))]+λ2(t)[βIq+r+p[θv1(s)+(1θ)v2(s)](ζ)βIr+p(λ(ζ)x(ζ))Iq+r[θv1(s)+(1θ)v2(s)](T)+Ir(λ(T)x(T))].

    Since F has convex values (SF,x is convex), therefore, θh1+(1θ)h2N(x).

    Step 2. N(x) maps bounded sets (balls) into bounded sets in C([0,T],R).

    For a positive number r, let Br={xC([0,T],R):xr} be a bounded ball in C([0,T],R). Then, for each hN(x),xBr, there exists vSF,x such that

    h(t)=Iq+rv(s)(t)Ir(λ(t)x(t))+λ1(t)[αIq+rνv(s)(η)αIrν(λ(η)x(η))Iq+r1v(s)(ξ)+Ir(λ(ξ)x(ξ))]+λ2(t)[βIq+r+pv(s)(ζ)βIr+p(λ(ζ)x(ζ))Iq+rv(s)(T)+Ir(λ(T)x(T))].

    In view of (H2), for each t[0,T], we have

    |h(t)|Iq+r|v(s)|(t)+ˉλ1[|α|Iq+rν|v(s)|(η)+Iq+r|f(s,x(s))|(ξ)]+ˉλ2[|β|Iq+r+p|v(s)|(ζ)+Iq+r|v(s)|(T)]+Ir|λ(t)x(t)|+ˉλ1[|α|Irν|λ(η)x(η)|+Ir|λ(ξ)x(ξ)|]+ˉλ2[|β|Ir+p|λ(ζ)x(ζ)|+Ir|λ(T)x(T)|]P(t)Q(x){Tq+rΓ(q+r+1)+ˉλ1[|α|ηq+rνΓ(q+rν+1)+ξq+rΓ(q+r+1)]+ˉλ2[|β|ζq+r+pΓ(q+r+p+1)+Tq+rΓ(q+r+1)]}+x{Ir|λ(T)|+ˉλ1[|α|Irν|λ(η)|+Ir|λ(ξ)|]+ˉλ2[|β|Ir+p|λ(ζ)|+Ir|λ(T)|]},

    which yields

    hΦPQ(r)+rΦλ.

    Step 3. N(x) maps bounded sets into equicontinuous sets of C([0,T],R).

    Let x be any element in Br and hN(x). Then there exists a function vSF,x such that, for each t[0,T] we have

    h(t)=Iq+rv(s)(t)Ir(λ(t)x(t))+λ1(t)[αIq+rνv(s)(η)αIrν(λ(η)x(η))Iq+r1v(s)(ξ)+Ir(λ(ξ)x(ξ))]+λ2(t)[βIq+r+pv(s)(ζ)βIr+p(λ(ζ)x(ζ))Iq+rv(s)(T)+Ir(λ(T)x(T))].

    Let τ1, τ2[0,T], τ1<τ2. Then

    |h(τ2)h(τ1)|Iq+r(|v(s)(ν2)v(s)(ν1)|)+|Γ||νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|α|ηq+rνΓ(q+rν+1)+ξq+rΓ(q+r+1)]+|νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|β|ζq+r+pΓ(q+r+p+1)+Tq+rΓ(q+r+1)]+|ν10[(ν2s)r1(ν1s)r1]Γ(r)λ(s)x(s)ds+ν2ν1(ν2s)r1Γ(r)λ(s)x(s)ds|+|Γ||νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|α|Irν|λ(η)|+Ir|λ(ξ)|]+|νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|β|Ir+p|λ(ζ)|+Ir|λ(T)|]PQ(r)Γ(q+r+1)[|tq+r2tq+r2|+2(t2t1)q+r]+|Γ||νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|α|ηq+rνΓ(q+rν+1)+ξq+rΓ(q+r+1)]+|νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|β|ζq+r+pΓ(q+r+p+1)+Tq+rΓ(q+r+1)]+|ν10[(ν2s)r1(ν1s)r1]Γ(r)λ(s)x(s)ds+ν2ν1(ν2s)r1Γ(r)λ(s)x(s)ds|+|Γ||νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|α|Irν|λ(η)|+Ir|λ(ξ)|]+|νq+r12νq+r11||Λ|Γ(q)Γ(q+r)[|β|Ir+p|λ(ζ)|+Ir|λ(T)|].

    The right hand of the above inequality tends to zero independently of xBr as ν1ν2.

    Combining the outcome of Steps 1–3 with Arzelá-Ascoli theorem, we deduce that N:C([0,T],R)P(C([0,T],R)) is completely continuous.

    In order to prove that the operator N is u.s.c., it is enough to establish that it has a closed graph by Proposition 1.2 in [11], as it is already shown to be completely continuous. This is done in the next step.

    Step 4. N has a closed graph.

    Let xnx, hnN(xn) and hnh. We need to show that hN(x). Now hnN(xn) implies that there exists vnSF,xn such that, for each t[0,T],

    hn(t)=Iq+rv(s)(t)Ir(λ(t)x(t))+λ1(t)[αIq+rνvn(s)(η)αIrν(λ(η)x(η))Iq+r1vn(s)(ξ)+Ir(λ(ξ)x(ξ))]+λ2(t)[βIq+r+pvn(s)(ζ)βIr+p(λ(ζ)x(ζ))Iq+rvn(s)(T)+Ir(λ(T)x(T))].

    Therefore, we must show that there exists vSF,x such that, for each t[0,T],

    h(t)=Iq+rv(s)(t)Ir(λ(t)x(t))+λ1(t)[αIq+rνv(s)(η)αIrν(λ(η)x(η))Iq+r1v(s)(ξ)+Ir(λ(ξ)x(ξ))]+λ2(t)[βIq+r+pv(s)(ζ)βIr+p(λ(ζ)x(ζ))Iq+rv(s)(T)+Ir(λ(T)x(T))].

    Consider the continuous linear operator Θ:L1([0,T],R)C([0,T],R) given by

    vΘ(v)(t)=Iq+rv(s)(t)Ir(λ(t)x(t))+λ1(t)[αIq+rνv(s)(η)αIrν(λ(η)x(η))Iq+r1v(s)(ξ)+Ir(λ(ξ)x(ξ))]+λ2(t)[βIq+r+pv(s)(ζ)βIr+p(λ(ζ)x(ζ))Iq+rv(s)(T)+Ir(λ(T)x(T))].

    Observe that hn(t)h(t)0 as n. So it follows from by a closed graph result obtained in [20] that ΘSF,x is a closed graph operator. Moreover, we have hnΘ(SF,xn). Since xnx, therefore we have

    h(t)=Iq+rv(s)(t)Ir(λ(t)x(t))+λ1(t)[αIq+rνv(s)(η)αIrν(λ(η)x(η))Iq+r1v(s)(ξ)+Ir(λ(ξ)x(ξ))]+λ2(t)[βIq+r+pv(s)(ζ)βIr+p(λ(ζ)x(ζ))Iq+rv(s)(T)+Ir(λ(T)x(T))],

    for some vSF,x

    Step 5. We show there exists an open set VC([0,T],R) with xθN(x) for any θ(0,1) and all xV.

    Let θ(0,1) and xθN(x). Then there exists vL1([0,T],R) with vSF,x such that, for t[0,T], we have

    x(t)=θIq+rv(s)(t)θIr(λ(t)x(t))+θλ1(t)[αIq+rνv(s)(η)αIrν(λ(η)x(η))Iq+r1v(s)(ξ)+θIr(λ(ξ)x(ξ))]+θλ2(t)[βIq+r+pv(s)(ζ)βIr+p(λ(ζ)x(ζ))Iq+rv(s)(T)+Ir(λ(T)x(T))].

    Using the computations done in Step 2, for each t[0,T], we get

    |x(t)|ΦPQ(x)+xΦλ,

    which can alternatively be written as

    (1Φλ)xΦPQ(x)1.

    In view of (A3), there exists M such that xM. Let us define a set

    V={xC([0,T],R):x<M}.

    The operator N:¯VP(C([0,T], R)) is a compact multi-valued map, u.s.c. with convex closed values. With the given choice of V, it is not possible to find xV satisfying xθN(x) for some θ(0,1). In consequence, we deduce by the nonlinear alternative of Leray-Schauder type [14] that the operator N has a fixed point x¯V, which corresponds to a solution of the problem (1.2). This completes the proof.

    Let (X,d) be a metric space induced from the normed space (X;). Consider Hd:P(X)×P(X)R{} defined by Hd(A,B)=max{supaAd(a,B),supbBd(A,b)}, where d(A,b)=infaAd(a;b) and d(a,B)=infbBd(a;b). Then (Pb,cl(X), Hd) is a metric space (see [17]), where Pb,cl(X)={YP(X):Y is bounded and closed},

    Now we present our existence result for the problem (1.2) with a non-convex valued right hand side by applying a fixed point theorem for multivalued maps due to Covitz and Nadler [10]: ``If N:XPcl(X) is a contraction, then FixN, where Pcl(X)={YP(X):Y is closed}".

    Theorem 4.2. Assume that the following conditions hold:

    (A4) F:[0,T]×RPcp(R) is such that F(,x):[0,T]Pcp(R) is measurable for each xR, where Pcp(R)={YP(R):Yis compact};

    (A5) Hd(F(t,x),F(t,ˉx))ϱ(t)|xˉx| for almost all t[0,T] and x,ˉxR with ϱC([0,T],R+) and d(0,F(t,0))ϱ(t) for almost all t[0,T].

    Then the nonlocal-terminal value problem (1.2) has at least one solution on [0,T] if

    Φϱ+Φλ<1,

    where Φ and Φλ are respectively given by (3.3) and (3.4).

    Proof. Once it is shown that the operator N:C([0,T],R)P(C([0,T],R)), defined in the beginning of the proof of Theorem 4.1, satisfies the assumptions of Covitz and Nadler fixed point theorem [10], we are done. We establish in two steps.

    Step Ⅰ. N(x) is nonempty and closed for every vSF,x.

    Observe that the set-valued map F(,x()) admits a measurable selection v:[0,T]R as it is measurable by the measurable selection theorem (e.g., [8, Theorem Ⅲ.6]). Moreover, by the assumption (A5), we have

    |v(t)|ϱ(t)(1+|x(t)|),

    that is, vL1([0,T],R) and hence F is integrably bounded. Therefore, SF,x.

    Next we show that N(x) is closed for each xC([0,T],R). Let {un}n0N(x) be such that unu as n in C([0,T],R). Then uC([0,T],R) and there exists vnSF,xn such that, for each t[0,T],

    un(t)=Iq+rvn(s)(t)Ir(λ(t)x(t))+λ1(t)[αIq+rνvn(s)(η)αIrν(λ(η)x(η))Iq+r1vn(s)(ξ)+Ir(λ(ξ)x(ξ))]+λ2(t)[βIq+r+pvn(s)(ζ)βIr+p(λ(ζ)x(ζ))Iq+rvn(s)(T)+Ir(λ(T)x(T))].

    As F has compact values, we pass onto a subsequence (if necessary) to obtain that v_n converges to v in L^1 ({[0, T]}, \mathbb{R}). Thus v \in S_{F, x} and for each t \in {[0, T]} , we have

    \begin{eqnarray*} &&u_n(t)\to v(t)\\ & = & I^{q+r}v(s)(t)- I^r (\lambda(t)x(t))\nonumber\\ &&+\lambda_1(t)\Big[\alpha I^{q+r-\nu}v(s)(\eta)- \alpha I^{r-\nu}(\lambda(\eta)x(\eta))-I^{q+r-1}v(s)(\xi)+ I^{r}(\lambda(\xi)x(\xi))\Big] \\ &&+\lambda_2(t)\Big[\beta I^{q+r+p}v(s)(\zeta)-\beta I^{r+p}(\lambda(\zeta)x(\zeta))-I^{q+r}v(s)(T)+ I^r (\lambda(T)x(T))\Big]. \end{eqnarray*}

    Hence u \in N(x).

    Step Ⅱ. We show that there exists 0 < \hat \theta < 1 ( \hat\theta = \Phi \|\varrho\|+\Phi_{\lambda} ) such that

    H_d(N(x),N(\bar{x}))\le \hat\theta \|x-\bar{x}\| \; \; \mbox{for each} \; \; x, \bar{x} \in C({[0,T]}, \mathbb{R}).

    Let x, \bar{x} \in C({[0, T]}, \mathbb{R}) and h_1 \in N(x) . Then there exists v_1(t) \in F(t, x(t)) such that, for each t \in {[0, T]} ,

    \begin{eqnarray*} h_1(t)& = & I^{q+r}v_1(s)(t)- I^r (\lambda(t)x(t))\nonumber\\ &&+\lambda_1(t)\Big[\alpha I^{q+r-\nu}v_1(s)(\eta)- \alpha I^{r-\nu}(\lambda(\eta)x(\eta))-I^{q+r-1}v_1(s)(\xi)+ I^{r}(\lambda(\xi)x(\xi))\Big] \\ &&+\lambda_2(t)\Big[\beta I^{q+r+p}v_1(s)(\zeta)-\beta I^{r+p}(\lambda(\zeta)x(\zeta))-I^{q+r}v_1(s)(T)+ I^r (\lambda(T)x(T))\Big]. \end{eqnarray*}

    By (A_5) , we have

    H_d(F(t,x), F(t,\bar{x}))\le \varrho(t)|x(t)-\bar{x}(t)|.

    So, there exists w(t) \in F(t, \bar{x}(t)) such that

    |v_1(t)-w|\le \varrho(t)|x(t)-\bar{x}(t)|, \; \; t \in {[0,T]}.

    Define {\mathcal U} : {[0, T]} \to \mathcal{P}(\mathbb{R}) by

    {\mathcal U}(t) = \{w \in \mathbb{R} : |v_1(t)-w|\le \varrho(t)|x(t)-\bar{x}(t)|\}.

    As the multivalued operator {\mathcal U}(t)\cap F(t, \bar{x}(t)) is measurable (see Proposition Ⅲ.4 [8]), there exists a function v_2(t) which is a measurable selection for {\mathcal U} . So v_2(t) \in F(t, \bar{x}(t)) and for each t \in {[0, T]} , we have |v_1(t)-v_2(t)|\le \varrho(t)|x(t)-\bar{x}(t)| .

    For each t \in {[0, T]} , let us define

    \begin{eqnarray*} h_2(t)& = & I^{q+r}v_2(s)(t)- I^r (\lambda(t)x(t))\nonumber\\ &&+\lambda_1(t)\Big[\alpha I^{q+r-\nu}v_2(s)(\eta)- \alpha I^{r-\nu}(\lambda(\eta)x(\eta))-I^{q+r-1}v_2(s)(\xi)+ I^{r}(\lambda(\xi)x(\xi))\Big] \\ &&+\lambda_2(t)\Big[\beta I^{q+r+p}v_2(s)(\zeta)-\beta I^{r+p}(\lambda(\zeta)x(\zeta))-I^{q+r}v_2(s)(T)+ I^r (\lambda(T)x(T))\Big]. \end{eqnarray*}

    In consequence, we get

    \begin{eqnarray*} \left|h_1(t)-h_2(t)\right| &\le& I^{q+r}|v_2(s)-v_1(s)|(t)+I^r (\lambda(t)|x(t)-\bar x(t)|)\nonumber\\ &&+\lambda_1(t)\Big[|\alpha| I^{q+r-\nu}|v_2(s)-v_1(s)|(\eta)+ |\alpha| I^{r-\nu}(\lambda(\eta)|x(\eta)-\bar x(\eta)|)\nonumber\\ &&+I^{q+r}|v_2(s)-v_1(s)|(\xi)+ I^{r-1}(\lambda(\xi)|x(\xi)-\bar x(\xi)|)\Big]\\ &&+\lambda_2(t)\Big[|\beta| I^{q+r+p}|v_2(s)-v_1(s)|(\zeta)+|\beta| I^{r+p}(\lambda(\zeta)|x(\zeta)-\bar x(\zeta)|)\\ &&+I^{q+r}|v_2(s)-v_1(s)|(T)+ I^r (\lambda(T)x(T))\Big]\\ &\le&(\Phi\|\varrho\|+\Phi_{\lambda})\|x-\overline{x}\|. \end{eqnarray*}

    Hence

    \| h_1-h_2\| \le (\Phi\|\varrho\|+\Phi_{\lambda})\|x-\overline{x}\|.

    Analogously, interchanging the roles of x and \overline{x} , we obtain

    H_d(N(x), N(\bar{x})) \le (\Phi\|\varrho\|+\Phi_{\lambda})\|x-\overline{x}\|.

    Therefore, N is a contraction in view of the given condition( (\Phi\|\varrho\|+\Phi_{\lambda}) < 1 ). Hence it follows by the conclusion of Covitz and Nadler fixed point theorem [10] that N has a fixed point x , which corresponds to a solution of (1.2). The proof is complete.

    Consider the following nonlocal-terminal value problem

    \begin{equation} \left\{\begin{array}{ll} & {}^{RL}D^{2/3}\Big({}^{C}D^{3/4}+t^2/2\Big)x(t) = f(t,x(t)), \; \; t\in [0,1],\\[0.3cm] & x(2/3) = 5\; ^C D^{1/2}x(1/3),\; \; \; x(1) = 4\; I^1x(4/5). \end{array}\right. \end{equation} (5.1)

    Here q = 2/3, r = 3/4, \xi = 2/3, \eta = 1/3, \nu = 1/2, \alpha = 5, \beta = 4, \zeta = 4/5, p = 1, T = 1. Using the given data in (2.2), we find that A\approx -5.739178, B \approx -1.616533, \Gamma \approx -1.915888, and \Lambda \approx 12.612155. From (3.5), we have \bar\lambda_1 \approx 0.360215 and \bar\lambda_2 \approx 0.576166. Inserting the given values in (3.3) and (3.4), we obtain \Phi \approx 2.538503 and \Phi_\lambda \approx 0.502890.

    In order to illustrate Theorem 3.1, we take

    \begin{equation} f(t, x) = (L/2)(\tan^{-1}x+|x|(1+|x|)^{-1}) +\sqrt{t^4+3}, \end{equation} (5.2)

    in (5.1), where L is to be fixed later. Obviously |f(t, x)-f(t, y)|\leq L |x-y|, \, t\in [0, 1], \, x, y \in {\mathbb R} and L\Phi+\Phi_{\lambda} < 1 is satisfied for L < 0.195828. Thus all the assumptions of Theorem 3.1 are satisfied. Hence, by the conclusion of Theorem 3.1, the problem (5.1) with f(t, x) given by (5.2) has a unique solution on [0, 1] .

    Also the conclusion of Theorem 3.2 applies to the problem (5.1) with f(t, x) given by (5.2) as |f(t, x)|\le [\pi+2 +\sqrt{t^4+3}] = \delta(t), where L = 4.

    For the illustration of Theorem 3.3, let us take

    \begin{equation} f(t, x) = \frac{e^{-t^2}}{\sqrt{t^4+100}}\Big(\sin x +\frac{1}{10}\Big), \end{equation} (5.3)

    and note that |f(t, x)| \le \phi(t) \psi(|x|), where \phi(t) = e^{-t^2}/\sqrt{t^4+100}, and \psi(|x|) = (|x|+1/10). So the assumption (H_3) holds. Moreover, there exists N > 0.1043536 satisfying (H_4) . Thus all the conditions of Theorem 3.3 are satisfied and consequently the problem (5.1) with f(t, x) given by (5.3) has at least one solution on [0, 1].

    Now we illustrate Theorem 4.2 by considering the following multivalued problem:

    \begin{equation} \left\{\begin{array}{ll} & {}^{RL}D^{2/3}\Big({}^{C}D^{3/4}+t^2/2\Big)x(t)\in F(t,x), \; \; t\in [0,1],\\[0.3cm] & x(2/3) = 5\; ^C D^{1/2}x(1/3),\; \; \; x(1) = 4\; I^1x(4/5), \end{array}\right. \end{equation} (5.4)

    where

    F(t,x) = \Bigg[\frac{e^{-t}}{19+t},\,\,\, \frac{x+\tan^{-1}(x)}{(t+4)^2}+\frac{1}{15}\Bigg].

    Clearly H_d(F(t, x), F(t, \bar{x}))\leq \varrho(t)|x-\bar{x}|, where \varrho(t) = \frac{2}{(t+4)^2}. Also d(0, F(t, 0))\le \varrho(t) for almost all t \in [0, 1] and \Phi \|\varrho\|+\Phi_{\lambda} \approx 0.820203 < 1. As the hypothesis of Theorem 4.2 is satisfied, therefore we conclude that the multivalued problem (5.4) has at least one solution on [0, 1].

    We have studied a new nonlocal-terminal value problem consisting of Langevin equation with variable coefficient involving both Riemann-Liouville and Caputo fractional derivatives, and equipped with nonlocal-terminal fractional integro-differential conditions. The results presented in this paper are new and enrich the existing literature on boundary value problems of Langevin equation.

    As an analogue of the problem (1.1), we interchange the role of Riemann-Liouville and Caputo fractional derivatives in the Langevin equation given by (1.1) and consider the following nonlocal boundary value problem:

    \begin{equation} \left\{\begin{array}{ll} & {}^{C}D^r({}^{RL}D^q +\lambda(t))x(t) = f(t,x(t)), \; \; t\in [0,T],\\[0.3cm] & x(0) = 0,\; \; \; x(T) = \beta\; I^px(\zeta),\; \; \zeta\in (0,T). \end{array}\right. \end{equation} (6.1)

    As argued in the proof of Lemma 2.6, the solution of the Langevin equation in (6.1) can be written as

    \begin{equation} x(t) = I^{q+r}f(t,x(t))- I^{q}(\lambda(t)x(t))+c_1 \frac{t^q}{\Gamma(q+1)}+c_2t^{q-1}. \end{equation} (6.2)

    Using the condition x(0) = 0 in (6.2) implies that c_2 = 0 . Inserting the value of c_2 in (6.2) and then using the resulting expression for x(t) in the condition x(T) = \beta\; I^px(\zeta) , we find that

    \begin{eqnarray} c_1 = \frac{1}{\Omega}\Big(\beta I^{q+r+p}f(\zeta,x(\zeta))- \beta I^{q+p}(\lambda(\zeta)x(\zeta)) \\-I^{r+q}f(T,x(T))+I^q(\lambda(T)x(T))\Big), \end{eqnarray} (6.3)

    where

    \Omega = \frac{T^q}{\Gamma(1+q)}-\frac{\beta\zeta^{q+p}}{\Gamma(q+p+1)}\ne 0.

    Thus the solution of the problem (6.1) is

    \begin{eqnarray*} \label{s-1} x(t)& = &I^{r+q}f(s,x(s))(t)-I^q(\lambda(t)x(t))\\ &&+\frac{t^q}{\Omega \Gamma(1+q)}\Big(\beta I^{q+r+p}f(\zeta,x(\zeta))- \beta I^{q+p}(\lambda(\zeta)x(\zeta)) \nonumber \\&&-I^{r+q}f(T,x(T))+I^q(\lambda(T)x(T))\Big). \end{eqnarray*}

    One can notice that the solution (6.2) becomes unbounded at t = 0 in view of the values of q \in (0, 1) , in contrast to the problem (1.1). So we impose the condition x(0) = 0 to ensure the boundedness of the solution of the Langevin equation in (6.1). It is equivalent to saying that the problem (1.1) is now well-posed one. The existence results for the problem (6.1), analogue to the ones for (1.1), can be obtained in a similar manner.

    Special cases. We can obtain several new results as special cases of the work presented in this paper by fixing the values of parameters involved in the problem at hand, which are listed below.

    ● By taking \alpha = 0, our results correspond to the ones for nonlocal-terminal fractional integral conditions: x(\xi) = 0, x(T) = \beta\; I^px(\zeta), \xi, \zeta\in (0, T) .

    ● The results of this paper reduce to the ones with boundary conditions of the form: x(\xi) = \alpha\; ^C D^{\nu}x(\eta), \; x(T) = 0, \; \xi, \eta \in (0, T) by fixing \beta = 0.

    ● Letting \alpha = 0, \beta = 0 in the results of this paper, we obtain the ones associated with the nonlocal-terminal conditions: x(\xi) = 0, x(T) = 0 .

    The authors thank the reviewers for their useful remarks on our paper.

    All authors declare no conflicts of interest in this paper.



    [1] L. V. Bertalanffy, Quantitative laws in metabolism and growth, Quarterly Rev. Biol., 32 (1957), 217–231.
    [2] M. Kühleitner, N. Brunner, W. Nowak, K. Renner-Martin, K. Scheicher, Best-fitting growth curves of the von Bertalanffy-Pütter type, Poultry Sci., 98 (2019), 3587–3592. https://doi.org/10.3382/ps/pez122 doi: 10.3382/ps/pez122
    [3] P. Román-Román, D. Romero, F. Torres-Ruiz, A diffusion process to model generalized von Bertalanffy growth patterns: fitting to real data, J. Theoret. Biol., 263 (2010), 59–69. https://doi.org/10.1016/j.jtbi.2009.12.009 doi: 10.1016/j.jtbi.2009.12.009
    [4] J. Calatayud, T. Caraballo, J. C. Cortés, M. Jornet, Mathematical methods for the randomized non-autonomous Bertalanffy model, Electron. J. Differ. Equat., 2020, 50.
    [5] M. P. Edwards, R. S. Anderssen, Symmetries and solutions of the non-autonomous von Bertalanffy equation, Commun. Nonlinear Sci. Numer. Simul., 22 (2015), 1062–1067. https://doi.org/10.1016/j.cnsns.2014.08.033 doi: 10.1016/j.cnsns.2014.08.033
    [6] J. Brzdęk, D. Popa, I. Rașa, B. Xu, Ulam stability of operators, Mathematical analysis and its applications, Academic Press, London, 2018.
    [7] M. Onitsuka, Hyers–Ulam stability of first order linear differential equations of Carathéodory type and its application, Appl. Math. Lett., 90 (2019), 61–68. https://doi.org/10.1016/j.aml.2018.10.013 doi: 10.1016/j.aml.2018.10.013
    [8] M. Onitsuka, T. Shoji, Hyers–Ulam stability of first-order homogeneous linear differential equations with a real-valued coefficient, Appl. Math. Lett., 63 (2017), 102–108. http://dx.doi.org/10.1016/j.aml.2016.07.020 doi: 10.1016/j.aml.2016.07.020
    [9] R. Fukutaka, M. Onitsuka, Best constant in Hyers–Ulam stability of first-order homogeneous linear differential equations with a periodic coefficient, J. Math. Anal. Appl., 473 (2019), 1432–1446. https://doi.org/10.1016/j.jmaa.2019.01.030 doi: 10.1016/j.jmaa.2019.01.030
    [10] R. Fukutaka, M. Onitsuka, A necessary and sufficient condition for Hyers–Ulam stability of first-order periodic linear differential equations, Appl. Math. Lett., 100 (2020), 106040. https://doi.org/10.1016/j.aml.2019.106040 doi: 10.1016/j.aml.2019.106040
    [11] D. Popa, I. Rașa, On the Hyers–Ulam stability of the linear differential equation, J. Math. Anal. Appl., 381 (2011), 530–537. https://doi.org/10.1016/j.jmaa.2011.02.051 doi: 10.1016/j.jmaa.2011.02.051
    [12] G. Wang, M. Zhou, L. Sun, Hyers–Ulam stability of linear differential equations of first order, Appl. Math. Lett., 21 (2008), 1024–1028. https://doi.org/10.1016/j.aml.2007.10.020 doi: 10.1016/j.aml.2007.10.020
    [13] A. Zada, O. Shah, R. Shah, Hyers–Ulam stability of non-autonomous systems in terms of boundedness of Cauchy problems, Appl. Math. Comput., 271 (2015), 512–518. https://doi.org/10.1016/j.amc.2015.09.040 doi: 10.1016/j.amc.2015.09.040
    [14] D. Dragičević, Hyers–Ulam stability for a class of perturbed Hill's equations, Results Math., 76 (2021). https://doi.org/10.1007/s00025-021-01442-1 doi: 10.1007/s00025-021-01442-1
    [15] R. Fukutaka, M. Onitsuka, Best constant for Ulam stability of Hill's equations, Bull. Sci. Math., 163 (2020), 102888. https://doi.org/10.1016/j.bulsci.2020.102888 doi: 10.1016/j.bulsci.2020.102888
    [16] M. Onitsuka, Hyers–Ulam stability for second order linear differential equations of Carathéodory type, J. Math. Inequal., 15 (2021), 1499–1518. https://doi.org/10.7153/jmi-2021-15-103 doi: 10.7153/jmi-2021-15-103
    [17] A. Akgül, A. Cordero, J. R. Torregrosa, A fractional Newton method with 2\alphath-order of convergence and its stability, Appl. Math. Lett., 98 (2019), 344–351. https://doi.org/10.1016/j.aml.2019.06.028 doi: 10.1016/j.aml.2019.06.028
    [18] N. Bouteraa, M. Inc, A. Akgül, Stability analysis of time-fractional differential equations with initial data, Math. Methods Appl. Sci., https://doi.org/10.1002/mma.7782
    [19] D. Dragičević, Hyers–Ulam stability for nonautonomous semilinear dynamics on bounded intervals, Mediterr. J. Math., 18 (2021), 71. https://doi.org/10.1007/s00009-021-01729-1 doi: 10.1007/s00009-021-01729-1
    [20] J. Huang, S-M. Jung, Y. Li, On Hyers–Ulam stability of nonlinear differential equations, Bull. Korean Math. Soc., 52 (2015), 685–697. https://doi.org/10.4134/BKMS.2015.52.2.685 doi: 10.4134/BKMS.2015.52.2.685
    [21] I. A. Rus, Ulam stability of ordinary differential equations, Stud. Univ. Babeș-Bolyai Math., 54 (2009), 125–133.
    [22] I. A. Rus, Ulam stabilities of ordinary differential equations in a Banach space, Carpathian J. Math., 26 (2010), 103–107.
    [23] M. Choubin, H. Javanshiri, A new approach to the Hyers–Ulam–Rassias stability of differential equations, Results Math., 76 (2021), 11. https://doi.org/10.1007/s00025-020-01318-w doi: 10.1007/s00025-020-01318-w
    [24] S.-M. Jung, A fixed point approach to the stability of differential equations y' = F(x, y), Bull. Malays. Math. Sci. Soc., 33 (2010), 47–56.
    [25] R. Murali, C. Park, A. Ponmana Selvan, Hyers–Ulam stability for an nth order differential equation using fixed point approach, J. Appl. Anal. Comput., 11 (2021), 614–631. https://doi.org/10.11948/20190093 doi: 10.11948/20190093
    [26] D. Popa, I. Rașa, A. Viorel, Approximate solutions of the logistic equation and Ulam stability, Appl. Math. Lett., 85 (2018), 64–69. https://doi.org/10.1016/j.aml.2018.05.018 doi: 10.1016/j.aml.2018.05.018
    [27] M. Onitsuka, Conditional Ulam stability and its application to the logistic model, Appl. Math. Lett., 122 (2021), 107565. https://doi.org/10.1016/j.aml.2021.107565 doi: 10.1016/j.aml.2021.107565
  • This article has been cited by:

    1. Bounmy Khaminsou, Chatthai Thaiprayoon, Jehad Alzabut, Weerawat Sudsutad, Nonlocal boundary value problems for integro-differential Langevin equation via the generalized Caputo proportional fractional derivative, 2020, 2020, 1687-2770, 10.1186/s13661-020-01473-7
    2. Mohamed I. Abbas, Maria Alessandra Ragusa, Solvability of Langevin equations with two Hadamard fractional derivatives via Mittag–Leffler functions, 2021, 0003-6811, 1, 10.1080/00036811.2020.1839645
    3. Choukri Derbazi, Hadda Hammouche, Caputo-Hadamard fractional differential equations with nonlocal fractional integro-differential boundary conditions via topological degree theory, 2020, 5, 2473-6988, 2694, 10.3934/math.2020174
    4. Shahram Rezapour, Bashir Ahmad, Sina Etemad, On the new fractional configurations of integro-differential Langevin boundary value problems, 2021, 60, 11100168, 4865, 10.1016/j.aej.2021.03.070
    5. Bashir Ahmad, Badrah Alghamdi, Ahmed Alsaedi, Sotiris K. Ntouyas, Existence results for Riemann-Liouville fractional integro-differential inclusions with fractional nonlocal integral boundary conditions, 2021, 6, 2473-6988, 7093, 10.3934/math.2021416
    6. John R. Graef, Kadda Maazouz, Moussa Daif Allah Zaak, A Comprehensive Study of the Langevin Boundary Value Problems with Variable Order Fractional Derivatives, 2024, 13, 2075-1680, 277, 10.3390/axioms13040277
    7. Fang Li, Ling Yang, Huiwen Wang, Analytic Solutions for Hilfer Type Fractional Langevin Equations with Variable Coefficients in a Weighted Space, 2024, 13, 2075-1680, 284, 10.3390/axioms13050284
    8. Keyu Zhang, Jiafa Xu, Solvability for a system of Hadamard-type hybrid fractional differential inclusions, 2023, 56, 2391-4661, 10.1515/dema-2022-0226
    9. Hamid Lmou, Omar Talhaoui, Ahmed Kajouni, Sina Etemad, Raaid Alubady, Existence theory on the Caputo-type fractional differential Langevin hybrid inclusion with variable coefficient, 2024, 2024, 1687-2770, 10.1186/s13661-024-01975-8
  • 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(2327) PDF downloads(81) Cited by(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog