Criteria for existence of solutions for a Liouville–Caputo boundary value problem via generalized Gronwall’s inequality

Abstract

In this research, we first investigate the existence of solutions for a new fractional boundary value problem in the Liouville–Caputo setting with mixed integro-derivative boundary conditions. To do this, Kuratowski’s measure of noncompactness and Sadovskii’s fixed point theorem are our tools to reach this aim. In the sequel, we discuss the continuous dependence of solutions on parameters by means of the generalized Gronwall inequality. Moreover, we consider an inclusion version of the given boundary problem in which we study its existence results by means of the endpoint theory. Finally, we prepare two simulative numerical examples to confirm the validity of the analytical findings.

1 Introduction

Fractional models have recently attracted the attention of many scientists and engineers of the modeling fields. This amount of attention can be a motivation for researchers and mathematicians to study the various new fractional models for the complex natural phenomena in the world around. By taking into account the importance and effectiveness of the mathematical modelings, researchers design various novel categories of fractional structures of initial and boundary value problems and study the existence, uniqueness, and stability of solutions by using some analytical and numerical techniques. Furthermore, many researchers endeavor to approximate the obtained solutions by utilizing some numerical and iterative algorithms. The most important to researchers today is the understanding of some qualitative properties of the solutions for such fractional dynamical systems. In this direction, numerous research articles on the existence and stability of solutions have been published, including [1–21].

As we know, the inequalities play an important role in all branches of the pure and applied mathematics, and we can use some of their properties to infer different desired identities and estimates. By following this path, one can observe different published research works in the recent years in which the researchers discuss the existence of solutions for the given boundary problems by combining some analytical methods and some well-known inequalities. For instance, we point out some papers, including [22–25]. Applications of the generalized Gronwall inequality can be found in some papers such as [26–29].

In 2016, Ahmad et al. [30] studied the existence of solutions for a mixed fractional initial value problem in the inclusion version involving Hadamard fractional derivative and Riemann–Liouville fractional integrals given by

$$ \textstyle\begin{cases} {}^{\mathcal{H}}\mathfrak{D}^{ \alpha } ( v(z) - \sum_{i=1}^{m} {}^{ \mathcal{RL}}\mathcal{I}^{\beta _{i}} h_{i}(z, v(z)) ) \in F (z, v(z) ) \quad (z\in J:= [1,T]), \\ v(1) = 0, \end{cases} $$

where ${}^{\mathcal{H}}\mathfrak{D}^{\alpha }$ is the Hadamard fractional derivative of order $\alpha \in (0,1]$ and ${}^{\mathcal{RL}}\mathcal{I}^{\phi }$ denotes the Riemann–Liouville fractional integral of order $0 < \phi \in \{ \beta _{1}, \beta _{2}, \dots, \beta _{m} \}$, $h_{i} \in \mathcal{C}(J\times \mathbb{R}, \mathbb{R})$ with $h_{i}(1,0)= 0\ (i=1,2,\dots, m)$, and finally $F: J \times \mathbb{R} \to \mathcal{P}(\mathbb{R})$ is a multifunction. The authors used the notions defined in the endpoint theory and proved the main existence results for the above initial inclusion problem with the aid of an approximate endpoint property [30].

Some existing structures in the works mentioned above motivate us to design a new fractional structure of the boundary value problem in the Liouville–Caputo framework (Lio–CapBvp) furnished with mixed integro-derivative conditions as

$$ \textstyle\begin{cases} {}^{\mathcal{C}}\mathfrak{D}_{0}^{ \xi } v(z) = \Phi _{*} (z, v(z) ) \\ v(0) = s_{1}^{*} {}^{\mathcal{C}}\mathfrak{D}_{0}^{1} v(0) + s_{2}^{*} {}^{\mathcal{C}}\mathfrak{D}_{0}^{1} v(1),\qquad {}^{\mathcal{C}} \mathfrak{D}_{0}^{\xi -1} v(1) = \int _{0}^{1} {}^{ \mathcal{C}}\mathfrak{D}_{0}^{\xi -1} v(\varsigma ) \,\mathrm{d} \varsigma \end{cases} $$

(1)

so that $z\in [0,1]$, $s_{1}^{*}, s_{2}^{*} \in \mathbb{R}^{+}$, ${}^{\mathcal{C}}\mathfrak{D}_{0}^{ \xi }$ stands for the Liouville–Caputo derivative of order $\xi \in (1,2)$, ${}^{\mathcal{C}}\mathfrak{D}_{0}^{1} = \frac{\mathrm{d}}{\mathrm{d}z}$ and $\Phi _{*} \in \mathcal{C}( [0,1] \times \mathbb{R}, \mathbb{R})$. To derive the existence criteria of solutions for the proposed fractional Lio–CapBvp (1), we invoke the notion of the well-known measure of noncompactness attributed to Kuratowski. In this direction, we construct a condensing operator and establish our existence result by utilizing a well-known fixed point result due to Sadovskii.

In the sequel, we investigate some estimates by utilizing the generalized Gronwall inequality. Moreover, we discuss an inclusion version of the given boundary problem in which the proof process is based on the approximate endpoint property and some properties of inequalities in relation to the Pompeiu–Hausdorff metric defined for multifunctions. It is notable that our proposed fractional Lio–CapBvp (1) involves new mixed integro-derivative boundary conditions simultaneously, and this formulated framework is novel at this moment.

An outline of the current research manuscript is prepared as follows: In Sect. 2, we assemble some auxiliary and key notions and theorems which are essential for the rest of the paper. In Sect. 3, an analytical criterion about the existence of solutions for the fractional Lio–CapBvp (1) is derived based on the fixed point result attributed to Sadovskii. In the following, we check the dependence of solutions by utilizing the generalized Gronwall inequality. In Sect. 4, we regard an inclusion version of the given structure (1) and find endpoints of this system by applying some inequalities in relation to multifunctions. In Sect. 5, we prepare two simulative numerical examples to confirm the validity of the obtained analytical findings.

2 Preliminaries

In this part of the article, we intend to recall some auxiliary and key concepts and theorems which are necessary in the rest of this manuscript. Take $\xi > 0$. The Riemann–Liouville integral (RL-integral) for a real-valued function v on $[a,b]$ in the fractional setting is given by

$$ {}^{\mathcal{RL}}\mathcal{I}_{0}^{ \xi } v(z)= \int _{0}^{z} \frac{(z - \varsigma )^{ \xi -1 }}{\Gamma ( \xi )} v( \varsigma ) \, \mathrm{d} \varsigma, $$

provided that the integral above is finite-valued (see [31, 32]). At this stage, let $n-1< \xi < n$ or $n = [ \xi ] +1$. The Liouville–Caputo derivative (Lio–Cap derivative) of a given function $v \in \mathcal{C}^{(n)}_{\mathbb{R}}([a,b])$ in the fractional setting is defined by

$$ {}^{\mathcal{C}}\mathfrak{D}_{0}^{ \xi } v(z)= \int _{0}^{z} \frac{(z - \varsigma )^{n- \xi -1}}{\Gamma (n- \xi )} v^{(n)}( \varsigma ) \,\mathrm{d} \varsigma $$

so that the right-hand side integral is finite-valued (see [31, 32]). In [33], one can see that the structure of the general series solution for the homogeneous differential equation ${}^{\mathcal{C}}\mathfrak{D}_{0}^{ \xi } v(z)=0$ is $v(z ) = \alpha _{*}^{(0)} + \alpha _{*}^{(1)} z + \alpha _{*}^{(2)} z^{2}+ \cdots + \alpha _{*}^{(n-1)} z^{n-1}$, and the following holds:

$$ {}^{\mathcal{RL}}\mathcal{I}_{0}^{ \xi } \bigl( {}^{\mathcal{C}} \mathfrak{D}_{0}^{ \xi } v(z) \bigr)= v(z) + \sum_{k=0}^{n-1} \alpha _{*}^{(k)} z^{k} = v(z)+ \alpha _{*}^{(0)} + \alpha _{*}^{(1)}z + \alpha _{*}^{(2)}z^{2}+ \cdots + \alpha _{*}^{(n-1)} z^{n-1}, $$

where $\alpha _{*}^{(0)},\ldots, \alpha _{*}^{(n-1)} \in \mathbb{R}$ and $n=[ \xi ]+1$. In the sequel, some important inequalities and analytical notions are reviewed. The measure of noncompactness Ω attributed to Kuratowski is introduced by

$$ \Omega (\mathfrak{B}):= \inf \Biggl\{ \epsilon > 0: \mathfrak{B} = \bigcup _{k=1}^{n} \mathfrak{B}_{k} \text{ and } \operatorname{diam} ( \mathfrak{B}_{k}) \leq \epsilon \text{ for } k=1,\dots,n \Biggr\} , $$

so that $\operatorname{diam}(\mathfrak{B}_{k})= \sup \{ | v - v' |: v,v'\in \mathfrak{B}_{k} \}$ and $\mathfrak{B}$ is assumed to be a bounded subset of the Banach space $\mathbb{V}$. Additionally, it is known that $0 \leq \Omega (\mathfrak{B}) \leq \operatorname{diam}(\mathfrak{B})< + \infty$ (see [34]).

Lemma 1

([34])

Let $\mathbb{V}$ be an arbitrary real Banach space and $\mathfrak{B}, \mathfrak{B}_{1}, \mathfrak{B}_{2}\subseteq \mathbb{V}$ be bounded subsets of $\mathbb{V}$. Then the following are valid:

• (Con1) $\mathfrak{B}$ is precompact iff $\Omega (\mathfrak{B})=0$;

• (Con2) $\Omega (\mathfrak{B}) = \Omega (\overline{\mathfrak{B}}) = \Omega ( \operatorname{cnvx} (\mathfrak{B}))$, where $\operatorname{cnvx} (\mathfrak{B})$ and $\overline{\mathfrak{B}}$ stand for the convex hull and the closure of $\mathfrak{B}$, respectively;

Moreover, the following inequalities hold:

• (Con3) if $\mathfrak{B}_{1}\subseteq \mathfrak{B}_{2}$, then $\Omega (\mathfrak{B}_{1})\leq \Omega (\mathfrak{B}_{2})$;

• (Con4) $\Omega (\mu + \mathfrak{B})\leq \Omega (\mathfrak{B})$ for each $\mu \in \mathbb{R}$;

• (Con5) $\Omega ( \mu \mathfrak{B})=| \mu | \Omega (\mathfrak{B})$ for all $\mu \in \mathbb{R}$;

