Systems of Hilfer–Hadamard Fractional Differential Equations with Nonlocal Coupled Boundary Conditions

Abstract

We study the existence and uniqueness of solutions for a system of Hilfer–Hadamard fractional differential equations. These equations are subject to coupled nonlocal boundary conditions that incorporate Riemann–Stieltjes integrals and a range of Hadamard fractional derivatives. To establish our key findings, we apply various fixed point theorems, notably including the Banach contraction mapping principle, the Krasnosel’skii fixed point theorem applied to the sum of two operators, the Schaefer fixed point theorem, and the Leray–Schauder nonlinear alternative.

1. Introduction

We examine the system of fractional differential equations

$$\left\{\begin{array}{c}{}^{HH}{D}_1^{\alpha ,\beta }u\left(t\right)=f(t,u\left(t\right),v\left(t\right)),t\in [1,T],\hfill \\ {}^{HH}{D}_1^{\gamma ,\delta }v\left(t\right)=g(t,u\left(t\right),v\left(t\right)),t\in [1,T],\hfill \end{array}\right.$$

(1)

subject to the nonlocal coupled boundary conditions

$$\left\{\begin{array}{c}{\displaystyle u\left(1\right)=0,{}^H{D}_1^{\varsigma }u\left(T\right)=\sum _{i=1}^m{\int }_1^T{}^H{D}_1^{{\varrho }_i}u\left(s\right)d{H}_i\left(s\right)+\sum _{i=1}^n{\int }_1^T{}^H{D}_1^{{\sigma }_i}v\left(s\right)d{K}_i\left(s\right),}\hfill \\ {\displaystyle v\left(1\right)=0,{}^H{D}_1^{\vartheta }v\left(T\right)=\sum _{i=1}^p{\int }_1^T{}^H{D}_1^{{\eta }_i}u\left(s\right)d{P}_i\left(s\right)+\sum _{i=1}^q{\int }_1^T{}^H{D}_1^{{\theta }_i}v\left(s\right)d{Q}_i\left(s\right),}\hfill \end{array}\right.$$

(2)

where $T>1$, $\alpha ,\gamma \in (1,2]$, $\beta ,\delta \in [0,1]$, $m,n,p,q\in \mathbb{N}$, $\varsigma ,\vartheta ,{\varrho }_i,{\sigma }_i,{\eta }_i,{\theta }_i\in [0,1]$, ${}^{HH}{D}_1^{{\kappa }_1,{\kappa }_2}$ denotes the Hilfer–Hadamard fractional derivative of order ${\kappa }_1$ and type ${\kappa }_2$ (for ${\kappa }_1=\alpha ,\gamma$ and ${\kappa }_2=\beta ,\delta$), ${}^H{D}_1^{\kappa }$ represents the Hadamard fractional derivative of order $\kappa$ (for $\kappa =\varsigma ,\vartheta ,{\varrho }_i,{\sigma }_j,{\eta }_k,{\theta }_{\iota },i=1,\dots ,m,$ $j=1,\dots ,n$, $k=1,\dots ,p$ and $\iota =1,\dots ,q$), the continuous functions f and g are defined on $[1,T]\times {\mathbb{R}}^2$, and the integrals from the boundary conditions (2) are Riemann–Stieltjes integrals with $H_i,K_j,P_k,Q_{\iota }$, $i=1,\dots ,m,$ $j=1,\dots ,n$, $k=1,\dots ,p$ and $\iota =1,\dots ,q$ functions of bounded variation.

In this paper, we present a variety of conditions for the functions f and g such that problem (1) and (2) has at least one solution. We will write our problem as an equivalent system of integral equations, and then we will associate it with an operator whose fixed points are our solutions. The proof of our primary outcomes involves the utilization of diverse fixed point theorems. Noteworthy among these theorems are the Banach contraction mapping principle, the Krasnosel’skii fixed point theorem applied to the sum of two operators, the Schaefer fixed point theorem, and the Leray–Schauder nonlinear alternative. The nonlocal boundary conditions (2) are general ones, and they include different particular cases. For example, if $\kappa =0$, for $\kappa =\varsigma ,\vartheta ,{\varrho }_i,{\sigma }_j,{\eta }_k,{\theta }_{\iota },i=1,\dots ,m,$ $j=1,\dots ,n$, $k=1,\dots ,p$ and $\iota =1,\dots ,q$, then the Hadamard derivative ${}^H{D}_1^{\kappa }z\left(t\right)$ coincides with $z\left(t\right)$. If one of the order of the Hadamard derivatives from the right-hand side of the relations from (2) is zero (for example, if ${\varrho }_1$ is zero), then the term ${\int }_1^T{}^H{D}_1^{{\varrho }_1}u\left(s\right)d{H}_1\left(s\right)$ becomes ${\int }_1^{T}u\left(s\right)d{H}_1\left(s\right)$, which contains the cases of the multi-point boundary conditions for the function u (if $H_1$ is a step function); a classical integral condition; a combination of them; or even a Hadamard fractional integral for a special form of $H_1$ (as we mentioned in [1]). If ${\varrho }_1\in (0,1]$ and $H_1$ is a step function, then ${\int }_1^T{}^H{D}_1^{{\varrho }_1}u\left(s\right)d{H}_1\left(s\right)={\sum }_{i=1}^{n_0}{}^H{D}_1^{{\varrho }_1}u\left({\xi }_i\right)$, which is a combination of the Hadamard fractional derivatives of function u in various points. If all functions $K_i,i=1,\dots ,n$ and $P_j,j=1,\dots ,p$ are constant functions, then the boundary conditions become uncoupled boundary conditions (where the Hadamard derivative of order $\varsigma$ of the function u in the point T is dependent only of the derivatives ${}^H{D}_1^{{\varrho }_i},i=1,\dots ,m$ of the function u, and the Hadamard derivative of order $\vartheta$ of the function v in the point T is dependent only of the derivatives ${}^H{D}_1^{{\theta }_i},i=1,\dots ,q$ of function v), and if $H_i,i=1,\dots ,m$ and $Q_j,j=1,\dots ,q$ are constant functions, then the boundary conditions become purely coupled boundary conditions (in which the Hadamard derivative of order $\varsigma$ of the function u in T is dependent only of the derivatives ${}^H{D}_1^{{\sigma }_i},i=1,\dots ,n$ of the function v, and the Hadamard derivative of order $\vartheta$ of the function v in T is dependent only of the derivatives ${}^H{D}_1^{{\eta }_i},i=1,\dots ,p$ of the function u).

Next, we will introduce some papers that are relevant to the issue posed by Equations (1) and (2). In [2], the authors investigated the existence and uniqueness of solutions for the Hilfer–Hadamard fractional differential equation with nonlocal boundary conditions

$$\left\{\begin{array}{c}{}^{HH}{D}_1^{\alpha ,\beta }x\left(t\right)=f(t,x\left(t\right)),t\in [1,T],\hfill \\ {\displaystyle x\left(1\right)=0,x\left(T\right)=\sum _{j=1}^m{\eta }_{j}x\left({\xi }_j\right)+\sum _{i=1}^n{\zeta }_i{}^H{I}_1^{{\varphi }_i}x\left({\theta }_i\right)+\sum _{k=1}^r{\lambda }_k{}^H{D}_1^{{\omega }_k}x\left({\mu }_k\right),}\hfill \end{array}\right.$$

(3)

where $\alpha \in (1,2]$, $\beta \in [0,1]$, ${\eta }_j,{\zeta }_i,{\lambda }_k\in \mathbb{R}$, $f:[1,T]\times \mathbb{R}\to \mathbb{R}$ is a continuous function, ${}^H{I}^{{\varphi }_i}$ is the Hadamard fractional integral operator of order ${\varphi }_i>0$, and ${\xi }_j,{\theta }_i,{\mu }_k\in (1,T)$ for $j=1,\dots ,m$, $i=1,\dots ,n$, $k=1,\dots ,r$. The multi-valued version of problem (3) is also studied. For the proof of the main results, they used differing fixed point theorems. In [3], the authors proved the existence of solutions for the system of sequential Hilfer–Hadamard fractional differential equations supplemented with boundary conditions

$$\left\{\begin{array}{c}{(}^{HH}{D}_1^{{\alpha }_1,{\beta }_1}+{\lambda }_1{}^{HH}{D}_1^{{\alpha }_1-1,{\beta }_1})u\left(t\right)=f(t,u\left(t\right),v\left(t\right)),t\in [1,e],\hfill \\ {(}^{HH}{D}_1^{{\alpha }_2,{\beta }_2}+{\lambda }_2{}^{HH}{D}_1^{{\alpha }_2-1,{\beta }_2})v\left(t\right)=g(t,u\left(t\right),v\left(t\right)),t\in [1,e],\hfill \\ u\left(1\right)=0,u\left(e\right)=A_1,v\left(1\right)=0,v\left(e\right)=A_2,\hfill \end{array}\right.$$

(4)

where ${\alpha }_1,{\alpha }_2\in (1,2]$, $A_1,A_2,{\lambda }_1,{\lambda }_2\in {\mathbb{R}}_{+}$, and $f,g:[1,e]\times {\mathbb{R}}^2\to \mathbb{R}$ are given continuous functions.

In paper [4], Hadamard defined a fractional derivative with a kernel involving a logarithmic function with an arbitrary exponent. In [5], Hilfer introduced a new fractional derivative (known as the Hilfer fractional derivative), which is a generalization of the Riemann-Liouville fractional derivative and the Caputo fractional derivative. Some applications of this new fractional derivative are presented in papers [6,7]. The Hilfer–Hadamard fractional derivative is an interpolation of the Hadamard fractional derivative, and it covers the cases of the Riemann–Liouville–Hadamard and Caputo–Hadamard fractional derivatives (see the definition in Section 2). The distinctive aspects of our presented challenge, (1) and (2), emerge from the exploration of a set of Hilfer–Hadamard fractional differential equations encompassing diverse orders and types. Additionally, the introduction of general nonlocal boundary conditions (2) contributes novelty, extending beyond numerous specific instances as previously observed. To the best of our knowledge, this issue represented by Equations (1) and (2) is a novel problem in the literature. Our theorems represent original contributions and make substantial advancements in the realm of coupled systems involving Hilfer–Hadamard fractional derivatives. Although the techniques employed in demonstrating our primary findings in Section 3 are conventional, their adaptation to address our problem (1) and (2) is innovative. For more recent investigations concerning Hadamard, Hilfer, and Hilfer–Hadamard fractional differential equations and their applications, we recommend the monograph [8] and the following papers: [1,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28].

The structure of the paper unfolds as follows: In Section 2, we offer definitions and properties related to fractional derivatives, along with a result regarding the existence of solutions for the linear boundary value problem linked to Equations (1) and (2). Moving on, Section 3 is dedicated to the core findings concerning the existence and uniqueness of solutions for problem (1) and (2). Subsequently, in Section 4, we provide illustrative examples that demonstrate the practical application of our theorems. Lastly, concluding insights for this paper can be found in Section 5.

2. Auxiliary Results

In this section, we present some definitions and properties of fractional derivatives and an existence result for the linear boundary value problem associated with (1) and (2).

Definition 1

(Hadamard fractional integral [29]). For a function $z:[a,\infty )\to \mathbb{R}$, ($a\ge 0$), the Hadamard fractional integral of order $p>0$ is defined by

$${\displaystyle {(}^H{I}_a^{p}z)\left(x\right)=\frac{1}{\mathsf{\Gamma }\left(p\right)}{\int }_a^x{\left({\displaystyle ln\frac{x}{t}}\right)}^{p-1}\frac{z\left(t\right)}{t}dt,x>a,}$$

(5)

and ${(}^H{I}_a^{0}z)\left(x\right)=z\left(x\right),x>a$.

Definition 2

(Hadamard fractional derivative [29]). For a function $z:[a,\infty )\to \mathbb{R}$, ($a\ge 0$), the Hadamard fractional derivative of order $p>0$ is defined by

$${\displaystyle {}^H{D}_a^{p}z\left(x\right)={\left({\displaystyle x\frac{d}{dx}}\right)}^n{}^H{I}_a^{n-p}z\left(x\right)=\frac{1}{\mathsf{\Gamma }(n-p)}{\left({\displaystyle x\frac{d}{dx}}\right)}^n{\int }_a^x{\left({\displaystyle ln\frac{x}{t}}\right)}^{n-p-1}\frac{z\left(t\right)}{t}dt,}$$

(6)

where $n-1<p<n$, ($n\in \mathbb{N}$). For $p=m\in \mathbb{N}$, ${}^H{D}_a^{m}z\left(x\right)=\left({\delta }^{m}z\right)\left(x\right),x>a$, where $\delta =x\frac{d}{dx}$ is the δ-derivative, and, for $p=0$, ${}^H{D}_a^{0}z\left(x\right)=z\left(x\right)$.

Lemma 1

([29]). If $\alpha ,\beta >0$, and $a>0$, then

$$\begin{array}{c}{\displaystyle \left({}^H{I}_a^{\alpha }{\left({\displaystyle ln\frac{t}{a}}\right)}^{\beta -1}\right)\left(x\right)=\frac{\mathsf{\Gamma }\left(\beta \right)}{\mathsf{\Gamma }(\beta +\alpha )}{\left({\displaystyle ln\frac{x}{a}}\right)}^{\beta +\alpha -1},}\hfill \\ {\displaystyle \left({}^H{D}_a^{\alpha }{\left({\displaystyle ln\frac{t}{a}}\right)}^{\beta -1}\right)\left(x\right)=\frac{\mathsf{\Gamma }\left(\beta \right)}{\mathsf{\Gamma }(\beta -\alpha )}{\left({\displaystyle ln\frac{x}{a}}\right)}^{\beta -\alpha -1}.}\hfill \end{array}$$

(7)

Definition 3

(Hilfer–Hadamard fractional derivative [2,7]). Let $z\in L^1(a,b)$ and $n-1<\alpha \le n$, ($n\in \mathbb{N}$), $0\le \beta \le 1$. The Hilfer–Hadamard fractional derivative of order α and type β for the function z is defined by

$$\begin{array}{c}\left({}^{HH}{D}_a^{\alpha ,\beta }z\right)\left(t\right)=\left({}^H{I}_a^{\beta (n-\alpha )}{\delta }^n{}^H{I}_a^{(n-\alpha )(1-\beta )}z\right)\left(t\right)\hfill \\ =\left({}^H{I}_a^{\beta (n-\alpha )}{\delta }^n{}^H{I}_a^{(n-\gamma )}z\right)\left(t\right)=\left({}^H{I}_a^{\beta (n-\alpha )}{}^H{D}_a^{\gamma }z\right)\left(t\right),\hfill \end{array}$$

(8)

where $\gamma =\alpha +n\beta -\alpha \beta$.

If $\beta =0$, the Hilfer–Hadamard fractional derivative ${}^{HH}{D}_a^{\alpha ,\beta }z$ coincides with the Hadamard fractional derivative ${}^H{D}_a^{\alpha }z$. If $\beta =1$, the fractional derivative ${}^{HH}{D}_a^{\alpha ,\beta }z$ coincides with the Caputo–Hadamard derivative, given by ${}^{CH}{D}_a^{\alpha }z\left(t\right)={(}^H{I}_a^{n-\alpha }{\delta }^{n}z)\left(t\right)$.

Theorem 1

([2]). Let $\alpha >0$, $n-1<\alpha \le n$, ($n\in \mathbb{N}$), $\beta \in [0,1]$, $\gamma =\alpha +n\beta -\alpha \beta$, and $0<a<b<\infty$. If $z\in L^1(a,b)$ and ${(}^H{I}_a^{n-\gamma }z)\left(t\right)\in A{C}_{\delta }^n[a,b]$, then the following relation holds

$$\begin{array}{c}{\displaystyle {}^H{I}_a^{\alpha }\left({}^{HH}{D}_a^{\alpha ,\beta }z\right)\left(t\right){=}^H{I}_a^{\gamma }\left({}^H{D}_a^{\gamma }z\right)\left(t\right)=z\left(t\right)-\sum _{i=0}^{n-1}\frac{\left({\delta }^{(n-i-1)}{(}^H{I}_a^{n-\gamma }z)\right)\left(a\right)}{\mathsf{\Gamma }(\gamma -i)}{\left({\displaystyle ln\frac{t}{a}}\right)}^{\gamma -i-1}.}\hfill \end{array}$$

(9)

We consider now the system of linear fractional differential equations

$$\left\{\begin{array}{c}{}^{HH}{D}_1^{\alpha ,\beta }u\left(t\right)=h\left(t\right),t\in [1,T],\hfill \\ {}^{HH}{D}_1^{\gamma ,\delta }v\left(t\right)=k\left(t\right),t\in [1,T],\hfill \end{array}\right.$$

(10)

subject to the boundary conditions (2), where $h,k\in C\left(\right[1,T],\mathbb{R})$. We denote by $\lambda =\alpha +(2-\alpha )\beta$, $\mu =\gamma +(2-\gamma )\delta$, and

$$\begin{array}{c}{\displaystyle a=\frac{\mathsf{\Gamma }\left(\lambda \right)}{\mathsf{\Gamma }(\lambda -\varsigma )}{(lnT)}^{\lambda -\varsigma -1}-\sum _{i=1}^m\frac{\mathsf{\Gamma }\left(\lambda \right)}{\mathsf{\Gamma }(\lambda -{\varrho }_i)}{\int }_1^T{(lns)}^{\lambda -{\varrho }_i-1}d{H}_i\left(s\right),}\hfill \\ {\displaystyle b=\sum _{i=1}^n\frac{\mathsf{\Gamma }\left(\mu \right)}{\mathsf{\Gamma }(\mu -{\sigma }_i)}{\int }_1^T{(lns)}^{\mu -{\sigma }_i-1}d{K}_i\left(s\right),}\hfill \\ {\displaystyle c=\sum _{i=1}^p\frac{\mathsf{\Gamma }\left(\lambda \right)}{\mathsf{\Gamma }(\lambda -{\eta }_i)}{\int }_1^T{(lns)}^{\lambda -{\eta }_i-1}d{P}_i\left(s\right),}\hfill \\ {\displaystyle d=\frac{\mathsf{\Gamma }\left(\mu \right)}{\mathsf{\Gamma }(\mu -\vartheta )}{(lnT)}^{\mu -\vartheta -1}-\sum _{i=1}^q\frac{\mathsf{\Gamma }\left(\mu \right)}{\mathsf{\Gamma }(\mu -{\theta }_i)}{\int }_1^T{(lns)}^{\mu -{\theta }_i-1}d{Q}_i\left(s\right),}\hfill \\ \Delta =ad-bc.\hfill \end{array}$$

(11)

Lemma 2.

We suppose that $a,b,c,d\in \mathbb{R}$, $\Delta \ne 0$, and $h,k\in C\left(\right[1,T],\mathbb{R})$. Then, the solution of problem (10) and (2) is given by

$$\begin{array}{c}{\displaystyle u\left(t\right){=}^H{I}_1^{\alpha }h\left(t\right)+\frac{{(lnt)}^{\lambda -1}}{\Delta }\left[-d{}^H{I}_1^{\alpha -\varsigma }h\left(T\right)+d(A_1\left(h\right)+A_2\left(k\right))\right.}\hfill \\ \left.-b{}^H{I}_1^{\gamma -\vartheta }k\left(T\right)+b(A_3\left(h\right)+A_4\left(k\right))\right],t\in [1,T],\hfill \\ {\displaystyle v\left(t\right){=}^H{I}_1^{\gamma }k\left(t\right)+\frac{{(lnt)}^{\mu -1}}{\Delta }\left[-a{}^H{I}_1^{\gamma -\vartheta }k\left(T\right)+a(A_3\left(h\right)+A_4\left(k\right))\right.}\hfill \\ \left.-c{}^H{I}_1^{\alpha -\varsigma }h\left(T\right)+c(A_1\left(h\right)+A_2\left(k\right))\right],t\in [1,T],\hfill \end{array}$$

(12)

where operators $A_i:C([1,T],\mathbb{R})\to \mathbb{R},i=1,\dots ,4$ are defined by

$$\begin{array}{c}{\displaystyle A_1\left(h\right)=\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}{\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}\frac{h\left(\tau \right)}{\tau }d\tau }\right)d{H}_i\left(s\right),}\hfill \\ {\displaystyle A_2\left(k\right)=\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}{\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}\frac{k\left(\tau \right)}{\tau }d\tau }\right)d{K}_i\left(s\right),}\hfill \\ {\displaystyle A_3\left(h\right)=\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}{\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}\frac{h\left(\tau \right)}{\tau }d\tau }\right)d{P}_i\left(s\right),}\hfill \\ {\displaystyle A_4\left(k\right)=\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}{\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}\frac{k\left(\tau \right)}{\tau }d\tau }\right)d{Q}_i\left(s\right).}\hfill \end{array}$$

