Open access peer-reviewed chapter - ONLINE FIRST

A Study of Nonlinear Boundary Value Problem

By Noureddine Bouteraa and Habib Djourdem

Reviewed: September 16th 2021Published: November 3rd 2021

DOI: 10.5772/intechopen.100491

Downloaded: 28


In this chapter, firstly we apply the iterative method to establish the existence of the positive solution for a type of nonlinear singular higher-order fractional differential equation with fractional multi-point boundary conditions. Explicit iterative sequences are given to approximate the solutions and the error estimations are also given. Secondly, we cover the multi-valued case of our problem. We investigate it for nonconvex compact valued multifunctions via a fixed point theorem for multivalued maps due to Covitz and Nadler. Two illustrative examples are presented at the end to illustrate the validity of our results.


  • Positive solution
  • Uniqueness
  • Iterative sequence
  • Green’s function
  • Fractional differential equation and inclusion
  • Existence
  • Nonlocal boundary value problem
  • Fixed point theorem

1. Introduction

In this chapter, we are interested in the existence of solutions for the nonlinear fractional boundary value problem (BVP)


We also cover the multi-valued case of problem


where D0+α,D0+βare the standard Riemann-Liouville fractional derivative of order


where D0+α,D0+βare the stantard Riemann-Liouville fractional derivative of order αn1n,β1n2forn3, the function fC01×RR, the multifunction F:01×R2Rare allowed to be singular at t=0and/or t=1and ajR+,j=1,2,,p,0<η1<η2<<ηp<1,forpN.

The first definition of fractional derivative was introduced at the end of the nineteenth century by Liouville and Riemann, but the concept of non-integer derivative and integral, as a generalization of the traditional integer order differential and integral calculus, was mentioned already in 1695 by Leibniz [1] and L’Hospital [2]. In fact, fractional derivatives provide an excellent tool for the description of memory and hereditary properties of various materials and processes. The mathematical modeling of systems and processes in the fields of physics, chemistry, aerodynamics, electrodynamics of complex medium, polymer rheology, Bode’s analysis of feedback amplifiers, capacitor theory, electrical circuits, electro-analytical chemistry, biology, control theory, fitting of experimental data, involves derivatives of fractional order. In consequence, the subject of fractional differential equations is gaining much importance and attention. For more details we refer the reader to [1, 2, 3, 4, 5, 6] and the references cited therein.

Boundary value problems for nonlinear differential equations arise in a variety of areas of applied mathematics, physics and variational problems of control theory. A point of central importance in the study of nonlinear boundary value problems is to understand how the properties of nonlinearity in a problem influence the nature of the solutions to the boundary value problems. The multi-point boundary conditions are important in various physical problems of applied science when the controllers at the end points of the interval (under consideration) dissipate or add energy according to the sensors located, at intermediate points, see [7, 8] and the references therein. We quote also that realistic problems arising from economics, optimal control, stochastic analysis can be modeled as differential inclusion. The study of fractional differential inclusions was initiated by El-Sayed and Ibrahim [9]. Also, recently, several qualitative results for fractional differential inclusion were obtained in [10, 11, 12, 13] and the references therein.

The techniques of nonlinear analysis, as the main method to deal with the problems of nonlinear differential equations (DEs), nonlinear fractional differential equations (FDEs), nonlinear partial differential equations (PDEs), nonlinear fractional partial differential equations (FPDEs), nonlinear stochastic fractional partial differential equations (SFPDEs), plays an essential role in the research of this field, such as establishing the existence, uniqueness and multiplicity of solutions (or positive solutions) and mild solutions for nonlinear of different kinds of FPDEs, FPDEs, SFPDEs, inclusion differential equations and inclusion fractional differential equations with various boundary conditions, by using different techniques (approaches). For more details, see [14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37] and the references therein. For example, iterative method is an important tool for solving linear and nonlinear Boundary Value Problems. It has been used in the research areas of mathematics and several branches of science and other fields. However, Many authors showed the existence of positive solutions for a class of boundary value problem at resonance case. Some recent devolopment for resonant case can be found in [38, 39]. Let us cited few papers. In [40], the authors studied the boundary value problems of the fractional order differential equation:


where 1<α2,0<η<1,0<a,β<1,fC01×R2Rand D0+α,D0+βare the stantard Riemann-Liouville fractional derivative of order α. They obtained the multiple positive solutions by the Leray-Schauder nonlinear alternative and the fixed point theorem on cones.