• (Con6) $\Omega (\mathfrak{B}_{1} + \mathfrak{B}_{2})\leq \Omega ( \mathfrak{B}_{1} )+\Omega (\mathfrak{B}_{2})$ so that $\mathfrak{B}_{1}+\mathfrak{B}_{2} =\{ v_{1} + v_{2}; v_{1}\in \mathfrak{B}_{1}, v_{2} \in \mathfrak{B}_{2}\}$;

• (Con7) $\Omega (\mathfrak{B}_{1}\cup \mathfrak{B}_{2})\leq \max \{ \Omega ( \mathfrak{B}_{1}), \Omega (\mathfrak{B}_{2}) \}$.

Lemma 2

([35])

Let $\mathbb{V}$ be a Banach space. For each bounded subset $\mathfrak{B}$ of $\mathbb{V}$, a countable set $\mathfrak{B}_{0}$ of $\mathfrak{B}$ exists so that $\Omega (\mathfrak{B})\leq 2 \Omega (\mathfrak{B}_{0})$.

Lemma 3

([34])

Let $\mathbb{V}$ be a Banach space. If $\mathfrak{B}\subseteq \mathcal{C}_{ \mathbb{V} }([a,b])$ is assumed to be a bounded and equicontinuous set, then $\Omega (\mathfrak{B}(z))$ is continuous on $[a,b]$ and $\Omega (\mathfrak{B})= \sup_{z\in [a,b]} \Omega (\mathfrak{B}(z))$.

Lemma 4

([34])

Let $\mathbb{V}$ be a Banach space. If $\mathfrak{B}= \{ v_{n} \}_{n \geq 1} \subseteq \mathcal{C}_{ \mathbb{V} }([a,b])$ is assumed to be a bounded and countable set, then $\Omega (\mathfrak{B}(z))$ is Lebesgue integrable on $[a,b]$ and the following inequality is valid:

$$ \Omega \biggl( \biggl\{ \int _{0}^{z} v_{n}( \varsigma ) \, \mathrm{d} \varsigma \biggr\} _{n\geq 1} \biggr) \leq 2 \int _{0}^{z} \Omega \bigl( \bigl\{ v_{n}( \varsigma ) \bigr\} _{n \geq 1} \bigr) \,\mathrm{d} \varsigma. $$

Definition 5

([34])

Let $\mathbb{V}$ be a Banach space and $\Phi _{*}: \mathbb{D}\subset \mathbb{V} \to \mathbb{V}$ be continuous and bounded. Then the operator $\Phi _{*}$ is called condensing if $\Omega (\Phi _{*} (\mathfrak{B})) <\Omega (\mathfrak{B})$ for any bounded closed set $\mathfrak{B} \subseteq \mathbb{D}$.

Theorem 6

([34], Sadovskii’s fixed point theorem)

Let $\mathbb{V}$ be a Banach space. Also assume that $\mathfrak{B}\subseteq \mathbb{V}$ is a bounded closed convex set and the continuous map $\Phi _{*}: \mathfrak{B}\rightarrow \mathfrak{B}$ is condensing. Then the map $\Phi _{*}$ has at least one fixed point in $\mathfrak{B}$.

Theorem 7

([36], The generalized Gronwall inequality)

Let r̆ be a nonnegative locally integrable function on the interval $[0, T]$ for some $T \leq \infty$ and h a nonnegative nondecreasing continuous mapping on $[0, T]$ satisfying $h(z)\leq M$ where M is a nonzero constant. In addition, let v̆ be a nonnegative locally integrable function on $[0, T]$ such that

$$ \breve{v}(z)\leq \breve{r}(z)+h(z) \int _{0}^{z}(z-\varsigma )^{q-1} \breve{v}(\varsigma ) \,\mathrm{d}\varsigma, $$

where $q > 0$. Then we have

$$ \breve{v}(z)\leq \breve{r}(z)+ \int _{0}^{z}\sum_{k=1}^{\infty } \biggl[ \frac{(h(z)\Gamma (q))^{k}}{\Gamma (kq)}(z- \varsigma )^{kq-1} \breve{r}(\varsigma ) \biggr] \,\mathrm{d}\varsigma \quad\bigl(z \in [0,T] \bigr). $$

Notation 8

In this manuscript, we denote the normed space by $(\mathbb{V},\Vert \cdot \Vert _{\mathbb{V} })$ and we introduce the notations ${\mathcal{P}}( \mathbb{V} )$, ${\mathcal{P}}_{cls}( \mathbb{V} )$, ${\mathcal{P}}_{bnd}( \mathbb{V} )$, ${\mathcal{P}}_{cmp}( \mathbb{V} )$, and ${\mathcal{P}}_{cvx}( \mathbb{V} )$ for the collections of all nonempty, closed, bounded, compact, and convex subsets in the space $\mathbb{V}$, respectively.

An element $v \in \mathbb{V}$ is said to be a fixed point for an abstract multivalued operator $\mathfrak{Q}_{*}: \mathbb{V} \to \mathcal{P}(\mathbb{V})$ if the inclusion $v\in \mathfrak{Q}_{*}(v)$ is satisfied [37]. At this phase, we denote the family of all existing fixed points of $\mathfrak{Q}_{*}$ by ${\mathfrak{FIX}} (\mathfrak{Q}_{*})$ [37]. Subsequently, the Pompeiu–Hausdorff metric $\mathbb{H}_{d_{\mathbb{V}}}: {\mathcal{P}}(\mathbb{V}) \times { \mathcal{P}}(\mathbb{V}) \to \mathbb{R} \cup \{\infty \}$ is defined by

$$ \mathbb{H}_{d_{\mathbb{V}}}(\mathcal{M}_{1}, \mathcal{M}_{2}) = \max \Bigl\{ \sup_{m_{1} \in \mathcal{M}_{1}}d_{\mathbb{V}}(m_{1}, \mathcal{M}_{2}), \sup_{m_{2} \in \mathcal{M}_{2}}d_{\mathbb{V}}( \mathcal{M}_{1},m_{2}) \Bigr\} $$

so that $d_{\mathbb{V}}(\mathcal{M}_{1},m_{2}) = \inf_{m_{1}\in \mathcal{M}_{1}}d_{\mathbb{V}}(m_{1},m_{2})$ and $d_{\mathbb{V}}(m_{1},\mathcal{M}_{2}) = \inf_{m_{2}\in \mathcal{M}_{2}}d_{\mathbb{V}}(m_{1},m_{2})$ (see [37]). A multifunction $\mathfrak{Q}_{*}: \mathbb{V} \to {\mathcal{P}}_{cls}(\mathbb{V})$ is called Lipschitz with real constant $\hat{\iota } >0$ if the inequality

$$ \mathbb{H}_{d_{\mathbb{V}}} \bigl(\mathfrak{Q}_{*}(v_{1}), \mathfrak{Q}_{*} (v_{2}) \bigr) \leq \hat{\iota } d_{\mathbb{V}}(v_{1}, v_{2}) $$

is valid for each $v_{1}, v_{2} \in \mathbb{V}$. It is notable that a Lipschitz map $\mathfrak{Q}_{*}$ is called a contraction if $\hat{\iota } \in (0,1)$ (see [37]). Also, we say that $\mathfrak{Q}_{*}$ possesses the complete continuity property if $\mathfrak{Q}_{*}(K)$ has the relative compactness property for any $K\in \mathcal{P}_{bnd}(\mathbb{V})$, and also $\mathfrak{Q}_{*}: [0,1] \to {\mathcal{P}}_{cls}(\mathbb{R})$ is measurable if $z \longmapsto d_{\mathbb{V}}( r, \mathfrak{Q}_{*}(z ))$ is measurable for any $r \in \mathbb{R}$ (see [37, 38]). Furthermore, $\mathfrak{Q}_{*}$ is called upper semicontinuous (u.s.c.) if for every $v \in \mathbb{V}$, the set $\mathfrak{Q}_{*}(v)$ belongs to $\mathcal{P}_{cls}(\mathbb{V})$ and for any open set $\mathbb{G}$ which contains $\mathfrak{Q}_{*}(v)$, a neighborhood $\mathcal{U}_{0}^{*}$ of v exists so that $\mathfrak{Q}_{*}(\mathcal{U}_{0}^{*})$ is contained in $\mathbb{G}$ [37].

It is an evident fact that $\mathfrak{Q}_{*}$ possesses convex values if $\mathfrak{Q}_{*}(v) \in \mathcal{P}_{cvx}(\mathbb{V})$ for any $v \in \mathbb{V}$. We denote the family of all existing selections of $\mathfrak{Q}_{*}$ at point $v \in \mathcal{C}_{\mathbb{R}}([0,1])$ by the following rule:

$$ ( \mathfrak{SEL})_{\mathfrak{Q}_{*},v}:= \bigl\{ \hbar \in \mathcal{L}_{ \mathbb{R}}^{1} \bigl([0,1] \bigr): \hbar (z)\in \mathfrak{Q}_{*} \bigl(z, v(z) \bigr) \bigr\} $$

for any $z \in [0,1]$ (a.e.) (see [37, 38]). It is important to pay attention to the issue that by assuming $\mathfrak{Q}_{*}$ as an arbitrary abstract multivalued operator, we deduce that for each $v \in \mathcal{C}_{\mathbb{V}}([0,1])$, the set $( \mathfrak{SEL})_{\mathfrak{Q}_{*},v}$ is nonempty if $\dim (\mathbb{V})$ is finite (see [37]).

Definition 9

([39])

An element $v \in \mathbb{V}$ is called an endpoint for the given multivalued operator $\mathfrak{Q}_{*}: \mathbb{V}\to \mathcal{P}(\mathbb{V})$ if we have $\mathfrak{Q}_{*}(v )=\{ v \}$.

The multivalued mapping $\mathfrak{Q}_{*}: \mathbb{V}\to \mathcal{P}(\mathbb{V})$ has an approximate endpoint property (approx-endpoint property) if we have $\inf_{v_{1} \in \mathbb{V}}\sup_{ v_{2} \in \mathfrak{Q}_{*}(v_{1} )}d_{\mathbb{V}}( v_{1},v_{2} )=0$ [39]. The next result is necessary for us in the rest of the article.

Theorem 10

([39], Endpoint theorem)

Assume that the metric space $(\mathbb{V},d_{\mathbb{V}})$ is complete and $\psi: [0, \infty ) \to [0, \infty )$ has an upper semicontinuity property provided that $\psi (z)< z$ and $\liminf_{z\to \infty }(z-\psi (z))>0$ for any $z>0$. Additionally, let us assume that $\mathfrak{Q}_{*}: \mathbb{V}\to \mathcal{P}_{cls,bnd}(\mathbb{V})$ is such that the inequality

$$ \mathbb{H}_{d_{\mathbb{V}}}(\mathfrak{Q}_{*} v_{1}, \mathfrak{Q}_{*}v_{2} )\leq \psi \bigl(d_{\mathbb{V}}(v_{1}, v_{2} ) \bigr) $$