(13)

Proof. 

We apply the integral operators ${}^H{I}_1^{\alpha }$ and ${}^H{I}_1^{\gamma }$, respectively, to equations of system (10). Then, the solutions of system (10) are given by

$$\left\{\begin{array}{c}u\left(t\right)=a_1{(lnt)}^{\lambda -1}+a_2{(lnt)}^{\lambda -2}{+}^H{I}_1^{\alpha }h\left(t\right),t\in [1,T],\hfill \\ v\left(t\right)=b_1{(lnt)}^{\mu -1}+b_2{(lnt)}^{\mu -2}{+}^H{I}_1^{\gamma }k\left(t\right),t\in [1,T],\hfill \end{array}\right.$$

(14)

where $a_i,b_i\in \mathbb{R},i=1,2$. Because $u\left(1\right)=v\left(1\right)=0$, we deduce that $a_2=b_2=0$. So, we obtain, for the solutions of (10), the formulas

$$\left\{\begin{array}{c}u\left(t\right)=a_1{(lnt)}^{\lambda -1}{+}^H{I}_1^{\alpha }h\left(t\right),t\in [1,T],\hfill \\ v\left(t\right)=b_1{(lnt)}^{\mu -1}{+}^H{I}_1^{\gamma }k\left(t\right),t\in [1,T].\hfill \end{array}\right.$$

(15)

For $\kappa =\varsigma ,{\varrho }_i,{\eta }_j$, $i=1,\dots ,m,j=1,\dots ,p$, we find

$$\begin{array}{c}{\displaystyle {}^H{D}_1^{\kappa }u\left(t\right)=a_1{}^H{D}_1^{\kappa }{(lnt)}^{\lambda -1}{+}^H{D}_1^{\kappa }{}^H{I}_1^{\alpha }h\left(t\right)=a_1\frac{\mathsf{\Gamma }\left(\lambda \right)}{\mathsf{\Gamma }(\lambda -\kappa )}{(lnt)}^{\lambda -\kappa -1}{+}^H{I}_1^{\alpha -\kappa }h\left(t\right)}\hfill \\ {\displaystyle =a_1\frac{\mathsf{\Gamma }\left(\lambda \right)}{\mathsf{\Gamma }(\lambda -\kappa )}{(lnt)}^{\lambda -\kappa -1}+\frac{1}{\mathsf{\Gamma }(\alpha -\kappa )}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\alpha -\kappa -1}\frac{h\left(s\right)}{s}ds,}\hfill \end{array}$$

(16)

and for $\tilde{\kappa }=\vartheta ,{\sigma }_i,{\theta }_j$, $i=1,\dots ,n,j=1,\dots ,q$, we obtain

$$\begin{array}{c}{\displaystyle {}^H{D}_1^{\tilde{\kappa }}v\left(t\right)=b_1{}^H{D}_1^{\tilde{\kappa }}{(lnt)}^{\mu -1}{+}^H{D}_1^{\tilde{\kappa }}{}^H{I}_1^{\gamma }k\left(t\right)=b_1\frac{\mathsf{\Gamma }\left(\mu \right)}{\mathsf{\Gamma }(\mu -\tilde{\kappa })}{(lnt)}^{\mu -\tilde{\kappa }-1}{+}^H{I}_1^{\gamma -\tilde{\kappa }}k\left(t\right)}\hfill \\ {\displaystyle =b_1\frac{\mathsf{\Gamma }\left(\mu \right)}{\mathsf{\Gamma }(\mu -\tilde{\kappa })}{(lnt)}^{\mu -\tilde{\kappa }-1}+\frac{1}{\mathsf{\Gamma }(\gamma -\tilde{\kappa })}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\gamma -\tilde{\kappa }-1}\frac{k\left(s\right)}{s}ds.}\hfill \end{array}$$

(17)

By applying the conditions ${}^H{D}_1^{\varsigma }u\left(T\right)={\sum }_{i=1}^m{\int }_1^T{}^H{D}_1^{{\varrho }_i}u\left(s\right)d{H}_i\left(s\right)+{\sum }_{i=1}^n{\int }_1^T{}^H{D}_1^{{\sigma }_i}v\left(s\right)d{K}_i\left(s\right)$, and ${}^H{D}_1^{\vartheta }v\left(T\right)={\sum }_{i=1}^p{\int }_1^T{}^H{D}_1^{{\eta }_i}u\left(s\right)d{P}_i\left(s\right)+{\sum }_{i=1}^q{\int }_1^T{}^H{D}_1^{{\theta }_i}v\left(s\right)d{Q}_i\left(s\right)$, we deduce

$$\begin{array}{c}{\displaystyle a_1\frac{\mathsf{\Gamma }\left(\lambda \right)}{\mathsf{\Gamma }(\lambda -\varsigma )}{(lnT)}^{\lambda -\varsigma -1}+\frac{1}{\mathsf{\Gamma }(\alpha -\varsigma )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}\frac{h\left(s\right)}{s}ds}\hfill \\ {\displaystyle =\sum _{i=1}^m{\int }_1^T\left[{\displaystyle a_1\frac{\mathsf{\Gamma }\left(\lambda \right)}{\mathsf{\Gamma }(\lambda -{\varrho }_i)}{(lns)}^{\lambda -{\varrho }_i-1}+\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}{\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}\frac{h\left(\tau \right)}{\tau }d\tau }\right]d{H}_i\left(s\right)}\hfill \\ {\displaystyle +\sum _{i=1}^n{\int }_1^T\left[{\displaystyle b_1\frac{\mathsf{\Gamma }\left(\mu \right)}{\mathsf{\Gamma }(\mu -{\sigma }_i)}{(lns)}^{\mu -{\sigma }_i-1}+\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}{\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}\frac{k\left(\tau \right)}{\tau }d\tau }\right]d{K}_i\left(s\right),}\hfill \\ {\displaystyle b_1\frac{\mathsf{\Gamma }\left(\mu \right)}{\mathsf{\Gamma }(\mu -\vartheta )}{(lnT)}^{\mu -\vartheta -1}+\frac{1}{\mathsf{\Gamma }(\gamma -\vartheta )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}\frac{k\left(s\right)}{s}ds}\hfill \\ {\displaystyle =\sum _{i=1}^p{\int }_1^T\left[{\displaystyle a_1\frac{\mathsf{\Gamma }\left(\lambda \right)}{\mathsf{\Gamma }(\lambda -{\eta }_i)}{(lns)}^{\lambda -{\eta }_i-1}+\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}{\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}\frac{h\left(\tau \right)}{\tau }d\tau }\right]d{P}_i\left(s\right)}\hfill \\ {\displaystyle +\sum _{i=1}^q{\int }_1^T\left[{\displaystyle b_1\frac{\mathsf{\Gamma }\left(\mu \right)}{\mathsf{\Gamma }(\mu -{\theta }_i)}{(lns)}^{\mu -{\theta }_i-1}+\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}{\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}\frac{k\left(\tau \right)}{\tau }d\tau }\right]d{Q}_i\left(s\right),}\hfill \end{array}$$

(18)

or

$$\begin{array}{c}{a}_1\left[{\displaystyle \frac{\mathsf{\Gamma }\left(\lambda \right)}{\mathsf{\Gamma }(\lambda -\varsigma )}{(lnT)}^{\lambda -\varsigma -1}-\sum _{i=1}^m\frac{\mathsf{\Gamma }\left(\lambda \right)}{\mathsf{\Gamma }(\lambda -{\varrho }_i)}{\int }_1^T{(lns)}^{\lambda -{\varrho }_i-1}d{H}_i\left(s\right)}\right]\hfill \\ {\displaystyle -b_1\sum _{i=1}^n\frac{\mathsf{\Gamma }\left(\mu \right)}{\mathsf{\Gamma }(\mu -{\sigma }_i)}{\int }_1^T{(lns)}^{\mu -{\sigma }_i-1}d{K}_i\left(s\right)}\hfill \\ ={-}^H{I}_1^{\alpha -\varsigma }h\left(T\right)+A_1\left(h\right)+A_2\left(k\right),\hfill \\ {\displaystyle -a_1\sum _{i=1}^p\frac{\mathsf{\Gamma }\left(\lambda \right)}{\mathsf{\Gamma }(\lambda -{\eta }_i)}{\int }_1^T{(lns)}^{\lambda -{\eta }_i-1}d{P}_i\left(s\right)}\hfill \\ +b_1\left[{\displaystyle \frac{\mathsf{\Gamma }\left(\mu \right)}{\mathsf{\Gamma }(\mu -\vartheta )}{(lnT)}^{\mu -\vartheta -1}-\sum _{i=1}^q\frac{\mathsf{\Gamma }\left(\mu \right)}{\mathsf{\Gamma }(\mu -{\theta }_i)}{\int }_1^T{(lns)}^{\mu -{\theta }_i-1}d{Q}_i\left(s\right)}\right]\hfill \\ ={-}^H{I}_1^{\gamma -\vartheta }k\left(T\right)+A_3\left(h\right)+A_4\left(k\right).\hfill \end{array}$$

(19)

The determinant of system (19) in the unknowns $a_1$ and $b_1$ is

$$\left|\begin{array}{cc}a& -b\\ -c& d\end{array}\right|=ad-bc,$$