In 2020Li et al. [41] consider the existence of a positive solution for the following BVP of nonlinear fractional differential equation with integral boundary conditions:


where 2<q3, 0<σ1, α,γ,δ0, and β>0satisfying 0<ρα+βγ+αδΓ2σ<βγ+δΓqΓqσ, f:01×0+0+and hii=12:010+are continuous. To obtain the existence results, the authors used the well-known GuoKrasnoselskiis fixed point theorem.

In 2017, Rezapour et al. [42] investigated a Caputo fractional inclusion with integral boundary condition for the following problem


where 1<α2,η,ν,β01,F:01×R×R×R2Ris a compact valued multifunction and cDαdenotes the Caputo fractional derivative of order α.

In 2018, Bouteraa and Benaicha [10] studied the existence of solutions for the Caputo fractional differential inclusion


subject to three-point boundary conditions


where 2<α3,1<p2,0<η<1,β,γR+,F:01×R×R2Ris a compact valued multifunction and cDαdenotes the Caputo fractional derivative of order α.

In 2019, Ahmad et al. [43] investigated the existence of solutions for the boundary value problem of coupled Caputo (Liouville-Caputo) type fractional differential inclusions:


subject to the coupled boundary conditions:


where CDα, CDβdenote the Caputo fractional derivatives of order αand βrespectively, F,G:0T×R×Rare given multivalued maps, PRis the family of all nonempty subsets of R, and νi,μi,i=1,2are real constants with νiμi1,i=1,2.

Inspired and motivated by the works mentioned above, we focus on the uniqueness of positive solutions for the nonlocal boundary value problem (1) with the iterative method and properties of ftu, explicit iterative sequences are given to approximate the solutions and the error estimations are also given. We also cover the multi-valued case of problem (2) when the right-hand side is nonconvex compact valued multifunctions via a fixed point theorem for multivalued maps due to Covitz and Nadler.

The chapter is organized as follows. In Section 2, we present some notations and lemmas that will be used to prove our main results of problem (1) and we discuss the uniqueness of problem (1). Finally, we give an example to illustrate our result. In Section 3, we introduce some definitions and preliminary results about essential properties of multifunction that will be used in the remainder of the chapter and we present existence results for the problem (2) when the right-hand side is a non-convex compact multifunction. We shall use the fixed point theorem for contraction multivalued maps due to Covitz and Nadler [44] to prove the uniqueness of solution of problem (2). Finally, we give an example to ascertain the main result.


2. Existence and uniqueness results for problem (2)

2.1 Preliminaries

In this section, we recall some definitions and facts which will be used in the later analysis. These details can be found in the recent literature; see [2, 4, 6, 45, 46, 47] and the references therein.

Let ACi01Rdenote the space of itimesdifferentiable functions u:01Rwhose ithderivative uiis absolutely continuous and αdonotes the integer part of number α.

Definition 2.1. Let α>0,n1<α<n,n=α+1and uACn0R.

The Caputo derivative of fractional order αfor the function u:0+Ris defined by


The Riemann-Liouville fractional derivative order αfor the function u:0+Ris defined by


provided that the right hand side is pointwise defined in 0and the function Γ:0R, defined by


is called Euler’s gamma function.

Definition 2.2. The Riemann-Liouville fractional integral of order α>0of a function u:0Ris given by


provided that the right hand side is pointwise defined in 0.

We recall in the following lemma some properties involving Riemann-Liouville fractional integral and Riemann-Liouville fractional derivative or Caputo fractional derivative which are need in Lemma 2.4.

Lemma 2.1. (([45], Prop.4.3), [46]) Letα,β0anduL101. Then the next formulas hold.

  1. DβIαut=Iαβut,

  2. DαIαut=ut,

  3. I0+αI0+βut=I0+α+βut.

  4. Ifβ>α>0, thenDαtβ1=Γβtβα1Γβα. whereDαandDβrepresents Riemann-Liouville’s or Caputo’s fractional derivative of orderαandβrespectively.

Lemma 2.2 [47]. Letα>0andyL101. Then, the general solution of the fractional differential equationD0+αut+yt=0,0<t<1is given by


wherec0,c1,,cn1are real constants andn=α+1.

Based on the previous Lemma 2.2, we will define the integral solution of our problem 1.