is satisfied for any $v_{1}, v_{2} \in \mathbb{V}$. Then a unique endpoint exists for $\mathfrak{Q}_{*}$ iff $\mathfrak{Q}_{*}$ possesses approx-endpoint property.

3 Criterion of existence for the fractional Lio–CapBvp (1)

In the present section, we derive a criterion of existence of solutions for the proposed Lio–CapBvp (1). In order to achieve this goal, we construct the space $\mathbb{V} = \{ v(z): v(z) \in \mathcal{C}_{ \mathbb{R} }([0,1]) \}$ furnished with the supremum norm $\Vert v \Vert _{\mathbb{V} } = \sup_{ z\in [0,1] } \vert v(z) \vert$. Then we can simply check that an ordered pair $( \mathbb{V}, \Vert \cdot \Vert _{\mathbb{V}} )$ is a Banach space. By following the contents of two previous sections, now we establish a new existence criterion and also discuss to the stability of the fractional Lio–CapBvp (1). To achieve this goal, we resort to the next structural lemma in relation to (1) at this moment.

Lemma 11

Let $\xi \in (1,2)$, $\xi -1 \in (0,1)$, $s_{1}^{*}, s_{2}^{*} \in \mathbb{R}^{+}$, and $g \in \mathcal{C}( [0,1], \mathbb{R})$. Then the solution of the given linear fractional Lio–CapBvp

$$ \textstyle\begin{cases} {}^{\mathcal{C}}\mathfrak{D}_{0}^{ \xi } v(z) = {g} (z),\quad (z \in [0,1]), \\ v(0) = s_{1}^{*} {}^{\mathcal{C}}\mathfrak{D}_{0}^{1} v(0) + s_{2}^{*} {}^{\mathcal{C}}\mathfrak{D}_{0}^{1} v(1), \qquad {}^{\mathcal{C}} \mathfrak{D}_{0}^{\xi -1} v(1) = \int _{0}^{1} {}^{ \mathcal{C}}\mathfrak{D}_{0}^{\xi -1} v(\varsigma ) \,\mathrm{d} \varsigma \end{cases} $$

(2)

is given by

$$\begin{aligned} v(z) = {}&\frac{1}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} {g}( \varsigma ) \,\mathrm{d} \varsigma + \frac{ s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1-\varsigma )^{\xi - 2} {g}( \varsigma ) \,\mathrm{d} \varsigma \\ &{}- \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} {g}(\varsigma ) \,\mathrm{d}\varsigma \\ &{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }{g}(\tau ) \,\mathrm{d}\tau \, \mathrm{d}\varsigma. \end{aligned}$$

(3)

Proof

At first, we assume that $v_{0}^{*}$ satisfies the fractional Lio–CapBvp (2). Then, in view of some properties of the Caputo and Riemann–Liouville fractional operators, two arbitrary real constants $\alpha _{*}^{(0)}$ and $\alpha _{*}^{(1)}$ exist so that the equation $v_{0}^{*}(z) = {}^{\mathcal{RL}}\mathcal{I}_{0}^{\xi }{g}(z) + \alpha _{*}^{(0)} + \alpha _{*}^{(1)}z$ is valid for any $z\in [0,1]$. Indeed, we get

$$ v_{0}^{*}(z) = \int _{0}^{z} \frac{(z-\varsigma )^{ \xi -1}}{ \Gamma ( \xi )} {g}( \varsigma ) \,\mathrm{d}\varsigma + \alpha _{*}^{(0)} + \alpha _{*}^{(1)}z. $$

(4)

In this case, we obtain ${}^{\mathcal{C}}\mathfrak{D}_{0}^{1} v_{0}^{*} (z) = \int _{0}^{z} \frac{(z-\varsigma )^{ \xi -2}}{ \Gamma ( \xi -1)} { g}(\varsigma ) \,\mathrm{d}\varsigma + \alpha _{*}^{(1)}$ and

$$\begin{aligned} & {}^{\mathcal{C}}\mathfrak{D}_{0}^{\xi -1} v_{0}^{*}(z) = \int _{0}^{z} {g}(\varsigma ) \,\mathrm{d} \varsigma + \alpha _{*}^{(1)} \frac{z^{2-\xi }}{\Gamma (3-\xi )}, \\ & \int _{0}^{1} {}^{\mathcal{C}} \mathfrak{D}_{0}^{\xi -1} v_{0}^{*}( \varsigma ) \,\mathrm{d}\varsigma = \int _{0}^{1} \int _{0}^{\varsigma }{g}(\tau ) \,\mathrm{d}\tau \, \mathrm{d}\varsigma + \alpha _{*}^{(1)} \frac{1}{\Gamma (4-\xi )}. \end{aligned}$$

Now, by utilizing the given integro-derivative boundary conditions in (2), the real constants $\alpha _{*}^{(0)}$ and $\alpha _{*}^{(1)}$ are obtained as follows:

$$\begin{aligned} \alpha _{*}^{(0)} = {}&\frac{s_{2}^{*} }{\Gamma (\xi -1)} \int _{0}^{1} (1-\varsigma )^{\xi - 2} {g}( \varsigma ) \,\mathrm{d}\varsigma - \frac{( s_{1}^{*} + s_{2}^{*} ) \Gamma ( 4-\xi ) }{ 2-\xi } \int _{0}^{1} {g}(\varsigma ) \,\mathrm{d}\varsigma \\ &{}+ \frac{( s_{1}^{*} + s_{2}^{*} ) \Gamma ( 4-\xi ) }{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }{g}(\tau ) \,\mathrm{d}\tau \, \mathrm{d}\varsigma \end{aligned}$$

and

$$ \alpha _{*}^{(1)} = - \frac{ \Gamma ( 4-\xi ) }{ 2-\xi } \int _{0}^{1} {g}( \varsigma ) \,\mathrm{d} \varsigma + \frac{ \Gamma ( 4-\xi ) }{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }{ g}(\tau ) \,\mathrm{d}\tau \, \mathrm{d}\varsigma. $$

Here, by inserting the obtained values of $\alpha _{*}^{(0)}$ and $\alpha _{*}^{(1)}$ into (4), one can observe that $v_{0}^{*}$ satisfies the integral equation (3). In other words,

$$\begin{aligned} v_{0}^{*}(z) ={}& \frac{1}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} {g}( \varsigma ) \,\mathrm{d} \varsigma + \frac{ s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1-\varsigma )^{\xi - 2} {g}( \varsigma ) \,\mathrm{d} \varsigma \\ &{}- \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} {g}(\varsigma ) \,\mathrm{d}\varsigma \\ &{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }{g}(\tau ) \,\mathrm{d}\tau \, \mathrm{d}\varsigma, \end{aligned}$$

and the proof is completed. □

In what follows, we invoke the notion of Kuratowski’s measure of noncompactness and establish a new criterion of solution existence for the relevant fractional Lio–CapBvp (1) based on some estimates and known existing properties of inequalities.

Theorem 12

Let $\Phi _{*}$ be a continuous real-valued function defined on $[0,1]\times \mathbb{V}$. Additionally, suppose a continuous function $\varphi: [0,1] \to \mathbb{R}^{+}$ exists such that the inequality

$$ \bigl\vert \Phi _{*} \bigl(z, v(z) \bigr) \bigr\vert \leq \varphi (z) $$

(5)

holds for any $z \in [0,1]$ and $v\in \mathbb{V}$. Furthermore, assume that a function $n_{ \Phi _{*}}: [0,1] \to \mathbb{R}^{+}$ exists so that the inequality

$$ \Omega \bigl( \Phi _{*} (z,\mathfrak{B} ) \bigr)\leq n_{ \Phi _{*}}(z)\Omega (\mathfrak{B}) $$

(6)

is valid for any bounded set $\mathfrak{B}\subset \mathbb{V}$ a.e. $z\in [0,1]$. Then, the given fractional Lio–CapBvp (1) possesses at least one solution on $[0,1]$ if

$$ \frac{ n_{ \Phi _{*}}^{*} }{\Gamma ( \xi + 1 )} + \frac{ s_{2}^{*} n_{ \Phi _{*}}^{*} }{ \Gamma (\xi ) } + \frac{ 3n_{ \Phi _{*}}^{*} ( s_{1}^{*} + s_{2}^{*} + 1)\Gamma (4-\xi )}{ 2(2-\xi ) } < \frac{1}{4}, $$

(7)

where $n_{ \Phi _{*}}^{*} = \sup_{z\in [0,1]} \vert n_{ \Phi _{*}}(z) \vert$.

Proof

To begin the current proof, in relation to Lio–CapBvp (1) and using Lemma 11, we define the mapping $\mathfrak{P}: \overline{ \mathfrak{E}_{\varepsilon }} \rightarrow \overline{ \mathfrak{E}_{\varepsilon }}$ as follows:

$$\begin{aligned} \mathfrak{P}(v) (z)={}& \frac{1}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} \Phi _{*} \bigl(\varsigma, v(\varsigma ) \bigr) \,\mathrm{d}\varsigma + \frac{ s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1-\varsigma )^{ \xi - 2} \Phi _{*} \bigl(\varsigma, v(\varsigma ) \bigr) \,\mathrm{d} \varsigma \\ &{}- \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \Phi _{*} \bigl(\varsigma, v(\varsigma ) \bigr) \,\mathrm{d}\varsigma \\ &{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }\Phi _{*} \bigl(\tau, v( \tau ) \bigr) \,\mathrm{d} \tau \,\mathrm{d}\varsigma, \end{aligned}$$

(8)

where $\overline{ \mathfrak{E}_{\varepsilon }}:= \{ v \in \mathbb{V}: \Vert v\Vert _{\mathbb{V} } \leq \varepsilon, \varepsilon \in \mathbb{R}^{+} \}$ is considered to be a convex, closed, and bounded subset of the Banach space $\mathbb{V}$. In this case, the suggested problem (1) is equivalent to the fixed point problem $v = \mathfrak{P}v$, and we have to prove that the operator $\mathfrak{P}$ has a fixed point, since the existence of a fixed point for the newly-introduced operator $\mathfrak{P}$ will imply the existence of a solution for the given fractional Lio–CapBvp (1).

To check all the hypotheses of Theorem 6, we first confirm the continuity of $\mathfrak{P}$ on $\overline{ \mathfrak{E}_{\varepsilon }}$. Let $\{ v_{n} \}_{n \geq 1}$ be a sequence contained in $\overline{ \mathfrak{E}_{\varepsilon }}$ so that $v_{n} \to v$ for each $v \in \overline{ \mathfrak{E}_{\varepsilon }}$. Due to the continuity of the function $\Phi _{*}$ on $[0,1]\times \mathbb{V}$, we obtain that $\lim_{n\to \infty }\Phi _{*} (z,v_{n}(z))=\Phi _{*} (z ,v(z ))$. In this case, by the help of the dominated convergence theorem of Lebesgue, we may write