(20)

that is, $\Delta$, given by (11), which is different than zero by the assumptions of this lemma. So, the solution of system (19) is unique, namely

$$\begin{array}{c}{\displaystyle a_1=-\frac{d}{\Delta }{}^H{I}_1^{\alpha -\varsigma }h\left(T\right)+\frac{d}{\Delta }(A_1\left(h\right)+A_2\left(k\right))-\frac{b}{\Delta }{}^H{I}_1^{\gamma -\vartheta }k\left(T\right)+\frac{b}{\Delta }(A_3\left(h\right)+A_4\left(k\right)),}\hfill \\ {\displaystyle b_1=-\frac{a}{\Delta }{}^H{I}_1^{\gamma -\vartheta }k\left(T\right)+\frac{a}{\Delta }(A_3\left(h\right)+A_4\left(k\right))-\frac{c}{\Delta }{}^H{I}_1^{\alpha -\varsigma }h\left(T\right)+\frac{c}{\Delta }(A_1\left(h\right)+A_2\left(k\right)).}\hfill \end{array}$$

(21)

By replacing the above formulas for $a_1$ and $b_1$ in (15), we obtain the solution of problems (10) and (2) given by (12). □

3. Existence Results

In this section, we will give the main existence and uniqueness theorems for the solutions of problem (1) and (2). By using Lemma 2, our problem (1) and (2) can be equivalently written as the following system of integral equations

$$\left\{\begin{array}{c}{\displaystyle u\left(t\right){=}^H{I}_1^{\alpha }{F}_{uv}\left(t\right)+\frac{{(lnt)}^{\lambda -1}}{\Delta }\left[-d{}^H{I}_1^{\alpha -\varsigma }{F}_{uv}\left(T\right)+d(A_1\left(F_{uv}\right)+A_2\left(G_{uv}\right))\right.}\hfill \\ \left.-b{}^H{I}_1^{\gamma -\vartheta }{G}_{uv}\left(T\right)+b(A_3\left(F_{uv}\right)+A_4\left(G_{uv}\right))\right],t\in [1,T],\hfill \\ {\displaystyle v\left(t\right){=}^H{I}_1^{\gamma }{G}_{uv}\left(t\right)+\frac{{(lnt)}^{\mu -1}}{\Delta }\left[-a{}^H{I}_1^{\gamma -\vartheta }{G}_{uv}\left(T\right)+a(A_3\left(F_{uv}\right)+A_4\left(G_{uv}\right))\right.}\hfill \\ \left.-c{}^H{I}_1^{\alpha -\varsigma }{F}_{uv}\left(T\right)+c(A_1\left(F_{uv}\right)+A_2\left(G_{uv}\right))\right],t\in [1,T],\hfill \end{array}\right.$$

(22)

where $F_{uv}\left(\tau \right)=f(\tau ,u\left(\tau \right),v\left(\tau \right)),G_{uv}\left(\tau \right)=g(\tau ,u\left(\tau \right),v\left(\tau \right)),$ $\tau \in [1,T]$.

We consider the Banach $\mathcal{X}=C\left(\right[1,T],\mathbb{R})$ with the supremum norm $\parallel u\parallel ={sup}_{t\in [1,T]}\left|u\left(t\right)\right|$ and the Banach space $\mathcal{Y}=\mathcal{X}\times \mathcal{X}$ with the norm ${\parallel (u,v)\parallel }_{\mathcal{Y}}=\parallel u\parallel +\parallel v\parallel$. We define the operator $\mathcal{A}:\mathcal{Y}\to \mathcal{Y}$, $\mathcal{A}(u,v)=({\mathcal{A}}_1(u,v),{\mathcal{A}}_2(u,v))$, with ${\mathcal{A}}_1,{\mathcal{A}}_2:\mathcal{Y}\to \mathcal{X}$ given by

$$\begin{array}{c}{\displaystyle {\mathcal{A}}_1(u,v)\left(t\right){=}^H{I}_1^{\alpha }{F}_{uv}\left(t\right)+\frac{{(lnt)}^{\lambda -1}}{\Delta }\left[-d{}^H{I}_1^{\alpha -\varsigma }{F}_{uv}\left(T\right)+d(A_1\left(F_{uv}\right)+A_2\left(G_{uv}\right))\right.}\hfill \\ \left.-b{}^H{I}_1^{\gamma -\vartheta }{G}_{uv}\left(T\right)+b(A_3\left(F_{uv}\right)+A_4\left(G_{uv}\right))\right],\hfill \\ {\displaystyle {\mathcal{A}}_2(u,v)\left(t\right){=}^H{I}_1^{\gamma }{G}_{uv}\left(t\right)+\frac{{(lnt)}^{\mu -1}}{\Delta }\left[-a{}^H{I}_1^{\gamma -\vartheta }{G}_{uv}\left(T\right)+a(A_3\left(F_{uv}\right)+A_4\left(G_{uv}\right))\right.}\hfill \\ \left.-c{}^H{I}_1^{\alpha -\varsigma }{F}_{uv}\left(T\right)+c(A_1\left(F_{uv}\right)+A_2\left(G_{uv}\right))\right],\hfill \end{array}$$

(23)

for all $t\in [1,T]$ and $(u,v)\in \mathcal{Y}$. We see that the solutions of problem (1) and (2) (or system (22) are the fixed points of operator $\mathcal{A}$. So, next, we will investigate the existence of the fixed points of this operator $\mathcal{A}$ in the space $\mathcal{Y}$.

We present now the basic assumptions that we will use in the next results.

(H1)

$\alpha ,\gamma \in (1,2]$; $\beta ,\delta \in [0,1]$; $m,n,p,q\in \mathbb{N}$; $\varsigma ,\vartheta ,{\varrho }_i,{\sigma }_j,{\eta }_k,{\theta }_{\iota }\in [0,1]$; $H_i,K_j,P_k,Q_{\iota }$ are bounded variation functions, for all $i=1,\dots ,m$, $j=1,\dots ,n$, $k=1,\dots ,p$, $\iota =1,\dots ,q$; $a,b,c,d\in \mathbb{R}$, and $\Delta \ne 0$ (given by (11)).

We also introduce the constants

$$\begin{array}{c}{\displaystyle {\mathsf{\Xi }}_1=\frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnT)}^{\alpha }+\frac{{(lnT)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \frac{\left|d\right|}{\mathsf{\Gamma }(\alpha -\varsigma +1)}{(lnT)}^{\alpha -\varsigma }}\right.}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\alpha -{\varrho }_i}d{H}_i\left(s\right)}\right|}\hfill \\ +\left.{\displaystyle \left|b\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\alpha -{\eta }_i}d{P}_i\left(s\right)}\right|}\right],\hfill \end{array}$$

$$\begin{array}{c}{\displaystyle {\mathsf{\Xi }}_2=\frac{{(lnT)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \left|d\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\gamma -{\sigma }_i}d{K}_i\left(s\right)}\right|}\right.}\hfill \\ +\left.{\displaystyle \frac{\left|b\right|}{\mathsf{\Gamma }(\gamma -\vartheta +1)}{(lnT)}^{\gamma -\vartheta }+\left|b\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\gamma -{\theta }_i}d{Q}_i\left(s\right)}\right|}\right],\hfill \\ {\displaystyle {\mathsf{\Xi }}_3=\frac{{(lnT)}^{\mu -1}}{|\Delta |}\left[{\displaystyle \left|a\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\alpha -{\eta }_i}d{P}_i\left(s\right)}\right|}\right.}\hfill \\ +\left.{\displaystyle \frac{\left|c\right|}{\mathsf{\Gamma }(\alpha -\varsigma +1)}{(lnT)}^{\alpha -\varsigma }+\left|c\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\alpha -{\varrho }_i}d{H}_i\left(s\right)}\right|}\right],\hfill \\ {\displaystyle {\mathsf{\Xi }}_4=\frac{1}{\mathsf{\Gamma }(\gamma +1)}{(lnT)}^{\gamma }+\frac{{(lnT)}^{\mu -1}}{|\Delta |}\left[{\displaystyle \frac{\left|a\right|}{\mathsf{\Gamma }(\gamma -\vartheta +1)}{(lnT)}^{\gamma -\vartheta }}\right.}\hfill \\ {\displaystyle +\left|a\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\gamma -{\theta }_i}d{Q}_i\left(s\right)}\right|}\hfill \\ \left.{\displaystyle +\left|c\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\gamma -{\sigma }_i}d{K}_i\left(s\right)}\right|}\right].\hfill \end{array}$$

(24)

Our first existence and uniqueness theorem for problem (1) and (2) is the following one, which is based on the Banach contraction mapping principle (see [30]).

Theorem 2.

We assume that assumption $\left(H1\right)$ holds. In addition, we suppose that the functions $f,g:[1,T]\times {\mathbb{R}}^2\to \mathbb{R}$ are continuous and satisfy the condition

(H2)

There exist $l_i>0$, $i=1,\dots ,4$ such that

$$\begin{array}{c}|f(t,x_1,y_1)-f(t,x_2,y_2)|\le l_1|x_1-x_2|+l_2|y_1-y_2|,\hfill \\ |g(t,x_1,y_1)-g(t,x_2,y_2)|\le l_3|x_1-x_2|+l_4|y_1-y_2|,\hfill \end{array}$$

(25)

for all $t\in [1,T]$ and $x_i,y_i\in \mathbb{R}$, $i=1,2$.

If

$$l_5({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)+l_6({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)<1,$$

(26)

where $l_5=max\{l_1,l_2\}$, $l_6=max\{l_3,l_4\}$, then the boundary value problem (1) and (2) has a unique solution $\left(u\right(t),v(t\left)\right),t\in [1,T].$

Proof. 

We will verify that operator $\mathcal{A}$ is a contraction in the space $\mathcal{Y}$. We denote this by ${\mathsf{\Lambda }}_1={sup}_{t\in [1,T]}\left|f(t,0,0)\right|$ and ${\mathsf{\Lambda }}_2={sup}_{t\in [1,T]}\left|g(t,0,0)\right|$. By using $\left(H2\right)$, we find

$$\begin{array}{c}|F_{uv}\left(t\right)|=|f(t,u\left(t\right),v\left(t\right))|\le |f(t,u\left(t\right),v\left(t\right))-f(t,0,0)|+|f(t,0,0)|\hfill \\ \le l_1\left|u\left(t\right)\right|+l_2\left|v\left(t\right)\right|+{\mathsf{\Lambda }}_1\le l_5\left(\right|u\left(t\right)|+\left|v\left(t\right)\right|)+{\mathsf{\Lambda }}_1,\hfill \\ |G_{uv}\left(t\right)|=|g(t,u\left(t\right),v\left(t\right))|\le |g(t,u\left(t\right),v\left(t\right))-g(t,0,0)|+|g(t,0,0)|\hfill \\ \le l_3\left|u\left(t\right)\right|+l_4\left|v\left(t\right)\right|+{\mathsf{\Lambda }}_2\le l_6\left(\right|u\left(t\right)|+\left|v\left(t\right)\right|)+{\mathsf{\Lambda }}_2,\hfill \end{array}$$

(27)

for all $t\in [1,T]$ and $(u,v)\in \mathcal{Y}$. We consider now the positive number

$${\displaystyle R\ge \frac{{\mathsf{\Lambda }}_1({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)+{\mathsf{\Lambda }}_2({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)}{1-l_5({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)-l_6({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)},}$$

(28)

and let the set $B_R={\{(u,v)\in \mathcal{Y},\parallel (u,v)\parallel }_{\mathcal{Y}}\le R\}$.

We will show firstly that $\mathcal{A}\left(B_R\right)\subset B_R$. Indeed, for this, let $(u,v)\in B_R$. Then, we obtain

$$\begin{array}{c}{\displaystyle |{\mathcal{A}}_1(u,v)\left(t\right)|\le \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\alpha -1}\left|F_{uv}\left(s\right)\right|\frac{ds}{s}}\hfill \\ {\displaystyle +\frac{{(lnt)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \frac{\left|d\right|}{\mathsf{\Gamma }(\alpha -\varsigma )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}\left|F_{uv}\left(s\right)\right|\frac{ds}{s}}\right.}\hfill \end{array}$$

$$\begin{array}{c}{\displaystyle +\left|d\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}\left|F_{uv}\left(\tau \right)\right|\frac{d\tau }{\tau }}\right)d{H}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}\left|G_{uv}\left(\tau \right)\right|\frac{d\tau }{\tau }}\right)d{K}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\frac{\left|b\right|}{\mathsf{\Gamma }(\gamma -\vartheta )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}\left|G_{uv}\left(s\right)\right|\frac{ds}{s}}\hfill \\ {\displaystyle +\left|b\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}\left|F_{uv}\left(\tau \right)\right|\frac{d\tau }{\tau }}\right)d{P}_i\left(s\right)}\right|}\hfill \\ +\left.{\displaystyle \left|b\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}\left|G_{uv}\left(\tau \right)\right|\frac{d\tau }{\tau }}\right)d{Q}_i\left(s\right)}\right|}\right]\hfill \\ {\displaystyle \le \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\alpha -1}[l_5\left(\right|u\left(s\right)|+\left|v\left(s\right)\right|)+{\mathsf{\Lambda }}_1]\frac{ds}{s}}\hfill \\ {\displaystyle +\frac{{(lnt)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \frac{\left|d\right|}{\mathsf{\Gamma }(\alpha -\varsigma )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}[l_5\left(\right|u\left(s\right)|+\left|v\left(s\right)\right|)+{\mathsf{\Lambda }}_1]\frac{ds}{s}}\right.}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}[l_5\left(\right|u\left(\tau \right)+\left|v\left(\tau \right)\right|)+{\mathsf{\Lambda }}_1]\frac{d\tau }{\tau }}\right)d{H}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}[l_6\left(\right|u\left(\tau \right)|+\left|v\left(\tau \right)\right|)+{\mathsf{\Lambda }}_2]\frac{d\tau }{\tau }}\right)d{K}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\frac{\left|b\right|}{\mathsf{\Gamma }(\gamma -\vartheta )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}[l_6\left(\right|u\left(s\right)|+\left|v\left(s\right)\right|)+{\mathsf{\Lambda }}_2]\frac{ds}{s}}\hfill \\ {\displaystyle +\left|b\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}[l_5\left(\right|u\left(\tau \right)|+\left|v\left(\tau \right)\right|)+{\mathsf{\Lambda }}_1]\frac{d\tau }{\tau }}\right)d{P}_i\left(s\right)}\right|}\hfill \\ +\left.{\displaystyle \left|b\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}[l_6\left(\right|u\left(\tau \right)|+\left|v\left(\tau \right)\right|)+{\mathsf{\Lambda }}_2]\frac{d\tau }{\tau }}\right)d{Q}_i\left(s\right)}\right|}\right].\hfill \end{array}$$

(29)

So, we find