Lemma 2.3. Letpj=1ajηjαβ101,αn1n,β1n2,n3andyC01. Then the solution of the fractional boundary value problem


is given by





Proof. By using Lemma 2.2, the solution of the equation D0+αut+yt=0is


where c1,c2,cnare arbitrary real constants.

From the boundary condition in (1), one can c2=c3=cn2=cn1=cn=0. Hence


By the last above equation and Lemma 2.1i, we get


this and by D0+βu1=j=1pajD0+βuηj, we have


Then, the unique solution of the problem (1) is given by


The proof is completed.

Lemma 2.4. Letj=1pajηjαβ101,αn1n,β1n2,n3. Then, the functionsgtsandhtsdefined by(6)and(7)have the following properties:

  1. The functionsgtsandhtsare continuous on01×01and for allt,s01


  • gtstα1Γαfor allt,s01.

  • gtstα1g1sfor allt,s01, where


    From the above properties, we deduce the following properties:

  • The functionGts0is continuous on01×01andGts>0for allt,s01.

  • maxt01Gts=G1s, for alls01, where

  • G1s=g1s+1dpj=1ajhηjs1sαβ1dΓα.

    Proof. It is easy to chek that i,v,viholds. So we prove that iiis true. Note that (6) and 01sαβ11. It follows that gtstα1Γαfor all t,s01. It remains to prove iii. We divide the proof into two cases and by (1), we have.

    Case1. When 0st1, we have


    Case2. When 0ts1, we have


    Hence gtstα1g1sfor all t,s01.

    2.2 Existence results

    First, for the uniqueness results of problem (1), we need the following assumptions.

    A1ftu1ftu2for any 0<t<1,0u1u2.

    A2For any r01, there exists a constant q01such that



    We shall consider the Banach space E=C01equipped with the norm u=max0t1utand let




    In view of Lemma 2.3, we define an operator Tas


    where Gtsis given by (5).

    By A1it is easy to see that the operator T:DC+01is increasing. Observe that the BVP (1) has a solution if and only if the operator Thas a fixed point.

    Obviously, from A1we obtain


    In what follows, we first prove T:DD. In fact, for any uD, there exist a positive constants 0<mu<1<Musuch that


    Then, from A1, ftunon-decreasing respect to uand A2, we can imply that for s01,q01


    From (11) and Lemma 2.4, we obtain




    Eqs. (12) and (13) and assumption A3imply that T:DD.

    Now, we are in the position to give the first main result of this chapter.

    Theorem 1.1 Suppose A1A3hold. Then problem (1) has a unique, nondecreasing solution uD, moreover, constructing successively the sequence of functions


    for any initial function h0tD, then hntmust converge to utuniformly on 01and the rate of convergence is


    where 0<θ<1, which depends on the initial function h0t.

    Proof. For any h0D, we let




    Since the operator Tis increasing, A1,A2and (16)(20) imply that there exist iterative sequences un,vnsatisfying


    In fact, from (19) and (20), we have




    Then, by (22)(24) and induction, the iterative sequences un,vnsatisfy


    Note that u0t=mMv0t, from A1, (10), (19) and (20), it can obtained by induction that


    where θ=mM.

    From (21) and (25) we know that


    and since 1θqnMh0t0,asn,this yields that there exists uDsuch that


    Moreover, from (26) and


    we have






    From A1,(19) and (20), we have


    This together with (27) and uniqueness of limit imply that usatisfy u=Tu,that is uDis a solution of BVP (1) and (2).

    From (19)(21) and A1, we obtain


    It follows from (26)(29) that




    Hence, (15) holds. Since h0tis arbitrary in Dwe know that utis the unique solution of the boundary value problem (1) in D.

    We construct an example to illustrate the applicability of the result presented.

    Example 2.1. Consider the following boundary value problem


    whereα=52,β=1,a1=22,η1=12andftu=u2316costtis increasing function with respect toufor allt01, so, assumptionA1satisfied.

    By simple calculation we haved=12212=12.

    For anyr01, there existsq=1201such that


    thus,ftusatisfiesA2and is singular att=0.

    On the other hand,


    so, assumptionA3is satisfied.

    Hence, all the assumptions of Theorem1.1are satisfied. Which implies that the boundary value30has an unique, nondecreasing solutionuD.


    3. Existence result for inclusion problem (2)

    We provide another result about the existence of solutions for the problem (2) by using the assumption of nonconvex compact values for multifunction. Our strategy to deal with this problem is based on the Covitz-Nadler theorem for the contraction multivalued maps [44] for lower semi-continuous maps with decomposable values.

    First, we will present notations, definitions and preliminary facts from multivalued analysis which are used throughout this chapter. For more details on the multivalued maps, see the book of Aubin and Cellina [48], Demling [49], Gorniewicz [50] and Hu and Papageorgiou [51], see also [44, 48, 49, 52, 53, 54].

    Here C01Rdenotes the Banach space of all continuous functions from 01into Rwith the norm u=suput:forallt01,L101R, the Banach space of measurable functions u:01Rwhich are Lebesgue integrable, normed by uL1=01utdt.

    Let Xdbe a metric space induced from the normed space X. We denote

    P0X=APX:Aϕ,PbX=AP0X:Ais bounded,PclX=AP0X:Ais closed,PcpX=AP0X:Ais compact,Pb,clX=AP0X:Ais closedandbounded,

    where PXis the family of all subsets of X.

    Definition 3.1. A multivalued mapG:XPX.

    1. Guis convex (closed) valued if Guis convex (closed) for all uX,

    2. is bounded on bounded sets if GB=uBGuis bounded in Xfor all BPbXi.e., supuBsupvvGu<,

    3. has a fixed point if there is uXsuch that uGu. The fixed point set of the multivalued operator Gwill be denote by Fix G.

    Definition 3.2. A multivalued map G:01PclRis said to be measurable if for every yRthe function


    is measurable.

    Definition 3.3. Let Ybe a nonempty closed subset of a Banach space Eand G:YPclEbe a multivalued operator with nonempty closed values.

    1. Gis said to be lower semi-continuous (l.s.c) if the set xX:GxUϕis open for any open set Uin E.

    2. Ghas a fixed point if there is xYsuch that xGx.

    For each uC01R, define the set of selection of Fby


    For PX=2X, consider the Pompeiu-Hausdorff metric (see [55]).

    Hd:2X×2X0given by


    where daB=infbBdaband dbA=infaAdab. Then Pb,clXHdis a metric space and PclXHdis a generalized metric space see [8].

    Definition 3.4. Let Abe a subset of 01×R. Ais LBmeasurable if Abelongs to the σalgebra generated by all sets of the J×D, where Jis Lebesgue measurable in 01and Dis Borel measurable in R.

    Definition 3.5. A subset Aof L101Ris decomposable if all u,vAand measurable J01=j, the function uχJ+vχj\JA, where χJstands for the caracteristic function of J.

    Definition 3.6. Let Ybe a separable metric space and N:YPL101Rbe a multivalued operator. We say Nhas property (BC) if Nis lower semi-continuous (l.s.c) and has nonempty closed and decomposable values.

    Let F:01×RPRbe a multivalued map with nonempty compact values. Define a multivalued operator


    by letting


    Definition 3.7. The operator Φis called the Niemytzki operator associated with F. We say Fis of the lower semi-continuous type (l.s.c type) if its associated Niemytzki operator Φhas (BC) property.

    Definition 3.8. A multivalued operator N:XPclXis called.

    1. ρLipschitz if and only if there exists ρ>0such that HdNuNvρduvfor each u,vX,

    2. a contraction if and only if it is ρLipschitz with ρ<1.

    Lemma 3.1. ([44] Covitz-Nadler). LetXdbe a complete metric space. IfN:XPclXis a contraction, then FixNϕ, where FixNis the fixed point of the operatorN.

    Definition 3.9. A measurable multivalued function F:01PXis said to be integrably bounded if there exists a function gL101Xsuch that, for all vFt,vgtfor a.e. t01.

    Let us introduce the following hypotheses.

    A4F:01×RPcpRbe a multivalued map verifying.

    1. tuFtuis LBmeasurable.

    2. uFtuis lower semi-continuous for a.e.t01.

    A5Fis integrably bounded, that is, there exists a function mL101R+such that Ftu=supv:vFtumtfor almost all t01.

    Lemma 3.2. [56] LetF:01×RPcpRbe a multivalued map. AssumeA4andA5hold. ThenFis of thel.s.c.type.

    Definition 3.10. A function uAC201Ris called a solution to the boundary value problem (2) if usatisfies the differential inclusion in (2)a.e.on 01and the conditions in (2).

    Finally, we state and prove the second main result of this Chapter. We prove the existence of solutions for the inclusion problem (2) with a nonconvex valued right hand side by applying a fixed point theorem for multivalued maps due to Covitz and Nadler. For investigation of the problem (2) we shall provide an application of the Lemma 3.4and the following Lemma.

    Lemma 3.3. ([13]) A multifunctionF:XCXis called a contraction whenever there existsγ01such thatHdNuNvγduvfor allu,vX.

    Now, we present second main result of this section.

    Theorem 1.2 Assume that the following hypothyses hold.

    H1F:J×RPcpRis an integrable bounded multifunction such that the map tFtuis measurable,

    H2HdFtu1Ftu2mtu1u2for almost all tJand u1,u2Rwith mL1JRand d0Ft0mtfor almost all tJ. Then the problem (2) has a solution provided that


    Proof. We transform problem (2) into a fixed point problem. Consider the operator N:C01PC01Rdefined by


    where Gtsdefined by (5). It is clear that fixed points of Nare solution of (2).

    We shall prove that Nfulfills the assumptions of Covitz-Nadler contraction principle.

    Note that, the multivalued map tFtutis measurable and closed for all uAC10(e.g., [52] Theorem III.6). Hence, it has a measurable selection and so the set SF,uis nonempty, so, Nuis nonempty for any uC0.

    First, we show that Nuis a closed subset of Xfor all uAC10R. Let uXand unn1be a sequence in Nuwith unu,asnin uC0. For each n, choose ynSF,usuch that


    Since Fhas compact values, we may pass onto a subsequence (if necessary) to obtain that ynconverges to yL101Rin L101R. In particular, ySF,uand for any t01, we have


    i.e., uNuand Nuis closed.

    Next, we show that Nis a contractive multifunction with constant l<1. Let u,vC01Rand h1Nu. Then there exist y1SF,usuch that


    By H2, we have


    for almost all tJ.

    So, there exists wSF,vsuch that


    for almost all tJ.

    Define the multifunction U:JPRby


    It is easy to chek that the multifunction V=UFvis measurable (e.g., [52] Theorem III.4).

    Thus, there exists a function y2twhich is measurable selection for V. So, y2SF,vand for each tJ, we have


    Now, consider h2Nuwhich is defined by


    and one can obtain




    Analogously, interchanging the roles of uand v, we obtain


    Since Nis a contraction, it follows by Lemma 3.1(by using the result of Covitz and Nadler) that Nhas a fixed point which is a solution to problem (2).

    We construct an example to illustrate the applicability of the result presented.

    Example 3.1. Consider the problem


    subject to the three-point boundary conditions


    whereα=52,β=1a1=12,a2=32,η1=116,η2=516. andFtut:01×R2Rmultivalued map given by





    we have


    which shows thatH2holds

    So, ifmt=t+12for allt01, then


    It can be easily found thatd=112116523251652=0,9176244637.



    Hence, all assumptions and conditions of Theorem1.2are satisfied. So, Theorem1.2implies that the inclusion problem(32)and(33)has at least one solution.


    4. Conclusions

    This chapter concerns the boundary value problem of a class of fractional differential equations involving the Riemann-Liouville fractional derivative with nonlocal boundary conditions. By using the properties of the Green’s function and the monotone iteration technique, one shows the existence of positive solutions and constructs two successively iterative sequences to approximate the solutions. In the multi-valued case, an existence result is proved by using fixed point theorem for contraction multivalued maps due to Covitz and Nadler. The results of the present chapter are significantly contribute to the existing literature on the topic.



    The authors want to thank the anonymous referee for the thorough reading of the manuscript and several suggestions that help us improve the presentation of the chapter.


    Conflict of interest

    The authors declare no conflict of interest.


    chapter PDF

    © 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

    How to cite and reference

    Link to this chapter Copy to clipboard

    Cite this chapter Copy to clipboard

    Noureddine Bouteraa and Habib Djourdem (November 3rd 2021). A Study of Nonlinear Boundary Value Problem [Online First], IntechOpen, DOI: 10.5772/intechopen.100491. Available from:

    chapter statistics

    28total chapter downloads

    More statistics for editors and authors

    Login to your personal dashboard for more detailed statistics on your publications.

    Access personal reporting

    We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

    More About Us