$$\begin{aligned} &\lim_{n\to \infty }(\mathfrak{P} v_{n}) (z) \\ &\quad= \frac{1}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} \lim _{n \to \infty } \Phi _{*} \bigl(\varsigma, v_{n}( \varsigma ) \bigr) \,\mathrm{d}\varsigma \\ &\qquad{}+ \frac{ s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1- \varsigma )^{\xi - 2} \lim _{n\to \infty } \Phi _{*} \bigl(\varsigma, v_{n}( \varsigma ) \bigr) \,\mathrm{d}\varsigma \\ &\qquad{}- \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \lim_{n\to \infty } \Phi _{*} \bigl(\varsigma, v_{n}(\varsigma ) \bigr) \,\mathrm{d} \varsigma \\ &\qquad{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }\lim_{n\to \infty } \Phi _{*} \bigl(\tau, v_{n}( \tau ) \bigr) \,\mathrm{d}\tau \, \mathrm{d}\varsigma \\ &\quad= \frac{1}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} \Phi _{*} \bigl(\varsigma, v(\varsigma ) \bigr) \,\mathrm{d} \varsigma + \frac{ s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1-\varsigma )^{\xi - 2} \Phi _{*} \bigl(\varsigma, v(\varsigma ) \bigr) \,\mathrm{d}\varsigma \\ &\qquad{}- \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \Phi _{*} \bigl(\varsigma, v(\varsigma ) \bigr) \,\mathrm{d}\varsigma \\ &\qquad{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }\Phi _{*} \bigl(\tau, v( \tau ) \bigr) \,\mathrm{d} \tau \,\mathrm{d}\varsigma \\ &\quad = (\mathfrak{P} v) (z) \end{aligned}$$

for any $z \in [0,1]$. Hence, we realize that $\lim_{n\to \infty } (\mathfrak{P} v_{n})(z) = (\mathfrak{P} v)(z)$, which confirms the continuity of $\mathfrak{P}$ on $\overline{ \mathfrak{E}_{\varepsilon }}$. In the sequel, we want to check the uniform boundedness of $\mathfrak{P}$ on $\overline{ \mathfrak{E}_{\varepsilon }}$. To reach this aim, consider an element v belonging to $\overline{ \mathfrak{E}_{\varepsilon }}$. In the light of the condition (5), we get

