The Existence Results of Solutions to the Nonlinear Coupled System of Hilfer Fractional Differential Equations and Inclusions

Abstract

This paper is dedicated to studying the existence results of solutions to the nonlinear coupled system of Hilfer fractional differential equations and inclusions, with multi-strip and multi-point mixed boundary conditions. Through tools such as the Leray-Schauder alternative and the nonlinear alternative of Leray-Schauder type, continuous and measurable selection theorems, along with Leray-Schauder degree theory, the main results can be obtained. The Hilfer fractional differential system has practical implications for specific physical phenomena. Examples are provided to clarify the application of our main results.

1. Introduction

Classical calculus began with Newton’s invention of forward flow in 1665, while fractional calculus was born in 1695 when Leibniz and L’Hospital discussed the significance of $\frac{1}{2}$ derivatives, see [1]. It can be said that fractional and integer calculus have almost the same research history span. However, fractional calculus has only been revitalized in the last few decades, thanks to the development of various fields that promote research in fractional calculus. The new study of fractional calculus can also offer novel ideas for solving challenging problems in various fields, see [2,3,4,5,6,7].

Although Riemann-Liouville and Caputo fractional derivatives are considered valuable tools for modeling many real-world problems, R. Hilfer found that the traditional fractional derivatives of Riemann-Liouville and Caputo could not meet the requirements for solving new problems during the study of fractional time evolution [8]. Therefore, in order to separate fractional integrals, a generalized definition of fractional derivatives is proposed based on the Riemann-Liouville integral by R. Hilfer, which is ${\hspace{0.17em}}^H{D}^{\alpha ,\beta }u\left(t\right)=I^{\beta (n-\alpha )}{D}^n{I}^{(1-\beta )(n-\alpha )}u\left(t\right),$ where $n-1<\alpha <n,$ $0\le \beta \le 1,$ $D=\frac{d}{dt}$. Many authors later called this definition as the Hilfer fractional derivative. The reader is referred to references [9] for the distinction between the same order $\alpha$ but with different values of $\beta$.

Initial value and boundary value problems involving the Hilfer fractional derivative have attracted a lot of research. In [10], K.M. Furati discussed the existence of solutions to a Hilfer fractional differential equation for the following initial value problem. Moreover, the stability of the solution to a weighted Cauchy-type problem is also analyzed.

$$\left\{\begin{array}{c}{\hspace{0.17em}}^H{D}_{a^{+}}^{\alpha ,\beta }y\left(x\right)=f(x,y),a<x,0<\alpha <1,0\le \beta \le 1,\hfill \\ I_{a^{+}}^{1-\gamma }y\left(a^{+}\right)=y_a,\gamma =\alpha +\beta -\alpha \beta ,\hfill \end{array}\right.$$

(1)

where ${\hspace{0.17em}}^H{D}^{\alpha ,\beta }$ is the Hilfer fractional derivative of order $\alpha$, and type $\beta$.

In [11], K. Dhawan investigated the coupled Hilfer fractional differential equations with nonlocal conditions. By applying the Leray-alternative Schauder’s and the Contraction principle, the author proved the existence and uniqueness of the solution. Furthermore, the Ulam stability of the solution was discussed for the defined problem.

$$\left\{\begin{array}{c}{\hspace{0.17em}}^H{D}^{{\alpha }_1,{\beta }_1}x\left(t\right)=f(t,x\left(t\right),y\left(t\right)),{\alpha }_1<v_1={\alpha }_1+{\beta }_1-{\alpha }_1{\beta }_1,t\in [0,T],\hfill \\ {\hspace{0.17em}}^H{D}^{{\alpha }_2,{\beta }_2}y\left(t\right)=f(t,x\left(t\right),y\left(t\right)),{\alpha }_2<v_2={\alpha }_2+{\beta }_2-{\alpha }_2{\beta }_2,t\in [0,T],\hfill \\ I_{0^{+}}^{1-v_1}x\left(0\right)=\sum _{i=1}^n{b}_{1i}x\left({\xi }_i\right),{\xi }_i\in [0,T],\hfill \\ I_{0^{+}}^{1-v_2}y\left(0\right)=\sum _{i=1}^n{b}_{2i}y\left({\eta }_i\right),{\eta }_i\in [0,T],\hfill \end{array}\right.$$

(2)

where ${\hspace{0.17em}}^H{D}^{{\alpha }_i,{\beta }_i}$ represents the Hilfer fractional derivative of order ${\alpha }_i$ and type ${\beta }_i$, $i=1,2$, ${\alpha }_i\in (0,1)$, ${\beta }_i\in [0,1]$. $I_{0^{+}}^{1-v_i}$ is the left-side Riemann-Liouville integral of order $1-v_i$, $i=1,2.$ Also, $f_1,f_2:$ $[0,T]\times R\times R\to R$ are given continuous nonlinear functions and $T>0$, $b_{11},b_{12}$ are real numbers, ${\xi }_i,{\eta }_i(i=1,2,\dots ,m)$ are prefixed points satisfying $0<{\xi }_1<{\xi }_2<\dots <T$ and $0<{\eta }_1<{\eta }_2<\dots <T$, respectively.

In [12], B. Ahmad studied the nonlinear generalized coupled fractional differential equations accompanied by nonlocal coupled multipoint Riemann-Stieltjes and generalized fractional integral boundary conditions using the Leray-Schauder alternative and Banach contraction mapping principle.

$$\left\{\begin{array}{c}{\hspace{0.17em}}^{\rho }{D}_{0^{+}}^{\alpha }u\left(t\right)=f(t,u\left(t\right),v\left(t\right){,}^{\rho }{D}_{0^{+}}^{{\gamma }_1}v\left(t\right)),t\in [0,T],\hfill \\ {\hspace{0.17em}}^{\rho }{D}_{0^{+}}^{\beta }v\left(t\right)=f(t,u\left(t\right){,}^{\rho }{D}_{0^{+}}^{{\gamma }_2}u\left(t\right),v\left(t\right)),t\in [0,T],\hfill \\ u\left(0\right)=v\left(0\right)=0,\hfill \\ {\int }_0^{T}u\left(s\right)d{H}_1\left(s\right)={\lambda }_1^{\rho }{I}_{0^{+}}^{{\delta }_1}v\left({\xi }_1\right)+\sum _{p=1}^m{a}_{p}v\left({\eta }_p\right),{\xi }_{1}a\in (0,T),\hfill \\ {\int }_0^{T}v\left(s\right)d{H}_2\left(s\right)={\lambda }_2^{\rho }{I}_{0^{+}}^{{\delta }_2}v\left({\xi }_2\right)+\sum _{p=1}^m{b}_{p}u\left({\eta }_p\right),{\xi }_2\in (0,T),\hfill \end{array}\right.$$

(3)

where ${\hspace{0.17em}}^{\rho }{D}_{0^{+}}^{\alpha }$ and ${\hspace{0.17em}}^{\rho }{D}_{0^{+}}^{\beta }$ are the generalized fractional derivative operators of order $1<\alpha ,\beta <2$, respectively, $0<{\gamma }_1,{\gamma }_2<1$, ${\hspace{0.17em}}^{\rho }{I}_{0^{+}}^{{\delta }_1}$ and ${\hspace{0.17em}}^{\rho }{I}_{0^{+}}^{{\delta }_2}$ are the generalized fractional integral operators of order ${\delta }_1,{\delta }_2>0$, respectively, ${\int }_0^{T}u\left(s\right)d{H}_i\left(s\right)(i=1,2)$ are the Riemann-Stieltjes integrals with respect to the functions $H_i:[0,T]\to \mathbb{R}$, $f,g\in C([0,T]\times {\mathbb{R}}^3,\mathbb{R})$, ${\lambda }_1,{\lambda }_2,a_p,b_p\in \mathbb{R}$ and ${\eta }_p\in (0,T),p=1,2,\dots ,m.$

It is worth noting that J. Pradeesh considered the existence of a mild solution for the following Hilfer fractional stochastic differential system in [13].

$$\left\{\begin{array}{c}{\hspace{0.17em}}^H{D}_{0^{+}}^{v,\mu }z\left(t\right)\in Az\left(t\right)+f(t,z\left(t\right))+G(t,z\left(t\right))\frac{dW\left(t\right)}{dt},t\in (0,C],\hfill \\ \left(I_{0^{+}}^{2-\gamma }z\right)\left(0\right)=z_0,{\left(I_{0^{+}}^{2-\gamma }z\right)}^{\prime }\left(0\right)=z_1,\hfill \end{array}\right.$$

(4)

where ${\hspace{0.17em}}^H{D}_{0^{+}}^{v,\mu }$ denotes the Hilfer fractional derivative of order $\mu$ and type v such that $1<\mu <2$ and $0\le v\le 1$, while $I_{0^{+}}^{2-\gamma }$ is the Riemann–Liouville integral operator with order $(2-\gamma )$, where $\gamma =\mu +v(2-\mu )$. $A:D\left(A\right)\subset H\to H$ is the infinitesimal generator of a strongly continuous cosine family ${C\left(t\right)}_{t\ge 0}$ of uniformly bounded linear operators and $z(·)$ takes values in a separable Hilbert space H with inner product $\langle ·,·\rangle$ and norm $\parallel ·\parallel$. The authors employed fractional calculus, multivalued analysis, sine and cosine operators, and Bohnenblust–Karlin’s fixed point theorem to investigate the existence of a mild solution for the Hilfer fractional stochastic differential system.

However, there are very few articles that have studied Hilfer fractional differential equations and inclusions simultaneously. Inspired by the aforementioned works, we are initiating a study on the existence of solutions to a coupled system of Hilfer fractional differential equations

$$\left\{\begin{array}{c}{\hspace{0.17em}}^H{D}^{{\alpha }_1,\beta }u\left(t\right)=f_1(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right)),t\in [0,1],\hfill \\ {\hspace{0.17em}}^H{D}^{{\alpha }_2,\beta }v\left(t\right)=f_2(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right)),t\in [0,1],\hfill \end{array}\right.$$

(5)

and inclusions

$$\left\{\begin{array}{c}{\hspace{0.17em}}^H{D}^{{\alpha }_1,\beta }u\left(t\right)\in F(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right)),t\in [0,1],\hfill \\ {\hspace{0.17em}}^H{D}^{{\alpha }_2,\beta }v\left(t\right)\in G(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right)),t\in [0,1],\hfill \end{array}\right.$$

(6)

subject to the coupled fractional integral and discrete mixed boundary conditions

$$\left\{\begin{array}{c}u\left(0\right)=0,v\left(0\right)=0,\hfill \\ u\left(1\right)=\sum _{i=1}^m{\lambda }_{1i}{I}^{l_{1i}}v\left({\xi }_i\right)+\sum _{j=1}^n{b}_{1j}v\left({\eta }_j\right),\hfill \\ v\left(1\right)=\sum _{i=1}^m{\lambda }_{2i}{I}^{l_{2i}}u\left({\xi }_i\right)+\sum _{j=1}^n{b}_{2j}u\left({\eta }_j\right),\hfill \end{array}\right.$$

(7)

where ${\hspace{0.17em}}^H{D}^{{\alpha }_k,\beta }$, ${\hspace{0.17em}}^H{D}^{{\gamma }_k,\beta }$ denote the Hilfer fractional derivative of order ${\alpha }_k$, ${\gamma }_k$ respectively and parameter $\beta$, $I^{l_{ki}}$ is the Riemann-Liouville fractional integral of order $l_{ki}$, $1<{\alpha }_k\le 2$, $0\le \beta \le 1,$ $0<{\gamma }_k<{\alpha }_k-1$, ${\lambda }_{ki}\ge 0$, $l_{ki}\ge 0$, $0\le {\xi }_i\le 1$, $b_{kj}\ge 0$, $0\le {\eta }_j\le 1$, for $k=1,2$; $i=1,2,\dots ,m$; $j=1,2,\dots ,n.$

In the research of fractional time evolution, a significant challenge arises when generalizing traditional equations of motion, as it involves deciding whether to utilize the Riemann-Liouville fractional derivative or the Caputo fractional derivative. R. Hilfer introduced the Hilfer fractional derivative during the study. When $\beta =0$, the Hilfer system can be simplified to the Riemann-Liouville system studied by authors such as X. Zhao, A. Guezane-Lakoud, T. Jankowski, and others. For more details, refer to [14,15,16]. When $\beta =1$, the Hilfer system can be reduced to the Caputo system discussed by investigators like Y. Zi, A. Alsaedi, J. Xie, B. Ahmad, and others. References can be found in [17,18,19,20]. We learned that models of fractional diffusion equations with Hilfer fractional derivatives are used in the context of glass relaxation and aquifer problems [21]. Also, fractional reaction-diffusion equations and space-time fractional diffusion equations involving the Hilfer fractional derivative are studied in [22,23]. And more importantly, differential inclusion, which can handle uncertainty problems, has been applied to dynamical systems and stochastic processes, such as control problems and sweeping processes [24,25]. Therefore, our studies on the nonlinear coupled Hilfer fractional differential inclusions have great practical applications in physical phenomena.

In the coupled Hilfer fractional differential system, there are $u\left(t\right)$, $v\left(t\right)$ that mutually influence each other. The nonlinear terms, $f_1$ and $f_2$, consist of $u\left(t\right)$, $v\left(t\right)$, ${\hspace{0.17em}}^H{D}^{{\gamma }_1,\beta }u\left(t\right)$, ${\hspace{0.17em}}^H{D}^{{\gamma }_2,\beta }v\left(t\right)$, where ${\gamma }_1,{\gamma }_2$ are less than ${\alpha }_1,{\alpha }_2$. The inclusion of ${\hspace{0.17em}}^H{D}^{{\gamma }_1,\beta }u\left(t\right)$, ${\hspace{0.17em}}^H{D}^{{\gamma }_2,\beta }v\left(t\right)$ in $f_1,f_2$ enhances the model’s capability to address real-world problems The coupled Hilfer fractional differential inclusions have the same structure. Moreover, the nonlocal boundary conditions consist of the Riemann-Liouville fractional integral and numerous discrete points. The value of the unknown function $u\left(t\right)$ at the right endpoint $t=1$ is equal to the sum of the values of the Riemann–Liouville fractional integral of the unknown function $v\left(t\right)$ on the subinterval $[0,{\xi }_i](i=1,2,\dots ,m)$ and the discrete values of the unknown function $v\left(t\right)$ at ${\eta }_j(j=1,2,\dots ,n)$.

In addition, there are now many well-established methods for studying fractional differential equations and inclusions, such as Guo-Krasnoselskii’s fixed-point theorem, the Banach contraction mapping principle, and monotone iteration techniques. We employed the Leray-Schauder alternative and the nonlinear alternative of Leray-Schauder type, continuous and measurable selection theorems, along with Leray-Schauder degree theory, to explore the existence of solutions for the Hilfer fractional differential equations and inclusions, respectively.

In fact, fractional derivatives have been greatly developed and applied, leading to the emergence of several mature definitions such as the Riemann-Liouville fractional derivative, the Hadamard fractional derivative, the Caputo-Katugampola fractional derivative, the Katugampola fractional derivative, and others. In contrast, the derivative under Hilfer’s definition requires more research efforts to promote its development. Whether transferring mature research techniques to Hilfer or developing new technical methods, the work is meaningful. Differential inclusions can be regarded as a collection of differential equations and inequalities. Moreover, this paper examines the coupled system of Hilfer differential equations and inclusions, which is of great significance for practical applications.

The rest of the paper is structured as follows: In Section 2, some fundamental concepts of fractional calculus and lemmas are presented. Section 3 is dedicated to presenting the main results, which are illustrated through examples. Section 4 contains a summary of previous work and future prospects.

2. Preliminaries

In this section, we present some basic definitions, lemmas, and auxiliary results for the proof that will be utilized in the next section.

Let $C\left(J\right)$, $AC\left(J\right)$ and $C^n\left(J\right)$ denote the spaces of continuous, absolutely continuous and n times continuously differentiable functions on $J:=[0,1]$, respectively. We denote by $L^p(0,1)$, $p\ge 1,$ the spaces of Lebesgue integrable functions on $(0,1)$.

Definition 1 

([10,26]). The Riemann-Liouville fractional integral of order $\alpha >0$ of a continuous function $u\in L^1(0,1)$ is defined by

$$I^{\alpha }u\left(t\right)=\frac{1}{\Gamma \left(\alpha \right)}{\int }_a^t{(t-s)}^{\alpha -1}u\left(s\right)ds,n-1<\alpha <n,$$

where $n=\left[\alpha \right]+1$, $\left[\alpha \right]$ denotes the integer part of real number $\alpha .$

Definition 2 

([10,26]). The Riemann-Liouville fractional derivative of order $\alpha >0$ of a continuous function $u\in AC(0,1)$ is defined by

$$D^{\alpha }u\left(t\right):=D^n{I}^{n-\alpha }u\left(t\right)=\frac{1}{\Gamma (n-\alpha )}{\left(\frac{d}{dt}\right)}^n{\int }_a^t{(t-s)}^{n-\alpha -1}u\left(s\right)ds,n-1<\alpha <n,$$

where $n=\left[\alpha \right]+1$.

Definition 3 

([27]). The Caputo fractional derivative of order $\alpha >0$ of a continuous function $u\in AC(0,1)$ is defined by

$${\hspace{0.17em}}^C{D}^{\alpha }u\left(t\right):=I^{n-\alpha }{D}^{n}u\left(t\right)=\frac{1}{\Gamma (n-\alpha )}{\int }_a^t{(t-s)}^{n-\alpha -1}{\left(\frac{d}{dt}\right)}^{n}u\left(s\right)ds,n-1<\alpha <n,$$

where $n=\left[\alpha \right]+1$.

Definition 4 

([10]). The Hilfer fractional derivative of order α and parameter β of a function $u\in AC(0,1)$ is defined by

$${\hspace{0.17em}}^H{D}^{\alpha ,\beta }u\left(t\right)=I^{\beta (n-\alpha )}{D}^n{I}^{(1-\beta )(n-\alpha )}u\left(t\right),$$

where $n-1<\alpha <n,$ $0\le \beta \le 1,$ $t>a,$ $D=\frac{d}{dt}.$

Remark 1. 

From Definitions 2 and 3, when $\beta =0$, the Hilfer fractional derivative corresponds to the Riemann-Liouville fractional derivative

$${\hspace{0.17em}}^H{D}^{\alpha ,0}u\left(t\right)=D^n{I}^{n-\alpha }u\left(t\right),$$

while when $\beta =1$, the Hilfer fractional derivative corresponds to the Caputo fractional derivative

$${\hspace{0.17em}}^H{D}^{\alpha ,1}u\left(t\right)=I^{n-\alpha }{D}^{n}u\left(t\right).$$

For convenience, let $\theta =\alpha +2\beta -\alpha \beta$, $1<\alpha \le 2$, we can know $1<\theta \le 2$ and $0\le 2-\theta <1$. We have the weighted spaces of continuous functions

$$C_{2-\theta }\left(J\right)=\{u:(0,1]\to R:t^{2-\theta }u\left(t\right)\in C\left(J\right)\},$$

(8)

with the norm

$${\parallel u\parallel }_{C_{2-\theta }}=\underset{t\in J}{sup}\left|t^{2-\theta }u\left(t\right)\right|,$$

(9)

and

$$C_{2-\theta }^{\alpha ,\beta }\left(J\right)=\{u\in C_{2-\theta }\left(J\right){,}^H{D}^{\alpha ,\beta }u\in C_{2-\theta }\left(J\right)\},$$

(10)

$$C_{2-\theta }^{\theta }\left(J\right)=\{u\in C_{2-\theta }\left(J\right),D^{\theta }u\in C_{2-\theta }\left(J\right)\}.$$

(11)

Lemma 1 

([10]). For $\alpha >0,I^{\alpha }$ maps $C\left(J\right)$ into $C\left(J\right).$

Lemma 2 

([10]). Let $\alpha >0,$ and $0\le 2-\theta <1$. Then $I^{\alpha }$ is bounded from $C_{2-\theta }\left(J\right)$ into $C_{2-\theta }\left(J\right).$

Lemma 3 

([10]). Let $\alpha >0,$ and $0\le 2-\theta <1$. If $2-\theta \le \alpha ,$ then $I^{\alpha }$ is bounded from $C_{2-\theta }\left(J\right)$ into $C\left(J\right).$

Lemma 4 

([27]). If $\alpha >0,b>0$, then

$$I^{\alpha }{t}^b=\frac{\Gamma (b+1)}{\Gamma (b+1+\alpha )}{t}^{b+\alpha },$$

$$D^{\alpha }{t}^b=\frac{\Gamma (b+1)}{\Gamma (b+1-\alpha )}{t}^{b-\alpha },$$

$${\hspace{0.17em}}^H{D}^{\alpha ,\beta }{t}^b=\frac{\Gamma (b+1)}{\Gamma (b+1-\alpha )}{t}^{b-\alpha }.$$

Lemma 5 

([10]). Let $0\le 2-\theta <1$ and $u\in C_{2-\theta }[0,1]$. Then

$$I^{\alpha }u\left(0\right):=\underset{t\to 0^{+}}{lim}{I}^{\alpha }u\left(t\right)=0,0\le 2-\theta <\alpha .$$

Lemma 6 

([10]). Let $\alpha >0,$ $0\le 2-\theta <1,$ and $u\in C_{2-\theta }\left(J\right).$ Then

$$D^{\alpha }{I}^{\alpha }u\left(t\right)=u\left(t\right),forallt\in (0,1].$$

Lemma 7 

([10]). For $t>0,$ we have

$$D_{0^{+}}^{\alpha }{t}^{\alpha -1}=0,1<\alpha \le 2.$$

Lemma 8 

([27]). Let $1<\alpha \le 2,$ $0\le 2-\theta <1.$ If $u\in C_{2-\theta }\left(J\right),$ and $I^{2-\alpha }u\in C_{2-\theta }^2\left(J\right),$ then

$$I^{\alpha }{D}^{\alpha }u\left(t\right)=u\left(t\right)-\frac{\left(I^{1-\alpha }u\right)\left(t\right)}{\Gamma \left(\alpha \right)}{t}^{\alpha -1}-\frac{\left(I^{2-\alpha }u\right)\left(t\right)}{\Gamma (\alpha -1)}{t}^{\alpha -2}.$$

for all $t\in (0,1].$

Lemma 9 

([27]). Let $1<\alpha \le 2,$ and $0\le \beta \le 1.$ If $u\in C_{2-\theta }\left(J\right),$ and $I^{1-\beta (2-\alpha )}u\in C_{2-\theta }^1\left(J\right),$ then ${\hspace{0.17em}}^H{D}^{\alpha ,\beta }{I}^{\alpha }u$ exists in $(0,1]$ and

$${\hspace{0.17em}}^H{D}^{\alpha ,\beta }{I}^{\alpha }u\left(t\right)=u\left(t\right),t\in (0,1].$$

Lemma 10 

([10]). Let $1<\alpha \le 2,$ and $0\le \beta \le 1.$ If $f\in C_{2-\theta }^{\theta }\left(J\right),$ then

$$I^{\theta }{D}^{\theta }u=I^{\alpha }{\hspace{0.17em}}^H{D}^{\alpha ,\beta }u,D^{\theta }{I}^{\alpha }u=D^{\beta (2-\alpha )}u.$$

Lemma 11 

(Leray–Schauder alternative [28]). Let X be a Banach space, $T:X\to X$ be a completely continuous operator. Let

$$\xi \left(T\right)=\{x\in X:X=\lambda T(x),0<\lambda <1\}.$$

Then, either the set $\xi \left(T\right)$ is unbounded, or T has at least one fixed point.

Lemma 12 

(Nonlinear alternative of Leray-Schauder type [29]). Let X be a Banach space, Ω be a closed convex subset of X, and U be an open subset of Ω with $0\in U$. Suppose that $F:\overline{U}\to P_{cp,cv}(\Omega )$ is an upper semicontinuous compact map. Then either (1) F has a fixed point in $\overline{U}$, or (2) there is $\partial U$ and $\lambda \in (0,1)$ such that $x\in \lambda F\left(x\right)$.

For each $z\in C\left(J\right)$, we defined the set of selections of F by

$$S_{F\circ z}=\{w\in L^1\left(J\right):w\left(t\right)\in F(t,z\left(t\right)),a.e.t\in J\}.$$

For convenience, we denote

$$\left\{\begin{array}{c}{l}_1=\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\theta }_2-1}ds+\sum _{j=1}^n{b}_{1j}{\eta }_j^{{\theta }_2-1},\hfill \\ l_2=\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\theta }_1-1}ds+\sum _{j=1}^n{b}_{2j}{\eta }_j^{{\theta }_1-1}.\hfill \end{array}\right.$$