$$\begin{array}{c}|{\mathcal{A}}_1(u,v)\left(t\right)|\le (l_5{\parallel (u,v)\parallel }_{\mathcal{Y}}+{\mathsf{\Lambda }}_1)\left\{{\displaystyle \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\alpha -1}\frac{ds}{s}}\right.\hfill \\ {\displaystyle +\frac{{(lnt)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \frac{\left|d\right|}{\mathsf{\Gamma }(\alpha -\varsigma )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}\frac{ds}{s}}\right.}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}\frac{d\tau }{\tau }}\right)d{H}_i\left(s\right)}\right|}\hfill \\ +\left.\left.{\displaystyle \left|b\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}\frac{d\tau }{\tau }}\right)d{P}_i\left(s\right)}\right|}\right]\right\}\hfill \\ {\displaystyle +(l_6{\parallel (u,v)\parallel }_{\mathcal{Y}}+{\mathsf{\Lambda }}_2)\frac{{(lnt)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \left|d\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}\frac{d\tau }{\tau }}\right)d{K}_i\left(s\right)}\right|}\right.}\hfill \\ {\displaystyle +\frac{\left|b\right|}{\mathsf{\Gamma }(\gamma -\vartheta )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}\frac{ds}{s}}\hfill \\ +\left.{\displaystyle \left|b\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}\frac{d\tau }{\tau }}\right)d{Q}_i\left(s\right)}\right|}\right]\hfill \\ =(l_5{\parallel (u,v)\parallel }_{\mathcal{Y}}+{\mathsf{\Lambda }}_1)\left\{{\displaystyle \frac{{(lnt)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\frac{{(lnt)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \frac{\left|d\right|}{\mathsf{\Gamma }(\alpha -\varsigma +1)}{(lnT)}^{\alpha -\varsigma }}\right.}\right.\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\alpha -{\varrho }_i}d{H}_i\left(s\right)}\right|}\hfill \\ +\left.\left.{\displaystyle \left|b\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\alpha -{\eta }_i}d{P}_i\left(s\right)}\right|}\right]\right\}\hfill \\ {\displaystyle +(l_6{\parallel (u,v)\parallel }_{\mathcal{Y}}+{\mathsf{\Lambda }}_2)\frac{{(lnt)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \left|d\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\gamma -{\sigma }_i}d{K}_i\left(s\right)}\right|}\right.}\hfill \\ {\displaystyle +\frac{\left|b\right|}{\mathsf{\Gamma }(\gamma -\vartheta +1)}{(lnT)}^{\gamma -\vartheta }+\left.{\displaystyle \left|b\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\gamma -{\theta }_i}d{Q}_i\left(s\right)}\right|}\right]}\hfill \\ \le (l_5{\parallel (u,v)\parallel }_{\mathcal{Y}}+{\mathsf{\Lambda }}_1){\mathsf{\Xi }}_1+(l_6{\parallel (u,v)\parallel }_{\mathcal{Y}}+{\mathsf{\Lambda }}_2){\mathsf{\Xi }}_2\hfill \\ \le R(l_5{\mathsf{\Xi }}_1+l_6{\mathsf{\Xi }}_2)+{\mathsf{\Lambda }}_1{\mathsf{\Xi }}_1+{\mathsf{\Lambda }}_2{\mathsf{\Xi }}_2,\forall t\in [1,T].\hfill \end{array}$$

(30)

In a similar manner, we obtain

$$\begin{array}{c}|{\mathcal{A}}_2(u,v)\left(t\right)|\le (l_5{\parallel (u,v)\parallel }_{\mathcal{Y}}+{\mathsf{\Lambda }}_1){\mathsf{\Xi }}_3+(l_6{\parallel (u,v)\parallel }_{\mathcal{Y}}+{\mathsf{\Lambda }}_2){\mathsf{\Xi }}_4\hfill \\ \le R(l_5{\mathsf{\Xi }}_3+l_6{\mathsf{\Xi }}_4)+{\mathsf{\Lambda }}_1{\mathsf{\Xi }}_3+{\mathsf{\Lambda }}_2{\mathsf{\Xi }}_4,\forall t\in [1,T].\hfill \end{array}$$

(31)

Then, by condition (26) and relations (30) and (31), we deduce

$$\begin{array}{c}{\parallel \mathcal{A}(u,v)\parallel }_{\mathcal{Y}}=\parallel {\mathcal{A}}_1(u,v)\parallel +\parallel {\mathcal{A}}_2(u,v)\parallel \le R[l_5({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)+l_6({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)]\hfill \\ +{\mathsf{\Lambda }}_1({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)+{\mathsf{\Lambda }}_2({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)\le R.\hfill \end{array}$$

(32)

So, $\mathcal{A}\left(B_R\right)\subset B_R$.

Next, we will prove that operator $\mathcal{A}$ is a contraction. For this, let $(u_1,v_1),(u_2,v_2)\in \mathcal{Y}$. Then, for any $t\in [1,T]$ we obtain

$$\begin{array}{c}|{\mathcal{A}}_1(u_1,v_1)\left(t\right)-{\mathcal{A}}_1(u_2,v_2)\left(t\right)|\hfill \\ {\displaystyle \le \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\alpha -1}|F_{u_1{v}_1}\left(s\right)-F_{u_2{v}_2}\left(s\right)|\frac{ds}{s}}\hfill \\ {\displaystyle +\frac{{(lnt)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \frac{\left|d\right|}{\mathsf{\Gamma }(\alpha -\varsigma )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}|F_{u_1{v}_1}\left(s\right)-F_{u_2{v}_2}\left(s\right)|\frac{ds}{s}}\right.}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}|F_{u_1{v}_1}\left(\tau \right)-F_{u_2{v}_2}\left(\tau \right)|\frac{d\tau }{\tau }}\right)d{H}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}|G_{u_1{v}_1}\left(\tau \right)-G_{u_2{v}_2}\left(\tau \right)|\frac{d\tau }{\tau }}\right)d{K}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\frac{\left|b\right|}{\mathsf{\Gamma }(\gamma -\vartheta )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}|G_{u_1{v}_1}\left(s\right)-G_{u_2{v}_2}\left(s\right)|\frac{ds}{s}}\hfill \\ {\displaystyle +\left|b\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}|F_{u_1{v}_1}\left(\tau \right)-F_{u_2{v}_2}\left(\tau \right)|\frac{d\tau }{\tau }}\right)d{P}_i\left(s\right)}\right|}\hfill \\ +\left.{\displaystyle \left|b\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}|G_{u_1{v}_1}\left(\tau \right)-G_{u_2{v}_2}\left(\tau \right)|\frac{d\tau }{\tau }}\right)d{Q}_i\left(s\right)}\right|}\right]\hfill \\ {\displaystyle \le l_5{\parallel (u_1,v_1)-(u_2,v_2)\parallel }_{\mathcal{Y}}\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\alpha -1}\frac{ds}{s}}\hfill \\ {\displaystyle +\frac{{(lnt)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \frac{\left|d\right|}{\mathsf{\Gamma }(\alpha -\varsigma )}{l}_5{\parallel (u_1,v_1)-(u_2,v_2)\parallel }_{\mathcal{Y}}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}\frac{ds}{s}}\right.}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}{l}_5{\parallel (u_1,v_1)-(u_2,v_2)\parallel }_{\mathcal{Y}}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}\frac{d\tau }{\tau }}\right)d{H}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}{l}_6{\parallel (u_1,v_1)-(u_2,v_2)\parallel }_{\mathcal{Y}}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}\frac{d\tau }{\tau }}\right)d{K}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\frac{\left|b\right|}{\mathsf{\Gamma }(\gamma -\vartheta )}{l}_6{\parallel (u_1,v_1)-(u_2,v_2)\parallel }_{\mathcal{Y}}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}\frac{ds}{s}}\hfill \\ {\displaystyle +\left|b\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}{l}_5{\parallel (u_1,v_1)-(u_2,v_2)\parallel }_{\mathcal{Y}}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}\frac{d\tau }{\tau }}\right)d{P}_i\left(s\right)}\right|}\hfill \\ +\left.{\displaystyle \left|b\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}{l}_6{\parallel (u_1,v_1)-(u_2,v_2)\parallel }_{\mathcal{Y}}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}\frac{d\tau }{\tau }}\right)d{Q}_i\left(s\right)}\right|}\right].\hfill \end{array}$$

(33)

Therefore, we find

$$\begin{array}{c}|{\mathcal{A}}_1(u_1,v_1)\left(t\right)-{\mathcal{A}}_1(u_2,v_2)\left(t\right)|\hfill \\ \le \parallel (u_1,v_1)-(u_2,v_2){\parallel }_{\mathcal{Y}}\left\{l_5\left\{{\displaystyle \frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnT)}^{\alpha }+\frac{{(lnT)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \frac{\left|d\right|}{\mathsf{\Gamma }(\alpha -\varsigma +1)}{(lnT)}^{\alpha -\varsigma }}\right.}\right.\right.\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\alpha -{\varrho }_i}d{H}_i\left(s\right)}\right|}\hfill \\ +\left.\left.{\displaystyle \left|b\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\alpha -{\eta }_i}d{P}_i\left(s\right)}\right|}\right]\right\}\hfill \\ {\displaystyle +l_6\frac{{(lnT)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \left|d\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\gamma -{\sigma }_i}d{K}_i\left(s\right)}\right|}\right.}\hfill \\ +\left.\left.{\displaystyle \frac{\left|b\right|}{\mathsf{\Gamma }(\gamma -\vartheta +1)}{(lnT)}^{\gamma -\vartheta }+\left|b\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\gamma -{\theta }_i}d{Q}_i\left(s\right)}\right|}\right]\right\}\hfill \\ =\parallel (u_1,v_1)-(u_2,v_2){\parallel }_{\mathcal{Y}}(l_5{\mathsf{\Xi }}_1+l_6{\mathsf{\Xi }}_2).\hfill \end{array}$$

(34)

In a similar manner, we obtain

$$\begin{array}{c}|{\mathcal{A}}_2(u_1,v_1)\left(t\right)-{\mathcal{A}}_2(u_2,v_2)\left(t\right)|\le \parallel (u_1,v_1)-(u_2,v_2){\parallel }_{\mathcal{Y}}(l_5{\mathsf{\Xi }}_3+l_6{\mathsf{\Xi }}_4),\forall t\in [1,T].\hfill \end{array}$$

(35)

Then, by relations (34) and (35), we deduce

$$\begin{array}{c}\parallel \mathcal{A}(u_1,v_1)-\mathcal{A}(u_2,v_2){\parallel }_{\mathcal{Y}}=\parallel {\mathcal{A}}_1(u_1,v_1)-{\mathcal{A}}_1(u_2,v_2)\parallel +\parallel {\mathcal{A}}_2(u_1,v_1)-{\mathcal{A}}_2(u_2,v_2)\parallel \hfill \\ \le [l_5({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)+l_6({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)]{\parallel (u_1,v_1)-(u_2,v_2)\parallel }_{\mathcal{Y}}.\hfill \end{array}$$

(36)

By (26), we conclude that operator $\mathcal{A}$ is a contraction. Therefore, operator $\mathcal{A}$ has a unique fixed point by the Banach contraction mapping principle. Hence, problem (1) and (2) has a unique solution $\left(u\right(t),v(t\left)\right),t\in [1,T]$. □

The next two results for the existence of solutions of problem (1) and (2) are based on the Krasnosel’skii fixed point theorem for the sum of two operators (see [31]).

Theorem 3.

We suppose that assumptions $\left(H1\right)$ and $\left(H2\right)$ hold. In addition, we assume that the functions $f,g:[1,T]\times {\mathbb{R}}^2\to \mathbb{R}$ are continuous and satisfy the following condition:

(H3)

There exist the continuous functions $\varphi ,\psi \in C([1,T],{\mathbb{R}}_{+})$, (${\mathbb{R}}_{+}=[0,\infty )$) such that

$$\left|f(t,x,y)\right|\le \varphi \left(t\right),\left|g(t,x,y)\right|\le \psi \left(t\right),\forall (t,x,y)\in [1,T]\times {\mathbb{R}}^2.$$

(37)

If

$$l_5\left[{\displaystyle {\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3-\frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnT)}^{\alpha }}\right]+l_6\left[{\displaystyle {\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4-\frac{1}{\mathsf{\Gamma }(\gamma +1)}{(lnT)}^{\gamma }}\right]<1,$$

(38)

then problem (1) and (2) has at least one solution $\left(u\right(t),v(t\left)\right),t\in [1,T]$.

Proof. 

We consider the number $r>0$ satisfying the condition

$$r\ge ({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)\parallel \varphi \parallel +({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)\parallel \psi \parallel ,$$

(39)

and the closed ball $B_r={\parallel (u,v)\in \mathcal{Y},\parallel (u,v)\parallel }_{\mathcal{Y}}\le r\}$. We will verify the assumptions of the Krasnosel’skii fixed point theorem for the sum of two operators. We split operator $\mathcal{A}$, defined on $B_r$, as $\mathcal{A}=\mathcal{B}+\mathcal{C}$, $\mathcal{B}=({\mathcal{B}}_1,{\mathcal{B}}_2)$, $\mathcal{C}=({\mathcal{C}}_1,{\mathcal{C}}_2)$, where ${\mathcal{B}}_i,{\mathcal{C}}_i,i=1,2$ are defined by

$$\begin{array}{c}{\displaystyle {\mathcal{B}}_1(u,v)\left(t\right)=\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\alpha -1}{F}_{uv}\left(s\right)\frac{ds}{s},}\hfill \\ {\mathcal{C}}_1(u,v)\left(t\right)={\mathcal{A}}_1(u,v)\left(t\right)-{\mathcal{B}}_1(u,v)\left(t\right),\hfill \\ {\displaystyle {\mathcal{B}}_2(u,v)\left(t\right)=\frac{1}{\mathsf{\Gamma }\left(\gamma \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\gamma -1}{G}_{uv}\left(s\right)\frac{ds}{s},}\hfill \\ {\mathcal{C}}_2(u,v)\left(t\right)={\mathcal{A}}_2(u,v)\left(t\right)-{\mathcal{B}}_2(u,v)\left(t\right),\hfill \end{array}$$

(40)

for all $t\in [1,T]$ and $(u,v)\in B_r$.

We will prove firstly that $\mathcal{B}(u_1,v_1)+\mathcal{C}(u_2,v_2)\in B_r$ for all $(u_1,v_1),(u_2,v_2)\in B_r$. For this, let $(u_1,v_1),(u_2,v_2)\in B_r$. Then, we obtain

$$\begin{array}{c}{\displaystyle |{\mathcal{B}}_1(u_1,v_1)\left(t\right)+{\mathcal{C}}_1(u_2,v_2)\left(t\right)|\le \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\alpha -1}\left|F_{u_1{v}_1}\left(s\right)\right|\frac{ds}{s}}\hfill \\ {\displaystyle +\frac{{(lnt)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \frac{\left|d\right|}{\mathsf{\Gamma }(\alpha -\varsigma )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}\left|F_{u_2{v}_2}\left(s\right)\right|\frac{ds}{s}}\right.}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}\left|F_{u_2{v}_2}\left(\tau \right)\right|\frac{d\tau }{\tau }}\right)d{H}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}\left|G_{u_2{v}_2}\left(\tau \right)\right|\frac{d\tau }{\tau }}\right)d{K}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\frac{\left|b\right|}{\mathsf{\Gamma }(\gamma -\vartheta )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}\left|G_{u_2{v}_2}\left(s\right)\right|\frac{ds}{s}}\hfill \\ {\displaystyle +\left|b\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}\left|F_{u_2{v}_2}\left(\tau \right)\right|\frac{d\tau }{\tau }}\right)d{P}_i\left(s\right)}\right|}\hfill \\ +\left.{\displaystyle \left|b\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}\left|G_{u_2{v}_2}\left(\tau \right)\right|\frac{d\tau }{\tau }}\right)d{Q}_i\left(s\right)}\right|}\right]\hfill \\ \le \parallel \varphi \parallel \left\{{\displaystyle \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\alpha -1}\frac{ds}{s}}\right.\hfill \\ {\displaystyle +\frac{{(lnt)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \frac{\left|d\right|}{\mathsf{\Gamma }(\alpha -\varsigma )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}\frac{ds}{s}}\right.}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}\frac{d\tau }{\tau }}\right)d{H}_i\left(s\right)}\right|}\hfill \\ +\left.\left.{\displaystyle \left|b\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}\frac{d\tau }{\tau }}\right)d{P}_i\left(s\right)}\right|}\right]\right\}\hfill \\ {\displaystyle +\parallel \psi \parallel \frac{{(lnt)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \left|d\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}\frac{d\tau }{\tau }}\right)d{K}_i\left(s\right)}\right|}\right.}\hfill \\ {\displaystyle +\frac{\left|b\right|}{\mathsf{\Gamma }(\gamma -\vartheta )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}\frac{ds}{s}}\hfill \\ +\left.{\displaystyle \left|b\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}\frac{d\tau }{\tau }}\right)d{Q}_i\left(s\right)}\right|}\right]\hfill \\ \le {\mathsf{\Xi }}_1\parallel \varphi \parallel +{\mathsf{\Xi }}_2\parallel \psi \parallel ,\forall t\in [1,T],\hfill \end{array}$$