$$\begin{aligned} \bigl\vert (\mathfrak{P} v) (z) \bigr\vert \leq {}& \frac{1}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} \bigl\vert \Phi _{*} \bigl(\varsigma, v(\varsigma ) \bigr) \bigr\vert \, \mathrm{d}\varsigma \\ &{}+ \frac{ s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1-\varsigma )^{ \xi - 2} \bigl\vert \Phi _{*} \bigl(\varsigma, v(\varsigma ) \bigr) \bigr\vert \, \mathrm{d}\varsigma \\ &{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \bigl\vert \Phi _{*} \bigl(\varsigma, v(\varsigma ) \bigr) \bigr\vert \,\mathrm{d}\varsigma \\ &{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma } \bigl\vert \Phi _{*} \bigl(\tau, v(\tau ) \bigr) \bigr\vert \,\mathrm{d}\tau \,\mathrm{d}\varsigma \\ \leq{}& \frac{z^{\xi }}{\Gamma ( \xi +1 )} \varphi (z) + \frac{ s_{2}^{*} }{ \Gamma ( \xi )} \varphi (z) \\ &{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \varphi (z) + \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2(2-\xi ) } \varphi (z) \end{aligned}$$

for any $z \in [0,1]$. As a consequence, the above estimate becomes $\Vert \mathfrak{P} v \Vert _{\mathbb{V}} \leq \hat{\Psi } \varphi ^{*} < \infty$, where

$$ \hat{\Psi } = \frac{1}{\Gamma ( \xi +1 )} + \frac{ s_{2}^{*} }{ \Gamma ( \xi )} + \frac{ 3( s_{1}^{*} + s_{2}^{*} + 1)\Gamma (4-\xi )}{ 2(2-\xi ) }. $$

(9)

Thus it is proved that $\mathfrak{P}(\overline{ \mathfrak{E}_{\varepsilon }})$ is a uniformly bounded set contained in $\mathbb{V}$. Now, we intend to verify that $\mathfrak{P}$ is equicontinuous. To confirm this, let $z_{1},z_{2} \in [0,1]$ with $z_{1}< z_{2}$ and $v \in \overline{ \mathfrak{E}_{\varepsilon }}$. Then, by assuming $\sup_{(z,v)\in [0,1]\times \overline{ \mathfrak{E}_{\varepsilon }}} \vert \Phi _{*} (z, v) \vert = \widetilde{\Phi }_{*} > 0$, we obtain

$$\begin{aligned} &\bigl\vert ( \mathfrak{P} v) (z_{2}) -(\mathfrak{P} v) (z_{1}) \bigr\vert \\ &\quad\leq \int _{0}^{z_{1}} \frac{ [ (z_{2}- \varsigma )^{ \xi -1} - (z_{1}- \varsigma )^{ \xi -1} ]}{\Gamma (\xi )} \bigl\vert \Phi _{*} \bigl(\varsigma, v(\varsigma ) \bigr) \bigr\vert \,\mathrm{d} \varsigma \\ &\qquad{}+ \int _{z_{1}}^{z_{2}} \frac{(z_{2}- \varsigma )^{ \xi -1}}{\Gamma (\xi )} \bigl\vert \Phi _{*} \bigl( \varsigma, v(\varsigma ) \bigr) \bigr\vert \,\mathrm{d} \varsigma \\ &\qquad{}+ \frac{ (z_{2} -z_{1})\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \bigl\vert \Phi _{*} \bigl(\varsigma, v(\varsigma ) \bigr) \bigr\vert \,\mathrm{d}\varsigma \\ &\qquad{}+ \frac{ (z_{2}-z_{1})\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma } \bigl\vert \Phi _{*} \bigl(\tau, v(\tau ) \bigr) \bigr\vert \,\mathrm{d}\tau \,\mathrm{d}\varsigma \\ &\quad\leq {\frac{\widetilde{\Phi }_{*}}{\Gamma (\xi +1)} \bigl( \bigl\vert z_{2}^{\xi }- z_{1}^{\xi } \bigr\vert + 2 (z_{1} - z_{2})^{\xi } \bigr)} \\ &\qquad{}+ {\frac{ (z_{2} -z_{1})\Gamma (4-\xi )\widetilde{\Phi }_{*} }{ 2-\xi } + \frac{ (z_{2}-z_{1})\Gamma (4-\xi )\widetilde{\Phi }_{*}}{ 2(2-\xi ) }. } \end{aligned}$$

When $z_{1}$ goes to $z_{2}$, the right-hand side of the above estimate approaches zero (independent of $v\in \overline{ \mathfrak{E}_{\varepsilon }}$). Thus, by assuming $z_{1} \to z_{2}$, we get $\Vert (\mathfrak{P} v)(z_{2}) -(\mathfrak{P} v)(z_{1}) \Vert _{ \mathbb{V} } \to 0$ and so $\mathfrak{P}$ is equicontinuous. Now, by the aid of the well-known theorem due to Arzela–Ascoli, $\mathfrak{P}$ is proved to be completely continuous, and so is compact on $\overline{ \mathfrak{E}_{\varepsilon }}$.

At this moment, we are going to verify that the operator $\mathfrak{P}$ is condensing on $\overline{ \mathfrak{E}_{\varepsilon }}$. By virtue of Lemma 2, we are sure that for each bounded subset $\mathfrak{B}\subset \overline{ \mathfrak{E}_{\varepsilon }}$, a countable set $\mathfrak{B}_{0}= \{v_{n}\}_{n \geq 1} \subset \mathfrak{B}$ exists such that the inequality $\Omega (\mathfrak{P}(\mathfrak{B})) \leq 2 \Omega (\mathfrak{P} ( \mathfrak{B}_{0}))$ is satisfied. Hence, by invoking Lemmas 1, 3, and 4, we have the following inequalities:

$$\begin{aligned} &\Omega \bigl(\mathfrak{P} \bigl(\mathfrak{B}(z) \bigr) \bigr) \\ &\quad\leq 2\Omega \bigl( \mathfrak{P} \bigl(\{v_{n} \}_{n\geq 1} \bigr) \bigr) \\ &\quad\leq \frac{2}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} \Omega \bigl( \Phi _{*} \bigl(\varsigma, \bigl\{ v_{n}( \varsigma ) \bigr\} _{n\geq 1} \bigr) \bigr) \,\mathrm{d}\varsigma \\ &\qquad{}+ \frac{ 2s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1- \varsigma )^{\xi - 2} \Omega \bigl( \Phi _{*} \bigl(\varsigma, \bigl\{ v_{n}( \varsigma ) \bigr\} _{n\geq 1} \bigr) \bigr) \,\mathrm{d}\varsigma \\ &\qquad{}+ \frac{ 2( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \Omega \bigl( \Phi _{*} \bigl(\varsigma, \bigl\{ v_{n}(\varsigma ) \bigr\} _{n\geq 1} \bigr) \bigr) \,\mathrm{d}\varsigma \\ &\qquad{}+ \frac{ 2( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }\Omega \bigl( \Phi _{*} \bigl(\tau, \bigl\{ v_{n}( \tau ) \bigr\} _{n\geq 1} \bigr) \bigr) \, \mathrm{d}\tau \,\mathrm{d}\varsigma \\ &\quad\leq \frac{4}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} n_{ \Phi _{*}}(\varsigma )\Omega \bigl( \bigl\{ v_{n}( \varsigma ) \bigr\} _{n\geq 1 } \bigr) \,\mathrm{d}\varsigma \\ &\qquad{}+ \frac{ 4s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1- \varsigma )^{\xi - 2} n_{ \Phi _{*}}(\varsigma )\Omega \bigl( \bigl\{ v_{n}( \varsigma ) \bigr\} _{n\geq 1 } \bigr) \,\mathrm{d}\varsigma \\ &\qquad{}+ \frac{ 4( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} n_{ \Phi _{*}}(\varsigma )\Omega \bigl( \bigl\{ v_{n}(\varsigma ) \bigr\} _{n \geq 1 } \bigr) \, \mathrm{d}\varsigma \\ &\qquad{}+ \frac{ 4( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }n_{ \Phi _{*}}(\tau )\Omega \bigl( \bigl\{ v_{n}( \tau ) \bigr\} _{n\geq 1 } \bigr) \,\mathrm{d}\tau \, \mathrm{d}\varsigma \\ &\quad\leq \frac{4 n_{ \Phi _{*}}^{*} \Omega (\mathfrak{B})}{\Gamma ( \xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} \, \mathrm{d}\varsigma + \frac{ 4s_{2}^{*} n_{ \Phi _{*}}^{*} \Omega (\mathfrak{B}) }{ \Gamma (\xi -1) } \int _{0}^{1} (1-\varsigma )^{\xi - 2} \, \mathrm{d}\varsigma \\ &\qquad{}+ \frac{ 4n_{ \Phi _{*}}^{*} \Omega (\mathfrak{B})( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \,\mathrm{d}\varsigma \\ &\qquad{}+ \frac{ 4n_{ \Phi _{*}}^{*} \Omega (\mathfrak{B})( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }\,\mathrm{d}\tau \,\mathrm{d}\varsigma \\ &\quad\leq \frac{4 n_{ \Phi _{*}}^{*} \Omega (\mathfrak{B})}{\Gamma ( \xi + 1 )} + \frac{ 4s_{2}^{*} n_{ \Phi _{*}}^{*} \Omega (\mathfrak{B}) }{ \Gamma (\xi ) } + \frac{ 12n_{ \Phi _{*}}^{*} \Omega (\mathfrak{B})( s_{1}^{*} + s_{2}^{*} + 1)\Gamma (4-\xi )}{ 2(2-\xi ) }. \end{aligned}$$

Hence,

$$\begin{aligned} \Omega \bigl(\mathfrak{P}(\mathfrak{B}) \bigr) &\leq 4 \biggl[ \frac{ n_{ \Phi _{*}}^{*} }{\Gamma ( \xi + 1 )} + \frac{ s_{2}^{*} n_{ \Phi _{*}}^{*} }{ \Gamma (\xi ) } + \frac{ 3n_{ \Phi _{*}}^{*} ( s_{1}^{*} + s_{2}^{*} + 1)\Gamma (4-\xi )}{ 2(2-\xi ) } \biggr] \Omega ( \mathfrak{B}). \end{aligned}$$

As a consequence, by taking into account the given condition (7), we get $\Omega (\mathfrak{P}(\mathfrak{B})) <\Omega (\mathfrak{B})$, and subsequently it is realized that $\mathfrak{P}$ is a condensing operator on $\overline{ \mathfrak{E}_{\varepsilon }}$. Eventually, by resorting to Theorem 6, it is deduced that the operator $\mathfrak{P}$ has at least one fixed point contained in $\overline{ \mathfrak{E}_{\varepsilon }}$ which is regarded as a solution for the given fractional Lio–CapBvp (1), and this ends the proof. □

At this moment, we are going to discuss the dependence of solutions for the given fractional Lio–CapBvp (1). Indeed, in this part we assume that the solution of the fractional Lio–CapBvp (1) depends on some parameters provided that the function $\Phi _{*}$ satisfies the conditions of Theorem 12 which guarantees the existence of solutions and, as we know, the continuous dependence of solutions on these parameters indicates the stability of solutions. So we investigate this property of the solutions of the fractional Lio–CapBvp (1) by making a small change in a parameter of the fractional Lio–CapBvp (1), i.e., its order. To do this, we will apply the generalized Gronwall inequality.

Theorem 13

Consider $\xi >0$ so that $1<\xi -\delta < \xi < 2$. Moreover, assume that $\Phi _{*}: [0,1] \times \mathbb{V} \to \mathbb{V}$ is continuous and a constant $\rho >0$ exists such that

$$ \bigl\vert \Phi _{*} \bigl(z,v(z) \bigr)- \Phi _{*} \bigl(z,v'(z) \bigr) \bigr\vert \leq \rho \bigl\vert v(z) - v'(z) \bigr\vert $$

for any $v,v'\in \mathbb{V}$ and $z\in [0,1]$. Additionally, let v be a solution of the fractional Lio–CapBvp (1) and v̄ be a solution of

$$ \textstyle\begin{cases} {}^{\mathcal{C}}\mathfrak{D}_{0}^{ \xi - \delta } \bar{v}(z) = \Phi _{*} (z, \bar{v}(z) ), \\ \bar{v}(0) = s_{1}^{*} {}^{\mathcal{C}}\mathfrak{D}_{0}^{1} \bar{v}(0) + s_{2}^{*} {}^{\mathcal{C}}\mathfrak{D}_{0}^{1} \bar{v}(1), \qquad {}^{ \mathcal{C}}\mathfrak{D}_{0}^{\xi - \delta -1} \bar{v}(1) = \int _{0}^{1} {}^{\mathcal{C}}\mathfrak{D}_{0}^{\xi - \delta -1} \bar{v}(\varsigma ) \,\mathrm{d}\varsigma. \end{cases} $$

(10)

Then, the following inequality is valid:

$$ \Vert v -\bar{v} \Vert _{\mathbb{V}} \leq \frac{\Delta + \Delta \sum_{k=1}^{\infty } \frac{\rho ^{k}}{\Gamma (k(\xi -\delta )+1)}}{1- \mho - \mho \sum_{k=1}^{\infty } \frac{\rho ^{k}}{\Gamma (k(\xi -\delta )+1)} } $$

provided that $\mho + \mho \sum_{k=1}^{\infty } \frac{\rho ^{k}}{\Gamma (k(\xi -\delta )+1)}<1$, where

$$\begin{aligned} \Delta ={}&\Vert \Phi _{*} \Vert _{\mathbb{V}} \sup_{z\in [0,1]} \biggl\vert \frac{z^{ \xi }}{\Gamma (\xi +1 )} - \frac{z^{ \xi -\delta }}{\Gamma (\xi -\delta +1)} \biggr\vert + \Vert \Phi _{*} \Vert _{\mathbb{V}} \sup _{z\in [0,1]} \biggl\vert \frac{s_{2}^{*} }{\Gamma (\xi )} - \frac{s_{2}^{*}}{\Gamma (\xi - \delta )} \biggr\vert \\ &{}+ \Vert \Phi _{*} \Vert _{\mathbb{V}} \sup _{z\in [0,1]} \biggl\vert \frac{(s_{1}^{*} + s_{2}^{*} +z)\Gamma (4-\xi )}{2-\xi } - \frac{(s_{1}^{*} + s_{2}^{*} +z) \Gamma (4- (\xi - \delta ))}{2- (\xi - \delta )} \biggr\vert \\ &{}+ \Vert \Phi _{*} \Vert _{\mathbb{V}} \sup _{z\in [0,1]} \biggl\vert \frac{(s_{1}^{*} + s_{2}^{*} +z)\Gamma (4-\xi )}{2(2-\xi )} - \frac{(s_{1}^{*} + s_{2}^{*} +z) \Gamma (4- (\xi - \delta ))}{2(2- (\xi - \delta ))} \biggr\vert \end{aligned}$$

(11)

and

$$ \mho = \biggl( \frac{s_{2}^{*}}{\Gamma (\xi -\delta )} + \frac{3(s_{1}^{*} + s_{2}^{*} + 1)\Gamma (4- (\xi - \delta ))}{2(2-(\xi -\delta )) } \biggr) \rho,$$

(12)

where $\Vert \Phi _{*} \Vert _{\mathbb{V}} = \sup_{z\in [0,1]} | \Phi _{*} ( z,v(z) ) |$.

Proof

Prior to deriving the desired inequality, in view of the above facts, it is known that the existence of solutions of both fractional Lio–CapBvps (1) and (10) is guaranteed by a similar argument as above and these solutions are given in the form of (8) and

$$\begin{aligned} \bar{v}(z) = {}&\frac{1}{\Gamma (\xi - \delta )} \int _{0}^{z} (z- \varsigma )^{ \xi - \delta -1} \Phi _{*} \bigl(\varsigma, \bar{v}( \varsigma ) \bigr) \,\mathrm{d} \varsigma \\ &{}+ \frac{ s_{2}^{*} }{ \Gamma (\xi - \delta -1) } \int _{0}^{1} (1-\varsigma )^{\xi - \delta - 2} \Phi _{*} \bigl(\varsigma, \bar{v}( \varsigma ) \bigr) \,\mathrm{d} \varsigma \\ &{}- \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4- (\xi - \delta ) )}{ 2- (\xi - \delta ) } \int _{0}^{1} \Phi _{*} \bigl(\varsigma, \bar{v}(\varsigma ) \bigr) \,\mathrm{d}\varsigma \\ &{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4- (\xi - \delta ))}{ 2-(\xi - \delta ) } \int _{0}^{1} \int _{0}^{\varsigma }\Phi _{*} \bigl(\tau, \bar{v}(\tau ) \bigr) \,\mathrm{d}\tau \,\mathrm{d}\varsigma, \end{aligned}$$

respectively. Then, an estimate for $v-\bar{v}$ is calculated as follows:

$$\begin{aligned} &\bigl\vert v(z) - \bar{v}(z) \bigr\vert \\ &\quad\leq \biggl\vert \int _{0}^{z} \biggl( \frac{1}{\Gamma (\xi )} (z- \varsigma )^{ \xi -1} - \frac{1}{\Gamma (\xi -\delta )} (z- \varsigma )^{ \xi -\delta -1} \biggr) \Phi _{*} \bigl( \varsigma,v(\varsigma ) \bigr) \,\mathrm{d} \varsigma \biggr\vert \\ &\qquad{}+ \int _{0}^{z} \frac{(z- \varsigma )^{ \xi -\delta -1}}{\Gamma (\xi -\delta )} \bigl\vert \Phi _{*} \bigl(\varsigma,v(\varsigma ) \bigr) - \Phi _{*} \bigl( \varsigma, \bar{v}(\varsigma ) \bigr) \bigr\vert \,\mathrm{d}\varsigma \\ &\qquad{}+ \biggl\vert \int _{0}^{1} \biggl( \frac{s_{2}^{*}(1- \varsigma )^{ \xi -2}}{\Gamma (\xi -1 )} - \frac{s_{2}^{*}(1- \varsigma )^{ \xi -\delta -2}}{\Gamma (\xi -\delta - 1 )} \biggr) \Phi _{*} \bigl(\varsigma,v(\varsigma ) \bigr) \,\mathrm{d}\varsigma \biggr\vert \\ &\qquad{} + \int _{0}^{1} \frac{s_{2}^{*}(1- \varsigma )^{ \xi -\delta -2}}{\Gamma (\xi -\delta - 1 )} \bigl\vert \Phi _{*} \bigl(\varsigma,v(\varsigma ) \bigr) - \Phi _{*} \bigl( \varsigma, \bar{v}(\varsigma ) \bigr) \bigr\vert \,\mathrm{d}\varsigma \\ &\qquad{}+ \biggl\vert \frac{(s_{1}^{*} + s_{2}^{*} +z)\Gamma (4-\xi )}{2-\xi } \int _{0}^{1} \Phi _{*} \bigl( \varsigma,v(\varsigma ) \bigr) \,\mathrm{d}\varsigma \\ &\qquad{}- \frac{(s_{1}^{*} + s_{2}^{*} +z)\Gamma (4- (\xi - \delta ) )}{2- (\xi - \delta ) } \int _{0}^{1} \Phi _{*} \bigl( \varsigma,v(\varsigma ) \bigr) \,\mathrm{d} \varsigma \biggr\vert \\ &\qquad{}+ \frac{(s_{1}^{*} + s_{2}^{*} +z)\Gamma (4- (\xi - \delta ) )}{2- (\xi - \delta ) } \int _{0}^{1} \bigl\vert \Phi _{*} \bigl(\varsigma,v(\varsigma ) \bigr) - \Phi _{*} \bigl( \varsigma, \bar{v}(\varsigma ) \bigr) \bigr\vert \,\mathrm{d}\varsigma \\ &\qquad{}+ \biggl\vert \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4- \xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }\Phi _{*} \bigl(\tau, v( \tau ) \bigr) \,\mathrm{d} \tau \,\mathrm{d}\varsigma \\ &\qquad{}- \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4- (\xi - \delta ))}{ 2-(\xi - \delta ) } \int _{0}^{1} \int _{0}^{\varsigma }\Phi _{*} \bigl(\tau, v( \tau ) \bigr) \,\mathrm{d}\tau \,\mathrm{d}\varsigma \biggr\vert \\ &\qquad{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4- (\xi - \delta ))}{ 2-(\xi - \delta ) } \int _{0}^{1} \int _{0}^{\varsigma } \bigl\vert \Phi _{*} \bigl(\tau,v(\tau ) \bigr) - \Phi _{*} \bigl(\tau,\bar{v}(\tau ) \bigr) \bigr\vert \,\mathrm{d}\tau \,\mathrm{d}\varsigma \\ &\quad\leq \Delta + \mho \Vert v-\bar{v} \Vert _{\mathbb{V}} + \int _{0}^{z} \frac{(z- \varsigma )^{ \xi -\delta -1}}{\Gamma (\xi - \delta )} \rho \bigl\vert v(\varsigma )- \bar{v}(\varsigma ) \bigr\vert \,\mathrm{d}\varsigma, \end{aligned}$$

where Δ and ℧ are introduced by (11) and (12). Hence with due attention to the generalized Gronwall inequality (Theorem 7) and by letting $\breve{v}(z)= \vert v(z) -\bar{v}(z) \vert$, $\breve{r}(z)=\Delta +\mho \Vert v- \bar{v} \Vert _{\mathbb{V}}$, and $h(z)=\frac{\rho }{\Gamma (\xi -\delta )}$, we reach

$$ \bigl\vert v(z) -\bar{v}(z) \bigr\vert \leq \Delta +\mho \Vert v- \bar{v} \Vert _{\mathbb{V}} + \int _{0}^{z} \sum_{k=1}^{\infty } \biggl[ \frac{\rho ^{k}(z- \varsigma )^{ k(\xi -\delta ) -1}}{\Gamma (k(\xi -\delta ))} \bigl(\Delta +\mho \Vert v-\bar{v} \Vert _{\mathbb{V}} \bigr) \biggr] \,\mathrm{d}\varsigma. $$

In conclusion, we obtain

$$ \Vert v -\bar{v} \Vert _{\mathbb{V}} \leq \frac{\Delta +\Delta \sum_{k=1}^{\infty } \frac{\rho ^{k}}{\Gamma (k(\xi -\delta )+1)}}{1-\mho - \mho \sum_{k=1}^{\infty } \frac{\rho ^{k}}{\Gamma (k(\xi -\delta )+1)} }, $$

and the latter inequality finishes the proof. □

Remark 1

From the above theorem, we deduce that if $\delta =0$, then $v=\bar{v}$ and so both considered solutions of the fractional Lio–CapBvp (1) are the same.

4 Inclusion version of the fractional Lio–CapBvp (1)

Here, we continue the above process to establish a criterion of existence of solutions for the inclusion version of the fractional Lio–CapBvp (1) given by

$$ \textstyle\begin{cases} {}^{\mathcal{C}}\mathfrak{D}_{0}^{ \xi } v(z) \in \mathfrak{Q}_{*} (z, v(z) ), \\ v(0) = s_{1}^{*} {}^{\mathcal{C}}\mathfrak{D}_{0}^{1} v(0) + s_{2}^{*} {}^{\mathcal{C}}\mathfrak{D}_{0}^{1} v(1), \qquad {}^{\mathcal{C}} \mathfrak{D}_{0}^{\xi -1} v(1) = \int _{0}^{1} {}^{ \mathcal{C}}\mathfrak{D}_{0}^{\xi -1} v(\varsigma ) \,\mathrm{d} \varsigma \end{cases} $$

(13)

so that $z\in [0,1]$, $s_{1}^{*}, s_{2}^{*} \in \mathbb{R}^{+}$, ${}^{\mathcal{C}}\mathfrak{D}_{0}^{ \xi }$ stands for the Cap-derivative of order $\xi \in (1,2)$, ${}^{\mathcal{C}}\mathfrak{D}_{0}^{1} = \frac{\,\mathrm{d}}{\,\mathrm{d}z}$ and $\mathfrak{Q}_{*}: [0,1] \times \mathbb{V} \to \mathcal{P}( \mathbb{V})$ is a multivalued mapping. To continue, we utilize the notion of the approx-endpoint property.

Definition 14

An absolutely continuous function $v:[0,1]\to \mathbb{R}$ is called a solution for the fractional inclusion Lio–CapBvp (13) whenever there is an $\hbar \in \mathcal{L}^{1}([0,1],\mathbb{R})$ such that $\hbar (z)\in \mathfrak{Q}_{*} (z, v(z) )$ for almost all $z \in [0,1]$ which satisfies integro-derivative boundary conditions

$$ v(0) = s_{1}^{*} {}^{\mathcal{C}}\mathfrak{D}_{0}^{1} v(0) + s_{2}^{*} {}^{\mathcal{C}}\mathfrak{D}_{0}^{1} v(1), \qquad{}^{\mathcal{C}} \mathfrak{D}_{0}^{\xi -1} v(1) = \int _{0}^{1} {}^{ \mathcal{C}} \mathfrak{D}_{0}^{\xi -1} v(\varsigma ) \,\mathrm{d} \varsigma $$

and we also have

$$\begin{aligned} v(z) = {}&\frac{1}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} \hbar (\varsigma ) \,\mathrm{d}\varsigma + \frac{ s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1- \varsigma )^{\xi - 2} \hbar (\varsigma ) \,\mathrm{d}\varsigma \\ &{}- \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \hbar (\varsigma ) \,\mathrm{d} \varsigma \\ &{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }\hbar (\tau ) \,\mathrm{d}\tau \, \mathrm{d} \varsigma \end{aligned}$$

for any $z \in [0,1]$.

For any $v\in \mathbb{V}$, the family of all existing selections of $\mathfrak{Q}_{*}$ is denoted by

$$ (\mathfrak{SEL})_{\mathfrak{Q}_{*},v}= \bigl\{ \hbar \in \mathcal{L}^{1} \bigl([0,1] \bigr): \hbar (z)\in \mathfrak{Q}_{*} \bigl(z, v(z) \bigr) \bigr\} $$

for almost all $z\in [0,1]$. Now, we define a multifunction $\mathfrak{F}: \mathbb{V}\to \mathcal{P}(\mathbb{V})$ by the following rule:

$$ \mathfrak{F}(v)= \bigl\{ \vartheta \in \mathbb{V}: \vartheta (z)= \kappa (z) \bigr\} , $$

(14)

where

$$\begin{aligned} \kappa (z) ={}&\frac{1}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} \hbar (\varsigma ) \,\mathrm{d}\varsigma + \frac{ s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1- \varsigma )^{\xi - 2} \hbar (\varsigma ) \,\mathrm{d}\varsigma \\ &{}- \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \hbar (\varsigma ) \,\mathrm{d} \varsigma \\ &{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }\hbar (\tau ) \,\mathrm{d}\tau \, \mathrm{d} \varsigma, \hbar \in (\mathfrak{SEL})_{\mathfrak{Q}_{*},v}. \end{aligned}$$

By utilizing the approx-endpoint property for the multivalued map $\mathfrak{F}$, we have the following:

Theorem 15

Regard $\mathfrak{Q}_{*}: [0,1]\times \mathbb{V}\to \mathcal{P}_{cmp}( \mathbb{V})$ as a multivalued map. Suppose that

($\mathfrak{C}_{1}$):

$\psi: [0, \infty ) \to [0, \infty )$ is increasing and u.s.c. with $\lim \inf_{z\to \infty }(z-\psi (z))>0$ and $\psi (z)< z$ for any $z>0$;

($\mathfrak{C}_{2}$):

the multifunction $\mathfrak{Q}_{*}: [0,1]\times \mathbb{V}\to \mathcal{P}_{cmp}( \mathbb{V})$ is integrable and bounded so that $\mathfrak{Q}_{*}(\cdot, v):[0,1]\to \mathcal{P}_{cmp}(\mathbb{V})$ is measurable for any $v\in \mathbb{V}$;

($\mathfrak{C}_{3}$):

a function $\varrho \in \mathcal{C}([0,1], [0, \infty ))$ exists such that

$$ \mathbb{H}_{d_{\mathbb{V}}} \bigl(\mathfrak{Q}_{*} \bigl(z, v_{1}(z) \bigr), \mathfrak{Q}_{*} \bigl(z, v_{2}(z) \bigr) \bigr)\leq \varrho (z)\psi \bigl( \bigl\vert v_{1}(z)-v_{2}(z) \bigr\vert \bigr) \frac{1}{ \hat{\mathbb{O}} } $$

for all $z\in [0,1]$ and $v_{1}, v_{2}\in \mathbb{V}$, where $\sup_{z\in [0,1]}|\varrho (z)|=\|\varrho \|$ and

$$ \hat{\mathbb{O}} = \biggl[ \frac{1}{ \Gamma ( \xi +1)} + \frac{ s_{2}^{*} }{ \Gamma ( \xi )} + \frac{ 3( s_{1}^{*} + s_{2}^{*} + 1)\Gamma (4-\xi )}{ 2(2-\xi ) } \biggr] \Vert \varrho \Vert ; $$

(15)

($\mathfrak{C}_{4}$):

the multifunction $\mathfrak{F}$ given by (14) has the approx-endpoint property.

Then a solution exists for the fractional inclusion Lio–CapBvp (13).

Proof

In the current proof, we try to confirm the existence of an endpoint for the multivalued map $\mathfrak{F}: \mathbb{V}\to \mathcal{P}(\mathbb{V})$ defined by (14). To proceed, we first need to check that $\mathfrak{F}(v)$ is closed for any $v\in \mathbb{V}$. In the light of the hypothesis ($\mathfrak{C}_{2}$), the mapping $z\mapsto \mathfrak{Q}_{*}(z, v(z))$ is a measurable and closed-valued multivalued map for each $v\in \mathbb{V}$. As a consequence, $\mathfrak{Q}_{*}$ has a measurable selection as $(\mathfrak{SEL})_{\mathfrak{Q}_{*},v}\neq \emptyset$. Now, we claim that $\mathfrak{F}(v)\subseteq \mathbb{V}$ is closed for all $v\in \mathbb{V}$. Consider a sequence $(v_{n})_{n\geq 1}$ contained in $\mathfrak{F}(v)$ with $v_{n} \to v^{*}$. For each n, an element $\hbar _{n} \in (\mathfrak{SEL})_{\mathfrak{Q}_{*},v}$ exists so that

$$\begin{aligned} v_{n}(z) = {}&\frac{1}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} \hbar _{n}(\varsigma ) \,\mathrm{d}\varsigma + \frac{ s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1- \varsigma )^{\xi - 2} \hbar _{n}(\varsigma ) \,\mathrm{d}\varsigma \\ &{}- \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \hbar _{n}(\varsigma ) \, \mathrm{d}\varsigma \\ &{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }\hbar _{n}(\tau ) \, \mathrm{d}\tau \,\mathrm{d} \varsigma \end{aligned}$$