(12)

In the forthcoming analysis, we always need to make the following assumptions:

(F₁)

$1<{\alpha }_k\le 2,0<\gamma <{\alpha }_k-1,l_{ki}>0,1<{\theta }_k\le 2,$ where ${\theta }_k={\alpha }_k+2\beta -{\alpha }_k\beta ,$ for $k=1,2$ and $i=1,2,\dots ,m$;

(F₂)

$0\le {\xi }_i,{\eta }_j\le 1,{\lambda }_{1i},{\lambda }_{2i}\ge 0,b_{1j},b_{2j}\ge 0,$ for $i=1,2,\dots ,m$, $j=1,2,\dots ,n;$

(F₃)

$1-l_1{l}_2>0,$ where $l_1,l_2$ are defined by (12).

Subject to BVP (5) and (7), we consider a corresponding linear differential system as follows and establish the expression of the corresponding Green’s functions.

Lemma 13. 

Assume that $\left({\mathrm{F}}_1\right)$–$\left({\mathrm{F}}_3\right)$ hold. For $h_1\in C_{2-{\theta }_1}\left(J\right),h_2\in C_{2-{\theta }_2}\left(J\right)$, the fractional differential system

$$\left\{\begin{array}{c}{\hspace{0.17em}}^H{D}^{{\alpha }_1,\beta }u\left(t\right)=h_1\left(t\right),\hfill \\ {\hspace{0.17em}}^H{D}^{{\alpha }_2,\beta }v\left(t\right)=h_2\left(t\right),\hfill \end{array}\right.$$

(13)

with boundary conditions (5) has an integral representation

$$\begin{array}{cc}\hfill & u\left(t\right)=I^{{\alpha }_1}{h}_1\left(t\right)+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ \hfill & [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{h}_1\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{h}_1\left(s\right)ds-I^{{\alpha }_2}{h}_2\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{h}_2\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{h}_2\left(s\right)ds-I^{{\alpha }_1}{h}_1\left(1\right)],\hfill \end{array}$$

(14)

$$\begin{array}{cc}\hfill & v\left(t\right)=I^{{\alpha }_2}{h}_2\left(t\right)+\frac{t^{{\theta }_2-1}}{(1-l_1{l}_2)}\hfill \\ \hfill & [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{h}_2\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{h}_2\left(s\right)ds-I^{{\alpha }_1}{h}_1\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{h}_1\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{h}_1\left(s\right)ds-I^{{\alpha }_2}{h}_2\left(1\right)].\hfill \end{array}$$

(15)

Proof. 

From Definition 1 and Lemma 10, using $I^{{\alpha }_1},I^{{\alpha }_2}$ to (13), we have:

$$\left\{\begin{array}{c}{I}^{{\theta }_1}{D}^{{\theta }_1}u\left(t\right)=I^{{\alpha }_1}{h}_1\left(t\right),{\theta }_1={\alpha }_1+\beta (2-{\alpha }_1),\hfill \\ I^{{\theta }_2}{D}^{{\theta }_2}v\left(t\right)=I^{{\alpha }_2}{h}_2\left(t\right),{\theta }_2={\alpha }_2+\beta (2-{\alpha }_2).\hfill \end{array}\right.$$

(16)

From Lemma 8, we obtain

$$\left\{\begin{array}{c}{I}^{{\theta }_1}{D}^{{\theta }_1}u\left(t\right)=u\left(t\right)-\frac{\left(I^{1-{\theta }_1}u\right)\left(t\right)}{\Gamma \left({\theta }_1\right)}{t}^{{\theta }_1-1}-\frac{\left(I^{2-{\theta }_1}u\right)\left(t\right)}{\Gamma ({\theta }_1-1)}{t}^{{\theta }_1-2},\hfill \\ I^{{\theta }_2}{D}^{{\theta }_2}v\left(t\right)=v\left(t\right)-\frac{\left(I^{1-{\theta }_2}u\right)\left(t\right)}{\Gamma \left({\theta }_2\right)}{t}^{{\theta }_2-1}-\frac{\left(I^{2-{\theta }_2}v\right)\left(t\right)}{\Gamma ({\theta }_2-1)}{t}^{{\theta }_2-2}.\hfill \end{array}\right.$$

(17)

Combining (16) and (17), we can simplify (13) to the following equivalent integral equations

$$\left\{\begin{array}{c}u\left(t\right)=I^{{\alpha }_1}{h}_1\left(t\right)+\frac{c_{11}}{\Gamma \left({\theta }_1\right)}{t}^{{\theta }_1-1}+\frac{c_{12}}{\Gamma ({\theta }_1-1)}{t}^{{\theta }_1-2},\hfill \\ v\left(t\right)=I^{{\alpha }_2}{h}_2\left(t\right)+\frac{c_{21}}{\Gamma \left({\theta }_2\right)}{t}^{{\theta }_2-1}+\frac{c_{22}}{\Gamma ({\theta }_2-1)}{t}^{{\theta }_2-2},\hfill \end{array}\right.$$

(18)

where $c_{11},c_{12},c_{21},c_{22}$ are constants.

From $u\left(0\right)=v\left(0\right)=0$, we obtain $c_{12}=c_{22}=0$, and we get

$$\left\{\begin{array}{c}u\left(t\right)=I^{{\alpha }_1}{h}_1\left(t\right)+\frac{c_{11}}{\Gamma \left({\theta }_1\right)}{t}^{{\theta }_1-1},\hfill \\ v\left(t\right)=I^{{\alpha }_2}{h}_2\left(t\right)+\frac{c_{21}}{\Gamma \left({\theta }_2\right)}{t}^{{\theta }_2-1}.\hfill \end{array}\right.$$

(19)

From the remaining conditions in (5), it can be inferred that

$$\left\{\begin{array}{c}u\left(1\right)=I^{{\alpha }_1}{h}_1\left(1\right)+\frac{c_{11}}{\Gamma \left({\theta }_1\right)}=\sum _{i=1}^m{\lambda }_{1i}{I}^{l_{1i}}v\left({\xi }_i\right)+\sum _{j=1}^n{b}_{1j}v\left({\eta }_j\right),\hfill \\ v\left(1\right)=I^{{\alpha }_2}{h}_2\left(1\right)+\frac{c_{21}}{\Gamma \left({\theta }_2\right)}=\sum _{i=1}^m{\lambda }_{2i}{I}^{l_{2i}}u\left({\xi }_i\right)+\sum _{j=1}^n{b}_{2j}u\left({\eta }_j\right).\hfill \end{array}\right.$$

(20)

$$\left\{\begin{array}{c}{c}_{11}=\Gamma \left({\theta }_1\right)\left(\sum _{i=1}^m{\lambda }_{1i}{I}^{l_{1i}}v\left({\xi }_i\right)+\sum _{j=1}^n{b}_{1j}v\left({\eta }_j\right)-I^{{\alpha }_1}{h}_1\left(1\right)\right),\hfill \\ c_{21}=\Gamma \left({\theta }_2\right)\left(\sum _{i=1}^m{\lambda }_{2i}{I}^{l_{2i}}u\left({\xi }_i\right)+\sum _{j=1}^n{b}_{2j}u\left({\eta }_j\right)-I^{{\alpha }_2}{h}_2\left(1\right)\right).\hfill \end{array}\right.$$

(21)

Further, we can reduce (18) to

$$\left\{\begin{array}{c}u\left(t\right)=I^{{\alpha }_1}{h}_1\left(t\right)+t^{{\theta }_1-1}\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}v\left(s\right)ds+\sum _{j=1}^n{b}_{1j}v\left({\eta }_j\right)-I^{{\alpha }_1}{h}_1\left(1\right)\right),\hfill \\ v\left(t\right)=I^{{\alpha }_2}{h}_2\left(t\right)+t^{{\theta }_2-1}\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}u\left(s\right)ds+\sum _{j=1}^n{b}_{2j}u\left({\eta }_j\right)-I^{{\alpha }_2}{h}_2\left(1\right)\right).\hfill \end{array}\right.$$

(22)

Then we can get

$$\begin{array}{cc}\hfill & \sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}v\left(s\right)ds+\sum _{j=1}^n{b}_{1j}v\left({\eta }_j\right)-I^{{\alpha }_1}{h}_1\left(1\right)\hfill \\ \hfill & =\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\theta }_2-1}ds+\sum _{j=1}^n{b}_{1j}{\eta }_j^{{\theta }_2-1}\right)\hfill \\ \hfill & ·\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}u\left(s\right)ds+\sum _{j=1}^n{b}_{2j}u\left({\eta }_j\right)-I^{{\alpha }_2}{h}_2\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{h}_2\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{h}_2\left(s\right)ds-I^{{\alpha }_1}{h}_1\left(1\right),\hfill \end{array}$$