(41)

and

$$\begin{array}{c}{\displaystyle |{\mathcal{B}}_2(u_1,v_1)\left(t\right)+{\mathcal{C}}_2(u_2,v_2)\left(t\right)|\le \frac{1}{\mathsf{\Gamma }\left(\gamma \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\gamma -1}\left|G_{u_1{v}_1}\left(s\right)\right|\frac{ds}{s}}\hfill \\ {\displaystyle +\frac{{(lnt)}^{\mu -1}}{|\Delta |}\left[{\displaystyle \frac{\left|a\right|}{\mathsf{\Gamma }(\gamma -\vartheta )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}\left|G_{u_2{v}_2}\left(s\right)\right|\frac{ds}{s}}\right.}\hfill \\ {\displaystyle +\left|a\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}\left|F_{u_2{v}_2}\left(\tau \right)\right|\frac{d\tau }{\tau }}\right)d{P}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\left|a\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}\left|G_{u_2{v}_2}\left(\tau \right)\right|\frac{d\tau }{\tau }}\right)d{Q}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\frac{\left|c\right|}{\mathsf{\Gamma }(\alpha -\varsigma )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}\left|F_{u_2{v}_2}\left(s\right)\right|\frac{ds}{s}}\hfill \\ {\displaystyle +\left|c\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}\left|F_{u_2{v}_2}\left(\tau \right)\right|\frac{d\tau }{\tau }}\right)d{H}_i\left(s\right)}\right|}\hfill \\ +\left.{\displaystyle \left|c\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}\left|G_{u_2{v}_2}\left(\tau \right)\right|\frac{d\tau }{\tau }}\right)d{K}_i\left(s\right)}\right|}\right]\hfill \\ {\displaystyle \le \parallel \varphi \parallel \frac{{(lnt)}^{\mu -1}}{|\Delta |}\left[{\displaystyle \left|a\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}\frac{d\tau }{\tau }}\right)d{P}_i\left(s\right)}\right|}\right.}\hfill \\ {\displaystyle +\frac{\left|c\right|}{\mathsf{\Gamma }(\alpha -\varsigma )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}\frac{ds}{s}}\hfill \\ +\left.{\displaystyle \left|c\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}\frac{d\tau }{\tau }}\right)d{H}_i\left(s\right)}\right|}\right]\hfill \\ +\parallel \psi \parallel \left\{{\displaystyle \frac{1}{\mathsf{\Gamma }\left(\gamma \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\gamma -1}\frac{ds}{s}}\right.\hfill \\ {\displaystyle +\frac{{(lnt)}^{\mu -1}}{|\Delta |}\left[{\displaystyle \frac{\left|a\right|}{\mathsf{\Gamma }(\gamma -\vartheta )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}\frac{ds}{s}}\right.}\hfill \\ {\displaystyle +\left|a\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}\frac{d\tau }{\tau }}\right)d{Q}_i\left(s\right)}\right|}\hfill \\ +\left.\left.{\displaystyle \left|c\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}\frac{d\tau }{\tau }}\right)d{K}_i\left(s\right)}\right|}\right]\right\}\hfill \\ \le {\mathsf{\Xi }}_3\parallel \varphi \parallel +{\mathsf{\Xi }}_4\parallel \psi \parallel ,\forall t\in [1,T].\hfill \end{array}$$

(42)

Then, by (41) and (42), we deduce

$$\begin{array}{c}\parallel \mathcal{B}(u_1,v_1)+\mathcal{C}(u_2,v_2){\parallel }_{\mathcal{Y}}=\parallel {\mathcal{B}}_1(u_1,v_1)+{\mathcal{C}}_1(u_2,v_2)\parallel +\parallel {\mathcal{B}}_2(u_1,v_1)+{\mathcal{C}}_2(u_2,v_2)\parallel \hfill \\ \le ({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)\parallel \varphi \parallel +({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)\parallel \psi \parallel \le r.\hfill \end{array}$$

(43)

Next, we will prove that operator $\mathcal{C}$ is a contraction mapping. Indeed, for all $(u_1,v_1),(u_2,v_2)\in B_r$, we find

$$\begin{array}{c}|{\mathcal{C}}_1(u_1,v_1)\left(t\right)-{\mathcal{C}}_1(u_2,v_2)\left(t\right)|\le \parallel (u_1,v_1)-(u_2,v_2){\parallel }_{\mathcal{Y}}\hfill \\ \times \left\{l_5\left[{\displaystyle {\mathsf{\Xi }}_1-\frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnT)}^{\alpha }}\right]+l_6{\mathsf{\Xi }}_2\right\},\forall t\in [1,T],\hfill \\ |{\mathcal{C}}_2(u_1,v_1)\left(t\right)-{\mathcal{C}}_2(u_2,v_2)\left(t\right)\parallel \le \parallel (u_1,v_1)-(u_2,v_2){\parallel }_{\mathcal{Y}}\hfill \\ \times \left\{l_5{\mathsf{\Xi }}_3+l_6\left[{\displaystyle {\mathsf{\Xi }}_4-\frac{1}{\mathsf{\Gamma }(\gamma +1)}{(lnT)}^{\gamma }}\right]\right\},\forall t\in [1,T].\hfill \end{array}$$

(44)

Therefore, by (44), we obtain

$$\begin{array}{c}\parallel \mathcal{C}(u_1,v_1)-\mathcal{C}(u_2,v_2){\parallel }_{\mathcal{Y}}\le {\parallel (u_1,v_1)-(u_2,v_2)\parallel }_{\mathcal{Y}}\hfill \\ \times \left\{l_5\left[{\displaystyle {\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3-\frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnT)}^{\alpha }}\right]+l_6\left[{\displaystyle {\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4-\frac{1}{\mathsf{\Gamma }(\gamma +1)}{(lnT)}^{\gamma }}\right]\right\}.\hfill \end{array}$$

(45)

By condition (38), we conclude that operator $\mathcal{C}$ is a contraction.

Operators ${\mathcal{B}}_1$, ${\mathcal{B}}_2$, and $\mathcal{B}$ are continuous by the continuity of functions f and g. In addition, $\mathcal{B}$ is uniformly bounded on $B_r$ because

$${\displaystyle \parallel {\mathcal{B}}_1(u,v)\parallel \le \frac{{(lnT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\parallel \varphi \parallel ,\parallel {\mathcal{B}}_2(u,v)\parallel \le \frac{{(lnT)}^{\gamma }}{\mathsf{\Gamma }(\gamma +1)}\parallel \psi \parallel ,\forall (u,v)\in B_r,}$$

(46)

and then

$${\displaystyle \parallel \mathcal{B}(u,v)\parallel \le \frac{{(lnT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\parallel \varphi \parallel +\frac{{(lnT)}^{\gamma }}{\mathsf{\Gamma }(\gamma +1)}\parallel \psi \parallel ,\forall (u,v)\in B_r.}$$

(47)

We finally prove that operator $\mathcal{B}$ is compact. Let $t_1,t_2\in [1,T]$, $t_1<t_2$. Then, for all $(u,v)\in B_r$, we find

$$\begin{array}{c}|{\mathcal{B}}_1(u,v)\left(t_2\right)-{\mathcal{B}}_1(u,v)\left(t_1\right)|\hfill \\ =\left|{\displaystyle \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{t_2}{\left({\displaystyle ln\frac{t_2}{s}}\right)}^{\alpha -1}{F}_{uv}\left(s\right)\frac{ds}{s}-\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{t_1}{\left({\displaystyle ln\frac{t_1}{s}}\right)}^{\alpha -1}{F}_{uv}\left(s\right)\frac{ds}{s}}\right|\hfill \\ {\displaystyle \le \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{t_1}\left[{\left({\displaystyle ln\frac{t_2}{s}}\right)}^{\alpha -1}-{\left({\displaystyle ln\frac{t_1}{s}}\right)}^{\alpha -1}\right]\left|F_{uv}\left(s\right)\right|\frac{ds}{s}}\hfill \\ {\displaystyle +\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_{t_1}^{t_2}{\left({\displaystyle ln\frac{t_2}{s}}\right)}^{\alpha -1}\left|F_{uv}\left(s\right)\right|\frac{ds}{s}}\hfill \\ \le \parallel \varphi \parallel \left\{{\displaystyle \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{t_1}\left[{\left({\displaystyle ln\frac{t_2}{s}}\right)}^{\alpha -1}-{\left({\displaystyle ln\frac{t_1}{s}}\right)}^{\alpha -1}\right]\frac{ds}{s}+\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_{t_1}^{t_2}{\left({\displaystyle ln\frac{t_2}{s}}\right)}^{\alpha -1}\frac{ds}{s}}\right\}\hfill \\ =\parallel \varphi \parallel \left\{{\displaystyle \frac{1}{\mathsf{\Gamma }(\alpha +1)}\left[{(ln{t}_2)}^{\alpha }-{\left({\displaystyle ln\frac{t_2}{t_1}}\right)}^{\alpha }-{(ln{t}_1)}^{\alpha }\right]+\frac{1}{\mathsf{\Gamma }(\alpha +1)}{\left({\displaystyle ln\frac{t_2}{t_1}}\right)}^{\alpha }}\right\}\hfill \\ {\displaystyle =\parallel \varphi \parallel \frac{1}{\mathsf{\Gamma }(\alpha +1)}\left[{(ln{t}_2)}^{\alpha }-{(ln{t}_1)}^{\alpha }\right],}\hfill \end{array}$$

(48)

which tends to zero as $t_2\to t_1$, independently of $(u,v)\in B_r$. We also have

$$\begin{array}{c}|{\mathcal{B}}_2(u,v)\left(t_2\right)-{\mathcal{B}}_2(u,v)\left(t_1\right)|\hfill \\ =\left|{\displaystyle \frac{1}{\mathsf{\Gamma }\left(\gamma \right)}{\int }_1^{t_2}{\left({\displaystyle ln\frac{t_2}{s}}\right)}^{\gamma -1}{G}_{uv}\left(s\right)\frac{ds}{s}-\frac{1}{\mathsf{\Gamma }\left(\gamma \right)}{\int }_1^{t_1}{\left({\displaystyle ln\frac{t_1}{s}}\right)}^{\gamma -1}{G}_{uv}\left(s\right)\frac{ds}{s}}\right|\hfill \\ {\displaystyle \le \frac{1}{\mathsf{\Gamma }\left(\gamma \right)}{\int }_1^{t_1}\left[{\left({\displaystyle ln\frac{t_2}{s}}\right)}^{\gamma -1}-{\left({\displaystyle ln\frac{t_1}{s}}\right)}^{\gamma -1}\right]\left|G_{uv}\left(s\right)\right|\frac{ds}{s}}\hfill \\ {\displaystyle +\frac{1}{\mathsf{\Gamma }\left(\gamma \right)}{\int }_{t_1}^{t_2}{\left({\displaystyle ln\frac{t_2}{s}}\right)}^{\gamma -1}\left|G_{uv}\left(s\right)\right|\frac{ds}{s}}\hfill \\ {\displaystyle \le \parallel \psi \parallel \frac{1}{\mathsf{\Gamma }(\gamma +1)}\left[{(ln{t}_2)}^{\gamma }-{(ln{t}_1)}^{\gamma }\right],}\hfill \end{array}$$

(49)

which tends to zero as $t_2\to t_1$, independently of $(u,v)\in B_r$.

Hence, by using (48) and (49), we obtain that operators ${\mathcal{B}}_1$, ${\mathcal{B}}_2$, and $\mathcal{B}$ are equicontinuous. By the Arzela–Ascoli theorem, we deduce that $\mathcal{B}$ is compact on $B_r$. Therefore, by the Krasnosel’skii fixed point theorem ([31]), we conclude that problem (1) and (2) has at least one solution $\left(u\right(t),v(t\left)\right),t\in [1,T]$. □

Theorem 4.

We suppose that assumption $\left(H1\right)$ holds and the functions $f,g:[1,T]\times {\mathbb{R}}^2\to \mathbb{R}$ are continuous and satisfy assumptions $\left(H2\right)$ and $\left(H3\right)$. If

$${\displaystyle l_5\frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnT)}^{\alpha }+l_6\frac{1}{\mathsf{\Gamma }(\gamma +1)}{(lnT)}^{\gamma }<1,}$$

(50)

then problem (1) and (2) has at least one solution $\left(u\right(t),v(t\left)\right),t\in [1,T]$.

Proof. 

As in the proof of Theorem 3, we consider the positive number $r\ge ({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)\parallel \varphi \parallel +({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)\parallel \psi \parallel$ and the closed ball $B_r$. We also split operator $\mathcal{A}$, defined on $B_r$, as $\mathcal{A}=\mathcal{B}+\mathcal{C}$, $\mathcal{B}=({\mathcal{B}}_1,{\mathcal{B}}_2)$, $\mathcal{C}=({\mathcal{C}}_1,{\mathcal{C}}_2)$, where ${\mathcal{B}}_i,{\mathcal{C}}_i,i=1,2$ are defined by (40).

For $(u_1,v_1),(u_2,v_2)\in B_r$, we obtain as in the proof of Theorem 3, that

$$\parallel \mathcal{B}(u_1,v_1)+\mathcal{C}(u_2,v_2){\parallel }_{\mathcal{Y}}\le r.$$

(51)

We will prove next that the operator $\mathcal{B}$ is a contraction. Indeed, we find

$$\begin{array}{c}{\displaystyle |{\mathcal{B}}_1(u_1,v_1)\left(t\right)-{\mathcal{B}}_1(u_2,v_2)\left(t\right)\parallel \le \parallel (u_1,v_1)-(u_2,v_2){\parallel }_{\mathcal{Y}}{l}_5\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\alpha -1}\frac{ds}{s}}\hfill \\ {\displaystyle =l_5\frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnt)}^{\alpha }{\parallel (u_1,v_1)-(u_2,v_2)\parallel }_{\mathcal{Y}}}\hfill \\ {\displaystyle \le l_5\frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnT)}^{\alpha }{\parallel (u_1,v_1)-(u_2,v_2)\parallel }_{\mathcal{Y}},\forall t\in [1,T],}\hfill \\ {\displaystyle |{\mathcal{B}}_2(u_1,v_1)\left(t\right)-{\mathcal{B}}_2(u_2,v_2)\left(t\right)\parallel \le \parallel (u_1,v_1)-(u_2,v_2){\parallel }_{\mathcal{Y}}{l}_6\frac{1}{\mathsf{\Gamma }\left(\gamma \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\gamma -1}\frac{ds}{s}}\hfill \\ {\displaystyle =l_6\frac{1}{\mathsf{\Gamma }(\gamma +1)}{(lnt)}^{\gamma }{\parallel (u_1,v_1)-(u_2,v_2)\parallel }_{\mathcal{Y}}}\hfill \\ {\displaystyle \le l_6\frac{1}{\mathsf{\Gamma }(\gamma +1)}{(lnT)}^{\gamma }{\parallel (u_1,v_1)-(u_2,v_2)\parallel }_{\mathcal{Y}},\forall t\in [1,T].}\hfill \end{array}$$

(52)

Then, we deduce

$$\begin{array}{c}\parallel \mathcal{B}(u_1,v_1)-\mathcal{B}(u_2,v_2){\parallel }_{\mathcal{Y}}\le \left[{\displaystyle l_5\frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnT)}^{\alpha }+l_6\frac{1}{\mathsf{\Gamma }(\gamma +1)}{(lnT)}^{\gamma }}\right]\hfill \\ \times \parallel (u_1,v_1)-(u_2,v_2){\parallel }_{\mathcal{Y}},\hfill \end{array}$$