for almost all $z\in [0,1]$. Since $\mathfrak{Q}_{*}$ is a compact multivalued mapping, we get a subsequence $\{\hbar _{n}\}_{n\geq 1}$ tending to $\hbar \in \mathcal{L}^{1}([0,1])$. Hence, we have $\hbar \in (\mathfrak{SEL})_{\mathfrak{Q}_{*},v}$ and

$$\begin{aligned} \lim_{n\to \infty }v_{n}(z) ={}& \frac{1}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} \hbar (\varsigma ) \,\mathrm{d}\varsigma + \frac{ s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1-\varsigma )^{ \xi - 2} \hbar ( \varsigma ) \,\mathrm{d}\varsigma \\ &{}- \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \hbar (\varsigma ) \,\mathrm{d} \varsigma \\ &{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }\hbar (\tau ) \,\mathrm{d}\tau \, \mathrm{d} \varsigma \\ ={}&v(z) \end{aligned}$$

for any $z\in [0,1]$. Thus $v \in \mathfrak{F}$ which confirms that $\mathfrak{F}$ is closed-valued. Therefore, we know that $\mathfrak{F}(v)$ is bounded for each $v \in \mathbb{V}$ since $\mathfrak{Q}_{*}$ is assumed to be compact. Now, we check whether the inequality