(23)

$$\begin{array}{cc}\hfill & \sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}u\left(s\right)ds+\sum _{j=1}^n{b}_{2j}u\left({\eta }_j\right)-I^{{\alpha }_2}{h}_2\left(1\right)\hfill \\ \hfill & =\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\theta }_1-1}ds+\sum _{j=1}^n{b}_{2j}{\eta }_j^{{\theta }_1-1}\right)\hfill \\ \hfill & ·\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}v\left(s\right)ds+\sum _{j=1}^n{b}_{1j}v\left({\eta }_j\right)-I^{{\alpha }_1}{h}_1\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{h}_1\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{h}_1\left(s\right)ds-I^{{\alpha }_2}{h}_2\left(1\right).\hfill \end{array}$$

(24)

Combining (23) and (24), it can be seen that

$$\begin{array}{cc}\hfill & \sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}v\left(s\right)ds+\sum _{j=1}^n{b}_{1j}v\left({\eta }_j\right)-I^{{\alpha }_1}{h}_1\left(1\right)\hfill \\ \hfill & =\frac{1}{1-l_1{l}_2}[l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{h}_1\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{h}_1\left(s\right)ds-I^{{\alpha }_2}{h}_2\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{h}_2\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{h}_2\left(s\right)ds-I^{{\alpha }_1}{h}_1\left(1\right)],\hfill \\ \hfill & \sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}u\left(s\right)ds+\sum _{j=1}^n{b}_{1j}u\left({\eta }_j\right)-I^{{\alpha }_2}{h}_2\left(1\right)\hfill \\ \hfill & =\frac{1}{1-l_1{l}_2}[l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{h}_2\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{h}_2\left(s\right)ds-I^{{\alpha }_1}{h}_1\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{h}_1\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{h}_1\left(s\right)ds-I^{{\alpha }_2}{h}_2\left(1\right)],\hfill \end{array}$$

(25)

where $l_k$ $(k=1,2)$ is defined by (12). From (22) and (25), we have

$$\begin{array}{cc}\hfill & u\left(t\right)=I^{{\alpha }_1}{h}_1\left(t\right)+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ \hfill & [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{h}_1\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{h}_1\left(s\right)ds-I^{{\alpha }_2}{h}_2\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{h}_2\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{h}_2\left(s\right)ds-I^{{\alpha }_1}{h}_1\left(1\right)],\hfill \end{array}$$

(26)

$$\begin{array}{cc}\hfill & v\left(t\right)=I^{{\alpha }_2}{h}_2\left(t\right)+\frac{t^{{\theta }_2-1}}{(1-l_1{l}_2)}\hfill \\ \hfill & [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{h}_2\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{h}_2\left(s\right)ds-I^{{\alpha }_1}{h}_1\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{h}_1\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{h}_1\left(s\right)ds-I^{{\alpha }_2}{h}_2\left(1\right)].\hfill \end{array}$$

(27)

This completes the proof of the lemma. □

Moreover, according to Lemma 4, and fractional order derivative of solution (11) can be expressed as

$$\begin{array}{cc}\hfill & {\hspace{0.17em}}^H{D}^{{\gamma }_1,\beta }u\left(t\right)=I^{{\alpha }_1-{\gamma }_1}{h}_1\left(t\right)+\frac{\Gamma \left({\theta }_1\right)t^{{\theta }_1-{\gamma }_1-1}}{\Gamma ({\theta }_1-{\gamma }_1)(1-l_1{l}_2)}\hfill \\ \hfill & [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{h}_1\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{h}_1\left(s\right)ds-I^{{\alpha }_2}{h}_2\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{h}_2\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{h}_2\left(s\right)ds-I^{{\alpha }_1}{h}_1\left(1\right)],\hfill \end{array}$$

(28)

$$\begin{array}{cc}\hfill & {\hspace{0.17em}}^H{D}^{{\gamma }_2,\beta }v\left(t\right)=I^{{\alpha }_2-{\gamma }_2}{h}_2\left(t\right)+\frac{\Gamma \left({\theta }_2\right)t^{{\theta }_2-{\gamma }_2-1}}{\Gamma ({\theta }_2-{\gamma }_2)(1-l_1{l}_2)}\hfill \\ \hfill & [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{h}_2\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{h}_2\left(s\right)ds-I^{{\alpha }_1}{h}_1\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{h}_1\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{h}_1\left(s\right)ds-I^{{\alpha }_2}{h}_2\left(1\right)].\hfill \end{array}$$

(29)

Lemma 14. 

Let $f_k:J\times R^4\to R$ be a function such that $f_k\in C_{2-{\theta }_k}\left(J\right),k=1,2$. If $u\in C_{2-{\theta }_1}^{{\theta }_1}\left(J\right),$ $v\in C_{2-{\theta }_2}^{{\theta }_2}\left(J\right),$ then $u,v$ satisfie (5) and (7) if and only if $u,v$ satisfies (30) and (31).