(53)

that is, by (50), operator $\mathcal{B}$ is a contraction.

Operators ${\mathcal{C}}_1,{\mathcal{C}}_2$, and $\mathcal{C}$ are continuous by the continuity of the functions f and g. Moreover, $\mathcal{C}$ is uniformly bounded on $B_r$ because we have

$$\begin{array}{c}{\displaystyle |{\mathcal{C}}_1(u,v)\left(t\right)|\le \frac{{(lnT)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \frac{\left|d\right|}{\mathsf{\Gamma }(\alpha -\varsigma )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}\parallel \varphi \parallel \frac{ds}{s}}\right.}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}\parallel \varphi \parallel \frac{d\tau }{\tau }}\right)d{H}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}\parallel \psi \parallel \frac{d\tau }{\tau }}\right)d{K}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\frac{\left|b\right|}{\mathsf{\Gamma }(\gamma -\vartheta )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}\parallel \psi \parallel \frac{ds}{s}}\hfill \\ {\displaystyle +\left|b\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}\parallel \varphi \parallel \frac{d\tau }{\tau }}\right)d{P}_i\left(s\right)}\right|}\hfill \\ +\left.{\displaystyle \left|b\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}\parallel \psi \parallel \frac{d\tau }{\tau }}\right)d{Q}_i\left(s\right)}\right|}\right]\hfill \\ {\displaystyle =\parallel \varphi \parallel \frac{{(lnT)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \frac{\left|d\right|}{\mathsf{\Gamma }(\alpha -\varsigma +1)}{(lnT)}^{\alpha -\varsigma }}\right.}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\alpha -{\varrho }_i}d{H}_i\left(s\right)}\right|}\hfill \\ +\left.{\displaystyle \left|b\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\alpha -{\eta }_i}d{P}_i\left(s\right)}\right|}\right]\hfill \\ {\displaystyle +\parallel \psi \parallel \frac{{(lnT)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \left|d\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\gamma -{\sigma }_i}d{K}_i\left(s\right)}\right|}\right.}\hfill \\ +\left.{\displaystyle \frac{\left|b\right|}{\mathsf{\Gamma }(\gamma -\vartheta +1)}{(lnT)}^{\gamma -\vartheta }+\left|b\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\gamma -{\theta }_i}d{Q}_i\left(s\right)}\right|}\right]\hfill \\ =\parallel \varphi \parallel \left[{\displaystyle {\mathsf{\Xi }}_1-\frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnT)}^{\alpha }}\right]+\parallel \psi \parallel {\mathsf{\Xi }}_2,\forall t\in [1,T],(u,v)\in B_r,\hfill \end{array}$$

(54)

and

$$\begin{array}{c}{\displaystyle |{\mathcal{C}}_2(u,v)\left(t\right)|\le \frac{{(lnT)}^{\mu -1}}{|\Delta |}\left[{\displaystyle \frac{\left|a\right|}{\mathsf{\Gamma }(\gamma -\vartheta )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}\parallel \psi \parallel \frac{ds}{s}}\right.}\hfill \\ {\displaystyle +\left|a\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}\parallel \varphi \parallel \frac{d\tau }{\tau }}\right)d{P}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\left|a\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}\parallel \psi \parallel \frac{d\tau }{\tau }}\right)d{Q}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\frac{\left|c\right|}{\mathsf{\Gamma }(\alpha -\varsigma )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}\parallel \varphi \parallel \frac{ds}{s}}\hfill \\ {\displaystyle +\left|c\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}\parallel \varphi \parallel \frac{d\tau }{\tau }}\right)d{H}_i\left(s\right)}\right|}\hfill \\ +\left.{\displaystyle \left|c\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}\parallel \psi \parallel \frac{d\tau }{\tau }}\right)d{K}_i\left(s\right)}\right|}\right]\hfill \\ {\displaystyle =\parallel \varphi \parallel \frac{{(lnT)}^{\mu -1}}{|\Delta |}\left[{\displaystyle \left|a\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\alpha -{\eta }_i}d{P}_i\left(s\right)}\right|}\right.}\hfill \\ +\left.{\displaystyle \frac{\left|c\right|}{\mathsf{\Gamma }(\alpha -\varsigma +1)}{(lnT)}^{\alpha -\varsigma }+\left|c\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\alpha -{\varrho }_i}d{H}_i\left(s\right)}\right|}\right]\hfill \\ {\displaystyle +\parallel \psi \parallel \frac{{(lnT)}^{\mu -1}}{|\Delta |}\left[{\displaystyle \frac{\left|a\right|}{\mathsf{\Gamma }(\gamma -\vartheta +1)}{(lnT)}^{\gamma -\vartheta }}\right.}\hfill \\ {\displaystyle +\left|a\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\gamma -{\theta }_i}d{Q}_i\left(s\right)}\right|}\hfill \\ +\left.{\displaystyle \left|c\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i+1)}\left|{\displaystyle {\int }_1^T{(lns)}^{\gamma -{\sigma }_i}d{K}_i\left(s\right)}\right|}\right]\hfill \\ =\parallel \varphi \parallel {\mathsf{\Xi }}_3+\parallel \psi \parallel \left[{\displaystyle {\mathsf{\Xi }}_4-\frac{1}{\mathsf{\Gamma }(\gamma +1)}{(lnT)}^{\gamma }}\right],\forall t\in [1,T],(u,v)\in B_r.\hfill \end{array}$$

(55)

Therefore, by (54) and (55), we obtain

$$\begin{array}{c}\parallel {\mathcal{C}}_1(u,v)\parallel \le \parallel \varphi \parallel \left[{\displaystyle {\mathsf{\Xi }}_1-\frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnT)}^{\alpha }}\right]+\parallel \psi \parallel {\mathsf{\Xi }}_2,\hfill \\ \parallel {\mathcal{C}}_2(u,v)\parallel \le \parallel \varphi \parallel {\mathsf{\Xi }}_3+\parallel \psi \parallel \left[{\displaystyle {\mathsf{\Xi }}_4-\frac{1}{\mathsf{\Gamma }(\gamma +1)}{(lnT)}^{\gamma }}\right],\hfill \end{array}$$

(56)

and then

$$\begin{array}{c}\parallel \mathcal{C}(u,v)\parallel \le \parallel \varphi \parallel \left[{\displaystyle {\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3-\frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnT)}^{\alpha }}\right]\hfill \\ +\parallel \psi \parallel \left[{\displaystyle {\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4-\frac{1}{\mathsf{\Gamma }(\gamma +1)}{(lnT)}^{\gamma }}\right],\forall (u,v)\in B_r.\hfill \end{array}$$

(57)

We finally prove that operator $\mathcal{C}$ is compact. Let $t_1,t_2\in [1,T]$, $t_1<t_2$. Then, for all $(u,v)\in B_r$, we find

$$\begin{array}{c}|{\mathcal{C}}_1(u,v)\left(t_2\right)-{\mathcal{C}}_1(u,v)\left(t_1\right)|\hfill \\ {\displaystyle \le \left\{\parallel \varphi \parallel \left[{\displaystyle {\mathsf{\Xi }}_1-\frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnT)}^{\alpha }}\right]+\parallel \psi \parallel {\mathsf{\Xi }}_2\right\}\frac{1}{{(lnT)}^{\lambda -1}}\left[{(ln{t}_2)}^{\lambda -1}-{(ln{t}_1)}^{\lambda -1}\right],}\hfill \end{array}$$

(58)

which tends to zero as $t_2\to t_1$, independently of $(u,v)\in B_r$. We also obtain

$$\begin{array}{c}|{\mathcal{C}}_2(u,v)\left(t_2\right)-{\mathcal{C}}_2(u,v)\left(t_1\right)|\hfill \\ {\displaystyle \le \left\{\parallel \varphi \parallel {\mathsf{\Xi }}_3+\parallel \psi \parallel \left[{\displaystyle {\mathsf{\Xi }}_4-\frac{1}{\mathsf{\Gamma }(\gamma +1)}{(lnT)}^{\gamma }}\right]\right\}\frac{1}{{(lnT)}^{\mu -1}}\left[{(ln{t}_2)}^{\mu -1}-{(ln{t}_1)}^{\mu -1}\right],}\hfill \end{array}$$

(59)

which tends to zero as $t_2\to t_1$, independently of $(u,v)\in B_r$.

So, by using inequalities (58) and (59), we obtain that operators ${\mathcal{C}}_1$, ${\mathcal{C}}_2$, and $\mathcal{C}$ are equicontinuous. By the Arzela–Ascoli theorem, we conclude that $\mathcal{C}$ is compact on $B_r$. Then, by applying the Krasnosel’skii fixed point theorem (see [31]), we deduce that problem (1) and (2) has at least one solution $\left(u\right(t),v(t\left)\right),t\in [1,T]$. □

Our next result is based on the Schaefer fixed point theorem (see [32]).

Theorem 5.

We assume that assumption $\left(H1\right)$ holds. In addition, we suppose that the functions $f,g:[1,T]\times {\left(\mathbb{R}\right)}^2\to \mathbb{R}$ are continuous and satisfy the following condition:

(H4)

There exist positive constants $M_1,M_2$ such that

$$\left|f(t,x,y)\right|\le M_1,\left|g(t,x,y)\right|\le M_2,\forall t\in [1,T],x,y\in \mathbb{R}.$$

(60)

Then, there exists at least one solution $\left(u\right(t),v(t\left)\right),t\in [1,T]$ for problem (1) and (2).

Proof. 

Firstly, we show that $\mathcal{A}$ is completely continuous. Operator $\mathcal{A}$ is continuous. Indeed, let $(u_n,v_n)\in \mathcal{Y}$, $n\in \mathbb{N}$, $(u_n,v_n)\to (u,v)$, as $n\to \infty$ in $\mathcal{Y}$. Then, for each $t\in [1,T]$, we obtain

$$\begin{array}{c}|{\mathcal{A}}_1(u_n,v_n)\left(t\right)-{\mathcal{A}}_1(u,v)\left(t\right)|\hfill \\ {\displaystyle \le \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\alpha -1}|F_{u_n{v}_n}\left(s\right)-F_{uv}\left(s\right)|\frac{ds}{s}}\hfill \\ {\displaystyle +\frac{{(lnt)}^{\lambda -1}}{|\Delta |}\left[{\displaystyle \frac{\left|d\right|}{\mathsf{\Gamma }(\alpha -\varsigma )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}|F_{u_n{v}_n}\left(s\right)-F_{uv}\left(s\right)|\frac{ds}{s}}\right.}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}|F_{u_n{v}_n}\left(\tau \right)-F_{uv}\left(\tau \right)|\frac{d\tau }{\tau }}\right)d{H}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\left|d\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}|G_{u_n{v}_n}\left(\tau \right)-G_{uv}\left(\tau \right)|\frac{d\tau }{\tau }}\right)d{K}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\frac{\left|b\right|}{\mathsf{\Gamma }(\gamma -\vartheta )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}|G_{u_n{v}_n}\left(s\right)-G_{uv}\left(s\right)|\frac{ds}{s}}\hfill \\ {\displaystyle +\left|b\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}|F_{u_n{v}_n}\left(\tau \right)-F_{uv}\left(\tau \right)|\frac{d\tau }{\tau }}\right)d{P}_i\left(s\right)}\right|}\hfill \\ +\left.{\displaystyle \left|b\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}|G_{u_n{v}_n}\left(\tau \right)-G_{uv}\left(\tau \right)|\frac{d\tau }{\tau }}\right)d{Q}_i\left(s\right)}\right|}\right],\hfill \end{array}$$

(61)

and