$$ \mathbb{H}_{d_{\mathbb{V}}} \bigl(\mathfrak{F}(v_{1}), \mathfrak{F}(v_{2}) \bigr) \leq \psi \bigl( \Vert v_{1} - v_{2} \Vert \bigr) $$

holds. Let $v_{1},v_{2}\in \mathbb{V}$ and $x_{1}\in \mathfrak{F}(v_{2})$. Select $\hbar _{1}\in (\mathfrak{SEL})_{\mathfrak{Q}_{*},v_{2}}$ so that

$$\begin{aligned} x_{1}(z)={}& \frac{1}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} \hbar _{1}(\varsigma ) \,\mathrm{d}\varsigma + \frac{ s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1- \varsigma )^{\xi - 2} \hbar _{1}(\varsigma ) \,\mathrm{d}\varsigma \\ &{}- \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \hbar _{1}(\varsigma ) \, \mathrm{d}\varsigma \\ &{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }\hbar _{1}(\tau ) \, \mathrm{d}\tau \,\mathrm{d} \varsigma \end{aligned}$$

for almost all $z\in [0,1]$. Since

$$ \mathbb{H}_{d_{\mathbb{V}}} \bigl(\mathfrak{Q}_{*} \bigl(z, v_{1}(z) \bigr), \mathfrak{Q}_{*} \bigl(z, v_{2}(z) \bigr) \bigr)\leq \varrho (z) \bigl(\psi \bigl(v_{1}(z)- v_{2}(z) \bigr) \bigr)\frac{1}{ \hat{\mathbb{O}} } $$

for any $z\in [0,1]$, we obtain that $\hbar ^{*}\in \mathfrak{Q}_{*}(z, v_{1}(z))$ exists such that

$$ \bigl\vert \hbar _{1}(z)- \hbar ^{*} \bigr\vert \leq \varrho (z) \bigl(\psi \bigl( v_{1}(z)- v_{2}(z) \bigr) \bigr)\frac{1}{ \hat{\mathbb{O}} } $$

for any $z\in [0,1]$. In the following, consider the multivalued map $\mathfrak{A}:[0,1]\to \mathcal{P}(\mathbb{V})$ introduced by

$$ \mathfrak{A}(z)= \biggl\{ \hbar ^{*}\in \mathbb{V}: \bigl\vert \hbar _{1}(z)- \hbar ^{*} \bigr\vert \leq \varrho (z) \bigl(\psi \bigl(v_{1}(z)- v_{2}(z) \bigr) \bigr) \frac{1}{ \hat{\mathbb{O}} } \biggr\} . $$

Since $\hbar _{1}$ and $\sigma =\varrho (\psi (v_{1} -v_{2} )) \frac{1}{ \hat{\mathbb{O}} }$ are measurable, we can choose $\hbar _{2}(z) \in \mathfrak{Q}_{*}(z, v_{1}(z))$ so that

$$ \bigl\vert \hbar _{1}(z)-\hbar _{2}(z) \bigr\vert \leq \varrho (z) \bigl(\psi \bigl(v_{1}(z)-v_{2}(z) \bigr) \bigr) \frac{1}{ \hat{\mathbb{O}} } $$

for all $z\in [0,1]$. Select $x_{2}\in \mathfrak{F}(v_{1})$ such that

$$\begin{aligned} x_{2}(z)={}& \frac{1}{\Gamma (\xi )} \int _{0}^{z} (z- \varsigma )^{ \xi -1} \hbar _{2}(\varsigma ) \,\mathrm{d}\varsigma + \frac{ s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1- \varsigma )^{\xi - 2} \hbar _{2}(\varsigma ) \,\mathrm{d}\varsigma \\ &{}- \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \hbar _{2}(\varsigma ) \, \mathrm{d}\varsigma \\ &{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma }\hbar _{2}(\tau ) \, \mathrm{d}\tau \,\mathrm{d} \varsigma \end{aligned}$$

for each $z\in [0,1]$. Hence, we get that

$$\begin{aligned} &\bigl\vert x_{1}(z) -x_{2}(z) \bigr\vert \\ &\quad\leq \frac{1}{ \Gamma ( \xi ) } \int _{0}^{z} (z-\varsigma )^{\xi -1} \bigl\vert \hbar _{1}( \varsigma )-\hbar (\varsigma ) \bigr\vert \, \mathrm{d}\varsigma \\ &\qquad{}+ \frac{ s_{2}^{*} }{ \Gamma (\xi -1) } \int _{0}^{1} (1- \varsigma )^{\xi - 2} \bigl\vert \hbar _{1}(\varsigma )-\hbar (\varsigma ) \bigr\vert \, \mathrm{d}\varsigma \\ &\qquad{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \bigl\vert \hbar _{1}( \varsigma )-\hbar (\varsigma ) \bigr\vert \,\mathrm{d}\varsigma \\ &\qquad{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + z)\Gamma (4-\xi )}{ 2-\xi } \int _{0}^{1} \int _{0}^{\varsigma } \bigl\vert \hbar _{1}( \tau )-\hbar (\tau ) \bigr\vert \,\mathrm{d} \tau \,\mathrm{d}\varsigma \\ &\quad\leq \frac{1}{ \Gamma ( \xi +1)} \Vert \varrho \Vert \psi \bigl( \Vert v_{1}- v_{2} \Vert \bigr)\frac{1}{ \hat{\mathbb{O}} } + \frac{ s_{2}^{*} }{ \Gamma (\xi ) } \Vert \varrho \Vert \psi \bigl( \Vert v_{1}- v_{2} \Vert \bigr)\frac{1}{ \hat{\mathbb{O}} } \\ &\qquad{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + 1)\Gamma (4-\xi )}{ 2-\xi } \Vert \varrho \Vert \psi \bigl( \Vert v_{1}- v_{2} \Vert \bigr) \frac{1}{ \hat{\mathbb{O}} } \\ &\qquad{}+ \frac{ ( s_{1}^{*} + s_{2}^{*} + 1)\Gamma (4-\xi )}{ 2(2-\xi ) } \Vert \varrho \Vert \psi \bigl( \Vert v_{1}- v_{2} \Vert \bigr) \frac{1}{ \hat{\mathbb{O}} } \\ &\quad= \biggl[ \frac{1}{ \Gamma ( \xi +1)} + \frac{ s_{2}^{*} }{ \Gamma ( \xi )} + \frac{ 3( s_{1}^{*} + s_{2}^{*} + 1)\Gamma (4-\xi )}{ 2(2-\xi ) } \biggr] \Vert \varrho \Vert \psi \bigl( \Vert v_{1} - v_{2} \Vert \bigr) \frac{1}{ \hat{\mathbb{O}} } \\ &\quad= \hat{\mathbb{O}} \psi \bigl( \Vert v_{1} - v_{2} \Vert \bigr) \frac{1}{ \hat{\mathbb{O}} }=\psi \bigl( \Vert v_{1}- v_{2} \Vert \bigr). \end{aligned}$$

This gives $\|x_{1}-x_{2}\|\leq \psi (\| v_{1} - v_{2} \| )$ and yields $\mathbb{H}_{d_{\mathbb{V}}}(\mathfrak{F}(v_{1}), \mathfrak{F}(v_{2})) \leq \psi (\|v_{1} - v_{2} \| )$ for any $v_{1}, v_{2} \in \mathbb{V}$. Also from ($\mathfrak{C}_{4}$), we understand that $\mathfrak{F}$ possesses the approx-endpoint property. Then Theorem 10 gives that $\mathfrak{F}$ has a unique endpoint. In other words, $v^{*} \in \mathbb{V}$ exists so that $\mathfrak{F}(v^{*} )=\{v^{*} \}$. As a consequence, $v^{*}$ is a solution of the fractional inclusion Lio–CapBvp (13) and the proof is finished. □