$$\begin{array}{cc}\hfill & u\left(t\right)=I^{{\alpha }_1}{f}_1(u,v)\left(t\right)+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ \hfill & [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_1(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_1(u,v)\left(s\right)ds-I^{{\alpha }_2}{f}_2(u,v)\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{f}_2(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{f}_2(u,v)\left(s\right)ds-I^{{\alpha }_1}{f}_1(u,v)\left(1\right)],\hfill \end{array}$$

(30)

$$\begin{array}{cc}\hfill & v\left(t\right)=I^{{\alpha }_2}{f}_2(u,v)\left(t\right)+\frac{t^{{\theta }_2-1}}{(1-l_1{l}_2)}\hfill \\ \hfill & [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{f}_2(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{f}_2(u,v)\left(s\right)ds-I^{{\alpha }_1}{f}_1(u,v)\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_1(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_1(u,v)\left(s\right)ds-I^{{\alpha }_2}{f}_2(u,v)\left(1\right)],\hfill \end{array}$$

(31)

where

$$\left\{\begin{array}{c}{f}_1(u,v)\left(s\right)\triangleq f_1(s,u\left(s\right),v\left(s\right){,}^H{D}^{{\gamma }_1,\beta }u\left(s\right){,}^H{D}^{{\gamma }_2,\beta }v\left(s\right)),\hfill \\ f_2(u,v)\left(s\right)\triangleq f_2(s,u\left(s\right),v\left(s\right){,}^H{D}^{{\gamma }_1,\beta }u\left(s\right){,}^H{D}^{{\gamma }_2,\beta }v\left(s\right)).\hfill \end{array}\right.$$

(32)

Proof. 

First we prove the necessity. Let $u\in C_{2-{\theta }_1}^{{\theta }_1}\left(J\right),$ $v\in C_{2-{\theta }_2}^{{\theta }_2}\left(J\right)$ be a solution of (5) and (7). We want to prove that $u,v$ are the solutions of the integral Equations (30) and (31). By the definition of $C_{2-\theta }^{\theta }\left(J\right),$ Lemma 3, we have

$$I^{2-{\theta }_1}u\in C\left(J\right),I^{2-{\theta }_2}v\in C\left(J\right),$$

(33)

and

$$\begin{array}{c}\hfill D^{{\theta }_1}u=D^2\left(I^{2-{\theta }_1}u\right)\in C_{2-{\theta }_1}\left(J\right),D^{{\theta }_2}v=D^2\left(I^{2-{\theta }_2}v\right)\in C_{2-{\theta }_2}\left(J\right).\end{array}$$

(34)

Thus by (11), we have

$$I^{2-{\theta }_1}u\in C_{2-{\theta }_1}^2\left(J\right),I^{2-{\theta }_2}v\in C_{2-{\theta }_2}^2\left(J\right).$$

(35)

We apply Lemma 8 to obtain

$$\begin{array}{cc}\hfill & I^{{\theta }_1}{D}^{{\theta }_1}u\left(t\right)=u\left(t\right)-\frac{\left(I^{1-{\theta }_1}u\right)\left(t\right)}{\Gamma \left({\theta }_1\right)}{t}^{{\theta }_1-1}-\frac{\left(I^{2-{\theta }_2}u\right)\left(t\right)}{\Gamma ({\theta }_1-1)}{t}^{{\theta }_1-2},\hfill \\ \hfill & I^{{\theta }_2}{D}^{{\theta }_2}v\left(t\right)=v\left(t\right)-\frac{\left(I^{1-{\theta }_2}v\right)\left(t\right)}{\Gamma \left({\theta }_2\right)}{t}^{{\theta }_2-1}-\frac{\left(I^{2-{\theta }_2}v\right)\left(t\right)}{\Gamma ({\theta }_2-1)}{t}^{{\theta }_2-2}.\hfill \end{array}$$

(36)

Since by $D^{{\theta }_1}u\in C_{2-{\theta }_1},D^{{\theta }_2}v\in C_{2-{\theta }_2},$ Lemma 10, we have

$$I^{{\theta }_1}{D}^{{\theta }_1}u=I^{{\alpha }_1}{\hspace{0.17em}}^H{D}^{{\alpha }_1,\beta }u=I^{{\alpha }_1}{f}_1,I^{{\theta }_2}{D}^{{\theta }_2}v=I^{{\alpha }_2}{\hspace{0.17em}}^H{D}^{{\alpha }_1,\beta }v=I^{{\alpha }_2}{f}_2.$$

(37)

Combining (36) and (37), we have

$$\left\{\begin{array}{c}u\left(t\right)=I^{{\alpha }_1}{f}_1(u,v)\left(t\right)+\frac{d_{11}}{\Gamma \left({\theta }_1\right)}{t}^{{\theta }_1-1}+\frac{d_{12}}{\Gamma ({\theta }_1-1)}{t}^{{\theta }_1-2},\hfill \\ v\left(t\right)=I^{{\alpha }_2}{f}_2(u,v)\left(t\right)+\frac{d_{21}}{\Gamma \left({\theta }_2\right)}{t}^{{\theta }_2-1}+\frac{d_{22}}{\Gamma ({\theta }_2-1)}{t}^{{\theta }_2-2},\hfill \end{array}\right.$$

(38)

where $d_{11},d_{12},d_{21},d_{22}$ are constants. According to boundary conditions (5) and calculation steps in Lemma 13, we obtain

$$\begin{array}{cc}\hfill & u\left(t\right)=I^{{\alpha }_1}{f}_1(u,v)\left(s\right)+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ \hfill & [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_1(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_1(u,v)\left(s\right)ds-I^{{\alpha }_2}{f}_2(u,v)\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{f}_2(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{f}_2(u,v)\left(s\right)ds-I^{{\alpha }_1}{f}_1(u,v)\left(1\right)],\hfill \end{array}$$

(39)

$$\begin{array}{cc}\hfill & v\left(t\right)=I^{{\alpha }_2}{f}_2(u,v)\left(s\right)+\frac{t^{{\theta }_2-1}}{(1-l_1{l}_2)}\hfill \\ \hfill & [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{f}_2(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{f}_2(u,v)\left(s\right)ds-I^{{\alpha }_1}{f}_1(u,v)\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_1(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_1(u,v)\left(s\right)ds-I^{{\alpha }_2}{f}_2(u,v)\left(1\right)],\hfill \end{array}$$

(40)

which is the integral Equations (30) and (31).

Now we prove the sufficiency. Let $u\in C_{2-{\theta }_1}^{{\theta }_1}\left(J\right),$ $v\in C_{2-{\theta }_2}^{{\theta }_2}\left(J\right)$ satisfy (30) and (31). Since ${\hspace{0.17em}}^H{D}^{{\alpha }_1,\beta }u=I^{\beta (2-{\alpha }_1)}{D}^{{\theta }_1}u,$ ${\hspace{0.17em}}^H{D}^{{\alpha }_2,\beta }v=I^{\beta (2-{\alpha }_2)}{D}^{\theta }v$, it follows from Lemma 2 that $C_{2-{\theta }_1}^{{\theta }_1}\left(J\right)\subset C_{2-{\theta }_1}^{{\alpha }_1,\beta }\left(J\right)$ and $C_{2-{\theta }_1}^{{\theta }_1}\left(J\right)\subset C_{2-{\theta }_1}^{{\alpha }_1,\beta }\left(J\right)$. We can get ${\hspace{0.17em}}^H{D}^{{\alpha }_1,\beta }u$ and ${\hspace{0.17em}}^H{D}^{{\alpha }_2,\beta }v$ exist. Applying the operator $D^{{\theta }_1},D^{{\theta }_2}$ to (30) and (31), it follows from Lemma 7, Lemma 10 that

$$D^{{\theta }_1}u=D^{\beta (2-{\alpha }_1)}{f}_1,D^{{\theta }_2}v=D^{\beta (2-{\alpha }_2)}{f}_2.$$

(41)

From (41) and $D^{{\theta }_1}u\in C_{2-{\theta }_1}\left(J\right),D^{{\theta }_2}v\in C_{2-{\theta }_2}\left(J\right),$ we have

$$D{I}^{1-\beta (2-{\alpha }_1)}{f}_1=D^{\beta (2-{\alpha }_1)}{f}_1\in C_{2-{\theta }_1}\left(J\right),$$

(42)

$$D{I}^{1-\beta (2-{\alpha }_2)}{f}_2=D^{\beta (2-{\alpha }_2)}{f}_2\in C_{2-{\theta }_2}\left(J\right).$$

(43)

Since $f_k\in C_{2-{\theta }_k},k=1,2,$ by Lemma 2, we have

$$I^{1-\beta (2-{\alpha }_1)}{f}_1\in C_{2-{\theta }_1}\left(J\right),I^{1-\beta (2-{\alpha }_2)}{f}_2\in C_{2-{\theta }_2}\left(J\right).$$

(44)

It follows from (11), (42)–(44) that

$$I^{1-\beta (2-{\alpha }_1)}{f}_1\in C_{2-{\theta }_1}^1\left(J\right),I^{1-\beta (2-{\alpha }_2)}{f}_2\in C_{2-{\theta }_2}^1\left(J\right).$$

(45)

Now by applying $I^{\beta (2-{\alpha }_1)},I^{\beta (2-{\alpha }_2)}$ to (41), using Lemma 8, we have

$${\hspace{0.17em}}^H{D}^{{\alpha }_1,\beta }u\left(t\right)=f_1+\frac{\left(I^{1-\beta (2-{\alpha }_1)}{f}_1\right)\left(0\right)}{\Gamma \left(\beta (2-{\alpha }_1)\right)}{t}^{\beta (2-{\alpha }_1)-1},$$

(46)

$${\hspace{0.17em}}^H{D}^{{\alpha }_2,\beta }v\left(t\right)=f_2+\frac{\left(I^{1-\beta (2-{\alpha }_2)}{f}_2\right)\left(0\right)}{\Gamma \left(\beta (2-{\alpha }_2)\right)}{t}^{\beta (2-{\alpha }_2)-1}.$$

(47)

According to ${\theta }_k={\alpha }_k+2\beta -{\alpha }_k\beta ,$ $1<{\alpha }_k\le 2$ and $0\le \beta \le 1$, we have $2-{\theta }_k<1-\beta (2-{\alpha }_k)$ for $k=1,2$. Lemma 5 implies that

$$\left(I^{1-\beta (2-{\alpha }_1)}{f}_1\right)\left(0\right)=0,\left(I^{1-\beta (2-{\alpha }_2)}{f}_2\right)\left(0\right)=0.$$

(48)

Hence the (46) and (47) reduce to

$${\hspace{0.17em}}^H{D}^{{\alpha }_1,\beta }u\left(t\right)=f_1,{\hspace{0.17em}}^H{D}^{{\alpha }_2,\beta }v\left(t\right)=f_2.$$

(49)

Now we show that the boundary conditions (5) also hold. From Lemma 5, (30) and (31), we have

$$u\left(0\right)=I^{{\alpha }_1}{f}_1\left(0\right)=0,v\left(0\right)=I^{{\alpha }_2}{f}_2\left(0\right)=0.$$

(50)

From (25), (30) and (31), we obtain

$$\begin{array}{cc}\hfill & u\left(1\right)=I^{{\alpha }_1}{f}_1(u,v)\left(1\right)+\frac{1}{(1-l_1{l}_2)}\hfill \\ \hfill & [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_1(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_1(u,v)\left(s\right)ds-I^{{\alpha }_2}{f}_2(u,v)\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{f}_2(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{f}_2(u,v)\left(s\right)ds-I^{{\alpha }_1}{f}_1(u,v)\left(1\right)]\hfill \\ \hfill & =I^{{\alpha }_1}{f}_1(u,v)\left(1\right)+\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}v\left(s\right)ds+\sum _{j=1}^n{b}_{1j}v\left({\eta }_j\right)-I^{{\alpha }_1}{f}_1(u,v)\left(1\right)\right)\hfill \\ \hfill & =\sum _{i=1}^m{\lambda }_{1i}{I}^{l_{1i}}v\left({\xi }_i\right)+\sum _{j=1}^n{b}_{1j}v\left({\eta }_j\right),\hfill \end{array}$$

(51)

$$\begin{array}{cc}\hfill & v\left(1\right)=I^{{\alpha }_2}{f}_2(u,v)\left(1\right)+\frac{1}{(1-l_1{l}_2)}\hfill \\ \hfill & [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{f}_2(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{f}_2(u,v)\left(s\right)ds-I^{{\alpha }_1}{f}_1(u,v)\left(1\right)\right)\hfill \\ \hfill & +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_1(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_1(u,v)\left(s\right)ds-I^{{\alpha }_2}{f}_2(u,v)\left(1\right)]\hfill \\ \hfill & =I^{{\alpha }_2}{f}_2(u,v)\left(1\right)+\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}u\left(s\right)ds+\sum _{j=1}^n{b}_{1j}u\left({\eta }_j\right)-I^{{\alpha }_2}{f}_2(u,v)\left(1\right)\right)\hfill \\ \hfill & =\sum _{i=1}^m{\lambda }_{2i}{I}^{l_{2i}}u\left({\xi }_i\right)+\sum _{j=1}^n{b}_{2j}u\left({\eta }_j\right).\hfill \end{array}$$

(52)

This completes the proof. □

3. Main Results

For computational convenience, we introduce the notations:

$$M_1=\frac{2-l_1{l}_2}{(1-l_1{l}_2)\Gamma ({\alpha }_1+1)}+\frac{l_1}{(1-l_1{l}_2)}\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)\Gamma ({\alpha }_1+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\alpha }_1}ds+\sum _{j=1}^n\frac{b_{2j}{\eta }_j^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}\right),$$

(53)

$$M_2=\frac{1}{(1-l_1{l}_2)}\left(\frac{l_1}{\Gamma ({\alpha }_2+1)}+\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)\Gamma ({\alpha }_2+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\alpha }_2}ds+\sum _{j=1}^n\frac{b_{1j}{\eta }_j^{{\alpha }_2}}{\Gamma ({\alpha }_2+1)}\right),$$

(54)

$$M_3=\frac{1}{\Gamma ({\alpha }_1-{\gamma }_1+1)}+\frac{\Gamma \left({\theta }_1\right)}{\Gamma ({\theta }_1-{\gamma }_1)(1-l_1{l}_2)}\left(\sum _{i=1}^m\frac{l_1{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)\Gamma ({\alpha }_1+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\alpha }_1}ds+\sum _{j=1}^n\frac{l_1{b}_{2j}{\eta }_j^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}+\frac{1}{\Gamma ({\alpha }_1+1)}\right),$$

(55)

$$M_4=\frac{\Gamma \left({\theta }_1\right)}{\Gamma ({\theta }_1-{\gamma }_1)(1-l_1{l}_2)}\left(\frac{l_1}{\Gamma ({\alpha }_2+1)}+\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)\Gamma ({\alpha }_2+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\alpha }_2}ds+\sum _{j=1}^n\frac{b_{1j}{\eta }_j^{{\alpha }_2}}{\Gamma ({\alpha }_2+1)}\right),$$

(56)

$$M_5=\frac{1}{(1-l_1{l}_2)}\left(\frac{l_2}{\Gamma ({\alpha }_1+1)}+\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)\Gamma ({\alpha }_1+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\alpha }_1}ds+\sum _{j=1}^n\frac{b_{2j}{\eta }_j^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}\right),$$

(57)

$$M_6=\frac{2-l_1{l}_2}{(1-l_1{l}_2)\Gamma ({\alpha }_2+1)}+\frac{l_2}{(1-l_1{l}_2)}\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)\Gamma ({\alpha }_2+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\alpha }_2}ds+\sum _{j=1}^n\frac{b_{1j}{\eta }_j^{{\alpha }_2}}{\Gamma ({\alpha }_2+1)}\right),$$

(58)

$$M_7=\frac{\Gamma \left({\theta }_2\right)}{\Gamma ({\theta }_2-{\gamma }_2)(1-l_1{l}_2)}\left(\frac{l_2}{\Gamma ({\alpha }_1+1)}+\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)\Gamma ({\alpha }_1+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\alpha }_1}ds+\sum _{j=1}^n\frac{b_{2j}{\eta }_j^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}\right),$$

(59)

$$M_8=\frac{1}{\Gamma ({\alpha }_2-{\gamma }_2+1)}+\frac{\Gamma \left({\theta }_2\right)}{\Gamma ({\theta }_2-{\gamma }_2)(1-l_1{l}_2)}\left(\sum _{i=1}^m\frac{{\lambda }_{1i}{l}_2}{\Gamma \left(l_{1i}\right)\Gamma ({\alpha }_2+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\alpha }_2}ds+\sum _{j=1}^n\frac{l_2{b}_{1j}{\eta }_j^{{\alpha }_2}}{\Gamma ({\alpha }_2+1)}+\frac{1}{\Gamma ({\alpha }_2+1)}\right).$$

(60)

Also, Let $C_1=\left\{u|u\in C_{2-{\theta }_1}{,}^H{D}^{{\gamma }_1,\beta }u\left(t\right)\in C_{2-{\theta }_1}\right\}$ be a Banach space endowed with the norm