$$\begin{array}{c}|{\mathcal{A}}_2(u_n,v_n)\left(t\right)-{\mathcal{A}}_2(u,v)\left(t\right)|\hfill \\ {\displaystyle \le \frac{1}{\mathsf{\Gamma }\left(\gamma \right)}{\int }_1^t{\left({\displaystyle ln\frac{t}{s}}\right)}^{\gamma -1}|G_{u_n{v}_n}\left(s\right)-G_{uv}\left(s\right)|\frac{ds}{s}}\hfill \\ {\displaystyle +\frac{{(lnt)}^{\mu -1}}{|\Delta |}\left[{\displaystyle \frac{\left|a\right|}{\mathsf{\Gamma }(\gamma -\vartheta )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\gamma -\vartheta -1}|G_{u_n{v}_n}\left(s\right)-G_{uv}\left(s\right)|\frac{ds}{s}}\right.}\hfill \\ {\displaystyle +\left|a\right|\sum _{i=1}^p\frac{1}{\mathsf{\Gamma }(\alpha -{\eta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\eta }_i-1}|F_{u_n{v}_n}\left(\tau \right)-F_{uv}\left(\tau \right)|\frac{d\tau }{\tau }}\right)d{P}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\left|a\right|\sum _{i=1}^q\frac{1}{\mathsf{\Gamma }(\gamma -{\theta }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\theta }_i-1}|G_{u_n{v}_n}\left(\tau \right)-G_{uv}\left(\tau \right)|\frac{d\tau }{\tau }}\right)d{Q}_i\left(s\right)}\right|}\hfill \\ {\displaystyle +\frac{\left|c\right|}{\mathsf{\Gamma }(\alpha -\varsigma )}{\int }_1^T{\left({\displaystyle ln\frac{T}{s}}\right)}^{\alpha -\varsigma -1}|F_{u_n{v}_n}\left(s\right)-F_{uv}\left(s\right)|\frac{ds}{s}}\hfill \\ {\displaystyle +\left|c\right|\sum _{i=1}^m\frac{1}{\mathsf{\Gamma }(\alpha -{\varrho }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\alpha -{\varrho }_i-1}|F_{u_n{v}_n}\left(\tau \right)-F_{uv}\left(\tau \right)|\frac{d\tau }{\tau }}\right)d{H}_i\left(s\right)}\right|}\hfill \\ +\left.{\displaystyle \left|c\right|\sum _{i=1}^n\frac{1}{\mathsf{\Gamma }(\gamma -{\sigma }_i)}\left|{\displaystyle {\int }_1^T\left({\displaystyle {\int }_1^s{\left({\displaystyle ln\frac{s}{\tau }}\right)}^{\gamma -{\sigma }_i-1}|G_{u_n{v}_n}\left(\tau \right)-G_{uv}\left(\tau \right)|\frac{d\tau }{\tau }}\right)d{K}_i\left(s\right)}\right|}\right].\hfill \end{array}$$

(62)

Because f and g are continuous, we find

$$\begin{array}{c}|F_{u_n{v}_n}\left(s\right)-F_{uv}\left(s\right)|=|f(s,u_n\left(s\right),v_n\left(s\right))-f(s,u\left(s\right),v\left(s\right))|\to 0,\mathrm{and}\hfill \\ |G_{u_n{v}_n}\left(s\right)-G_{uv}\left(s\right)|=|g(s,u_n\left(s\right),v_n\left(s\right))-g(s,u\left(s\right),v\left(s\right))|\to 0,\hfill \end{array}$$

(63)

as $n\to \infty$, for all $s\in [1,T]$. So, by relations (61)–(63), we deduce

$$\parallel {\mathcal{A}}_1(u_n,v_n)-{\mathcal{A}}_1(u,v)\parallel \to 0,\parallel {\mathcal{A}}_2(u_n,v_n)-{\mathcal{A}}_2(u,v)\parallel \to 0,\mathrm{as}n\to \infty ,$$

(64)

and then $\parallel \mathcal{A}(u_n,v_n)-\mathcal{A}(u,v){\parallel }_{\mathcal{Y}}\to 0$, as $n\to \infty$; i.e., $\mathcal{A}$ is a continuous operator.

We prove now that $\mathcal{A}$ maps bounded sets into bounded sets in $\mathcal{Y}$. For $R>0$, let $B_R={\{(u,v)\in \mathcal{Y},\parallel (u,v)\parallel }_{\mathcal{Y}}\le R\}$. Then, by using (60) and similar computations to those in the first part of the proof of Theorem 2, we obtain

$$|{\mathcal{A}}_1(u,v)\left(t\right)|\le M_1{\mathsf{\Xi }}_1+M_2{\mathsf{\Xi }}_2,\left|{\mathcal{A}}_2(u,v)\left(t\right)\right|\le M_1{\mathsf{\Xi }}_3+M_2{\mathsf{\Xi }}_4,$$

(65)

for all $t\in [1,T]$ and $(u,,v)\in B_R$. Then, by (65), we conclude