5 Numerical examples

Here, we prepare two numerical examples to confirm the validity of the analytical findings. The first example illustrates Theorem 12.

Example 1

With due attention to the proposed fractional Lio–CapBvp (1), we formulate the following structure of the fractional Lio–CapBvp as follows:

$$ \textstyle\begin{cases} {}^{\mathcal{C}}\mathfrak{D}_{0}^{ 1.5 } v(z) = \frac{3 + 2e^{-z}}{10\text{,}000} \sin (v(z)) \\ v(0) = 0.3 {}^{\mathcal{C}}\mathfrak{D}_{0}^{1} v(0) + 0.4 {}^{ \mathcal{C}}\mathfrak{D}_{0}^{1} v(1), \qquad {}^{\mathcal{C}} \mathfrak{D}_{0}^{0.5} v(1) = \int _{0}^{1} {}^{ \mathcal{C}}\mathfrak{D}_{0}^{0.5} v(\varsigma ) \,\mathrm{d} \varsigma \end{cases} $$

(16)

so that $\xi = 1.5$, $\xi -1 = 0.5$, $s_{1}^{*} = 0.3$, $s_{2}^{*} = 0.4$, and $z\in [0,1]$. In addition, we introduce a continuous function $\Phi _{*}: [0,1]\times \mathbb{R} \to \mathbb{R}$ given by

$$ \Phi _{*} \bigl(z, v(z) \bigr) = \frac{3 + 2e^{-z}}{10\text{,}000} \sin \bigl(v(z) \bigr). $$

Then, for any $v\in \mathbb{R}$, one can write

$$ \bigl\vert \Phi _{*} \bigl(z, v(z) \bigr) \bigr\vert \leq \frac{3 + 2e^{-z}}{10\text{,}000} \bigl\vert \sin \bigl(v(z) \bigr) \bigr\vert \leq \frac{3 + 2e^{-z}}{10\text{,}000} = \varphi (z), $$

where a continuous function $\varphi: [0,1] \to \mathbb{R}^{+}$ is defined by $\varphi (z) = \frac{3 + 2e^{-z}}{10\text{,}000}$. On the other hand, for every $v_{1}, v_{2} \in \mathbb{R}$, we have

$$\begin{aligned} \bigl\vert \Phi _{*} \bigl(z, v_{1}(z) \bigr) - \Phi _{*} \bigl(z, v_{2}(z) \bigr) \bigr\vert &\leq \frac{3 + 2e^{-z}}{10\text{,}000} \bigl\vert \sin \bigl(v_{1}(z) \bigr)- \sin \bigl(v_{2}(z) \bigr) \bigr\vert \\ &\leq \frac{3 + 2e^{-z}}{10\text{,}000} \bigl\vert v_{1}(z)- v_{2}(z) \bigr\vert . \end{aligned}$$

As a consequence, for any bounded set $\mathfrak{B}\subset \mathbb{R}$, we may write

$$ \Omega \bigl( \Phi _{*}(z,\mathfrak{B}) \bigr) \leq \frac{3 + 2e^{-z}}{10\text{,}000} \Omega (\mathfrak{B} ):= n_{ \Phi _{*}}(z)\Omega ( \mathfrak{B} ) $$

so that $n_{\Phi _{*}}^{*} = \sup_{z\in [0,1]} \vert n_{\Phi _{*}} \vert = 0.0005$. Then in view of above values, since

$$ \frac{ n_{ \Phi _{*}}^{*} }{\Gamma ( \xi + 1 )} + \frac{ s_{2}^{*} n_{ \Phi _{*}}^{*} }{ \Gamma (\xi ) } + \frac{ 3n_{ \Phi _{*}}^{*} ( s_{1}^{*} + s_{2}^{*} + 1)\Gamma (4-\xi )}{ 2(2-\xi ) } \simeq 0.00399152 < \frac{1}{4}, $$

the inequality (7) holds. In conclusion, all the hypotheses of Theorem 12 are valid and thus this result implies that there exists at least one solution for the fractional Lio–CapBvp (16).

The next example illustrates Theorem 15.

Example 2

With due attention to (13), we design the fractional inclusion Cap–Bvp in the form

$$ \textstyle\begin{cases} {}^{\mathcal{C}}\mathfrak{D}_{0}^{ 1.5 } v(z) \in [0, \frac{e^{-2z} \vert \arctan (v(z)) \vert }{12(50 + z)(1+ \vert \arctan (v(z)) \vert )} ], \\ v(0) = 0.3 {}^{\mathcal{C}}\mathfrak{D}_{0}^{1} v(0) + 0.4 {}^{ \mathcal{C}}\mathfrak{D}_{0}^{1} v(1), \qquad {}^{\mathcal{C}} \mathfrak{D}_{0}^{0.5} v(1) = \int _{0}^{1} {}^{ \mathcal{C}}\mathfrak{D}_{0}^{0.5} v(\varsigma ) \,\mathrm{d} \varsigma \end{cases} $$

(17)

with the same values $\xi = 1.5$, $\xi -1 = 0.5$, $s_{1}^{*} = 0.3$, $s_{2}^{*} = 0.4$ and $z\in [0,1]$. At first, we introduce the Banach space $\mathbb{V}=\{ v(z): v(z)\in \mathcal{C}_{\mathbb{R}}([0, 1])\}$ furnished with $\| v \|_{\mathbb{V}}=\sup_{z\in [0, 1]}| v(z)|$. In addition, we define a multifunction $\mathfrak{Q}_{*}: [0, 1]\times \mathbb{V}\to \mathcal{P}(\mathbb{V})$ by

$$ \mathfrak{Q}_{*} \bigl(z, v(z) \bigr)= \biggl[0, \frac{e^{-2z} \vert \arctan (v(z)) \vert }{12(50 + z)(1+ \vert \arctan (v(z)) \vert )} \biggr] $$

for all $z\in [0, 1]$. In the next stage, we regard an increasing u.s.c. mapping $\psi: [0, \infty )\to [0, \infty )$ by $\psi (z)=\frac{z}{3}$ for all $z>0$. It is obvious that $\lim \inf_{z\to \infty }(z-\psi (z))>0$ and $\psi (z)< z$ for all $z>0$. At this moment, for any $v_{1}, v_{2}\in \mathbb{V}$, we have

$$\begin{aligned} &\mathbb{H}_{d_{\mathbb{V}}} \bigl(\mathfrak{Q}_{*} \bigl(z, v_{1}(z) \bigr), \mathfrak{Q}_{*} \bigl(z, v_{2}(z) \bigr) \bigr) \\ &\quad\leq \frac{e^{-2z}}{12(50 +z)} \bigl( \bigl\vert \arctan \bigl(v_{1}(z) \bigr) - \arctan \bigl(v_{2}(z) \bigr) \bigr\vert \bigr) \\ &\quad\leq \frac{e^{-2z}}{12(50 +z)} \bigl( \bigl\vert v_{1}(z)- v_{2}(z) \bigr\vert \bigr) \\ &\quad= \frac{e^{-2z}}{4(50 +z)} \psi \bigl( \bigl\vert v_{1}(z)- v_{2}(z) \bigr\vert \bigr) \\ &\quad\leq \varrho (z)\psi \bigl( \bigl\vert v_{1}(z)- v_{2}(z) \bigr\vert \bigr) \frac{1}{\hat{\mathbb{O}}}, \end{aligned}$$

where

$$ \hat{\mathbb{O}} = \biggl[ \frac{1}{ \Gamma ( \xi +1)} + \frac{ s_{2}^{*} }{ \Gamma ( \xi )} + \frac{ 3( s_{1}^{*} + s_{2}^{*} + 1)\Gamma (4-\xi )}{ 2(2-\xi ) } \biggr] \Vert \varrho \Vert \simeq 0.0399152785, $$

and we obtain $\varrho \in \mathcal{C}([0, 1], [0, \infty ))$ given by $\varrho (z)= \frac{e^{-2z}}{4(50 +z)}$ for all z. In this phase, we get $\|\varrho \|=\sup_{z\in [0, 1]} \vert \varrho (z) \vert = 0.005$. Finally, we introduce the multifunction $\mathfrak{F}: \mathbb{V}\to \mathcal{P}(\mathbb{V})$ by

$$ \mathfrak{F}(v)= \bigl\{ \vartheta \in \mathbb{V}: \text{ there exists } \hbar \in ( \mathfrak{SEL})_{\mathfrak{Q}_{*},v} \text{ s.t. } \vartheta (z)= \kappa (z), \forall z \in [0, 1] \bigr\} , $$

where

$$\begin{aligned} \kappa (z) ={}& \frac{1}{\Gamma (1.5 )} \int _{0}^{z} (z- \varsigma )^{ 1.5 - 1} \hbar (\varsigma ) \,\mathrm{d}\varsigma + \frac{ 0.4 }{ \Gamma (0.5) } \int _{0}^{1} (1- \varsigma )^{1.5 - 2} \hbar (\varsigma ) \,\mathrm{d}\varsigma \\ &{}- \frac{ ( 0.7 + z)\Gamma (2.5 )}{ 0.5 } \int _{0}^{1} \hbar ( \varsigma ) \,\mathrm{d} \varsigma + \frac{ ( 0.7 + z)\Gamma (2.5)}{ 0.5 } \int _{0}^{1} \int _{0}^{\varsigma }\hbar (\tau ) \,\mathrm{d}\tau \, \mathrm{d}\varsigma. \end{aligned}$$

Thus we observe that all conditions of Theorem 15 are valid. Hence the fractional inclusion Lio–CapBvp (17) has a solution.

6 Conclusion

In this paper, we designed a new Liouville–Caputo fractional boundary value problem with mixed integro-derivative boundary conditions. To obtain some criteria establishing the existence of solutions for the proposed problem, we used condensing operators and proved the main result with the help of the measure of noncompactness due to Kuratowski. Next, the continuous dependence of solutions was checked by utilizing the generalized Gronwall inequality. In the next step, we considered an inclusion version of the suggested boundary value problem in which we derived existence results based on approximate endpoint property for the defined multifunction. In the final step, we prepared two illustrative examples to show the validity of the results. In spite of many papers published in this field in which the authors use some standard fixed point theorems for deriving the existence results, we here utilized two theoretical techniques contained in measure theory and that of endpoint simultaneously. This defined boundary value problem is an instance of abstract boundary value problem in which one can extend boundary conditions to multipoint multiorder multistrip conditions in the future.