$${\parallel u\parallel }_{C_1}={\parallel u\parallel }_{C_{2-{\theta }_1}}+{{\parallel }^H{D}^{{\gamma }_1,\beta }u\parallel }_{C_{2-{\theta }_1}},$$

(61)

which defined by (9) and $C_2=\left\{v|v\in C_{2-{\theta }_2}{,}^H{D}^{{\gamma }_2,\beta }v\left(t\right)\in C_{2-{\theta }_2}\right\}$ be a Banach space endowed with the norm

$${\parallel v\parallel }_{C_2}={\parallel v\parallel }_{C_{2-{\theta }_2}}+{{\parallel }^H{D}^{{\gamma }_2,\beta }v\parallel }_{C_{2-{\theta }_2}},$$

(62)

which defined by (9). We can have the poiduct space $\mathcal{C}:=C_1\times C_2$ with the norm

$${\parallel (u,v)\parallel }_{\mathcal{C}}={\parallel u\parallel }_{C_1}+{\parallel v\parallel }_{C_2}.$$

(63)

3.1. The Existence Results of the Coupled Hilfer Fractional Differential System

Theorem 1. 

Let $f_k:J\times R^4\to R$ be a function such that $f_k\in C_{2-{\theta }_k}\left(J\right),k=1,2$. Assume that there exist real constants $a_i,b_i\ge 0$ for $i=1,2,3,4$ and $a_i,b_i>0$ such that, for any $x_i\in R(i=1,2,3,4),$ we have

$$\begin{array}{c}\left|f_1(t,x_1,x_2,x_3,x_4)\right|\le a_0+a_1|x_1|+a_2|x_2|+a_3|x_3|+a_4\left|x_4\right|,\hfill \\ \left|f_2(t,x_1,x_2,x_3,x_4)\right|\le b_0+b_1|x_1|+b_2|x_2|+b_3|x_3|+b_4\left|x_4\right|.\hfill \end{array}$$

(64)

If $(a_1+a_3)(M_1+M_3+M_5+M_7)+(b_1+b_3)(M_2+M_4+M_6+M_8)<1$ and $(a_2+a_4)\hspace{0.17em}\hspace{0.17em}(M_1+M_3+M_5+M_7)+(b_2+b_4)(M_2+M_4+M_6+M_8)<1$, where $M_i(i=1,2,3,4,5,6,7,8)$ are defined by (53)–(60), then (5) and (7) has at least one solution on J.

Proof. 

We define the operator $T:\mathcal{C}\to \mathcal{C}$ by

$$\begin{array}{c}\hfill T(u,v)\left(t\right)=\left(\begin{array}{c}{T}_1(u,v)\left(t\right)\\ T_2(u,v)\left(t\right)\end{array}\right),\end{array}$$

(65)

where

$$\begin{array}{c}{T}_1(u,v)\left(t\right)=I^{{\alpha }_1}{f}_1(u,v)\left(t\right)+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_1(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_1(u,v)\left(s\right)ds-I^{{\alpha }_2}{f}_2(u,v)\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{f}_2(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{f}_2(u,v)\left(s\right)ds-I^{{\alpha }_1}{f}_1(u,v)\left(1\right)],\hfill \end{array}$$

(66)

$$\begin{array}{c}{T}_2(u,v)\left(t\right)=I^{{\alpha }_2}{f}_2(u,v)\left(t\right)+\frac{t^{{\theta }_2-1}}{(1-l_1{l}_2)}\hfill \\ [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{f}_2(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{f}_2(u,v)\left(s\right)ds-I^{{\alpha }_1}{f}_1(u,v)\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_1(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_1(u,v)\left(s\right)ds-I^{{\alpha }_2}{f}_2(u,v)\left(1\right)].\hfill \end{array}$$

(67)

Also, according to (28) and (29), it is easy to see that

$$\begin{array}{c}{\hspace{0.17em}}^H{D}^{{\gamma }_1,\beta }{T}_1(u,v)\left(t\right)=I^{{\alpha }_1-{\gamma }_1}{f}_1(u,v)\left(t\right)+\frac{\Gamma \left({\theta }_1\right)t^{{\theta }_1-{\gamma }_1-1}}{\Gamma ({\theta }_1-{\gamma }_1)(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_1(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_1(u,v)\left(s\right)ds-I^{{\alpha }_2}{f}_2(u,v)\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{f}_2(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{f}_2(u,v)\left(s\right)ds-I^{{\alpha }_1}{f}_1(u,v)\left(1\right)],\hfill \end{array}$$

(68)

$$\begin{array}{c}{\hspace{0.17em}}^H{D}^{{\gamma }_2,\beta }{T}_2(u,v)\left(t\right)=I^{{\alpha }_2-{\gamma }_2}{f}_2(u,v)\left(t\right)+\frac{\Gamma \left({\theta }_2\right)t^{{\theta }_2-{\gamma }_2-1}}{\Gamma ({\theta }_2-{\gamma }_2)(1-l_1{l}_2)}\hfill \\ [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{f}_2(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{f}_2(u,v)\left(s\right)ds-I^{{\alpha }_1}{f}_1(u,v)\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_1(u,v)\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_1(u,v)\left(s\right)ds-I^{{\alpha }_2}{f}_2(u,v)\left(1\right)].\hfill \end{array}$$

(69)

We will show that the operator $T:\mathcal{C}\to \mathcal{C}$ is completely continuous.

(a) The operator $T:\mathcal{C}\to \mathcal{C}$ is continuous.

By the continuity of the functions $f_1$ and $f_2$, the operator T is continuous.

(b) The operator $T:\mathcal{C}\to \mathcal{C}$ is uniformly bounded. Let $\Omega$ be any bounded subset of $\mathcal{C}$. There exist positive constants $L_1,L_2$ such that

$$|f_1(u,v)\left(t\right)|\le L_1,\left|f_2(u,v)\left(t\right)\right|\le L_2.$$

(70)

For any $(u,v)\in \Omega$, we get

$$\begin{array}{c}|t^{2-{\theta }_1}{T}_1(u,v))\left(t\right)|\le I^{{\alpha }_1}|f_1(u,v)|\left(t\right)+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}|f_1(u,v)|\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}|f_1(u,v)\left(s\right)|ds+I^{{\alpha }_2}|f_2(u,v)|\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}|f_2(u,v)|\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}|f_2(u,v)|\left(s\right)ds+I^{{\alpha }_1}|f_1(u,v)|\left(1\right)]\hfill \\ \le \frac{t^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}{L}_1+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}{L}_1}{\Gamma \left(l_{2i}\right)\Gamma ({\alpha }_1+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\alpha }_1}ds+\sum _{j=1}^n\frac{b_{2j}{\eta }_j^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}{L}_1+\frac{L_2}{\Gamma ({\alpha }_2+1)}\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}{L}_2}{\Gamma \left(l_{1i}\right)\Gamma ({\alpha }_2+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\alpha }_2}ds+\sum _{j=1}^n\frac{b_{1j}{\eta }_j^{{\alpha }_2}}{\Gamma ({\alpha }_2+1)}{L}_2+\frac{L_1}{\Gamma ({\alpha }_1+1)}]\hfill \\ \le L_1{M}_1+L_2{M}_2,\hfill \end{array}$$

(71)

which yields

$$\parallel T_1(u,v){\parallel }_{C_{2-{\theta }_1}}\le L_1{M}_1+L_2{M}_2,$$

(72)

$$\begin{array}{c}|t^{2-{\theta }_1}{\hspace{0.17em}}^H{D}^{{\gamma }_1,\beta }{T}_1(u,v)\left(t\right)|\le I^{{\alpha }_1-{\gamma }_1}|f_1(u,v)|\left(t\right)+\frac{\Gamma \left({\theta }_1\right)t^{{\theta }_1-{\gamma }_1-1}}{\Gamma ({\theta }_1-{\gamma }_1)(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}|f_1(u,v)|\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}|f_1(u,v)|\left(s\right)ds+I^{{\alpha }_2}|f_2(u,v)|\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}|f_2(u,v)|\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}|f_2(u,v)|\left(s\right)ds+I^{{\alpha }_1}|f_1(u,v)|\left(1\right)]\hfill \\ \le \frac{t^{{\alpha }_1-{\gamma }_1}}{\Gamma ({\alpha }_1-{\gamma }_1+1)}{L}_1+\frac{\Gamma \left({\theta }_1\right)t^{{\theta }_1-{\gamma }_1-1}}{\Gamma ({\theta }_1-{\gamma }_1)(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}{L}_1}{\Gamma \left(l_{2i}\right)\Gamma ({\alpha }_1+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\alpha }_1}ds+\sum _{j=1}^n\frac{b_{2j}{\eta }_j^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}{L}_1+\frac{L_2}{\Gamma ({\alpha }_2+1)}\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}{L}_2}{\Gamma \left(l_{1i}\right)\Gamma ({\alpha }_2+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\alpha }_2}ds+\sum _{j=1}^n\frac{b_{1j}{\eta }_j^{{\alpha }_2}}{\Gamma ({\alpha }_2+1)}{L}_2+\frac{L_1}{\Gamma ({\alpha }_1+1)}]\hfill \\ \le L_1{M}_3+L_2{M}_4,\hfill \end{array}$$

(73)

which yields

$$\parallel {\hspace{0.17em}}^H{D}^{{\gamma }_1,\beta }{T}_1{(u,v)\parallel }_{C_{2-{\theta }_1}}\le L_1{M}_3+L_2{M}_4.$$

(74)

So we can get

$$\parallel T_1(u,v){\parallel }_{C_1}\le L_1(M_1+M_3)+L_2(M_2+M_4).$$

(75)

Similarly, we obtain that

$$\parallel T_2{(u,v)\parallel }_{C_{2-{\theta }_2}}\le L_1{M}_5+L_2{M}_6,$$

(76)

$${\parallel }^H{D}^{{\gamma }_2,\beta }{T}_2{(u,v)\parallel }_{C_{2-{\theta }_2}}\le L_1{M}_7+L_2{M}_8,$$

(77)

and

$$\parallel T_2{(u,v)\parallel }_{C_2}\le L_1(M_5+M_7)+L_2(M_6+M_8).$$

(78)

Hence, for $(u,v)\in \Omega$, $T_1,T_2$ is uniformly bounded. Thus it follows from the above inequalities that the set $T\Omega$ is uniformly bounded.

(c) The operator $T\Omega$ is equicontinous.

For any $(u,v)\in \Omega$ and $t_1,t_2\in J$ such that $t_1<t_2$, we have

$$\begin{array}{c}|t_2^{2-{\theta }_1}{T}_1(u,v)\left(t_2\right)-t_1^{2-{\theta }_1}{T}_1(u,v)\left(t_1\right)|\hfill \\ \le \frac{1}{\Gamma \left({\alpha }_1\right)}\left|{\int }_0^{t_1}\left[t_2^{2-{\theta }_1}{(t_2-s)}^{{\alpha }_1-1}-t_1^{2-{\theta }_1}{(t_1-s)}^{{\alpha }_1-1}\right]f_1(u,v)\left(s\right)ds+{\int }_{t_1}^{t_2}{t}_2^{2-{\theta }_1}{(t_2-s)}^{{\alpha }_1-1}{f}_1(u,v)\left(s\right)ds\right|+\frac{t_2-t_1}{(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}|f_1(u,v)|\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}|f_1(u,v)\left(s\right)|ds+I^{{\alpha }_2}|f_2(u,v)|\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}|f_2(u,v)|\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}|f_2(u,v)|\left(s\right)ds+I^{{\alpha }_1}|f_1(u,v)|\left(1\right)]\hfill \\ \le \frac{t_2^{{\alpha }_1}-t_1^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}{L}_1+\frac{t_2-t_1}{(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}{L}_1}{\Gamma \left(l_{2i}\right)\Gamma ({\alpha }_1+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\alpha }_1}ds+\sum _{j=1}^n\frac{b_{2j}{\eta }_j^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}{L}_1+\frac{L_2}{\Gamma ({\alpha }_2+1)}\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}{L}_2}{\Gamma \left(l_{1i}\right)\Gamma ({\alpha }_2+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\alpha }_2}ds+\sum _{j=1}^n\frac{b_{1j}{\eta }_j^{{\alpha }_2}}{\Gamma ({\alpha }_2+1)}{L}_2+\frac{L_1}{\Gamma ({\alpha }_1+1)}]\hfill \\ \to 0as{t}_2\to t_1,\hfill \end{array}$$

(79)

and

$$\begin{array}{c}|t_2^{2-{\theta }_1}{\hspace{0.17em}}^H{D}^{{\gamma }_1,\beta }{T}_1(u,v)\left(t_2\right)-t_1^{2-{\theta }_1}{\hspace{0.17em}}^H{D}^{{\gamma }_1,\beta }{T}_1(u,v)\left(t_1\right)|\hfill \\ \le \frac{1}{\Gamma ({\alpha }_1-{\gamma }_1)}\left|{\int }_0^{t_1}\left[t_2^{2-{\theta }_1}{(t_2-s)}^{{\alpha }_1-{\gamma }_1-1}-t_1^{2-{\theta }_1}{(t_1-s)}^{{\alpha }_1-{\gamma }_1-1}\right]f_1(u,v)\left(s\right)ds+{\int }_{t_1}^{t_2}{t}_2^{2-{\theta }_1}{(t_2-s)}^{{\alpha }_1-{\gamma }_1-1}{f}_1(u,v)\left(s\right)ds\right|\hfill \\ +\frac{\Gamma \left({\theta }_1\right)(t_2^{1-{\gamma }_1}-t_{!}^{1-{\gamma }_1})}{\Gamma ({\theta }_1-{\gamma }_1)(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}|f_1(u,v)|\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}|f_1(u,v)|\left(s\right)ds+I^{{\alpha }_2}|f_2(u,v)|\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}|f_2(u,v)|\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}|f_2(u,v)|\left(s\right)ds+I^{{\alpha }_1}|f_1(u,v)|\left(1\right)]\hfill \\ \le \frac{t_2^{{\alpha }_1-{\gamma }_1}-t_1^{{\alpha }_1-{\gamma }_1}}{\Gamma ({\alpha }_1-{\gamma }_1+1)}{L}_1+\frac{\Gamma \left({\theta }_1\right)(t_2^{1-{\gamma }_1}-t_{!}^{1-{\gamma }_1})}{\Gamma ({\theta }_1-{\gamma }_1)(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}{L}_1}{\Gamma \left(l_{2i}\right)\Gamma ({\alpha }_1+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\alpha }_1}ds+\sum _{j=1}^n\frac{b_{2j}{\eta }_j^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}{L}_1+\frac{L_2}{\Gamma ({\alpha }_2+1)}\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}{L}_2}{\Gamma \left(l_{1i}\right)\Gamma ({\alpha }_2+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\alpha }_2}ds+\sum _{j=1}^n\frac{b_{1j}{\eta }_j^{{\alpha }_2}}{\Gamma ({\alpha }_2+1)}{L}_2+\frac{L_1}{\Gamma ({\alpha }_1+1)}]\hfill \\ \to 0as{t}_2\to t_1.\hfill \end{array}$$