$${\parallel \mathcal{A}(u,v)\parallel }_{\mathcal{Y}}\le M_1({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)+M_2({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4),\forall (u,v)\in B_R;$$

(66)

i.e., $\mathcal{A}\left(B_R\right)$ is bounded.

In the following, we will prove that $\mathcal{A}$ maps bounded sets into equicontinuous sets. For this, let $t_1,t_2\in [1,T]$, $t_1<t_2$, and $(u,v)\in B_R$. Then, by using similar computations to those in the proofs of Theorems 3 and 4, we find

$$\begin{array}{c}|{\mathcal{A}}_1(u,v)\left(t_2\right)-{\mathcal{A}}_1(u,v)\left(t_1\right)|\hfill \\ \le |{\mathcal{B}}_1(u,v)\left(t_2\right)-{\mathcal{B}}_1(u,v)\left(t_1\right)|+|{\mathcal{C}}_1(u,v)\left(t_2\right)-{\mathcal{C}}_1(u,v)\left(t_1\right)|\hfill \\ {\displaystyle \le M_1\frac{1}{\mathsf{\Gamma }(\alpha +1)}\left[{(ln{t}_2)}^{\alpha }-{(ln{t}_1)}^{\alpha }\right]}\hfill \\ +\left\{M_1\left[{\displaystyle {\mathsf{\Xi }}_1-\frac{1}{\mathsf{\Gamma }(\alpha +1)}{(lnT)}^{\alpha }}\right]+M_2{\mathsf{\Xi }}_2\right\}\hfill \\ {\displaystyle \times \frac{1}{{(lnT)}^{\lambda -1}}\left[{(ln{t}_2)}^{\lambda -1}-{(ln{t}_1)}^{\lambda -1}\right]\to 0,\mathrm{as}{t}_2\to t_1,}\hfill \end{array}$$

(67)

independently of $(u,v)\in B_R$, and

$$\begin{array}{c}|{\mathcal{A}}_2(u,v)\left(t_2\right)-{\mathcal{A}}_2(u,v)\left(t_1\right)|\hfill \\ \le |{\mathcal{B}}_2(u,v)\left(t_2\right)-{\mathcal{B}}_2(u,v)\left(t_1\right)|+|{\mathcal{C}}_2(u,v)\left(t_2\right)-{\mathcal{C}}_2(u,v)\left(t_1\right)|\hfill \\ {\displaystyle \le M_2\frac{1}{\mathsf{\Gamma }(\gamma +1)}\left[{(ln{t}_2)}^{\gamma }-{(ln{t}_1)}^{\gamma }\right]}\hfill \\ +\left\{M_1{\mathsf{\Xi }}_3+M_2\left[{\displaystyle {\mathsf{\Xi }}_4-\frac{1}{\mathsf{\Gamma }(\gamma +1)}{(lnT)}^{\gamma }}\right]\right\}\hfill \\ {\displaystyle \times \frac{1}{{(lnT)}^{\mu -1}}\left[{(ln{t}_2)}^{\mu -1}-{(ln{t}_1)}^{\mu -1}\right]\to 0,\mathrm{as}{t}_2\to t_1,}\hfill \end{array}$$

(68)

independently of $(u,v)\in B_R$.

Therefore, by using relations (67) and (68), we obtain that operators ${\mathcal{A}}_1$ and ${\mathcal{A}}_2$ are equicontinuous, and so, $\mathcal{A}$ is equicontinuous. So, the operator $\mathcal{A}:\mathcal{Y}\to \mathcal{Y}$ is completely continuous, by using the Arzela–Ascoli theorem.

Finally, we show that set $\mathcal{U}=\left\{\right(u,v)\in \mathcal{Y},(u,v)=\omega \mathcal{A}(u,v),0\le \omega \le 1\}$ is bounded. Let $(u,v)\in \mathcal{U}$, i.e., there exists $\omega \in [0,1]$ such that $(u,v)=\omega \mathcal{A}(u,v)$ or $u\left(t\right)=\omega {\mathcal{A}}_1(u,v)\left(t\right)$ and $v\left(t\right)=\omega {\mathcal{A}}_2(u,v)\left(t\right)$ for all $t\in [1,T]$. Then, by $\left(H4\right)$, we obtain in a similar manner as that used in the first part of this proof that

$$\begin{array}{c}\left|u\left(t\right)\right|=\omega |{\mathcal{A}}_1(u,v)\left(t\right)|\le |{\mathcal{A}}_1(u,v)\left(t\right)|\le M_1{\mathsf{\Xi }}_1+M_2{\mathsf{\Xi }}_2,\forall t\in [1,T],\hfill \\ \left|v\left(t\right)\right|=\omega |{\mathcal{A}}_2(u,v)\left(t\right)|\le |{\mathcal{A}}_2(u,v)\left(t\right)|\le M_1{\mathsf{\Xi }}_3+M_2{\mathsf{\Xi }}_4,\forall t\in [1,T],\hfill \end{array}$$

(69)

and then

$${\parallel (u,v)\parallel }_{\mathcal{Y}}=\parallel u\parallel +\parallel v\parallel \le M_1({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)+M_2({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4).$$

(70)

This shows that the set $\mathcal{U}$ is bounded. Therefore, by the Schaefer fixed point theorem (see [32]), we deduce that operator $\mathcal{A}$ has at least one fixed point. Hence, problem (1) and (2) has at least one solution. □

In our last existence result, we will use the Leray–Schauder nonlinear alternative (see [33]).

Theorem 6.

We suppose that assumption $\left(H1\right)$ holds. Moreover, we assume that the functions $f,g:[1,T]\times {\mathbb{R}}^2\to \mathbb{R}$ are continuous, and the following conditions are satisfied:

(H5)

There exist the functions $p_1,p_2\in C([1,T],{\mathbb{R}}_{+})$ and the functions ${\phi }_1,{\phi }_2\in C({\mathbb{R}}_{+}\times {\mathbb{R}}_{+},{\mathbb{R}}_{+})$ nondecreasing in each of both variables such that

$$\left|f(t,x,y)\right|\le p_1\left(t\right){\phi }_1\left(\right|x|,|y\left|\right),\left|g(t,x,y)\right|\le p_2\left(t\right){\phi }_2\left(\right|x|,|y\left|\right),\forall t\in [1,T],x,y\in \mathbb{R}.$$

(71)

(H6)

There exists a positive constant L such that

$${\displaystyle \frac{L}{\parallel p_1\parallel {\phi }_1(L,L)({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)+\parallel p_2\parallel {\phi }_2(L,L)({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)}>1.}$$

(72)

Then, the fractional boundary value problem (1) and (2) has at least one solution $\left(u\right(t),v(t\left)\right),t\in [1,T]$.

Proof. 

We define the set $\mathcal{V}=\{(u,v)\in \mathcal{Y},\parallel (u,v){\parallel }_{\mathcal{Y}}<L\}$, where L is the constant given by (72). The operator $\mathcal{A}:\overline{\mathcal{V}}\to \mathcal{Y}$ is completely continuous.

We assume that there exist $(u,v)\in \partial \mathcal{V}$ such that $(u,v)=\nu \mathcal{A}(u,v)$ for some $\nu \in (0,1)$. Then, we find

$$\begin{array}{c}\left|u\left(t\right)\right|=\nu |{\mathcal{A}}_1(u,v)\left(t\right)|\le |{\mathcal{A}}_1(u,v)\left(t\right)|\le \parallel p_1\parallel {\phi }_1(\parallel u\parallel ,\parallel v\parallel ){\mathsf{\Xi }}_1+\parallel p_2\parallel {\phi }_2(\parallel u\parallel ,\parallel v\parallel ){\mathsf{\Xi }}_2,\hfill \\ \left|v\left(t\right)\right|=\nu |{\mathcal{A}}_2(u,v)\left(t\right)|\le |{\mathcal{A}}_2(u,v)\left(t\right)|\le \parallel p_1\parallel {\phi }_1(\parallel u\parallel ,\parallel v\parallel ){\mathsf{\Xi }}_3+\parallel p_2\parallel {\phi }_2(\parallel u\parallel ,\parallel v\parallel ){\mathsf{\Xi }}_4,\hfill \end{array}$$

(73)

for all $t\in [1,T]$, and, therefore,

$$\begin{array}{c}{\parallel (u,v)\parallel }_{\mathcal{Y}}=\parallel u\parallel +\parallel v\parallel \le \parallel p_1\parallel {\phi }_1(\parallel u\parallel ,\parallel v\parallel )({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)+\parallel p_2\parallel {\phi }_2(\parallel u\parallel ,\parallel v\parallel )({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)\hfill \\ \le \parallel p_1\parallel {\phi }_1{(\parallel (u,v)\parallel }_{\mathcal{Y}}{,\parallel (u,v)\parallel }_{\mathcal{Y}})({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)+\parallel p_2\parallel {\phi }_2(\parallel (u,v){\parallel }_{\mathcal{Y}},\parallel (u,v){\parallel }_{\mathcal{Y}})({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4).\hfill \end{array}$$

(74)

So, we obtain

$${\displaystyle \frac{L}{\parallel p_1\parallel {\phi }_1(L,L)({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)+\parallel p_2\parallel {\phi }_2(L,L)({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)}\le 1,}$$

(75)

which, based on (72), is a contradiction.

We deduce that there is no $(u,v)\in \partial \mathcal{V}$ such that $(u,v)=\nu \mathcal{A}(u,v)$ for some $\nu \in (0,1)$. Therefore, by the Leray–Schauder nonlinear alternative (see [33]), we conclude that $\mathcal{A}$ has a fixed point $(u,v)\in \overline{\mathcal{V}}$, which is a solution of problem (1) and (2). □

4. Examples

In this section, we will present some examples that illustrate our theorems obtained in Section 3.

We consider $T=9$, $\alpha =\frac{3}{2}$, $\gamma =\frac{4}{3}$, $\beta =\frac{2}{5}$, $\delta =\frac{1}{7}$, $\varsigma =\frac{1}{4}$, $\vartheta =0$, $m=2$, $n=2$, $p=2$, $q=2$, ${\varrho }_1=0$, ${\varrho }_2=\frac{2}{11}$, ${\sigma }_1=\frac{3}{8}$, ${\sigma }_2=\frac{5}{12}$, ${\eta }_1=0$, ${\eta }_2=\frac{1}{2}$, ${\theta }_1=0$, and ${\theta }_2=\frac{4}{9}$.

In addition, we introduce the following functions

$$\begin{array}{l}{H}_1\left(s\right)=\left\{\begin{array}{c}{\displaystyle -\frac{513}{22},s\in \left[{\displaystyle 1,\frac{19}{6}}\right);}\hfill \\ {\displaystyle \frac{15}{22},s\in \left[{\displaystyle \frac{19}{6},\frac{11}{2}}\right);}\hfill \\ {\displaystyle \frac{8}{27}{(s-2)}^{9/2}+\frac{15}{22}-\frac{7^{9/2}}{27·2^{3/2}},s\in \left[{\displaystyle \frac{11}{2},\frac{25}{3}}\right);}\hfill \\ {\displaystyle \frac{8·{19}^{9/2}}{3^{15/2}}+\frac{15}{22}-\frac{7^{9/2}}{27·2^{3/2}},s\in \left[\frac{25}{3},9\right];}\hfill \end{array}\right.\\ H_2\left(s\right)=\left\{\begin{array}{c}{\displaystyle \frac{35}{12},s\in \left[1,\frac{54}{13}\right);}\hfill \\ {\displaystyle \frac{8}{3},s\in \left[\frac{54}{13},9\right];}\hfill \end{array}\right.\\ K_1\left(s\right)=-\frac{67}{72}{(s-1)}^3,s\in [1,9];\\ K_2\left(s\right)=\left\{\begin{array}{c}2,s\in [1,7);\hfill \\ {\displaystyle \frac{97}{38},s\in [7,9];}\hfill \end{array}\right.\\ P_1\left(s\right)=\left\{\begin{array}{c}2s+15,s\in \left[{\displaystyle 1,\frac{72}{11}}\right);\hfill \\ {\displaystyle 2s+\frac{49}{6},s\in \left[{\displaystyle \frac{72}{11},9}\right];}\hfill \end{array}\right.\\ P_2\left(s\right)=\left\{\begin{array}{c}{\displaystyle \frac{1}{2},s\in \left[{\displaystyle 1,\frac{33}{8}}\right);}\hfill \\ {\displaystyle \frac{23}{26},s\in \left[{\displaystyle \frac{33}{8},9}\right];}\hfill \end{array}\right.\\ Q_1\left(s\right)=\left\{\begin{array}{c}17,s\in [1,2);\hfill \\ 14,s\in [2,7);\hfill \\ {\displaystyle \frac{5}{68}{(s-3)}^{17/5}+14-\frac{5}{17}·2^{24/5},s\in [7,9];}\hfill \end{array}\right.\\ Q_2\left(s\right)=\left\{\begin{array}{c}{\displaystyle \frac{102}{25}{\left(s-\frac{1}{2}\right)}^{-25/6},s\in \left[{\displaystyle 1,\frac{53}{8}}\right);}\hfill \\ {\displaystyle \frac{51·2^{27/2}}{25·{49}^{25/6}},s\in \left[{\displaystyle \frac{53}{8},9}\right].}\hfill \end{array}\right.\end{array}$$

(76)

We consider the system of fractional differential equations

$$\left\{\begin{array}{c}{}^{HH}{D}_1^{\frac{3}{2},\frac{2}{5}}u\left(t\right)=f(t,u\left(t\right),v\left(t\right)),t\in [1,T],\hfill \\ {}^{HH}{D}_1^{\frac{4}{3},\frac{1}{7}}v\left(t\right)=g(t,u\left(t\right),v\left(t\right)),t\in [1,T],\hfill \end{array}\right.$$

(77)

subject to the nonlocal coupled boundary conditions

$$\left\{\begin{array}{c}{\displaystyle u\left(1\right)=0,{}^H{D}_1^{\frac{1}{4}}u\left(9\right)=24u\left({\displaystyle \frac{19}{6}}\right)+\frac{4}{3}{\int }_{\frac{11}{2}}^{\frac{25}{3}}{(s-2)}^{\frac{7}{2}}u\left(s\right)ds-\frac{1}{4}{}^H{D}_1^{\frac{2}{11}}u\left({\displaystyle \frac{54}{13}}\right)}\hfill \\ {\displaystyle -\frac{67}{24}{\int }_1^9{(s-1)}^2{}^H{D}_1^{\frac{3}{8}}v\left(s\right)ds+\frac{21}{38}{}^H{D}_1^{\frac{5}{12}}v\left(\frac{7}{3}\right),}\hfill \\ {\displaystyle v\left(1\right)=0,v\left(9\right)=2{\int }_1^{9}u\left(s\right)ds-\frac{41}{6}u\left({\displaystyle \frac{72}{11}}\right)+\frac{5}{13}{}^H{D}_1^{1/2}u\left({\displaystyle \frac{33}{8}}\right)}\hfill \\ {\displaystyle -3v\left(2\right)+\frac{1}{4}{\int }_7^9{(s-3)}^{\frac{12}{5}}v\left(s\right)ds-17{\int }_1^{\frac{53}{8}}{\left(s-\frac{1}{2}\right)}^{-\frac{31}{6}}{}^H{D}_1^{\frac{4}{9}}v\left(s\right)ds.}\hfill \end{array}\right.$$

(78)

After some computations, using the Mathematica program, we obtain $\lambda =\frac{17}{10}$, $\mu =\frac{10}{7}$, $a\approx -1831.69757626$, $b\approx -449.04583604$, $c\approx 10.3109466$, $d\approx 37.64671349$, and $\Delta \approx -64327.3062105\ne 0$. So, assumption $\left(H1\right)$ is satisfied.

In addition, we find ${\mathsf{\Xi }}_1\approx 4.99459858$, ${\mathsf{\Xi }}_2\approx 1.73975237$, ${\mathsf{\Xi }}_3\approx 0.93885308$, and ${\mathsf{\Xi }}_4\approx 5.31130101$.

Example 1. We consider the functions

$$\begin{array}{c}{\displaystyle f(t,x,y)=\frac{1}{27}{e}^{-{(t-1)}^3}\sqrt{x^2+1}-\frac{2}{21}{e}^{-3t^2}arctany+\frac{cos(4t+1)}{\sqrt[3]{t^2+1}},}\hfill \\ {\displaystyle g(t,x,y)=\frac{1}{12(t^2+1)}·\frac{\left|x\right|}{1+\left|x\right|}+\frac{3}{41(t+2)}{cos}^{2}y-t^3+7,}\hfill \end{array}$$

(79)

for all $t\in [1,9]$ and $x,y\in \mathbb{R}$.

We have

$$\begin{array}{c}{\displaystyle |f(t,x_1,y_1)-f(t,x_2,y_2)|\le \frac{1}{27}|x_1-x_2|+\frac{2}{21}|y_1-y_2|,}\hfill \\ {\displaystyle |g(t,x_1,y_1)-g(t,x_2,y_2)|\le \frac{1}{24}|x_1-x_2|+\frac{2}{41}|y_1-y_2|,}\hfill \end{array}$$

(80)

for all $t\in [1,9]$ and $x_i,y_i\in \mathbb{R}$, $i=1,2$. So, we find $l_1=\frac{1}{27}$, $l_2=\frac{2}{21}$, $l_3=\frac{1}{24}$, $l_4=\frac{2}{41}$ (from assumption $\left(H2\right)$), and then $l_5=\frac{2}{21}$ and $l_6=\frac{2}{41}$. Because $l_5({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)+l_6({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)\approx 0.909<1$, then condition (26) is satisfied. Therefore, by Theorem 2, we conclude that the boundary value problem (77) and (78) with the nonlinearities (79) has a unique solution $\left(u\right(t),v(t\left)\right),t\in [1,9]$.

Example 2. We consider the functions

$$\begin{array}{c}{\displaystyle f(t,x,y)=\frac{1}{t^2+12}·\frac{x}{x^2+1}+\frac{(t^7+1)e^{-2t+3}}{15(t^3+31)}siny-\frac{cos(4t^2+9)}{t+6},}\hfill \\ {\displaystyle g(t,x,y)=\frac{2t+1}{7\sqrt{t^4+18}}{e}^{-x^2}-\frac{4t^2+1}{3(t^5+42)}·\frac{2y^4+1}{y^4+3}+\frac{e^{-t+6}}{5t+4},}\hfill \end{array}$$

(81)

for all $t\in [1,9]$ and $x,y\in \mathbb{R}$.

We have the following inequalities

$$\begin{array}{c}{\displaystyle \left|f(t,x,y)\right|\le \frac{1}{2(t^2+12)}+\frac{(t^7+1)e^{-2t+3}}{15(t^3+31)}+\frac{1}{t+6}=:\varphi \left(t\right),}\hfill \\ {\displaystyle \left|g(t,x,y)\right|\le \frac{2t+1}{7\sqrt{t^4+18}}+\frac{2(4t^2+1)}{3(t^5+42)}+\frac{e^{-t+6}}{5t+4}=:\psi \left(t\right),}\hfill \end{array}$$

(82)

for all $t\in [1,9]$ and $x,y\in \mathbb{R}$, and

$$\begin{array}{c}|f(t,x_1,y_1)-f(t,x_2,y_2)|\le 0.077|x_1-x_2|+0.1265|y_1-y_2|,\hfill \\ |g(t,x_1,y_1)-g(t,x_2,y_2)|\le 0.1057|x_1-x_2|+0.1041|y_1-y_2|,\hfill \end{array}$$

(83)

for all $t\in [1,9]$ and $x_i,y_i\in \mathbb{R},i=1,2$. So $l_1=0.077$, $l_2=0.1265$, $l_3=0.1057$, $l_4=0.1041$, and so $l_5=0.1265$ and $l_6=0.1057$. So, assumptions $\left(H2\right)$ and $\left(H3\right)$ are satisfied. In addition, we obtain $z_1:={\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3-\frac{1}{\mathsf{\Gamma }(5/2)}{(ln9)}^{3/2}\approx 3.48339869$ and $z_2:={\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4-\frac{1}{\mathsf{\Gamma }(7/3)}{(ln9)}^{4/3}\approx 4.65193119$. Then, we find $l_5{z}_1+l_6{z}_2\approx 0.932359<1$, i.e., condition (38) is also satisfied. By Theorem 3, we deduce that problem (77) and (78) with the nonlinearities (81) has at least one solution $\left(u\right(t),v(t\left)\right),t\in [1,9]$.

Example 3. We consider the functions

$$\begin{array}{c}{\displaystyle f(t,x,y)=\frac{e^{-3t+1}}{\sqrt[3]{t^4+2}}sin(x+y)+\frac{1}{5}\sqrt{t^2+13},}\hfill \\ {\displaystyle g(t,x,y)=\frac{t}{t^2+1}·\frac{3x^2+1}{x^2+4}-(t^3-2t)e^{-y^2}+arctan\frac{t^2}{t^6+2},}\hfill \end{array}$$

(84)

for all $t\in [1,9]$ and $x,y\in \mathbb{R}$. We obtain $\left|f\right(t,x,y\left)\right|\le 2$ and $\left|g(t,x,y)\right|\le 712.5+\frac{\pi }{2}$, for all $t\in [1,9]$ and $x,y\in \mathbb{R}$. So, $M_1=2$ and $M_2=712.5+\frac{\pi }{2}$ (from assumption $\left(H4\right)$). By Theorem 5, we conclude that the boundary value problem (77) and (78) with the nonlinearities (84) has at least one solution $\left(u\right(t),v(t\left)\right),t\in [1,9]$.

Example 4. We consider the functions

$$\begin{array}{c}{\displaystyle f(t,x,y)=\frac{3}{t^2+81}\left({\displaystyle \frac{x}{28}cos(x^3+y)-\frac{y}{50}arctan\frac{x}{y^2+1}+\frac{1}{2}}\right),}\hfill \\ {\displaystyle g(t,x,y)=\frac{2t}{t^3+93}\left({\displaystyle \frac{5x}{23}sin(x^2-4y)+\frac{y}{34}arctan\frac{x+1}{y^4+2}-\frac{7}{8}}\right),}\hfill \end{array}$$

(85)

for all $t\in [1,9]$ and $x,y\in \mathbb{R}$. We have the inequalities

$$\begin{array}{c}{\displaystyle \left|f(t,x,y)\right|\le \frac{3}{t^2+81}\left({\displaystyle \frac{\left|x\right|}{28}+\frac{\pi \left|y\right|}{100}+\frac{1}{2}}\right),\left|g(t,x,y)\right|\le \frac{2t}{t^3+93}\left({\displaystyle \frac{5\left|x\right|}{23}+\frac{\pi \left|y\right|}{34}+\frac{7}{8}}\right),}\hfill \end{array}$$

(86)

for all $t\in [1,9]$ and $x,y\in \mathbb{R}$. So, we find $p_1\left(t\right)=\frac{3}{t^2+81}$, $p_2\left(t\right)=\frac{2t}{t^3+93}$, ${\phi }_1(u,v)=\frac{u}{28}+\frac{\pi v}{100}+\frac{1}{2}$, ${\phi }_2(u,v)=\frac{5u}{23}+\frac{\pi v}{34}+\frac{7}{8}$, for all $t\in [1,9]$ and $u,v\in {\mathbb{R}}_{+}$, and then assumption $\left(H5\right)$ is satisfied. In addition, we obtain $\parallel p_1\parallel \approx 0.03658537$ and $\parallel p_2\parallel \approx 0.05155531$. The condition from assumption $\left(H6\right)$ becomes

$${\displaystyle L>\frac{\frac{1}{2}({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)\parallel p_1\parallel +\frac{7}{8}({\mathsf{\Xi }}_3+{\mathsf{\Xi }}_4)\parallel p_2\parallel }{1-(\frac{1}{28}+\frac{\pi }{100})({\mathsf{\Xi }}_1+{\mathsf{\Xi }}_3)\parallel p_1\parallel -(\frac{5}{23}+\frac{\pi }{34})({\mathsf{\Xi }}_2+{\mathsf{\Xi }}_4)\parallel p_2\parallel }\approx 0.48878552,}$$

(87)

So, if $L\ge 0.4888$, then assumption $\left(H6\right)$ is also satisfied. Therefore, by Theorem 6, we deduce the existence of at least one solution $\left(u\right(t),v(t\left)\right),t\in [1,9]$ for problem (77) and (78) with the nonlinearities (85).

5. Conclusions

In this paper, we investigated the existence and uniqueness of solutions for a system of fractional differential equations denoted as (1). These equations are subject to nonlocal boundary conditions as specified in (2). System (1) encompasses Hilfer–Hadamard fractional derivatives that vary in orders and types, while the conditions (2) are nonlocal, featuring a combination of Riemann–Stieltjes integrals and Hadamard derivatives with varying orders. It is worth noting that these conditions are general ones, encompassing scenarios that range from uncoupled boundary conditions (in the event that all functions $K_i$ for $i=1,\dots ,n$ and $P_j$ for $j=1,\dots ,p$ are constants) to more complex cases that generalize multi-point boundary conditions, classical integral conditions, and various combinations thereof. In Section 2, we have provided an existence theorem for the linear fractional differential problem associated with (1) and (2). In Section 3, we have presented our primary findings, supported by rigorous proofs in which we have employed various fixed point theorems. These theorems include the Banach contraction mapping principle (applied to prove Theorem 2), the Krasnosel’skii fixed point theorem for the sum of two operators (utilized in proving Theorems 3 and 4), the Schaefer fixed point theorem (employed for Theorem 5), and the Leray–Schauder nonlinear alternative (used to establish Theorem 6). Finally, in Section 4, we have provided several illustrative examples to elucidate the implications of our main existence results. Going forward, our aim is to investigate different sets of fractional equations, which include fractional derivatives of various kinds and are subject to diverse boundary conditions.