(80)

Therefore the set $T_1\Omega$ is equicontinous for all $(u,v)\in \Omega$. Similarly, we can get the set $T_2\Omega$ is equicontinous for all $(u,v)\in \Omega$. As a consequence, the set $T\Omega$ is equicontinous for all $(u,v)\in \Omega$. By applying the Arzelá–Ascoli theorem, the set $T\Omega$ is relative compact which implies that the operator T is completely continuous.

Lastly, we shall show that the set $\xi =\left\{\right(u,v)\in \mathcal{C}:(u,v)=\lambda T(u,v),0\le \lambda \le 1\}$ is bounded. Let any $(u,v)\in \xi ,$ then $(u,v)=\lambda T(u,v).$ For any $t\in J,$ we have

$$\begin{array}{c}\hfill u\left(t\right)=\lambda T_1(u,v)\left(t\right),\\ \hfill v\left(t\right)=\lambda T_2(u,v)\left(t\right).\end{array}$$

(81)

Then, we get

$$\begin{array}{c}{\parallel u\parallel }_{C_1}\le \left[a_0+(a_1+a_3){\parallel u\parallel }_{C_1}+(a_2+a_4){\parallel v\parallel }_{C_2}\right](M_1+M_3)+\left[b_0+(b_1+b_3){\parallel u\parallel }_{C_1}+(b_2+b_4){\parallel v\parallel }_{C_2}\right](M_2+M_4),\hfill \\ {\parallel v\parallel }_{C_2}\le \left[a_0+(a_1+a_3){\parallel u\parallel }_{C_1}+(a_2+a_4){\parallel v\parallel }_{C_2}\right](M_5+M_7)+\left[b_0+(b_1+b_3){\parallel u\parallel }_{C_1}+(b_2+b_4){\parallel v\parallel }_{C_2}\right](M_6+M_8),\hfill \end{array}$$

(82)

which imply that

$$\begin{array}{cc}\hfill {\parallel u\parallel }_{C_1}+{\parallel v\parallel }_{C_2}& \le (M_1+M_3+M_5+M_7)a_0+(M_2+M_4+M_6+M_8)b_0\hfill \\ & +[(a_1+a_3)(M_1+M_3+M_5+M_7)+(b_1+b_3)(M_2+M_4+M_6+M_8)]{\parallel u\parallel }_{C_1}\hfill \\ & +[(a_2+a_4)(M_1+M_3+M_5+M_7)+(b_2+b_4)(M_2+M_4+M_6+M_8)]{\parallel v\parallel }_{C_2}.\hfill \end{array}$$

(83)

Thus, we obtain

$${\parallel (u,v)\parallel }_{\mathcal{C}}\le \frac{(M_1+M_3+M_5+M_7)a_0+(M_2+M_4+M_6+M_8)b_0}{M^{*}},$$

(84)

where $M^{*}=min\{1-(a_1+a_3)(M_1+M_3+M_5+M_7)-(b_1+b_3)(M_2+M_4+M_6+M_8),\hspace{0.17em}\hspace{0.17em}1-(a_2+a_4)(M_1+M_3+M_5+M_7)-(b_2+b_4)(M_2+M_4+M_6+M_8)\},$ which shows that the set $\xi$ is bounded. Therefore, by applying Lemma 11, the operator T has at least one fixed point. Therefore, we deduce that problem (5) and (7) has at least one solution on J.

The proof is completed. □

3.2. The Existence Results of the Coupled Hilfer Fractional Differential Inclusions

X is a real (or complex) separable Banach space with a norm $\parallel ·\parallel$, defined by $\parallel u\parallel ={sup}_{t\in J}\left|u\left(t\right)\right|$, $\mathcal{P}\left(X\right)$ is the family of all nonempty subsets of X. For a normed space $(X,\parallel ·\parallel )$, let Y be a subset of X. We denote

(i)

$\mathcal{P}$(X) = {Y ⊆ X: $Y\ne \varnothing$};

(ii)

${\mathcal{P}}_{cl}$(X) = {$Y\in \mathcal{P}\left(X\right)$: Y closed};

(iii)

${\mathcal{P}}_b$(X) = {$Y\in \mathcal{P}\left(X\right)$: Y bounded};

(iv)

${\mathcal{P}}_{cp}$(X) = {$Y\in \mathcal{P}\left(X\right)$: Y compact};

(v)

${\mathcal{P}}_{cv}$(X) = {$Y\in \mathcal{P}\left(X\right)$: Y convex};

(vi)

${\mathcal{P}}_{cp,cv}$(X) = {$Y\in \mathcal{P}\left(X\right)$: Y compact, convex};

$F,G$: $J\times X^4\to \mathcal{P}\left(X\right),$ $i=1,2,$ are given multivalued maps. When $F,G$ are convex valued, to complete our result we need the following assumtions:

(H1)

$F,G:J\times X^4\to {\mathcal{P}}_{cp,cv}\left(X\right)$ are $L^1$-Carathéodory multivalued maps;

(H2)

There exist $m_1,m_2\in C(J,X^{+})$ and ${\varphi }_1,{\varphi }_2,{\psi }_1,{\psi }_2,{\varrho }_1,{\varrho }_2,{\rho }_1,{\rho }_2:[0,\infty )\to (0,\infty )$ continuous, nondecreasing such that

$$∥F(t,u,v,w,z)∥=sup\left\{\right|f|:f\in F(t,u,v,w,z)\}\le m_1\left(t\right)({\varphi }_1\left(\right|u\left|\right)+{\psi }_1\left(\right|v\left|\right)+{\varrho }_1\left(\right|w\left|\right)+{\rho }_1\left(\right|z\left|\right)),$$

(85)

and

$$∥G(t,u,v,w,z)∥=sup\left\{\right|g|:g\in F(t,u,v,w,z)\}\le m_2\left(t\right)({\varphi }_2\left(\right|u\left|\right)+{\psi }_2\left(\right|v\left|\right)+{\varrho }_2\left(\right|w\left|\right)+{\rho }_2\left(\right|z\left|\right)),$$

(86)

for $u,v,w,z\in X$ and a.e. $t\in J$.

Lemma 15. 

[28] Let X be a Banach space. Let $F:J\times X^4\to {\mathcal{P}}_{cp,cv}\left(X\right)$ be an $L^1$-Carathéodory multivalued map and T be a linear continuous mapping from $L^1(J,X)$ to $C(J,X).$ Then the operator

$$T\circ S_F:C(J,X)\to {\mathcal{P}}_{cp,cv}\left(C(J,X)\right),y\mapsto (T\circ S_F)\left(y\right)=T(S_F,y),$$

(87)

is a closed graph operator in $C(J,X)\times C(J,X)$.

Theorem 2. 

Assume that (H1), (H2) are satisfied and there exists $K>0$ such that

$$\begin{array}{cc}\hfill K& >(M_1+M_3+M_5+M_7)∥m_1∥({\varphi }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\psi }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\varrho }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\rho }_1(\parallel (u,v){\parallel }_{\mathcal{C}}\left)\right)\hfill \\ & +(M_2+M_4+M_6+M_8)∥m_2∥({\varphi }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\psi }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\varrho }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\rho }_2(\parallel (u,v){\parallel }_{\mathcal{C}}\left)\right),\hfill \end{array}$$

where $M_i(i=1,2,3,4,5,6,7,8)$ are defined by (53)–(60), then the coupled system (6) and (7) has at least one solution on J.

Proof. 

For each $(u,v)\in \mathcal{C}$, define the sets of selections of $F,G$ by

$$S_{F,(u,v)}=\left\{f\left(t\right)\in L^1(J,E):f\left(t\right)\in F(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right))fora.e.t\in J\right\},$$

(88)

and

$$S_{G,(u,v)}=\left\{g\left(t\right)\in L^1(J,E):g\left(t\right)\in G(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right))fora.e.t\in J\right\}.$$

(89)

Define the multivalued operators $N_1:C_1\to \mathcal{P}\left(C_1\right)$ and $N_2:C_2\to \mathcal{P}\left(C_2\right)$ by

$$N_1(u,v)=\left\{h_1\in C_1:thereexistf\in S_{F,(u,v)}andg\in S_{G,(u,v)}suchthat{h}_1(u,v)\left(t\right)=A_1(t,u,v)\right\},$$

(90)

and

$$N_2(u,v)=\left\{h_2\in C_2:thereexistf\in S_{F,(u,v)}andg\in S_{G,(u,v)}suchthat{h}_2(u,v)\left(t\right)=A_2(t,u,v)\right\},$$

(91)

where

$$\begin{array}{c}{A}_1(t,u,v)=I^{{\alpha }_1}f\left(t\right)+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}f\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}f\left(s\right)ds-I^{{\alpha }_2}g\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}g\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}g\left(s\right)ds-I^{{\alpha }_1}f\left(1\right)],\hfill \end{array}$$

(92)

$$\begin{array}{c}{A}_2(t,u,v)=I^{{\alpha }_2}g\left(t\right)+\frac{t^{{\theta }_2-1}}{(1-l_1{l}_2)}\hfill \\ [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}g\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}g\left(s\right)ds-I^{{\alpha }_1}f\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}f\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}f\left(s\right)ds-I^{{\alpha }_2}g\left(1\right)].\hfill \end{array}$$

(93)

Consider the continuous operator $N:\mathcal{C}\to \mathcal{P}\left(\mathcal{C}\right)$ defined by

$$N(u,v)=\{(h_1,h_2)\in \mathcal{C}:h_1\in N_1(u,v),h_2\in N_2(u,v)\}.$$

(94)

Clearly, the fixed points of N are solutions of the system (6) and (7).

Step 1. $N(u,v)$ is convex valued.

Suppose $(h_i,\widetilde{h_i})\in (N_1,N_2)(i=1,2).$ Then there exist $f_i\in S_{F,(u,v)}$, $g_i\in S_{G,(u,v)}\hspace{0.17em}\hspace{0.17em}(i=1,2)$ such that for any $t\in J,t=1,2,$ we have

$$\begin{array}{c}{h}_i\left(t\right)=I^{{\alpha }_1}{f}_i\left(t\right)+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_i\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_i\left(s\right)ds-I^{{\alpha }_2}{g}_i\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{g}_i\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{g}_i\left(s\right)ds-I^{{\alpha }_1}{f}_i\left(1\right)],\hfill \end{array}$$

(95)

$$\begin{array}{c}\widetilde{h_i}\left(t\right)=I^{{\alpha }_2}{g}_i\left(t\right)+\frac{t^{{\theta }_2-1}}{(1-l_1{l}_2)}\hfill \\ [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{g}_i\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{g}_i\left(s\right)ds-I^{{\alpha }_1}{f}_i\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_i\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_i\left(s\right)ds-I^{{\alpha }_2}{g}_i\left(1\right)].\hfill \end{array}$$

(96)

Let $0\le \lambda \le 1$. Then, for any $t\in J$, we have

$$\begin{array}{c}(\lambda h_1+(1-\lambda )h_2)\left(t\right)=I^{{\alpha }_1}(\lambda f_1+(1-\lambda )f_2)\left(t\right)+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ \left[l_1\right(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}(\lambda f_1+(1-\lambda )f_2)\left(\tau \right)d\tau ds\hfill \\ +\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}(\lambda f_1+(1-\lambda )f_2)\left(s\right)ds-I^{{\alpha }_2}(\lambda g_1+(1-\lambda )g_2)\left(1\right))\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}(\lambda g_1+(1-\lambda )g_2)\left(\tau \right)d\tau ds\hfill \\ +\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}(\lambda g_1+(1-\lambda )g_2)\left(s\right)ds-I^{{\alpha }_1}(\lambda f_1+(1-\lambda )f_2)\left(1\right)],\hfill \end{array}$$

(97)

and

$$\begin{array}{c}(\lambda {\tilde{h}}_1+(1-\lambda ){\tilde{h}}_2)\left(t\right)=I^{{\alpha }_2}(\lambda g_1+(1-\lambda )g_2)\left(t\right)+\frac{t^{{\theta }_2-1}}{(1-l_1{l}_2)}\hfill \\ \left[l_2\right(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}(\lambda g_1+(1-\lambda )g_2)\left(\tau \right)d\tau ds\hfill \\ +\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}(\lambda g_1+(1-\lambda )g_2)\left(s\right)ds-I^{{\alpha }_1}(\lambda f_1+(1-\lambda )f_2)\left(1\right))\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}(\lambda f_1+(1-\lambda )f_2)\left(\tau \right)d\tau ds\hfill \\ +\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}(\lambda f_1+(1-\lambda )f_2)\left(s\right)ds-I^{{\alpha }_2}(\lambda g_1+(1-\lambda )g_2)\left(1\right)].\hfill \end{array}$$

(98)

Since F and G are convex valued, we infer that $S_{F,(u,v)}$ and $S_{G,(u,v)}$ are convex. Obviously, $\lambda h_1+(1-\lambda )h_2\in N_1$, $\lambda {\tilde{h}}_1+(1-\lambda ){\tilde{h}}_2\in N_2$. Therefore, $\lambda (h_1,{\tilde{h}}_1)+(1-\lambda )(h_2,{\tilde{h}}_2)\in N.$

Step 2. N maps bounded sets into bounded sets in $\mathcal{C}$.

Let $r>0$, $B_r=\{(u,v)\in \mathcal{C}:{∥(u,v)∥}_{\mathcal{C}}\le r\}$ be a bounded subset of $\mathcal{C}$, $(h_1,h_2)\in N(u,v)$ and $(u,v)\in B_r.$ Then there exist $f\in S_{F,(u,v)}$ and $g\in S_{G,(u,v)}$ such that for any $t\in J$,

$$\begin{array}{c}{h}_1(u,v)\left(t\right)=I^{{\alpha }_1}{f}_1\left(t\right)+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_1\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_1\left(s\right)ds-I^{{\alpha }_2}{f}_2\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{f}_2\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{f}_2\left(s\right)ds-I^{{\alpha }_1}{f}_1\left(1\right)],\hfill \end{array}$$

(99)

$$\begin{array}{c}{h}_2(u,v)\left(t\right)=I^{{\alpha }_2}{f}_2\left(t\right)+\frac{t^{{\theta }_2-1}}{(1-l_1{l}_2)}\hfill \\ [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{f}_2\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{f}_2\left(s\right)ds-I^{{\alpha }_1}{f}_1\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_1\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_1\left(s\right)ds-I^{{\alpha }_2}{f}_2\left(1\right)].\hfill \end{array}$$

(100)

We have

$$\begin{array}{c}|t^{2-{\theta }_1}{h}_1(u,v)|\le \frac{1}{\Gamma ({\alpha }_1+1)}∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+\frac{1}{(1-l_1{l}_2)}\hfill \\ \left[l_1\right(\sum _{i=1}^m\frac{{\lambda }_{2i}∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))}{\Gamma \left(l_{2i}\right)\Gamma ({\alpha }_1+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\alpha }_1}ds\hfill \\ +\sum _{j=1}^n\frac{b_{2j}{\eta }_j^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+\frac{∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right))}{\Gamma ({\alpha }_2+1)})\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right))}{\Gamma \left(l_{1i}\right)\Gamma ({\alpha }_2+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\alpha }_2}ds\hfill \\ +\sum _{j=1}^n\frac{b_{1j}{\eta }_j^{{\alpha }_2}}{\Gamma ({\alpha }_2+1)}∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right))+\frac{∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))}{\Gamma ({\alpha }_1+1)}]\hfill \\ \le M_1∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+M_2∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right)),\hfill \end{array}$$

(101)

which yields

$$\parallel h_1(u,v){\parallel }_{C_{2-{\theta }_1}}\le M_1∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+M_2∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right)),$$

(102)

$$\begin{array}{c}|t^{2-{\theta }_1}{\hspace{0.17em}}^H{D}^{{\gamma }_1,\beta }{h}_1(u,v)\left(t\right)|\le \frac{1}{\Gamma ({\alpha }_1-{\gamma }_1+1)}{M}_1∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+\frac{\Gamma \left({\theta }_1\right)}{\Gamma ({\theta }_1-{\gamma }_1)(1-l_1{l}_2)}\hfill \\ \left[l_1\right(\sum _{i=1}^m\frac{{\lambda }_{2i}{M}_1∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))}{\Gamma \left(l_{2i}\right)\Gamma ({\alpha }_1+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\alpha }_1}ds\hfill \\ +\sum _{j=1}^n\frac{b_{2j}{\eta }_j^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}{M}_1∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+\frac{∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right))}{\Gamma ({\alpha }_2+1)})\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right))}{\Gamma \left(l_{1i}\right)\Gamma ({\alpha }_2+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\alpha }_2}ds\hfill \\ +\sum _{j=1}^n\frac{b_{1j}{\eta }_j^{{\alpha }_2}}{\Gamma ({\alpha }_2+1)}∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right))+\frac{M_1∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))}{\Gamma ({\alpha }_1+1)}]\hfill \\ \le M_3∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+M_4∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right)),\hfill \end{array}$$

(103)

which yields

$${\parallel }^H{D}^{{\gamma }_1,\beta }{h}_1{(u,v)\parallel }_{C_{2-{\theta }_1}}\le M_3∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+M_4∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right)).$$

(104)

Thus

$$\parallel h_1(u,v){\parallel }_{C_1}\le (M_1+M_3)∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+(M_2+M_4)∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right)).$$

(105)

In a similar manner, we have

$$\parallel h_2{(u,v)\parallel }_{C_{2-{\theta }_2}}\le M_5∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+M_6∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right)),$$

(106)

$${\parallel }^H{D}^{{\gamma }_2,\beta }{h}_2{(u,v)\parallel }_{C_{2-{\theta }_2}}\le M_7∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+M_8∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right)),$$

(107)

and

$$\parallel h_2{(u,v)\parallel }_{C_2}\le (M_5+M_7)∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+(M_6+M_8)∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right)).$$

(108)

Hence we have

$$\begin{array}{cc}\hfill \parallel (h_1,h_2){\parallel }_{\mathcal{C}}& =\parallel h_1(u,v){\parallel }_{C_1}+{\parallel h_2(u,v)\left(t\right)\parallel }_{C_2}\hfill \\ & \le (M_1+M_3+M_5+M_7)∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))\hfill \\ & +(M_2+M_4+M_6+M_8)∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right)).\hfill \end{array}$$

(109)

Step 3. N maps bounded sets into equicontinuous sets in $\mathcal{C}$.

Let $B_r$ be a bounded set of $\mathcal{C}$ as in step 2. Let $0\le t_1\le t_2\le 1$ and $(u,v)\in B_r.$

$$\begin{array}{c}|t_2^{2-{\theta }_1}{h}_1(u,v)\left(t_2\right)-t_1^{2-{\theta }_1}{h}_1(u,v)\left(t_1\right)|\le \frac{t_2^{{\alpha }_1}-t_1^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+\frac{t_2-t_1}{(1-l_1{l}_2)}\hfill \\ \left[l_1\right(\sum _{i=1}^m\frac{{\lambda }_{2i}∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))}{\Gamma \left(l_{2i}\right)\Gamma ({\alpha }_1+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\alpha }_1}ds\hfill \\ +\sum _{j=1}^n\frac{b_{2j}{\eta }_j^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+\frac{∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right))}{\Gamma ({\alpha }_2+1)})\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right))}{\Gamma \left(l_{1i}\right)\Gamma ({\alpha }_2+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\alpha }_2}ds\hfill \\ +\sum _{j=1}^n\frac{b_{1j}{\eta }_j^{{\alpha }_2}}{\Gamma ({\alpha }_2+1)}∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right))+\frac{∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))}{\Gamma ({\alpha }_1+1)}]\hfill \\ \to 0as{t}_2\to t_1,\hfill \end{array}$$

(110)

$$\begin{array}{c}|t_2^{2-{\theta }_1}{\hspace{0.17em}}^H{D}^{{\gamma }_1,\beta }{h}_1(u,v)\left(t_2\right)-t_1^{2-{\theta }_1}{\hspace{0.17em}}^H{D}^{{\gamma }_1,\beta }{h}_1(u,v)\left(t_1\right)|\le \frac{t_2^{{\alpha }_1-{\gamma }_1}-t_1^{{\alpha }_1-{\gamma }_1}}{\Gamma ({\alpha }_1-{\gamma }_1+1)}∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))\hfill \\ +\frac{\Gamma \left({\theta }_1\right)(t_2^{1-{\gamma }_1}-t_1^{1-{\gamma }_1})}{\Gamma ({\theta }_1-{\gamma }_1)(1-l_1{l}_2)}\left[l_1\right(\sum _{i=1}^m\frac{{\lambda }_{2i}∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))}{\Gamma \left(l_{2i}\right)\Gamma ({\alpha }_1+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{s}^{{\alpha }_1}ds\hfill \\ +\sum _{j=1}^n\frac{b_{2j}{\eta }_j^{{\alpha }_1}}{\Gamma ({\alpha }_1+1)}∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))+\frac{∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right))}{\Gamma ({\alpha }_2+1)})\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right))}{\Gamma \left(l_{1i}\right)\Gamma ({\alpha }_2+1)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{s}^{{\alpha }_2}ds\hfill \\ +\sum _{j=1}^n\frac{b_{1j}{\eta }_j^{{\alpha }_2}}{\Gamma ({\alpha }_2+1)}∥m_2∥({\varphi }_2\left(r\right)+{\psi }_2\left(r\right)+{\varrho }_2\left(r\right)+{\rho }_2\left(r\right))+\frac{∥m_1∥({\varphi }_1\left(r\right)+{\psi }_1\left(r\right)+{\varrho }_1\left(r\right)+{\rho }_1\left(r\right))}{\Gamma ({\alpha }_1+1)}]\hfill \\ \to 0as{t}_2\to t_1.\hfill \end{array}$$

(111)

Analogously, we can obtain

$$\parallel h_2(u,v)\left(t_2\right)-h_2(u,v)\left(t_1\right)\parallel \to 0as{t}_2\to t_1,$$

(112)

$${\parallel }^H{D}^{{\gamma }_2,\beta }{h}_2(u,v)\left(t_2\right)-{\hspace{0.17em}}^H{D}^{{\gamma }_2,\beta }{h}_2(u,v)\left(t_1\right)\parallel \to 0as{t}_2\to t_1.$$

(113)

Therefore, the operator $N(u,v)$ is equicontinuous. By the Arzelá-Ascoli theorem, we infer that the operator $N(u,v)$ is completely continuous.

Step 4. N has a closed graph.

Let $(u_n,v_n)\to (u_{*},v_{*}),$ $(h_n,{\overline{h}}_n)\in N(u_n,v_n)$ and $(h_n,{\overline{h}}_n)\to (h_{*},{\overline{h}}_{*})$, we need to proof $(h_{*},{\overline{h}}_{*})\in N(u_{*},v_{*})$. $(h_n,{\overline{h}}_n)\in N(u_n,v_n)$ implies that there exist $f_n\in S_{F,(u_n,v_n)}$ and $g_n\in S_{G,(u_n,v_n)}$ such that for all $t\in J$,

$$\begin{array}{c}{h}_n(u_n,v_n)\left(t\right)=I^{{\alpha }_1}{f}_n\left(t\right)+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_n\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_n\left(s\right)ds-I^{{\alpha }_2}{g}_n\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{g}_n\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{g}_n\left(s\right)ds-I^{{\alpha }_1}{f}_n\left(1\right)],\hfill \end{array}$$

(114)

$$\begin{array}{c}{\overline{h}}_n(u_n,v_n)\left(t\right)=I^{{\alpha }_2}{g}_n\left(t\right)+\frac{t^{{\theta }_2-1}}{(1-l_1{l}_2)}\hfill \\ [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{g}_n\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{g}_n\left(s\right)ds-I^{{\alpha }_1}{f}_n\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_n\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_n\left(s\right)ds-I^{{\alpha }_2}{g}_n\left(1\right)].\hfill \end{array}$$

(115)

Let us consider the continuous linear operators ${\Phi }_1,{\Phi }_2:$ $L^1(J,C)\to C(J,C)$ given by

$$\begin{array}{c}{\Phi }_1(u,v)\left(t\right)=I^{{\alpha }_1}f\left(t\right)+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}f\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}f\left(s\right)ds-I^{{\alpha }_2}g\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}g\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}g\left(s\right)ds-I^{{\alpha }_1}f\left(1\right)],\hfill \end{array}$$

(116)

$$\begin{array}{c}{\Phi }_2(u,v)\left(t\right)=I^{{\alpha }_2}g\left(t\right)+\frac{t^{{\theta }_2-1}}{(1-l_1{l}_2)}\hfill \\ [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}g\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}g\left(s\right)ds-I^{{\alpha }_1}f\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}f\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}f\left(s\right)ds-I^{{\alpha }_2}g\left(1\right)].\hfill \end{array}$$

(117)

From Lemma 15, we know that $({\Phi }_1,{\Phi }_2)\circ (S_F,S_G)$ is a closed graph operator. Moreover, we get $(h_n,{\overline{h}}_n)\in ({\Phi }_1,{\Phi }_2)\circ (S_{F,(u_n,v_n)},S_{G,(u_n,v_n)})$ for all n. Since $(u_n,v_n)\to (u_{*},v_{*})$, $(h_n,{\overline{h}}_n)\to (h_{*},{\overline{h}}_{*})$, it follows the existence of $f_{*}\in S_{F,(u_{*},v_{*})}$ and $g_{*}\in S_{G,(u_{*},v_{*})}$ such that

$$\begin{array}{c}{h}_{*}(u_{*},v_{*})\left(t\right)=I^{{\alpha }_1}{f}_{*}\left(t\right)+\frac{t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_{*}\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_{*}\left(s\right)ds-I^{{\alpha }_2}{g}_{*}\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{g}_{*}\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{g}_{*}\left(s\right)ds-I^{{\alpha }_1}{f}_{*}\left(1\right)],\hfill \end{array}$$

(118)

$$\begin{array}{c}{\overline{h}}_{*}(u_{*},v_{*})\left(t\right)=I^{{\alpha }_2}{g}_{*}\left(t\right)+\frac{t^{{\theta }_2-1}}{(1-l_1{l}_2)}\hfill \\ [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}{g}_{*}\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}{g}_{*}\left(s\right)ds-I^{{\alpha }_1}{f}_{*}\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}{f}_{*}\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}{f}_{*}\left(s\right)ds-I^{{\alpha }_2}{g}_{*}\left(1\right)],\hfill \end{array}$$

(119)

that is, $(h_{*},{\overline{h}}_{*})\in N(u_{*},v_{*}).$

Step 5. A priori bounds on solutions.

Let $(u,v)\in \lambda N(u,v)$ for some $\lambda \in (0,1).$ Then there exist $f\in S_{F,(u,v)}$ and $g\in S_{G,(u,v)}$ such that for all $t\in J$,

$$\begin{array}{c}u\left(t\right)=\lambda I^{{\alpha }_1}f\left(t\right)+\frac{\lambda t^{{\theta }_1-1}}{(1-l_1{l}_2)}\hfill \\ [l_1\left(\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}f\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}f\left(s\right)ds-I^{{\alpha }_2}g\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}g\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}g\left(s\right)ds-I^{{\alpha }_1}f\left(1\right)],\hfill \end{array}$$

(120)

and

$$\begin{array}{c}v\left(t\right)=\lambda I^{{\alpha }_2}g\left(t\right)+\frac{\lambda t^{{\theta }_2-1}}{(1-l_1{l}_2)}\hfill \\ [l_2\left(\sum _{i=1}^m\frac{{\lambda }_{1i}}{\Gamma \left(l_{1i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{1i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_2-1}g\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{1j}}{\Gamma \left({\alpha }_2\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_2-1}g\left(s\right)ds-I^{{\alpha }_1}f\left(1\right)\right)\hfill \\ +\sum _{i=1}^m\frac{{\lambda }_{2i}}{\Gamma \left(l_{2i}\right)}{\int }_0^{{\xi }_i}{\left({\xi }_i-s\right)}^{(l_{2i}-1)}{\int }_0^s{(s-\tau )}^{{\alpha }_1-1}f\left(\tau \right)d\tau ds+\sum _{j=1}^n\frac{b_{2j}}{\Gamma \left({\alpha }_1\right)}{\int }_0^{{\eta }_j}{({\eta }_j-s)}^{{\alpha }_1-1}f\left(s\right)ds-I^{{\alpha }_2}g\left(1\right)].\hfill \end{array}$$

(121)

With the same arguments as in Step 2 of our proof, for each $(u,v)\in \mathcal{C},$ we obtain

$$\begin{array}{cc}\hfill {\parallel u\parallel }_{C_1}& \le (M_1+M_3)∥m_1∥({\varphi }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\psi }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\varrho }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\rho }_1(\parallel (u,v){\parallel }_{\mathcal{C}}\left)\right)\hfill \\ & +(M_2+M_4)∥m_2∥({\varphi }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\psi }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\varrho }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\rho }_2(\parallel (u,v){\parallel }_{\mathcal{C}}\left)\right),\hfill \end{array}$$

(122)

and

$$\begin{array}{cc}\hfill {\parallel v\parallel }_{C_2}& \le (M_5+M_7)∥m_1∥({\varphi }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\psi }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\varrho }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\rho }_1(\parallel (u,v){\parallel }_{\mathcal{C}}\left)\right)\hfill \\ & +(M_6+M_8)∥m_2∥({\varphi }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\psi }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\varrho }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\rho }_2(\parallel (u,v){\parallel }_{\mathcal{C}}\left)\right).\hfill \end{array}$$

(123)

Hence we have

$$\begin{array}{cc}\hfill {\parallel (u,v)\parallel }_{\mathcal{C}}& ={\parallel u\parallel }_{C_1}+{\parallel v\parallel }_{C_2}\hfill \\ & \le (M_1+M_3+M_5+M_7)∥m_1∥({\varphi }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\psi }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\varrho }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\rho }_1(\parallel (u,v){\parallel }_{\mathcal{C}}\left)\right)\hfill \\ & +(M_2+M_4+M_6+M_8)∥m_2∥({\varphi }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\psi }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\varrho }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\rho }_2(\parallel (u,v){\parallel }_{\mathcal{C}}\left)\right).\hfill \end{array}$$

(124)

Now we set $U=\{(u,v)\in \mathcal{C}:{∥(u,v)∥}_{\mathcal{C}}<K\}$. Clearly, U is an open subset of $\mathcal{C}$ and $(0,0)\in U.$ As a consequence of Steps 1–4, together with the Arzelá-Ascoli theorem, we can conclude that $N:\overline{U}\to {\mathcal{P}}_{cp,cv}\left(C_1\right)\times {\mathcal{P}}_{cp,cv}\left(C_2\right)$ is upper semicontinuous and completely continuous. From the choice of U, there is no $(u,v)\in \partial U$ such that $(u,v)\in \lambda N(u,v)$ for some $\lambda \in (0,1)$. Therefore, by Lemma 12, we deduce that N has a fixed point $(u,v)\in \overline{U},$ which is a solution of the coupled system (6) and (7).

This completes the proof. □

Example 1. 

For $t\in J$, consider the following fractional differential system:

$$\left\{\begin{array}{c}{\hspace{0.17em}}^H{D}^{2,0.5}u\left(t\right)=f_1(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right)),\hfill \\ {\hspace{0.17em}}^H{D}^{2,0.5}v\left(t\right)=f_2(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right)),\hfill \end{array}\right.$$

(125)

with the coupled integral and discrete mixed boundary conditions:

$$\left\{\begin{array}{c}\hfill u\left(0\right)=0,u\left(1\right)=\sum _{i=1}^2{\lambda }_{1i}{I}^{l_{1i}}v\left({\xi }_i\right)+\sum _{j=1}^2{b}_{1j}v\left({\eta }_j\right),\\ \hfill v\left(0\right)=0,v\left(1\right)=\sum _{i=1}^2{\lambda }_{2i}{I}^{l_{2i}}u\left({\xi }_i\right)+\sum _{j=1}^2{b}_{2j}u\left({\eta }_j\right),\end{array}\right.$$

(126)

where ${\lambda }_{11}=0.45$, ${\lambda }_{21}=0.37$, $l_{11}=1.1$, $l_{21}=1.5$, ${\xi }_1=0.32$, $b_{11}=0.47$, $b_{21}=0.19$, ${\eta }_1=0.65$, ${\lambda }_{12}=0.55$, ${\lambda }_{22}=0.25$, $l_{12}=1.7$, $l_{22}=1.6$, ${\xi }_2=0.75$, $b_{12}=0.36$, $b_{22}=0.28$, ${\eta }_2=0.45$, the nonlinear functions $f_1$ and $f_2$ are defined by

$$\begin{array}{c}{f}_1=0.7t^2+0.09tu\left(t\right)+0.01t^{2}v\left(t\right)+0.08t{\hspace{0.17em}}^H{D}^{0.1,0.5}u\left(t\right)+0.06t^2{\hspace{0.17em}}^H{D}^{0.1,0.5}v\left(t\right),\hfill \\ f_2=0.5t+0.04t^{2}u\left(t\right)+0.07tv\left(t\right)+0.06t^2{\hspace{0.17em}}^H{D}^{0.1,0.5}u\left(t\right)+0.09t{\hspace{0.17em}}^H{D}^{0.1,0.5}v\left(t\right),\hfill \end{array}$$

(127)

By calculation we get $l_1=0.5468591,$ $l_2=0.2877841,$ $M_1=1.1425206,$ $M_2=0.5024695,$ $M_3=1.1921522,$ $M_4=0.5224447,$ $M_5=0.2606166,$ $M_6=1.1446027,$ $M_7=0.2709772,$ and $M_8=1.2174674.$

From (127), we can get

$$\begin{array}{c}|f_1(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right))|\le 0.7+0.09u\left(t\right)+0.01v\left(t\right)+0.08{\hspace{0.17em}}^H{D}^{0.1,0.5}u\left(t\right)+0.06{\hspace{0.17em}}^H{D}^{0.1,0.5}v\left(t\right),\hfill \\ |f_2(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right))|\le 0.5+0.04u\left(t\right)+0.07v\left(t\right)+0.06{\hspace{0.17em}}^H{D}^{0.1,0.5}u\left(t\right)+0.09{\hspace{0.17em}}^H{D}^{0.1,0.5}v\left(t\right).\hfill \end{array}$$

(128)

Let $a_0=0.7,$ $a_1=0.09,$ $a_2=0.01,$ $a_3=0.08,$ $a_4=0.06,$ $b_0=0.5,$ $b_1=0.04,$ $b_2=0.07,$ $b_3=0.06$ and $b_4=0.09,$ Theorem 3.1 can be applied to problem (3.73) and (3.74). We find $(a_1+a_3)(M_1+M_3+M_5+M_7)+(b_1+b_3)(M_2+M_4+M_6+M_8)=0.8260<1$ and $(a_2+a_4)(M_1+M_3+M_5+M_7)+(b_2+b_4)(M_2+M_4+M_6+M_8)=0.7426<1$. Therefore, the conclusion of Theorem 3.1 implies that problem (3.73) and (3.74) have at least one solution $(u,v)$ on J.

Example 2. 

For $t\in J$, consider the following fractional differential inclusions:

$$\left\{\begin{array}{c}{\hspace{0.17em}}^H{D}^{2,0.5}u\left(t\right)\in F(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right)),\hfill \\ {\hspace{0.17em}}^H{D}^{2,0.5}v\left(t\right)\in G(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right)),\hfill \end{array}\right.$$

(129)

with the coupled integral and discrete mixed boundary conditions:

$$\left\{\begin{array}{c}\hfill u\left(0\right)=0,u\left(1\right)=\sum _{i=1}^2{\lambda }_{1i}{I}^{l_{1i}}v\left({\xi }_i\right)+\sum _{j=1}^2{b}_{1j}v\left({\eta }_j\right),\\ \hfill v\left(0\right)=0,v\left(1\right)=\sum _{i=1}^2{\lambda }_{2i}{I}^{l_{2i}}u\left({\xi }_i\right)+\sum _{j=1}^2{b}_{2j}u\left({\eta }_j\right),\end{array}\right.$$

(130)

where ${\lambda }_{11}=0.45$, ${\lambda }_{21}=0.37$, $l_{11}=1.1$, $l_{21}=1.5$, ${\xi }_1=0.32$, $b_{11}=0.47$, $b_{21}=0.19$, ${\eta }_1=0.65$, ${\lambda }_{12}=0.55$, ${\lambda }_{22}=0.25$, $l_{12}=1.7$, $l_{22}=1.6$, ${\xi }_2=0.75$, $b_{12}=0.36$, $b_{22}=0.28$, ${\eta }_2=0.45$, and $F,G$: $J\times X^4\to \mathcal{P}\left(X\right),$ $i=1,2,$ $F,G:J\times X^4\to \mathcal{P}\left(X\right)$ are multivalued maps given by

$$\begin{array}{c}F(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right))=\hfill \\ \left\{f\in X:0\le f\le 0.7t^2+0.09tu\left(t\right)+0.01t^{2}v\left(t\right)+0.08t{\hspace{0.17em}}^H{D}^{0.1,0.5}u\left(t\right)+0.06t^2{\hspace{0.17em}}^H{D}^{0.1,0.5}v\left(t\right)\right\},\hfill \\ G(t,u\left(t\right),v\left(t\right){,}^H{D}^{{\gamma }_1,\beta }u\left(t\right){,}^H{D}^{{\gamma }_2,\beta }v\left(t\right))=\hfill \\ \left\{g\in X:0\le g\le 0.5t+0.04t^{2}u\left(t\right)+0.07tv\left(t\right)+0.06t^2{\hspace{0.17em}}^H{D}^{0.1,0.5}u\left(t\right)+0.09t{\hspace{0.17em}}^H{D}^{0.1,0.5}v\left(t\right)\right\}.\hfill \end{array}$$

(131)

From (131), we know that $F,G$ are $L^1$-Carathéodory and have convex values satisfying

$$\begin{array}{c}\parallel F(t,u,v,w,z)\parallel =sup\left\{\left|f\right|:f\in F(t,u,v,w,z)\right\}\le 5,foreach(t,u,v,w,z)\in J\times X^4,\hfill \\ \parallel G(t,u,v,w,z)\parallel =sup\left\{\left|g\right|:g\in G(t,u,v,w,z)\right\}\le 6,foreach(t,u,v,w,z)\in J\times X^4,\hfill \end{array}$$

(132)

with $m_1\left(t\right)=m_2\left(t\right)={\varphi }_1\left(\right|u\left|\right)={\varphi }_2\left(\right|u\left|\right)={\psi }_1\left(\right|v\left|\right)={\psi }_2\left(\right|v\left|\right)={\varrho }_1\left(\right|w\left|\right)\equiv 1,\hspace{0.17em}\hspace{0.17em}{\rho }_1\left(\right|z\left|\right)={\varrho }_2\left(\right|w\left|\right)={\rho }_2\left(\right|z\left|\right)\equiv 2.$

Let K be any number satisfying

$$\begin{array}{cc}\hfill K& >(M_1+M_3+M_5+M_7)∥m_1∥({\varphi }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\psi }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\varrho }_1(\parallel (u,v){\parallel }_{\mathcal{C}})+{\rho }_1(\parallel (u,v){\parallel }_{\mathcal{C}}\left)\right)\hfill \\ & +(M_2+M_4+M_6+M_8)∥m_2∥({\varphi }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\psi }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\varrho }_2(\parallel (u,v){\parallel }_{\mathcal{C}})+{\rho }_2(\parallel (u,v){\parallel }_{\mathcal{C}}\left)\right)\hfill \\ & =34.6532.\hfill \end{array}$$

(133)

Clearly, all the conditions of Theorem 2 are satisfied. So, there exists at least one solution to problem (129) and (130) on J.

4. Conclusions

In this paper, we study the nonlinear coupled system of Hilfer fractional differential equations and inclusions with multi-strip and multi-point mixed boundary conditions. The existence results can be derived using tools such as the Leray-Schauder alternative, the Arzelá-Ascoli theorem, etc. It is very significant to study Hilfer fractional differential equations and inclusions for biological models and physical phenomena. However, for the Hilfer fractional differential inclusion, we assume a convex function. After that, we will continue studying the conclusions in the non-convex case.