Tissue distribution patterns of solubilized metals from internalized tungsten alloy in the F344 rat

Vernieda B. Vergara; Christy A. Emond; John F. Kalinich; Vernieda B. Vergara; Christy A. Emond; John F. Kalinich

doi:10.3934/environsci.2016.2.290

AIMS Environmental Science

2016, Volume 3, Issue 2: 290-304. doi: 10.3934/environsci.2016.2.290

Previous Article Next Article

Research article Special Issues

Tissue distribution patterns of solubilized metals from internalized tungsten alloy in the F344 rat

1.
Internal Contamination and Metal Toxicity Program, Armed Forces Radiobiology Research Institute, Uniformed Services University, 8901 Wisconsin Avenue, Bethesda, Maryland 20889-5603, USA
2.
Current address: Department of Veterans Affairs, Manchester Veterans Affairs Medical Center, 718 Smyth Road, Manchester, New Hampshire 03104, USA

Received: 01 April 2016 Accepted: 30 May 2016 Published: 01 June 2016

Because of its unique physical and chemical properties, tungsten has been increasingly utilized in a variety of civilian and military applications. This expanded use also raises the risk of human exposure through internalization by various routes. In most cases the toxicological and carcinogenic properties of these tungsten-based compounds are not known nor are the dissolution biokinetics and ultimate fate of the associated metals. Using a laboratory rodent model system designed to assess the health effects of embedded metals, and a tungsten alloy comprised of tungsten (91.1%), nickel (6.0%), and cobalt (2.9%), we investigated the tissue distribution patterns of the metals over a six month period. Despite its perceived insolubility, tungsten rapidly solubilized from the implanted metal fragments, as did nickel and cobalt. All three metals distributed systemically over time with extremely elevated levels of all three metals found in kidney, liver, and spleen. Unexpectedly, tungsten was found to cross the blood-brain and blood-testis barriers and localize in those tissues. These results, along with recent reports suggesting that tungsten is a tumor promoter, raises serious concerns as to the long-term health effects of exposure to tungsten and tungsten-based compounds.

Keywords:

tungsten alloy (WA),
tungsten (W),
nickel (Ni),
cobalt (Co),
wound,
contamination,
rat,
inductively coupled-plasma mass spectrometry (ICP-MS)

Citation: Vernieda B. Vergara, Christy A. Emond, John F. Kalinich. Tissue distribution patterns of solubilized metals from internalized tungsten alloy in the F344 rat[J]. AIMS Environmental Science, 2016, 3(2): 290-304. doi: 10.3934/environsci.2016.2.290

Related Papers:

[1]	Bashir Ahmad, Manal Alnahdi, Sotiris K. Ntouyas, Ahmed Alsaedi . On a mixed nonlinear boundary value problem with the right Caputo fractional derivative and multipoint closed boundary conditions. AIMS Mathematics, 2023, 8(5): 11709-11726. doi: 10.3934/math.2023593
[2]	Donny Passary, Sotiris K. Ntouyas, Jessada Tariboon . Hilfer fractional quantum system with Riemann-Liouville fractional derivatives and integrals in boundary conditions. AIMS Mathematics, 2024, 9(1): 218-239. doi: 10.3934/math.2024013
[3]	Lakhlifa Sadek, Tania A Lazǎr . On Hilfer cotangent fractional derivative and a particular class of fractional problems. AIMS Mathematics, 2023, 8(12): 28334-28352. doi: 10.3934/math.20231450
[4]	Sunisa Theswan, Sotiris K. Ntouyas, Jessada Tariboon . Coupled systems of $ \psi $-Hilfer generalized proportional fractional nonlocal mixed boundary value problems. AIMS Mathematics, 2023, 8(9): 22009-22036. doi: 10.3934/math.20231122
[5]	Ugyen Samdrup Tshering, Ekkarath Thailert, Sotiris K. Ntouyas . Existence and stability results for a coupled system of Hilfer-Hadamard sequential fractional differential equations with multi-point fractional integral boundary conditions. AIMS Mathematics, 2024, 9(9): 25849-25878. doi: 10.3934/math.20241263
[6]	Weerawat Sudsutad, Wicharn Lewkeeratiyutkul, Chatthai Thaiprayoon, Jutarat Kongson . Existence and stability results for impulsive $ (k, \psi) $-Hilfer fractional double integro-differential equation with mixed nonlocal conditions. AIMS Mathematics, 2023, 8(9): 20437-20476. doi: 10.3934/math.20231042
[7]	Ishfaq Mallah, Idris Ahmed, Ali Akgul, Fahd Jarad, Subhash Alha . On $ \psi $-Hilfer generalized proportional fractional operators. AIMS Mathematics, 2022, 7(1): 82-103. doi: 10.3934/math.2022005
[8]	Md. Asaduzzaman, Md. Zulfikar Ali . Existence of positive solution to the boundary value problems for coupled system of nonlinear fractional differential equations. AIMS Mathematics, 2019, 4(3): 880-895. doi: 10.3934/math.2019.3.880
[9]	Saima Rashid, Abdulaziz Garba Ahmad, Fahd Jarad, Ateq Alsaadi . Nonlinear fractional differential equations and their existence via fixed point theory concerning to Hilfer generalized proportional fractional derivative. AIMS Mathematics, 2023, 8(1): 382-403. doi: 10.3934/math.2023018
[10]	Ahmed Alsaedi, Fawziah M. Alotaibi, Bashir Ahmad . Analysis of nonlinear coupled Caputo fractional differential equations with boundary conditions in terms of sum and difference of the governing functions. AIMS Mathematics, 2022, 7(5): 8314-8329. doi: 10.3934/math.2022463

Abstract

1. Introduction

In recent years, the researchers and modelers have shown a keen interest in the topic of fractional differential equations. In fact, such equations appear in the mathematical models of several real-world phenomena occurring in pure, applied and technical sciences, for instance, see the books ^[1,2,3]. Unlike the classical derivative, there do exist many definitions of fractional derivatives and integrals. In ^[4], Hilfer proposed an important definition of fractional derivative (known as Hilfer fractional derivative), which represents both Riemann-Liouville and Caputo fractional derivatives under suitable choice of parameters. Several authors studied initial value problems involving Hilfer fractional derivatives, for example, see ^[5,6,7,8,9]. Some interesting results on boundary value problems involving Hilfer fractional differential equations can be found in the literature. For example, we refer the reader to works on nonlocal Hilfer problems ^[10,11], Hilfer Langevin equations ^[12,13], Hilfer Katugampola operators ^[14], Hilfer Erdelyi-Kober operators ^[15], Hilfer inclusion problems ^[16], Hilfer stochastic differential equations ^[17], $ \psi $-Hilfer problems ^[18], $ \psi $-Hilfer coupled systems ^[19], delay Hilfer fractional differential equations ^[20], Hilfer equations with variable coefficients ^[21], Hilfer sequential fractional differential equations ^[22,23], Hilfer approximate controllability ^[24] and Hilfer-Hadamard boundary value problems ^[25]. A variety of recent results on boundary value problems and coupled systems of Hilfer fractional differential equations and inclusions can be found in the survey paper ^[26].

In ^[27], the authors introduced and developed the existence and uniqueness of solutions for a new class of coupled systems of Hilfer-type fractional differential equations with nonlocal integral boundary conditions of the form

$\begin{equation} \left\{\begin{array}{ll} ^HD^{\alpha , \beta}x(t) = f(t, x(t), y(t)), \; \; \; t\in [a, b], \\ ^HD^{\alpha_1 , \beta_1}y(t) = g(t, x(t), y(t)), \; \; t\in [a, b], \\ x(a) = 0, \quad x(b) = \sum\limits_{i = 1}^{m}\theta_i I^{\varphi_i} y(\xi_i), \; \; y(a) = 0, \quad y(b) = \sum\limits_{j = 1}^{n}\zeta_j I^{\psi_j} x(z_j), \end{array}\right. \end{equation}$ $

(1.1)

where $ ^HD^{\alpha, \beta} $, $ ^HD^{\alpha_1, \beta_1} $ are the Hilfer fractional derivatives of orders $ \alpha $, $ \alpha_1 $, $ 1 < \alpha, \alpha_1 < 2 $, and parameters $ \beta $, $ \beta_1 $, $ 0\leq\beta, \beta_1\leq1 $, respectively, and $ I^{\varphi_i} $, $ I^{\psi_j} $ are the Riemann-Liouville fractional integrals of order $ \varphi_i > 0 $ and $ \psi_j > 0 $, respectively, the points $ \xi_i, z_j \in (a, b), a\geq 0, $ $ f, g: [a, b]\times {\mathbb R}\times {\mathbb R}\to {\mathbb R} $ are continuous functions and $ \theta_i $, $ \zeta_j\in \mathbb{R} $, $ i = 1, 2, \ldots, m $, $ j = 1, 2, \ldots, n $ are given real constants.

Recently, in ^[28], the authors studied a coupled system of $ \psi $-Hilfer fractional order Langevin equations with nonlocal integral boundary conditions given by

$\begin{equation} \left\{ \begin{array}{c} {}^{H}\mathfrak{D}_{a^+}^{\alpha_{1}, \beta_{1}; \psi }\left({}^{H}\mathfrak{D}_{a^+}^{p_{1}, q_{1}; \psi } + \lambda_{1}\right)x(t) = f(t, x(t), y(t)), \quad t\in J: = [a, b], \\ {}^{H}\mathfrak{D}_{a^+}^{\alpha_{2}, \beta_{2}; \psi }\left({}^{H}\mathfrak{D}_{a^+}^{p_{2}, q_{2}; \psi } + \lambda_{2}\right)y(t) = g(t, x(t), y(t)), \quad t\in J: = [a, b], \\ [0.2cm] x(a) = 0, \; \; x(b) = \sum\limits_{i = 1}^{m}\eta_{i}\mathcal{I}_{a^+}^{\delta_{i}; \psi }y(\theta_{i}), \; \; y(a) = 0, \; \; y(b) = \sum\limits_{j = 1}^{n}\mu_{j}\mathcal{I}_{a^+}^{\kappa_{j}; \psi }x(\xi_{j}), \end{array} \right. \end{equation}$ $

(1.2)

where $ {}^{H}\mathfrak{D}_{a^+}^{u, v; \psi } $ is $ \psi $-Hilfer fractional derivatives of order $ u \in \{\alpha_{1}, \alpha_{2}, p_{1}, p_{2}\} $ with $ 0 < u \leq 1 $ and $ v \in \{\beta_{1}, \beta_{2}, q_{1}, q_{2}\} $ with $ 0 \leq v \leq 1 $, $ \mathcal{I}_{a^+}^{w; \psi} $ is $ \psi $-Riemann-Liouville fractional integral of order $ w = \{\delta_{i}, \kappa_{j}\}, $ $ w > 0 $, the points $ \theta_i $, $ \xi_j \in (a, b) $, $ i = 1, 2, \ldots, m $, $ j = 1, 2, \ldots, n $, $ \lambda_{1} $, $ \lambda_{2} \in\mathbb{R} $, $ f $, $ g \in \mathcal{C}([a, b]\times \mathbb{R}^{2}, \mathbb{R}) $ and $ b > a \geq 0 $.

The objective of the present paper is to investigate the existence and uniqueness of solutions for a new class of coupled systems of Langevin type Hilfer fractional differential equations of different orders involving non-integral and autonomous type Riemann-Liouville mixed integral nonlinearities complemented with nonlocal coupled multi-point and Riemann-Liouville integral boundary conditions. This work is motivated by ^[27] and ^[28]. In precise terms, we consider the following problem:

$\begin{equation} \left\{ \begin{array}{ll} ^HD^{\alpha_{1}, \beta_{1}}(^HD^{\alpha_{2}, \beta_{2}}+\lambda_{1} )x(t) = I^{\zeta_1}_{a^+}g_1(x(t), y(t))+f_{1}(t, x(t), y(t)), \; t \in [a, b], \\ ^HD^{\alpha_{3}, \beta_{3}}(^HD^{\alpha_{4}, \beta_{4}}+\lambda_{2} )y(t) = I^{\zeta_2}_{b^-}g_2(x(t), y(t))+f_{2}(t, x(t), y(t)), \; t \in [a, b], \\ x(a) = 0, \, \, x(b) = \sum\limits_{i = 1}^{m}\mu_{i}y(\eta_{i})+\sum\limits_{k = 1}^{n}\nu_{k}I_{a^+}^{q_{k}}y(\xi_{k}), \; q_{k} > 0, \; \\ y(a) = 0, \, \, y(b) = \sum\limits_{i = 1}^{m}\delta_{i}x(\eta_{i})+\sum\limits_{k = 1}^{n}\theta_{k}I_{a^+}^{p_{k}}x(\xi_{k}), \; p_k > 0, \end{array} \right. \end{equation}$ $

(1.3)

where $ ^HD^{\alpha_j, \beta_j} $ represents Hilfer fractional derivative operator of order $ \alpha_j \in (0, 1) $ with parameter $ \beta_j \in [0, 1], $ $ j = 1, 2, 3, 4 $, $ \lambda_{1}, \lambda_{2}, \mu_{i}, \nu_{k}, \delta_{i} $ and $ \theta_{k}, i = 1, 2, ..., m, k = 1, 2, ..., n $ are constants, $ a < \eta_{i}, \xi_{k} < b, $ where $ a\geq 0 $ and $ m, n \in \mathbb{N}, $ $ I_{a^+}^{\zeta_1}, I_{a^+}^{q_{k}}, I_{a^+}^{p_{k}} $ denote the left Riemann-Liouville fractional integral operators of orders $ \zeta_1 > 0, q_{k} > 0, p_{k} > 0 $ respectively, while $ I^{\zeta_2}_{b^-} $ denotes the right Riemann-Liouville fractional integral operator of order $ \zeta_2 > 0, $ and $ f_{1}, f_{2}: [a, b]\times\mathbb{R}\times\mathbb{R}\rightarrow \mathbb{R}, $ $ g_{1}, g_{2}: \mathbb{R}\times\mathbb{R}\rightarrow \mathbb{R} $ are given continuous functions.

Note that problem (1.3) is more general than problem (1.2), since it contains non-integral as well as Riemann-Liouville mixed integral nonlinearities and nonlocal coupled multi-point and Riemann-Liouville integral boundary conditions.

The rest of the paper is organized as follows. In Section 2, we present some necessary material related to our study and prove an auxiliary lemma to define the solution for the problem at hand. Section 3 contains the main results which rely on Banach contraction mapping principle, Krasnosel'ski${\rm{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over i} }}$'s fixed point theorem and Leray-Schauder alternative. In Section 4, we construct examples for the illustration of the results obtained in Section 3.

2. Preliminaries

We begin this section with some basic concepts used in our study.

Definition 2.1. (^[3]) The left and right Riemann–Liouville fractional integrals of order $ \omega > 0 $ for a continuous function $ g, $ existing almost everywhere on $ [a, b], $ are respectively defined by

$ I^{\omega}_{a+} g(t) = \int_{a}^t \frac{(t-s)^{\omega-1}}{\Gamma (\omega)}g(s)ds \, \, \, \, \;{{and}}\; \, \, \, \, I^{\omega}_{b-} g(t) = \int_t^b \frac{(s-t)^{\omega-1}}{\Gamma (\omega)}g(s)ds. $

For the sake of simplicity, we write $ I^{\omega}_{a+} $ and $ I^{\omega}_{b-} $ as $ I^{\omega}_{a} $ and $ I^{\omega}_{b} $ respectively.

Definition 2.2. (^[4]) For $ n-1 < \alpha < n, \; 0\leq\beta\leq 1, $ the Hilfer fractional derivative of order $ \alpha $ and parameter $ \beta $ for a continuous function $ g $ is defined by

$ ^HD^{\alpha, \beta}g(t) = I_{a}^{\beta(n-\alpha)}D^{n}I_{a}^{(1-\beta)(n-\alpha)}g(t), D = \dfrac{d}{dt}, $

where

$ I_{a}^\omega g\left( t \right) = {\rm{ }}\frac{1}{{\Gamma \left( \omega \right)}}\int\limits_a^t {\left( {t - s} \right)^{\omega - 1} g\left( s \right)} ds, \; a\ge 0, $

with $ \omega \in \{\beta(n-\alpha), (1-\beta)(n-\alpha)\}. $

Lemma 2.1. (^[16]) Let $ h \in L(a, b), \; n-1 < \gamma_{1} \leq n, \; n \in \mathbb{N}, \; 0\leq \gamma_{2} \leq1 $ and $ I_a^{(n-\gamma_{1})(1-\gamma_{2})}h \in AC^{k}[a, b] $. Then

$\begin{equation*} I_a^{\gamma_{1}}( ^HD^{\gamma_{1}, \gamma_{2}}h)(t) = h(t)-\sum\limits_{k = 0}^{n-1}\dfrac{(t-a)^{k-(n-\gamma_{1})(1-\gamma_{2})}}{\Gamma(k-(n-\gamma_{1})(1-\gamma_{2})+1)} \lim\limits_{t\longrightarrow a^{+}}\dfrac{d^{k}}{dt^{k}}(I_a^{(1-\gamma_{2})(n-\gamma_{1})}h)(t). \end{equation*}$ $

In the following lemma, we solve the linear variant of the problem (1.3).

Lemma 2.2. Let $ h_{1}, h_{2}:[a, b]\rightarrow \mathbb{R} $ be continuous functions and $ \Delta \ne 0. $ Then the unique solution of the following coupled system:

$\begin{equation} \left\{ \begin{array}{ll} ^HD^{\alpha_{1}, \beta_{1}}(^HD^{\alpha_{2}, \beta_{2}}+\lambda_{1} )x(t) = h_{1}(t), \; \; \; \; \; \; t \in [a, b], \\ ^HD^{\alpha_{3}, \beta_{3}}(^HD^{\alpha_{4}, \beta_{4}}+\lambda_{2} )y(t) = h_{2}(t), \; \; \; \; \; \; t \in [a, b], \\ x(a) = 0, \, \, x(b) = \sum\limits_{i = 1}^{m}\mu_{i}y(\eta_{i})+\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}}y(\xi_{k}), \; \; q_{k} > 0, \; \\ y(a) = 0, \, \, y(b) = \sum\limits_{i = 1}^{m}\delta_{i}x(\eta_{i})+\sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}}x(\xi_{k}), \; \; \; p_{k} > 0, \end{array} \right. \end{equation}$ $

(2.1)

is given by

$\begin{eqnarray} x(t)& = &I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(t)-\lambda_{1} I_{a}^{\alpha_{2}}x(t)+\dfrac{(t-a)^{\alpha_{2}+\epsilon_{1}-1}}{\Delta\Gamma(\alpha_{2}+\epsilon_{1})}\Big\{\Omega_{2}\Big(\lambda_{1} I_{a}^{\alpha_{2}}x(b)-I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(b) \\ &&-\lambda_{2} \sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{4}}y(\eta_{i}) +\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(\eta_{i})+\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}} h_{2}(\xi_{k})\\ &&-\lambda_{2} \sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{4}}y(\xi_{k})\Big)+\Omega_{4}\Big(\lambda_{2} I_{a}^{\alpha_{4}}y(b)-I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(b)+\sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(\eta_{i})\\ &&-\lambda_{1} \sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{2}}x(\eta_{i}) +\sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}h_{1}(\xi_{k})-\lambda_{1} \sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{2}}x(\xi_{k})\Big) \Big\}, \end{eqnarray}$ $

(2.2)

$\begin{eqnarray} y(t)& = &I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(t)-\lambda_{2} I_{a}^{\alpha_{4}}y(t)+\dfrac{(t-a)^{\alpha_{4}+\epsilon_{3}-1}}{\Delta\Gamma(\alpha_{4}+\epsilon_{3})}\Big\{\Omega_{3}\Big(\lambda_{1} I_{a}^{\alpha_{2}}x(b)-I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(b)\\ &&-\lambda_{2} \sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{4}}y(\eta_{i}) +\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(\eta_{i}) +\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}} h_{2}(\xi_{k})\\ &&-\lambda_{2} \sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{4}}y(\xi_{k})\Big)+\Omega_{1}\Big(-I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(b)+\lambda_{2} I_{a}^{\alpha_{4}}y(b) +\sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(\eta_{i})\\ &&-\lambda_{1} \sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{2}}x(\eta_{i}) +\sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}h_{1}(\xi_{k}) -\lambda_{1} \sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{2}}x(\xi_{k}) \Big)\Big\}, \end{eqnarray}$ $

(2.3)

where $ \Delta, \Omega_{i}, \; i = 1, 2, 3, 4 $ are given by

$\begin{eqnarray} \Omega_{1} & = &\dfrac{(b-a)^{\alpha_{2}+\epsilon_{1}-1}}{\Gamma(\alpha_{2}+\epsilon_{1})}, \; \; \Omega_{2} = \dfrac{(b-a)^{\alpha_{4}+\epsilon_{3}-1}}{\Gamma(\alpha_{4}+\epsilon_{3})}, \\ \Omega_{3}& = &\sum\limits_{i = 1}^{m}\delta_{i} \dfrac{(\eta_{i}-a)^{\alpha_{2}+\epsilon_{1}-1}}{\Gamma(\alpha_{2}+\epsilon_{1})}+ \sum\limits_{k = 1}^{n}\theta_{k}\dfrac{(\xi_{k}-a)^{p_{k}+\alpha_{2}+\epsilon_{1}-1}}{\Gamma(p_{k}+\alpha_{2}+\epsilon_{1})}, \\ \Omega_{4}& = &\sum\limits_{i = 1}^{m}\mu_{i}\dfrac{(\eta_{i}-a)^{\alpha_{4}+\epsilon_{3}-1}}{\Gamma(\alpha_{4}+\epsilon_{3})}+ \sum\limits_{k = 1}^{n}\nu_{k}\dfrac{(\xi_{k}-a)^{q_{k}+\alpha_{4}+\epsilon_{3}-1}}{\Gamma(q_{k}+\alpha_{4}+\epsilon_{3})}, \\ \Delta & = & \Omega_{1}\Omega_{2}-\Omega_{3}\Omega_{4}, \end{eqnarray}$ $

(2.4)

and $ \epsilon_{i} = \alpha_{i}+\beta_{i}-\alpha_{i}\beta_{i}, \; i = 1, 2, 3, 4 $.

Proof. Applying the integral operators $ I_{a}^{\alpha_{1}} $ and $ I_{a}^{\alpha_{3}} $ on the first and second Hilfer fractional differential equations in (2.1) respectively and using Lemma 2.1, we obtain

$\begin{equation} (^HD^{\alpha_{2}, \beta_{2}}+\lambda_{1} )x(t)-\dfrac{c_{0}(t-a)^{\epsilon_{1}-1}}{\Gamma(\epsilon_{1})} = I_{a}^{\alpha_{1}}h_{1}(t), \end{equation}$ $

(2.5)

$\begin{equation} (^HD^{\alpha_{4}, \beta_{4}}+\lambda_{2} )y(t)-\dfrac{d_{0}(t-a)^{\epsilon_{3}-1}}{\Gamma(\epsilon_{3})} = I_{a}^{\alpha_{3}}h_{2}(t). \end{equation}$ $

(2.6)

Now operating $ I_{a}^{\alpha_{2}} $ and $ I_{a}^{\alpha_{4}} $ respectively to the Eqs (2.5) and (2.6), we get

$\begin{eqnarray} x(t)+\lambda_{1} I_{a}^{\alpha_{2}}x(t)-\dfrac{c_{1}(t-a)^{\epsilon_{2}-1}}{\Gamma(\epsilon_{2})}-\dfrac{c_{0}(t-a)^{\alpha_{2}+\epsilon_{1}-1}}{\Gamma(\alpha_{2}+\epsilon_{1})} = I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(t), \end{eqnarray}$ $

(2.7)

$\begin{eqnarray} y(t)+\lambda_{2} I_{a}^{\alpha_{4}}y(t)-\dfrac{d_{1}(t-a)^{\epsilon_{4}-1}}{\Gamma(\epsilon_{4})}-\dfrac{d_{0}(t-a)^{\alpha_{4}+\epsilon_{3}-1}}{\Gamma(\alpha_{4}+\epsilon_{3})} = I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(t). \end{eqnarray}$ $

(2.8)

Using the conditions $ x(a) = 0 $ and $ y(a) = 0 $ in (2.7) and (2.8) respectively, we find that $ c_{1} = d_{1} = 0. $ Thus we have

$\begin{eqnarray} x(t) = I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(t)-\lambda_{1} I_{a}^{\alpha_{2}}x(t)+\dfrac{c_{0}(t-a)^{\alpha_{2}+\epsilon_{1}-1}}{\Gamma(\alpha_{2}+\epsilon_{1})}, \end{eqnarray}$ $

(2.9)

$\begin{eqnarray} y(t) = I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(t)-\lambda_{2} I_{a}^{\alpha_{4}}y(t)+\dfrac{d_{0}(t-a)^{\alpha_{4}+\epsilon_{3}-1}}{\Gamma(\alpha_{4}+\epsilon_{3})}. \end{eqnarray}$ $

(2.10)

Inserting (2.9) and (2.10) in the condition $ x(b) = \sum\limits_{i = 1}^{m}\mu_{i}y(\eta_{i})+\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}}y(\xi_{k}), $ we find that

$\begin{eqnarray*} I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(b)-\lambda_{1} I_{a}^{\alpha_{2}}x(b)+\dfrac{c_{0}(b-a)^{\alpha_{2}+\epsilon_{1}-1}}{\Gamma(\alpha_{2}+\epsilon_{1})} = \sum\limits_{i = 1}^{m}\mu_{i} \Big\{I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(\eta_{i}) -\lambda_{2} I_{a}^{\alpha_{4}}y(\eta_{i})\\ +\dfrac{d_{0}(\eta_{i}-a)^{\alpha_{4}+\epsilon_{3}-1}}{\Gamma(\alpha_{4}+\epsilon_{3})}\Big\} +\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}} \Big\{I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(\xi_{k})-\lambda_{2} I_{a}^{\alpha_{4}}y(\xi_{k})+\dfrac{d_{0}(\xi_{k}-a)^{\alpha_{4}+\epsilon_{3}-1}}{\Gamma(\alpha_{4}+\epsilon_{3})}\Big\}, \end{eqnarray*}$ $

which can alternatively be written as

$\begin{eqnarray} &&c_{0}\dfrac{(b-a)^{\alpha_{2}+\epsilon_{1}-1}}{\Gamma(\alpha_{2}+\epsilon_{1})}- d_{0}\Big\{\sum\limits_{i = 1}^{m}\mu_{i}\dfrac{(\eta_{i}-a)^{\alpha_{4}+\epsilon_{3}-1}}{\Gamma(\alpha_{4}+\epsilon_{3})}+ \sum\limits_{k = 1}^{n}\nu_{k}\dfrac{(\xi_{k}-a)^{q_{k}+\alpha_{4}+\epsilon_{3}-1}}{\Gamma(q_{k}+\alpha_{4}+\epsilon_{3})}\Big \} \\ & = &\lambda_{1} I_{a}^{\alpha_{2}}x(b)-I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(b)-\lambda_{2} \sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{4}}y(\eta_{i}) +\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(\eta_{i})\\ &&+\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}} h_{2}(\xi_{k})-\lambda_{2} \sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{4}}y(\xi_{k}). \end{eqnarray}$ $

(2.11)

In a similar manner, making use of (2.9) and (2.10) in the condition: $ y(b) = \sum_{i = 1}^{m}\delta_{i}x(\eta_{i})+\sum_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}}x(\xi_{k}), $ leads to

$\begin{eqnarray} &&-c_{0}\Big\{\sum\limits_{i = 1}^{m}\delta_{i} \dfrac{(\eta_{i}-a)^{\alpha_{2}+\epsilon_{1}-1}}{\Gamma(\alpha_{2}+\epsilon_{1})}+ \sum\limits_{k = 1}^{n}\theta_{k}\dfrac{(\xi_{k}-a)^{p_{k}+\alpha_{2}+\epsilon_{1}-1}}{\Gamma(p_{k}+\alpha_{2}+\epsilon_{1})}\Big\} +d_{0}\dfrac{(b-a)^{\alpha_{4}+\epsilon_{3}-1}}{\Gamma(\alpha_{4}+\epsilon_{3})}\\ && = \lambda_{2} I_{a}^{\alpha_{4}}y(b) -I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(b)+\sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(\eta_{i})-\lambda_{1} \sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{2}}x(\eta_{i}) \\ &&+\sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}h_{1}(\xi_{k})-\lambda_{1} \sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{2}}x(\xi_{k}). \end{eqnarray}$ $

(2.12)

Making use of the notation in (2.4), we can write (2.11) and (2.12) as

$\begin{eqnarray*} \Omega_{1} c_{0} -\Omega_{4}d_{0}& = &\lambda_{1} I_{a}^{\alpha_{2}}x(b)-I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(b)-\lambda_{2} \sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{4}}y(\eta_{i}) +\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(\eta_{i})\\ &&+\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}} h_{2}(\xi_{k})-\lambda_{2} \sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{4}}y(\xi_{k}), \\ -\Omega_{3}c_{0} +\Omega_{2}d_{0}& = &\lambda_{2} I_{a}^{\alpha_{4}}y(b)-I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(b) +\sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(\eta_{i})-\lambda_{1} \sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{2}}x(\eta_{i}) \\ &&+\sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}h_{1}(\xi_{k})-\lambda_{1} \sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{2}}x(\xi_{k}), \end{eqnarray*}$ $

which, on solving for $ c_{0} $ and $ d_{0}, $ yields

$\begin{eqnarray*} c_{0}& = &\dfrac{1}{\Delta}\Big\{\Omega_{2}\Big(\lambda_{1} I_{a}^{\alpha_{2}}x(b)-I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(b)-\lambda_{2} \sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{4}}y(\eta_{i}) +\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(\eta_{i})\\ &&+\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}} h_{2}(\xi_{k})-\lambda_{2} \sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{4}}y(\xi_{k})\Big)+\Omega_{4}\Big(\lambda_{2} I_{a}^{\alpha_{4}}y(b)-I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(b)\\ &&+\sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(\eta_{i})-\lambda_{1} \sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{2}}x(\eta_{i}) +\sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}h_{1}(\xi_{k})\\ &&-\lambda_{1} \sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{2}}x(\xi_{k})\Big) \Big\}, \\ d_{0}& = &\dfrac{1}{\Delta}\Big\{\Omega_{3}\Big(\lambda_{1} I_{a}^{\alpha_{2}}x(b)-I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(b)-\lambda_{2} \sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{4}}y(\eta_{i}) +\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(\eta_{i})\\ &&+\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}} h_{2}(\xi_{k})-\lambda_{2} \sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{4}}y(\xi_{k})\Big)+\Omega_{1}\Big(\lambda_{2} I_{a}^{\alpha_{4}}y(b) -I_{a}^{\alpha_{3}+\alpha_{4}}h_{2}(b)\\&&+\sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{1}+\alpha_{2}}h_{1}(\eta_{i})-\lambda_{1} \sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{2}}x(\eta_{i}) +\sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}h_{1}(\xi_{k})\\&&-\lambda_{1} \sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{2}}x(\xi_{k})\Big)\Big\}. \end{eqnarray*}$ $

Substituting the values of $ c_0 $ and $ d_0 $ in (2.9) and (2.10) respectively together with (2.4), we get the solution (2.2) and (2.3). By direct computation, one can obtain the converse of this lemma. The proof is finished.

3. Existence and uniqueness results

Let $ \mathcal{X} = C([a, b], \mathbb{R}) $ denote the Banach space of all continuous functions from $ [a, b] $ to $ \mathbb{R} $ with the norm $ \|x\| = \sup_{t \in [a, b]}|x(t)|. $ Then the product space $ (\mathcal{X} \times \mathcal{X}, \|\cdot\|) $ is also a Banach space endowed with the norm $ \|(x, y)\| = \|x\|+\|y\| $ for $ (x, y) \in \mathcal{X} \times \mathcal{X} $.

In view of Lemma 2.2, we introduce an operator $ \mathcal{T}:\mathcal{X}\times \mathcal{X} \rightarrow \mathcal{X} \times \mathcal{X} $ as

$\begin{equation*} \mathcal{T}(x, y)(t) = \left( \begin{array}{ll} \mathcal{T}_{1}(x, y)(t) \\ \mathcal{T}_{2}(x, y)(t) \end{array} \right), \end{equation*}$ $

where

$\begin{eqnarray} \mathcal{T}_{1}(x, y)(t)& = & I_{a}^{\alpha_{1}+\alpha_{2}+\zeta_1} g_1(x(t), y(t)) + I_{a}^{\alpha_{1}+\alpha_{2}}f_{1}(t, x(t), y(t))-\lambda_{1} I_{a}^{\alpha_{2}}x(t)\\ && +\dfrac{(t-a)^{\alpha_{2}+\epsilon_{1}-1}}{\Delta\Gamma(\alpha_{2}+\epsilon_{1})}\times \Big\{\Omega_{2}\Big(\lambda_{1} I_{a}^{\alpha_{2}}x(b) -I_{a}^{\alpha_{1}+\alpha_{2}+\zeta_1} g_1(x(b), y(b))\\ && - I_{a}^{\alpha_{1}+\alpha_{2}}f_{1}(b, x(b), y(b)) -\lambda_{2} \sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{4}}y(\eta_{i})\\ &&+\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}g_2(x(\eta_{i}), y(\eta_{i}))) +\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}f_{2}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))\\ && +\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}g_2(x(\xi_{k}), y(\xi_{k}))) +\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}} f_{2}(\xi_{k}, x(\xi_{k}), y(\xi_{k})) \\ && -\lambda_{2} \sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{4}}y(\xi_{k})\Big)+\Omega_{4}\Big(\lambda_{2} I_{a}^{\alpha_{4}}y(b) -I_{a}^{\alpha_{3}+\alpha_{4}}f_{2}(b, x(b), y(b))\\ &&+\sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{1}+\alpha_{2}+\zeta_1} g_1(x(\eta_{i}), y(\eta_{i}))+\sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{1}+\alpha_{2}}f_{1}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))\\ && -\lambda_{1} \sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{2}}x(\eta_{i}) +\sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_1}g_{1}(x(\xi_{k}), y(\xi_{k}))\\ &&+\sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}f_{1}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))-\lambda_{1} \sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{2}}x(\xi_{k}) \Big)\Big\}, \end{eqnarray}$ $

(3.1)

and

$\begin{eqnarray} \mathcal{T}_{2}(x, y)(t)& = & I_{a}^{\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}g_2(x(t), y(t)))+ I_{a}^{\alpha_{3}+\alpha_{4}}f_{2}(t, x(t), y(t)) -\lambda_{2} I_{a}^{\alpha_{4}}y(t)\\ && +\dfrac{(t-a)^{\alpha_{4}+\epsilon_{3}-1}}{\Delta\Gamma(\alpha_{4}+\epsilon_{3})}\times \Big\{\Omega_{3}\Big(\lambda_{1} I_{a}^{\alpha_{2}}x(b) -I_{a}^{\alpha_{1}+\alpha_{2}+\zeta_1}g_{1}(x(b), y(b)) \\ &&-I_{a}^{\alpha_{1}+\alpha_{2}}f_{1}(b, x(b), y(b))-\lambda_{2} \sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{4}}y(\eta_{i}) \\ && +\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}g_2(x(\eta_{i}), y(\eta_{i}))) +\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}f_{2}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))\\ && +\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}g_2(x(\xi_{k}), y(\xi_{k}))) +\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}} f_{2}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))\\ &&-\lambda_{2} \sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{4}}y(\xi_{k})\Big)+\Omega_{1}\Big(\lambda_{2} I_{a}^{\alpha_{4}}y(b) -I_{a}^{\alpha_{3}+\alpha_{4}}f_{2}(b, x(b), y(b))\\ &&+\sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{1}+\alpha_{2}+\zeta_1}g_{1}(x(\eta_{i}), y(\eta_{i})) +\sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{1}+\alpha_{2}}f_{1}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))\\ && -\lambda_{1} \sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{2}}x(\eta_{i})+\sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_1}g_{1}(x(\xi_{k}), y(\xi_{k}))\\ && +\sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}f_{1}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))-\lambda_{1} \sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{2}}x(\xi_{k}) \Big)\Big\}. \end{eqnarray}$ $

(3.2)

For computational facilitation, we set

$\begin{eqnarray} \sigma_{1}& = & \dfrac{(b-a)^{\alpha_{1}+\alpha_{2}}}{\Gamma (\alpha_{1}+\alpha_{2}+1)}+\dfrac{\Omega_{1}}{|\Delta|} \Big\{\dfrac{\Omega_{2}(b-a)^{\alpha_{1}+\alpha_{2}}}{\Gamma (\alpha_{1}+\alpha_{2}+1)}+|\Omega_{4}|\Big( \sum\limits_{i = 1}^{m}\dfrac{|\delta_{i}|(\eta_{i}-a)^{\alpha_{1}+\alpha_{2}}}{\Gamma(\alpha_{1}+\alpha_{2}+1)} \\&& +\sum\limits_{k = 1}^{n}\dfrac{|\theta_{k}|(\xi_{k}-a)^{p_{k}+\alpha_{1}+\alpha_{2}}}{\Gamma(p_{k}+\alpha_{1}+\alpha_{2}+1)}\Big)\Big\}, \end{eqnarray}$ $

(3.3)

$\begin{eqnarray} \sigma_{2}& = &\dfrac{\Omega_{1}}{|\Delta|}\Big\{\Omega_{2}\Big( \sum\limits_{i = 1}^{m}|\dfrac{\mu_{i}|(\eta_{i}-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma(\alpha_{3}+\alpha_{4}+1)} + \sum\limits_{k = 1}^{n}\dfrac{|\nu_{k}|(\xi_{k}-a)^{q_{k}+\alpha_{3}+\alpha_{4}}}{\Gamma(q_{k}+\alpha_{3}+\alpha_{4}+1)}\Big) \\&& + \dfrac{|\Omega_{4}|(b-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma (\alpha_{3}+\alpha_{4}+1)} \Big\}, \end{eqnarray}$ $

(3.4)

$\begin{eqnarray} \sigma_{3}& = &\dfrac{\Omega_{2}}{|\Delta|}\Big\{\dfrac{|\Omega_{3}|(b-a)^{\alpha_{1}+\alpha_{2}}}{\Gamma (\alpha_{1}+\alpha_{2}+1)}+\Omega_{1}\Big( \sum\limits_{i = 1}^{m}\dfrac{|\delta_{i}|(\eta_{i}-a)^{\alpha_{1}+\alpha_{2}}}{\Gamma(\alpha_{1}+\alpha_{2}+1)}\\ &&+ \sum\limits_{k = 1}^{n}\dfrac{|\theta_{k}|(\xi_{k}-a)^{p_{k}+\alpha_{1}+\alpha_{2}}}{\Gamma(p_{k}+\alpha_{1}+\alpha_{2}+1)}\Big)\Big\}, \end{eqnarray}$ $

(3.5)

$\begin{eqnarray} \sigma_{4}& = &\dfrac{(b-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma (\alpha_{3}+\alpha_{4}+1)} + \dfrac{\Omega_{2}}{| \Delta|} \Big\{|\Omega_{3} |\Big( \sum\limits_{i = 1}^{m} \dfrac{|\mu_{i}|(\eta_{i}-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma(\alpha_{3}+\alpha_{4}+1)}\\ &&+ \sum\limits_{k = 1}^{n}\dfrac{|\nu_{k}|(\xi_{k}-a)^{q_{k}+\alpha_{3}+\alpha_{4}}}{\Gamma(q_{k}+\alpha_{3}+\alpha_{4}+1)}\Big)+\dfrac{\Omega_{1}(b-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma (\alpha_{3}+\alpha_{4}+1)}\Big \}, \end{eqnarray}$ $

(3.6)

$\begin{eqnarray} \sigma_{5} & = & \dfrac{1}{|\Delta|} \Big\{ \dfrac{|\lambda_{1}|(|\Delta|+\Omega_{2}\Omega_{1})(b-a)^{\alpha_{2}}}{\Gamma(\alpha_{2}+1)}+ \dfrac{\Omega_{1}|\lambda_{2}\Omega_{4}|(b-a)^{\alpha_{4}}}{\Gamma(\alpha_{4}+1)} \\ &&+ \Omega_{1}|\lambda_{1}\Omega_{4}| \Big(\sum\limits_{i = 1}^{m}\dfrac{|\delta_{i}|(\eta_{i}-a)^{\alpha_{2}}}{\Gamma(\alpha_{2}+1)} + \sum\limits_{k = 1}^{n}\dfrac{|\theta_{k}|(\xi_{k}-a)^{p_{k}+\alpha_{2}}}{\Gamma(p_{k}+\alpha_{2}+1)}\Big) \\ &&+ |\lambda_{2}|\Omega_{2}\Omega_{1}\Big( \sum\limits_{k = 1}^{n}\dfrac{|\nu_{k}|(\xi_{k}-a)^{q_{k}+\alpha_{4}}}{\Gamma(q_{k}+\alpha_{4}+1)}+\sum\limits_{i = 1}^{m}\dfrac{|\mu_{i}|(\eta_{i}-a)^{\alpha_{4}}}{ \Gamma(\alpha_{4}+1)}\Big)\Big\}, \end{eqnarray}$ $

(3.7)

$\begin{eqnarray} \sigma_{6}& = &\dfrac{1}{|\Delta|} \Big\{ \dfrac{\Omega_{2}|\lambda_{1}\Omega_{3}|(b-a)^{\alpha_{2}}}{\Gamma(\alpha_{2}+1)}+ \dfrac{|\lambda_{2}|(|\Delta|+\Omega_{1}\Omega_{2})(b-a)^{\alpha_{4}}}{\Gamma(\alpha_{4}+1)}\\ &&+ |\lambda_{1}|\Omega_{1}\Omega_{2}\Big( \sum\limits_{i = 1}^{m}\dfrac{|\delta_{i}|(\eta_{i}-a)^{\alpha_{2}}}{\Gamma(\alpha_{2}+1)}+\sum\limits_{k = 1}^{n}\dfrac{|\theta_{k}|(\xi_{k}-a)^{p_{k}+\alpha_{2}}}{\Gamma(p_{k}+\alpha_{2}+1)}\Big)\\ &&+ \Omega_{2}|\lambda_{2}\Omega_{3}| \Big(\sum\limits_{i = 1}^{m}\dfrac{|\mu_{i}|(\eta_{i}-a)^{\alpha_{4}}}{\Gamma(\alpha_{4}+1)} +\sum\limits_{k = 1}^{n}\dfrac{|\nu_{k}|(\xi_{k}-a)^{q_{k}+\alpha_{4}}}{\Gamma(q_{k}+\alpha_{4}+1)}\Big)\Big\}, \end{eqnarray}$ $

(3.8)

$\begin{eqnarray} \sigma_{7} & = & \dfrac{(b-a)^{\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma (\alpha_{1}+\alpha_{2}+\zeta_{1}+1)}+\dfrac{\Omega_{1}}{|\Delta|} \Big\{\dfrac{\Omega_{2}(b-a)^{\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma (\alpha_{1}+\alpha_{2}+\zeta_{1}+1)}\\ &&+|\Omega_{4}|\Big(\sum\limits_{i = 1}^{m}\dfrac{|\delta_{i}|(\eta_{i}-a)^{\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma(\alpha_{1}+\alpha_{2}+\zeta_{1}+1)}+ \sum\limits_{k = 1}^{n}\dfrac{|\theta_{k}|(\xi_{k}-a)^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma(p_{k}+\alpha_{1}+\alpha_{2}+\zeta_{1}+1)}\Big)\Big\}, \end{eqnarray}$ $

(3.9)

$\begin{eqnarray} \sigma_{8}& = &\dfrac{\Omega_{1}\Omega_{2}(b-a)^{\zeta_{2}}}{|\Delta|\Gamma(\zeta_{2}+1)}\Big(\sum\limits_{i = 1}^{m}\dfrac{|\mu_{i}|(\eta_{i}-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma(\alpha_{3}+\alpha_{4}+1)}+ \sum\limits_{k = 1}^{n}\dfrac{|\nu_{k}|(\xi_{k}-a)^{q_{k}+\alpha_{3}+\alpha_{4}}}{\Gamma(q_{k}+\alpha_{3}+\alpha_{4}+1)}\Big), \end{eqnarray}$ $

(3.10)

$\begin{eqnarray} \sigma_{9}& = & \dfrac{\Omega_{2}}{|\Delta|}\Big\{\dfrac{|\Omega_{3}|(b-a)^{\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma (\alpha_{1}+\alpha_{2}+\zeta_{1}+1)}+\Omega_{1}\Big( \sum\limits_{i = 1}^{m}\dfrac{|\delta_{i}|(\eta_{i}-a)^{\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma(\alpha_{1}+\alpha_{2}+\zeta_{1}+1)}\\ &&+ \sum\limits_{k = 1}^{n}\dfrac{|\theta_{k}|(\xi_{k}-a)^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma(p_{k}+\alpha_{1}+\alpha_{2}+\zeta_{1}+1)}\Big)\Big\}, \end{eqnarray}$ $

(3.11)

$\begin{eqnarray} \sigma_{10}& = & \dfrac{(b-a)^{\zeta_{2}}}{\Gamma(\zeta_{2}+1)}\Big(\dfrac{(b-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma (\alpha_{3}+\alpha_{4}+1)} + \dfrac{\Omega_{2}|\Omega_{3} |}{| \Delta|} \Big( \sum\limits_{i = 1}^{m}\dfrac{|\mu_{i}|(\eta_{i}-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma(\alpha_{3}+\alpha_{4}+1)}\\ &&+ \sum\limits_{k = 1}^{n}\dfrac{|\nu_{k}|(\xi_{k}-a)^{q_{k}+\alpha_{3}+\alpha_{4}}}{\Gamma(q_{k}+\alpha_{3}+\alpha_{4}+1)}\Big)\Big). \end{eqnarray}$ $

(3.12)

In the sequel, we suppose that $ f_{1}, f_{2} :[a, b]\times {\mathbb R}\times {\mathbb R}\to {\mathbb R} $ and $ g_{1}, g_{2}: \mathbb{R}\times\mathbb{R}\rightarrow \mathbb{R} $ are continuous functions satisfying the following assumptions:

$ (H_{1}) $ $ \forall (x_{1}, y_{1}), (x_{2}, y_{2}) \in \mathbb{R}^2, $ there exist positive real constants $ K_{i}, $ i = 1, 2, such that

$\begin{eqnarray*} |f_{1}(t, x_{1}, y_{1})-f_{1}(t, x_{2}, y_{2})|\leq K_{1}(|x_{1}-x_{2}|+|y_{1}-y_{2}|), \\ |f_{2}(t, x_{1}, y_{1})-f_{2}(t, x_{2}, y_{2})| \leq K_{2}(|x_{1}-x_{2}|+|y_{1}-y_{2}|); \end{eqnarray*}$ $

$ (H_{2}) \; \forall (x_{1}, y_{1}), (x_{2}, y_{2}) \in \mathbb{R}^2, $ there exist positive real constants $ L_{i}, $ i = 1, 2, such that

$\begin{eqnarray*} |g_{1}(x_{1}, y_{1})-g_{1}(x_{2}, y_{2})|\leq L_{1}(|x_{1}-x_{2}|+|y_{1}-y_{2}|), \\ |g_{2}(x_{1}, y_{1})-g_{2}(x_{2}, y_{2})| \leq L_{2}(|x_{1}-x_{2}|+|y_{1}-y_{2}|); \end{eqnarray*}$ $

$ (H_{3}) $ We can find real constants $ u_{k}, v_{k}, \omega_{k}, \tau_{k}\geqslant 0, k = 0, 1, 2 $ with $ u_{0}, v_{0}, \omega_{0}, \tau_{0}\neq 0 $ such that

$\begin{eqnarray*} |f_{1}(t, x, y)|&\leq& u_{0}+u_{1}|x|+u_{2}|y|, |f_{2}(t, x, y)|\leq v_{0}+v_{1}|x|+v_{2}|y|, \\ |g_{1}(x, y)| &\leq& \omega_{0}+\omega_{1}|x|+\omega_{2}|y|, |g_{2}(x, y)|\leq \tau_{0}+\tau_{1}|x|+\tau_{2}|y|; \end{eqnarray*}$ $

$ (H_{4}) $ There exist nonnegative functions $ \phi _{1}, \phi _{2} \in C([a, b], {\mathbb R}^{+}), $ and positive constants $ \Lambda_{1}, \Lambda_{2} $ such that

$ |f_{1}(t, x, y)|\le \phi_{1}(t) $, $ |f_{2}(t, x, y)|\le \phi_{2}(t) $, $ |g_1(x, y)|\leqslant \Lambda_{1} $, $ |g_2(x, y)|\leqslant \Lambda_{2} $ for all $ (t, x, y)\in [a, b] \times {\mathbb R}\times {\mathbb R}. $

Now we present our first main result dealing with the uniqueness of solutions for the system (1.3), which relies on Banach contraction mapping principle ^[29].

Theorem 3.1. Assume that conditions ($ H_{1}) $ and $ (H_{2} $) hold. Then the system (1.3) has a unique solution on $ [a, b] $ provided that

$\begin{equation} (\sigma_{1}+\sigma_{3})K_{1}+(\sigma_{4}+\sigma_{2})K_{2}+(\sigma_{9}+\sigma_{7})L_{1}+(\sigma_{10}+\sigma_{8})L_{2} +\sigma_{5}+\sigma_{6} < 1, \end{equation}$ $

(3.13)

where $ \sigma_1, \dots, \sigma_{10} $ are given in (3.3)–(3.12).

Proof. Let us fix $ \sup_{t\in[a, b]}|f_{i}(t, 0, 0)| = M_{i} < \infty, |g_{i}(0, 0)| = 0, \; i = 1, 2. $ In order to satisfy the hypotheses of Banach contraction mapping principle, we first show that $ \mathcal{T}B_{\rho}\subset B_{\rho}, $ where $ B_{\rho} $ is a closed bounded ball $ B_{\rho}\subset \mathcal{X} \times \mathcal{X} $ defined by

$\begin{equation*} B_{\rho} = \{(x, y) \in \mathcal{X} \times \mathcal{X} :\|(x, y)\|\leq \rho \}, \end{equation*}$ $

with

$\begin{equation} \rho\geq \dfrac{M_{1}(\sigma_{1}+\sigma_{3})+M_{2}(\sigma_{2}+\sigma_{4})}{1-[K_{1}(\sigma_{1}+\sigma_{3})+K_{2}(\sigma_{2}+\sigma_{4})+L_{1}(\sigma_{7}+\sigma_{9})+L_{2}(\sigma_{8}+\sigma_{10})+\sigma_{5}+\sigma_{6}]}. \end{equation}$ $

(3.14)

For an arbitrary element $ (x, y) \in B_{\rho} $ and for each $ t \in [a, b], $ we have

$\begin{eqnarray*} |\mathcal{T}_{1}(x, y)(t)| &\leq& I_a^{\alpha_{1}+\alpha_{2}+\zeta_1}(|g_1(x(t), y(t)) -g_{1}(0, 0)|+|g_{1}(0, 0)|)\\\nonumber && + I_{a}^{\alpha_{1}+\alpha_{2}}(|f_{1}(t, x(t), y(t))-f_{1}(t, 0, 0)|+|f_{1}(t, 0, 0)|)\\\nonumber && +|\lambda_{1}| I_{a}^{\alpha_{2}}|x(t)| +\dfrac{(b-a)^{\alpha_{2}+\epsilon_{1}-1}}{|\Delta|\Gamma(\alpha_{2}+\epsilon_{1})} \Big\{\Omega_{2} \Big(|\lambda_{1}| I_{a}^{\alpha_{2}}|x(b)|\\\nonumber && +I_a^{\alpha_{1}+\alpha_{2}+\zeta_1} (|g_1(x(b), y(b)) -g_{1}(0, 0)|+|g_{1}(0, 0)|)\\ \nonumber && +I_{a}^{\alpha_{1}+\alpha_{2}}(|f_{1}(b, x(b), y(b))-f_{1}(b, 0, 0)|+|f_{1}(b, 0, 0)|) +|\lambda_{2}| \sum\limits_{i = 1}^{m}|\mu_{i}| I_{a}^{\alpha_{4}}|y(\eta_{i})|\\\nonumber &&+\sum\limits_{i = 1}^{m}|\mu_{i}| I_{a}^{\alpha_{3}+\alpha_{4}}I^{\zeta_2}_{b}(|g_2(x(\eta_{i}), y(\eta_{i})) -g_{2}(0, 0)|+|g_{2}(0, 0)|)\\\nonumber &&+\sum\limits_{i = 1}^{m}|\mu_{i}| I_{a}^{\alpha_{3}+\alpha_{4}}(|f_{2}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))-f_{2}(\eta_{i}, 0, 0)|+|f_{2}(\eta_{i}, 0, 0)|)\\\nonumber && +\sum\limits_{k = 1}^{n}|\nu_{k}|I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}}I^{\zeta_2}_{b}(|g_2(x(\xi_{k}), y(\xi_{k})) -g_{2}(0, 0)|+|g_{2}(0, 0)|)\\\nonumber &&+\sum\limits_{k = 1}^{n}|\nu_{k}|I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}}(| f_{2}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))-f_{2}(\xi_{k}, 0, 0)|+|f_{2}(\xi_{k}, 0, 0)|) \\\nonumber &&+|\lambda_{2} |\sum\limits_{k = 1}^{n}|\nu_{k}|I_{a}^{q_{k}+\alpha_{4}}|y(\xi_{k})|\Big) +\Omega_{4}\Big(|\lambda_{2} |I_{a}^{\alpha_{4}}|y(b)|\\ \nonumber &&+I_{a}^{\alpha_{3}+\alpha_{4}}(| f_{2}(b, x(b), y(b))-f_{2}(b, 0, 0)|+|f_{2}(b, 0, 0)|) \\\nonumber &&+\sum\limits_{i = 1}^{m}|\delta_{i}|I_a^{\alpha_{1}+\alpha_{2}+\zeta_1} (|g_1(x(\eta_{i}), y(\eta_{i})) -g_{1}(0, 0)|+|g_{1}(0, 0)|)\\ \nonumber &&+\sum\limits_{i = 1}^{m}|\delta_{i}|I_{a}^{\alpha_{1}+\alpha_{2}}(|f_{1}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))-f_{1}(\eta_{i}, 0, 0)|+|f_{1}(\eta_{i}, 0, 0)|) \\\nonumber && +\sum\limits_{k = 1}^{n}|\theta_{k}|I_a^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_1}(|g_1(x(\xi_{k}), y(\eta_{i})) -g_{1}((0, 0)|+|g_{1}(0, 0)|)\\\nonumber && +\sum\limits_{k = 1}^{n}|\theta_{k}|I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}(|f_{1}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))-f_{1}(\xi_{k}, 0, 0)|+|f_{1}(\xi_{k}, 0, 0)|)\\\nonumber && +|\lambda_{1}| \sum\limits_{i = 1}^{m}|\delta_{i}|I_{a}^{\alpha_{2}}|x(\eta_{i})|+|\lambda_{1}| \sum\limits_{k = 1}^{n}|\theta_{k}|I_{a}^{p_{k}+\alpha_{2}}|x(\xi_{k})| \Big)\Big\}\\ &\le& \sigma_{1}[K_1(\|x\|+\|y\|)+M_1]+ \sigma_{2}[K_2(\|x\|+\|y\|)+M_2]+\sigma_{7}L_{1}(\|x\|+\|y\|)\\ \nonumber &&+\sigma_{8}L_{2}(\|x\|+\|y\|)+\sigma_{5}(\|x\|+\|y\|). \end{eqnarray*}$ $

In a similar manner, one can find that

$\begin{eqnarray*} \|\mathcal{T}_{2}(x, y)\|& \leq & \sigma_{3}[K_1(\|x\|+\|y\|)+M_1]+ \sigma_{4}[K_2(\|x\|+\|y\|)+M_2]\\ &&+ \sigma_{9}L_{1}(\|x\|+\|y\|)+ \sigma_{10}L_{2}(\|x\|+\|y\|)+\sigma_{6} (\|x\|+\|y\|). \end{eqnarray*}$ $

Adding the last two inequalities and using (3.14), we obtain

$\begin{eqnarray*} \|\mathcal{T}(x, y)\| &\leq& ( \sigma_{1}+\sigma_{3})[K_1(\|x\|+\|y\|)+M_1]+ (\sigma_{2}+\sigma_{4})[K_2(\|x\|+\|y\|)+M_2]\\ &&+(\sigma_{7}+\sigma_{9})L_{1}(\|x\|+\|y\|)+(\sigma_{8}+\sigma_{10})L_{2}(\|x\|+\|y\|)\\ &&+(\sigma_{5}+\sigma_{6})(\|x\|+\|y\|)\\ \nonumber &\leq& \rho, \end{eqnarray*}$ $

which shows that $ \mathcal{T}B_{\rho}\subset B_{\rho} $.

Next, we show that $ \mathcal{T} $ is a contraction on $ \mathcal{X}\times \mathcal{X}. $ For that, let $ (x, y), (x_1, y_1) \in \mathcal{X} \times \mathcal{X}. $ Then we have

$\begin{eqnarray*} &&|\mathcal{T}_{1}(x, y)(t)-\mathcal{T}_{1}(x_1, y_1)(t)|\\ &\leq& I_a^{\alpha_{1}+\alpha_{2}+\zeta_1} |g_1(x(t), y(t))-g_1(x_{1}(t), y_{1}(t))|+I_{a}^{\alpha_{1}+\alpha_{2}}|f_{1}(t, x(t), y(t))-f_{1}(t, x_1(t), y_1(t))| \\&& +|\lambda_{1}| I_{a}^{\alpha_{2}}|x(t)-x_1(t)| +\dfrac{\Omega_{1}}{|\Delta|}\Big\{\Omega_{2} \Big(I_{a}^{\alpha_{1}+\alpha_{2}}|f_{1}(b, x(b), y(b))-f_{1}(b, x_1(b), y_1(b))|\\&& +|\lambda_{1}| I_{a}^{\alpha_{2}}|x(b)-x_1(b)|+I_a^{\alpha_{1}+\alpha_{2}+\zeta_1} |g_1(x(b), y(b))-g_1(x_{1}(b), y_{1}(b))| \\\nonumber && +|\lambda_{2}| \sum\limits_{i = 1}^{m}|\mu_{i}| I_{a}^{\alpha_{4}}|y(\eta_{i})-y_1(\eta_{i})|\\\nonumber &&+\sum\limits_{i = 1}^{m}|\mu_{i}| I_{a}^{\alpha_{3}+\alpha_{4}}I^{\zeta_2}_{b}(|g_2(x(\eta_{i}), y(\eta_{i}))-g_2(x_{1}(\eta_{i}), y_{1}(\eta_{i}))|)\\\nonumber &&+\sum\limits_{i = 1}^{m}|\mu_{i}| I_{a}^{\alpha_{3}+\alpha_{4}}|f_{2}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))-f_{2}(\eta_{i}, x_1(\eta_{i}), y_1(\eta_{i}))| \\&&+\sum\limits_{k = 1}^{n}|\nu_{k}|I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}} |f_{2}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))-f_{2}(\xi_{k}, x_1(\xi_{k}), y_1(\xi_{k}))|\\\nonumber && +\sum\limits_{k = 1}^{n}|\nu_{k}|I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}}I^{\zeta_2}_{b}(|g_2(x(\xi_{k}), y(\xi_{k}))-g_2 (x_1(\xi_{k}), y_1(\xi_{k}))|) \\\nonumber &&+|\lambda_{2}| \sum\limits_{k = 1}^{n}|\nu_{k}|I_{a}^{q_{k}+\alpha_{4}}|y(\xi_{k})-y_1(\xi_{k})|\Big)\\ &&+|\Omega_{4}|\Big( I_{a}^{\alpha_{3}+\alpha_{4}}|f_{2}(b, x(b), y(b))-f_{2}(b, x_1(b), y_1(b))|+ |\lambda_{2} |I_{a}^{\alpha_{4}}|y(b)-y_1(b)|\\&& +\sum\limits_{i = 1}^{m}|\delta_{i}|I_{a}^{\alpha_{1}+\alpha_{2}}|f_{1}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))-f_{1}(\eta_{i}, x_1(\eta_{i}), y_1(\eta_{i}))|\\\nonumber &&+\sum\limits_{i = 1}^{m}|\delta_{i}|I_a^{\alpha_{1}+\alpha_{2}+\zeta_1}| g_1(x(\eta_{i}), y(\eta_{i}))- g_1(x_{1}(\eta_{i}), y_{1}(\eta_{i}))|\\\nonumber &&+\sum\limits_{k = 1}^{n}|\theta_{k}|I_a^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_1}|g_{1}(x(\xi_{k}), y(\xi_{k}))-g_{1}(x_{1}(\xi_{k}), y_{1}(\xi_{k}))|\\\nonumber && +|\lambda_{1}| \sum\limits_{k = 1}^{n}|\theta_{k}|I_{a}^{p_{k}+\alpha_{2}}|x(\xi_{k})-x_1(\xi_{k})| \\&& +\sum\limits_{k = 1}^{n}|\theta_{k}|I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}|f_{1}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))-f_{1}(\xi_{k}, x_1(\xi_{k}), y_1(\xi_{k}))|\\&&+|\lambda_{1}| \sum\limits_{i = 1}^{m}|\delta_{i}|I_{a}^{\alpha_{2}}|x(\eta_{i})-x_{1}(\eta_{i})|\Big)\Big\}, \end{eqnarray*}$ $

which leads to

$\begin{eqnarray*} \|\mathcal{T}_{1}(x, y)-\mathcal{T}_{1}(x_1, y_1)\| &\leq & [\sigma_{1}K_{1}+\sigma_{2}K_{2}+\sigma_{7}L_{1}+\sigma_{8}L_{2}+\sigma_{5}] (\|x-x_1\|+\|y-y_1\|). \end{eqnarray*}$ $

Similarly one can obtain

$\begin{eqnarray*} \|\mathcal{T}_{2}(x, y)-\mathcal{T}_{2}(x_1, y_1)\|\leq [\sigma_{3}K_{1}+\sigma_{4}K_{2}+\sigma_{9}L_{1}+\sigma_{10}L_{2}+\sigma_{6}] (\|x-x_1\|+\|y-y_1\|). \end{eqnarray*}$ $

It follows from the last two inequalities that

$\begin{eqnarray*} \|\mathcal{T}(x, y)-\mathcal{T}(x_1, y_1)\| &\leq &[(\sigma_{1}+\sigma_{3})K_{1}+(\sigma_{4}+\sigma_{2})K_{2}+(\sigma_{9}+\sigma_{7})L_{1}+(\sigma_{10}+\sigma_{8})L_{2} \\&&+ \sigma_{5}+\sigma_{6}] [\|x-x_1\|+\|y-y_1\|]. \end{eqnarray*}$ $

Since $ (\sigma_{1}+\sigma_{3})K_{1}+(\sigma_{4}+\sigma_{2})K_{2}+(\sigma_{9}+\sigma_{7})L_{1}+(\sigma_{10}+\sigma_{8})L_{2} +\sigma_{5}+\sigma_{6} < 1 $ by (3.13), therefore $ \mathcal{T} $ is a contraction and hence by Banach's contraction mapping principle, the operator $ \mathcal{T} $ has a unique fixed point. In consequence, the problem (1.3) has a unique solution on $ [a, b]. $ The proof is completed.

The next existence result is based on Krasnosel'ski${\rm{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over i} }}$'s fixed point theorem.

Lemma 3.1. (Krasnosel'ski${{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over i} }}$'s fixed point theorem). ^[30] Let $ B $ be a closed, convex, bounded and nonempty subset of a Banach space $ X. $ Let $ E_1 $ and $ E_2 $ be the operators such that (i) $ E_1x+E_2y \in B $ whenever $ x, y \in B; $ (ii) $ E_1 $ is compact and continuous; (iii) $ E_2 $ is a contractionmapping. Then there exists $ z \in B $ such that $ z = E_1z+E_2z. $

Theorem 3.2. Assume that $ (H_1), $ $ (H_2) $ and $ (H_{4}) $ hold and

$\begin{equation} \sigma_{5}+\sigma_{6} < 1, \end{equation}$ $

(3.15)

where $ \sigma_{5} $ and $ \sigma_{6} $ are given by (3.7) and (3.8) respectively. Then the problem (1.3) has at least one solution on $ [a, b]. $

Proof. Let us split the operators $ \mathcal{T}_{1} $ and $ \mathcal{T}_{2} $ defined by (3.1) and (3.2) respectively into four operators as follows

$ \mathcal{T}_{1}(x, y)(t) = \mathcal{T}_{1, 1}(x, y)(t)+\mathcal{T}_{1, 2}(x, y)(t), \, \, \mathcal{T}_{2}(x, y)(t) = \mathcal{T}_{2, 1}(x, y)(t)+\mathcal{T}_{2, 2}(x, y)(t), $

where

$\begin{eqnarray*} \mathcal{T}_{1, 1}(x, y)(t)& = & I_a^{\alpha_{1}+\alpha_{2}+\zeta_1} g_1(x(t), y(t)) + I_{a}^{\alpha_{1}+\alpha_{2}}f_{1}(t, x(t), y(t))\\ \nonumber && +\dfrac{(t-a)^{\alpha_{2}+\epsilon_{1}-1}}{\Delta\Gamma(\alpha_{2}+\epsilon_{1})}\times \Big\{\Omega_{2}\Big( -I_a^{\alpha_{1}+\alpha_{2}+\zeta_1} g_1(x(b), y(b)) - I_{a}^{\alpha_{1}+\alpha_{2}}f_{1}(b, x(b), y(b))\\&& +\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}g_2(x(\eta_{i}), y(\eta_{i}))) +\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}f_{2}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))\\&& +\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}g_2(x(\xi_{k}), y(\xi_{k}))) +\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}} f_{2}(\xi_{k}, x(\xi_{k}), y(\xi_{k})) \Big)\\&& +\Omega_{4}\Big(-I_{a}^{\alpha_{3}+\alpha_{4}}f_{2}(b, x(b), y(b)) +\sum\limits_{i = 1}^{m}\delta_{i}I_a^{\alpha_{1}+\alpha_{2}+\zeta_1} g_1(x(\eta_{i}), y(\eta_{i}))\\&& +\sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{1}+\alpha_{2}}f_{1}(\eta_{i}, x(\eta_{i}), y(\eta_{i})) +\sum\limits_{k = 1}^{n}\theta_{k}I_a^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_1}g_{1}(x(\xi_{k}), y(\xi_{k}))\\&& +\sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}f_{1}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))\Big)\Big\}, \\ \mathcal{T}_{1, 2}(x, y)(t)& = & -\lambda_{1} I_{a}^{\alpha_{2}}x(t)+\dfrac{(t-a)^{\alpha_{2}+\epsilon_{1}-1}}{\Delta\Gamma(\alpha_{2}+\epsilon_{1})} \Big\{\Omega_{2}\Big(-\lambda_{1} I_{a}^{\alpha_{2}}x(b) - \lambda_{2} \sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{4}}y(\eta_{i})\\&& - \lambda_{2} \sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{4}}y(\xi_{k})\Big)+\Omega_{4}\Big(\lambda_{2} I_{a}^{\alpha_{4}}y(b)-\lambda_{1} \sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{2}}x(\eta_{i})\\& & - \lambda_{1} \sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{2}}x(\xi_{k}) \Big)\Big\}, \end{eqnarray*}$ $

$\begin{eqnarray*} \mathcal{T}_{2, 1}(x, y)(t)& = &I_{a}^{\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}g_2(x(t), y(t)))+ I_{a}^{\alpha_{3}+\alpha_{4}}f_{2}(t, x(t), y(t)) \\ \nonumber && +\dfrac{(t-a)^{\alpha_{4}+\epsilon_{3}-1}}{\Delta\Gamma(\alpha_{4}+\epsilon_{3})}\times \Big\{\Omega_{3}\Big(-I_a^{\alpha_{1}+\alpha_{2}+\zeta_1}g_{1}(x(b), y(b)) -I_{a}^{\alpha_{1}+\alpha_{2}}f_{1}(b, x(b), y(b))\\&& +\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}g_2(x(\eta_{i}), y(\eta_{i}))) +\sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{3}+\alpha_{4}}f_{2}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))\\&& +\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}g_2(x(\xi_{k}), y(\xi_{k}))) +\sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}} f_{2}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))\Big)\\&& +\Omega_{1}\Big( -I_{a}^{\alpha_{3}+\alpha_{4}}f_{2}(b, x(b), y(b)) +\sum\limits_{i = 1}^{m}\delta_{i}I_a^{\alpha_{1}+\alpha_{2}+\zeta_1}g_{1}(x(\eta_{i}), y(\eta_{i})) \\ \nonumber &&+\sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{1}+\alpha_{2}}f_{1}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))+\sum\limits_{k = 1}^{n}\theta_{k}I_a^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_1}g_{1}(x(\xi_{k}), y(\xi_{k}))\\ \nonumber &&+\sum\limits_{k = 1}^{n}\theta_{k}I_a^{p_{k}+\alpha_{1}+\alpha_{2}}f_{1}(\xi_{k}, x(\xi_{k}), y(\xi_{k})) \Big)\Big\}, \\ \mathcal{T}_{2, 2}(x, y)(t)& = & -\lambda_{2} I_{a}^{\alpha_{4}}y(t) +\dfrac{(t-a)^{\alpha_{4}+\epsilon_{3}-1}}{\Delta\Gamma(\alpha_{4}+\epsilon_{3})}\times \Big\{\Omega_{3}\Big(\lambda_{1} I_{a}^{\alpha_{2}}x(b) -\lambda_{2} \sum\limits_{i = 1}^{m}\mu_{i} I_{a}^{\alpha_{4}}y(\eta_{i}) \\&& -\lambda_{2} \sum\limits_{k = 1}^{n}\nu_{k}I_{a}^{q_{k}+\alpha_{4}}y(\xi_{k})\Big)+\Omega_{1}\Big(\lambda_{2} I_{a}^{\alpha_{4}}y(b) -\lambda_{1} \sum\limits_{i = 1}^{m}\delta_{i}I_{a}^{\alpha_{2}}x(\eta_{i})\\ \nonumber &&-\lambda_{1} \sum\limits_{k = 1}^{n}\theta_{k}I_{a}^{p_{k}+\alpha_{2}}x(\xi_{k}) \Big)\Big\}. \end{eqnarray*}$ $

Now we verify the hypotheses of Krasnosel'ski${\rm{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over i} }}$'s fixed point theorem (Lemma 3.1) in three steps.

(ⅰ) In this step, it will be shown that $ \mathcal{T}_{1}(x, y)+\mathcal{T}_{2}(\widehat{x}, \widehat{y}) \in B_{r} $ for all $ (x, y), (\widehat{x}, \widehat{y}) \in B_{r}, $ where $ B_{r}\subset \mathcal{X}\times \mathcal{X} $ is a bounded closed ball with radius

$ r \geqslant \dfrac{ (\sigma_{1}+\sigma_{3})\|\phi_{1}\|+(\sigma_{2}+\sigma_{4})\|\phi_{2}\|+(\sigma_{7}+\sigma_{9})\Lambda_{1}+(\sigma_{8}+\sigma_{10})\Lambda_{2}}{1-\sigma_{5}-\sigma_{6}}. $

As in the proof of Theorem 3.1, we can find that

$\begin{eqnarray*} |\mathcal{T}_{1, 1}(x, y)(t)+\mathcal{T}_{1, 2}(x, y)(t)| \leq \sigma_{1}\|\phi_{1}\|+\sigma_{2}\|\phi_{2}\|+ \sigma_{7}\Lambda_{1}+ \sigma_{8}\Lambda_{2}+\sigma_{5}r, \end{eqnarray*}$ $

and

$\begin{eqnarray*} |\mathcal{T}_{2, 1}(\widehat{x}, \widehat{y})(t)+\mathcal{T}_{2, 2}(\widehat{x}, \widehat{y})(t)| \leq \sigma_{3}\|\phi_{1}\|+\sigma_{4}\|\phi_{2}\|+ \sigma_{9}\Lambda_{1}+ \sigma_{10}\Lambda_{2}+\sigma_{6}r, \end{eqnarray*}$ $

which lead to the inequality

$\begin{eqnarray*} \| \mathcal{T}_{1}(x, y)+\mathcal{T}_{2}(\widehat{x}, \widehat{y})\| &\leq & (\sigma_{1}+\sigma_{3})\|\phi_{1}\|+(\sigma_{2}+\sigma_{4})\|\phi_{2}\|+(\sigma_{7}+\sigma_{9})\Lambda_{1} \\&&+ (\sigma_{8}+\sigma_{10})\Lambda_{2}+ (\sigma_{5}+\sigma_{6})r \leq r. \end{eqnarray*}$ $

Thus $ \mathcal{T}_{1}(x, y)+\mathcal{T}_{2}(\widehat{x}, \widehat{y}) \in B_{r} $.

(ⅱ) Here we establish that $ (\mathcal{T}_{1, 2}, \mathcal{T}_{2, 2}) $ is a contraction mapping. Let $ (x, y), (\widehat{x}, \widehat{y}) \in B_{r}. $ Then it is easy to find that

$ |\mathcal{T}_{1, 2}(x, y)(t)-\mathcal{T}_{1, 2}(\widehat{x}, \widehat{y})(t)| \leq \sigma_{5}[\|x-\widehat{x}\|+\|y-\widehat{y}\|], $

$ |\mathcal{T}_{2, 2}(x, y)(t)-\mathcal{T}_{2, 2}(\widehat{x}, \widehat{y})| \leq \sigma_{6}[\|x-\widehat{x}\|+\|y-\widehat{y}\|]. $

Consequently, we get

$\begin{eqnarray*} \|(\mathcal{T}_{1, 2}, \mathcal{T}_{2, 2})(x, y)-(\mathcal{T}_{1, 2}, \mathcal{T}_{2, 2})(\widehat{x}, \widehat{y})\| \leq (\sigma_{5}+\sigma_{6}) [\|x-\widehat{x}\|+\|y-\widehat{y}\|], \end{eqnarray*}$ $

which, by (3.15), implies that $ (\mathcal{T}_{1, 2}, \mathcal{T}_{2, 2}) $ is a contraction.

(ⅲ) We show that $ (\mathcal{T}_{1, 1}, \mathcal{T}_{2, 1}) $ is compact and continuous.

Continuity of $ (\mathcal{T}_{1, 1}, \mathcal{T}_{2, 1}) $ is obvious. For $ (x, y) \in B_{r}, $ we have

$\begin{eqnarray*} |\mathcal{T}_{1, 1}(x, y)(t)| &\leq& I_a^{\alpha_{1}+\alpha_{2}+\zeta_1}(|g_1(x(t), y(t))|) + I_{a}^{\alpha_{1}+\alpha_{2}}(|f_{1}(t, x(t), y(t))|)\\\nonumber && +\dfrac{(b-a)^{\alpha_{2}+\epsilon_{1}-1}}{|\Delta|\Gamma(\alpha_{2}+\epsilon_{1})} \Big\{\Omega_{2} \Big( I_a^{\alpha_{1}+\alpha_{2}+\zeta_1} (|g_1(x(b), y(b)) |) +I_{a}^{\alpha_{1}+\alpha_{2}}(|f_{1}(b, x(b), y(b))|)\\&& +\sum\limits_{i = 1}^{m}|\mu_{i}| I_{a}^{\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}(|g_2(x(\eta_{i}), y(\eta_{i})) |) +\sum\limits_{i = 1}^{m}|\mu_{i}| I_{a}^{\alpha_{3}+\alpha_{4}}(|f_{2}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))|)\\&& +\sum\limits_{k = 1}^{n}|\nu_{k}|I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}(|g_2(x(\xi_{k}), y(\xi_{k}))|)+\sum\limits_{k = 1}^{n}|\nu_{k}|I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}}(| f_{2}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))|) \Big) \\&&+\Omega_{4}\Big(I_{a}^{\alpha_{3}+\alpha_{4}}(| f_{2}(b, x(b), y(b))|) +\sum\limits_{i = 1}^{m}\delta_{i}I_a^{\alpha_{1}+\alpha_{2}+\zeta_1} (|g_1(x(\eta_{i}), y(\eta_{i})) |)\\&& +\sum\limits_{i = 1}^{m}|\delta_{i}|I_{a}^{\alpha_{1}+\alpha_{2}}(|f_{1}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))|) +\sum\limits_{k = 1}^{n}\theta_{k}I_a^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_1}(|g_1(x(\xi_{k}), y(\eta_{i})) |)\\&& +\sum\limits_{k = 1}^{n}|\theta_{k}|I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}(|f_{1}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))|) \Big)\Big\}\\ &\le& \sigma_{1}\|\phi_{1}\|+\sigma_{2}\|\phi_{2}\|+ \sigma_{7}\Lambda_{1}+ \sigma_{8}\Lambda_{2}. \end{eqnarray*}$ $

In a similar manner, we can get $ |\mathcal{T}_{2, 1}(x, y)(t)| \leq \sigma_{3}\|\phi_{1}\|+\sigma_{4}\|\phi_{2}\|+ \sigma_{9}\Lambda_{1}+ \sigma_{10}\Lambda_{2}. $ Thus

$ \|(\mathcal{T}_{1, 1}, \mathcal{T}_{2, 1})(x, y)\|\leq (\sigma_{1}+\sigma_{3})\|\phi_{1}\|+(\sigma_{2}+\sigma_{4})\|\phi_{2}\|+(\sigma_{7}+\sigma_{9})\Lambda_{1} + (\sigma_{8}+\sigma_{10})\Lambda_{2}, $

which means that $ (\mathcal{T}_{1, 1}, \mathcal{T}_{2, 1}) $ is uniformly bounded on $ B_{r}. $

In order to show the equicontinuity of $ (\mathcal{T}_{1, 1}, \mathcal{T}_{2, 1}), $ we take $ t_{1}, t_{2} \in [a, b] $ with $ t_{1} < t_{2}. $ Then, for arbitrary $ (x, y) \in B_{r}, $ we obtain

$\begin{eqnarray*} &&| \mathcal{T}_{1, 1}(x, y)(t_{2})-\mathcal{T}_{1, 1}(x, y)(t_{1})|\\ &\leq& \Big|\int_{a}^{t_{2}}\dfrac{(t_{2}-s)^{\alpha_{1}+\alpha_{2}-1}}{\Gamma(\alpha_{1}+\alpha_{2})}f_{1}(s, x(s), y(s))ds -\int_{a}^{t_{1}}\dfrac{(t_{1}-s)^{\alpha_{1}+\alpha_{2}-1}}{\Gamma(\alpha_{1}+\alpha_{2})}f_{1}(s, x(s), y(s))ds\Big| \\&& + \Big|\int_{a}^{t_{2}}\dfrac{(t_{2}-s)^{\alpha_{1}+\alpha_{2}+\zeta_1-1}}{\Gamma(\alpha_{1}+\alpha_{2}+\zeta_1)}g_1(x(s), y(s))ds -\int_{a}^{t_{1}}\dfrac{(t_{1}-s)^{\alpha_{1}+\alpha_{2}+\zeta_1-1}}{\Gamma(\alpha_{1}+\alpha_{2})+\zeta_1}g_1(x(s), y(s))ds\Big| \\&&+\dfrac{ |(t_{2}-a)^{\alpha_{2}+\epsilon_{1}-1}-(t_{1}-a)^{\alpha_{2}+\epsilon_{1}-1}|}{|\Delta| \Gamma(\alpha_{2}+\epsilon_{1})} \Big\{\Omega_{2} \Big( I_a^{\alpha_{1}+\alpha_{2}+\zeta_1} (|g_1(x(b), y(b)) |)\\&& +I_{a}^{\alpha_{1}+\alpha_{2}}(|f_{1}(b, x(b), y(b))|) +\sum\limits_{i = 1}^{m}|\mu_{i}| I_{a}^{\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}(|g_2(x(\eta_{i}), y(\eta_{i})) |)\\&& +\sum\limits_{i = 1}^{m}|\mu_{i}| I_{a}^{\alpha_{3}+\alpha_{4}}(|f_{2}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))|) +\sum\limits_{k = 1}^{n}|\nu_{k}|I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}(|g_2(x(\xi_{k}), y(\xi_{k}))|)\\&& +\sum\limits_{k = 1}^{n}|\nu_{k}|I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}}(| f_{2}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))|) \Big) +\Omega_{4}\Big(I_{a}^{\alpha_{3}+\alpha_{4}}(| f_{2}(b, x(b), y(b))|)\\&& +\sum\limits_{i = 1}^{m}|\delta_{i}|I_a^{\alpha_{1}+\alpha_{2}+\zeta_1} (|g_1(x(\eta_{i}), y(\eta_{i})) |) +\sum\limits_{i = 1}^{m}|\delta_{i}|I_{a}^{\alpha_{1}+\alpha_{2}}(|f_{1}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))|)\\&& +\sum\limits_{k = 1}^{n}|\theta_{k}|I_a^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_1}(|g_1(x(\xi_{k}), y(\eta_{i})) |) +\sum\limits_{k = 1}^{n}|\theta_{k}|I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}(|f_{1}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))|) \Big)\Big\}\\ & \leqslant & \dfrac{\|\phi_{1}\|}{\Gamma(\alpha_{1}+\alpha_{2}+1)}\Big\{2(t_{2}-t_1)^{\alpha_{1}+\alpha_{2}}+|(t_{2}-a)^{\alpha_{1}+\alpha_{2}}-(t_{1}-a)^{\alpha_{1}+\alpha_{2}}|\Big\}\\&& +\dfrac{\Lambda_{1}}{\Gamma(\alpha_{1}+\alpha_{2}+\zeta_1+1)}\Big\{2(t_{2}-t_1)^{\alpha_{1}+\alpha_{2}+\zeta_1}+|(t_{2}-a)^{\alpha_{1}+\alpha_{2}+\zeta_1}-(t_{1}-a)^{\alpha_{1}+\alpha_{2}+\zeta_1}|\Big\}\\&& +\dfrac{ |(t_{2}-a)^{\alpha_{2}+\epsilon_{1}-1}-(t_{1}-a)^{\alpha_{2}+\epsilon_{1}-1}|}{|\Delta| \Gamma(\alpha_{2}+\epsilon_{1})} \Big\{\Omega_{2}\Big(\dfrac{\|\phi_{1}\|(b-a)^{\alpha_{1}+\alpha_{2}}}{\Gamma(\alpha_{1}+\alpha_{2}+1)}\\&& +\dfrac{\Lambda_{1}(b-a)^{\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma(\alpha_{1}+\alpha_{2}+\zeta_{1}+1)} +\sum\limits_{i = 1}^{m} \dfrac{|\mu_{i}| \|\phi_{2}\|(\eta_{i}-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma(\alpha_{3}+\alpha_{4}+1)}\\&& +\sum\limits_{i = 1}^{m} \dfrac{|\mu_{i}| \Lambda_{2}(b-a)^{\zeta_{2}}(\eta_{i}-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma(\alpha_{3}+\alpha_{4}+1)\Gamma(\zeta_{2}+1)} +\sum\limits_{k = 1}^{n}\dfrac{|\nu_{k}| \|\phi_{2}\|(\xi_{k}-a)^{q_{k}+\alpha_{3}+\alpha_{4}}}{\Gamma(q_{k}+\alpha_{3}+\alpha_{4}+1)}\Big)\\&& +\sum\limits_{k = 1}^{n}\dfrac{|\nu_{k}| \Lambda_{2}(b-a)^{\zeta_{2}}(\xi_{k}-a)^{q_{k}+\alpha_{3}+\alpha_{4}}}{\Gamma(\zeta_{2}+1)\Gamma(q_{k}+\alpha_{3}+\alpha_{4}+1)}\Big)\\&& +|\Omega_{4}|\Big( \dfrac{\|\phi_{2}\|(b-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma(\alpha_{3}+\alpha_{4}+1)} +\sum\limits_{i = 1}^{m}\dfrac{|\delta_{i}| \|\phi_{1}\|(\eta_{i}-a)^{\alpha_{1}+\alpha_{2}}}{\Gamma(\alpha_{1}+\alpha_{2}+1)}\\&& +\sum\limits_{i = 1}^{m}\dfrac{|\delta_{i}| \Lambda_{1}(\eta_{i}-a)^{\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma(\alpha_{1}+\alpha_{2}+\zeta_{1}+1)} +\sum\limits_{k = 1}^{n}\dfrac{|\theta_{k}| \|\phi_{1}\|(\xi_{k}-a)^{p_{k}+\alpha_{1}+\alpha_{2}}}{\Gamma(p_{k}+\alpha_{1}+\alpha_{2}+1)} \\&& +\sum\limits_{k = 1}^{n}\dfrac{|\theta_{k}| \Lambda_{1}(\xi_{k}-a)^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma(p_{k}+\alpha_{1}+\alpha_{2}+\zeta_{1}+1)}\Big)\Big\}\to 0, \end{eqnarray*}$ $

as $ t_{2}\rightarrow t_{1} $ independently of $ (x, y) \in B_{r}. $ Also

$\begin{eqnarray*} &&| \mathcal{T}_{2, 1}(x, y)(t_{2})-\mathcal{T}_{2, 1}(x, y)(t_{1})|\\ &\leq& \dfrac{\|\phi_{2}\|}{\Gamma(\alpha_{3}+\alpha_{4}+1)}\Big\{2(t_{2}-t_1)^{\alpha_{3}+\alpha_{4}}+|(t_{2}-a)^{\alpha_{3}+\alpha_{4}}-(t_{1}-a)^{\alpha_{3}+\alpha_{4}}|\Big\}\\&& +\dfrac{\Lambda_{2}(b-a)^{\zeta_{2}}}{\Gamma(\zeta_{2}+1)\Gamma(\alpha_{3}+\alpha_{4}+1)}\Big\{2(t_{2}-t_1)^{\alpha_{3}+\alpha_{4}}+|(t_{2}-a)^{\alpha_{3}+\alpha_{4}}-(t_{1}-a)^{\alpha_{3}+\alpha_{4}}|\Big\}\\&& +\dfrac{ |(t_{2}-a)^{\alpha_{4}+\epsilon_{3}-1}-(t_{1}-a)^{\alpha_{4}+\epsilon_{3}-1}|}{|\Delta| \Gamma(\alpha_{4}+\epsilon_{3})} \Big\{|\Omega_{3}|\Big(\dfrac{\|\phi_{1}\|(b-a)^{\alpha_{1}+\alpha_{2}}}{\Gamma(\alpha_{1}+\alpha_{2}+1)}\\&& +\dfrac{\Lambda_{1}(b-a)^{\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma(\alpha_{1}+\alpha_{2}+\zeta_{1}+1)} +\sum\limits_{i = 1}^{m} \dfrac{|\mu_{i}| \|\phi_{2}\|(\eta_{i}-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma(\alpha_{3}+\alpha_{4}+1)}\\&& +\sum\limits_{i = 1}^{m} \dfrac{|\mu_{i}| \Lambda_{2}(b-a)^{\zeta_{2}}(\eta_{i}-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma(\zeta_{2}+1)\Gamma(\alpha_{3}+\alpha_{4}+1)} +\sum\limits_{k = 1}^{n}\dfrac{|\nu_{k}| \|\phi_{2}\|(\xi_{k}-a)^{q_{k}+\alpha_{3}+\alpha_{4}}}{\Gamma(q_{k}+\alpha_{3}+\alpha_{4}+1)}\\&& +\sum\limits_{k = 1}^{n}\dfrac{|\nu_{k}| \Lambda_{2}(b-a)^{\zeta_{2}}(\xi_{k}-a)^{q_{k}+\alpha_{3}+\alpha_{4}}}{\Gamma(\zeta_{2}+1)\Gamma(q_{k}+\alpha_{3}+\alpha_{4}+1)}\Big) +\Omega_{1}\Big( \dfrac{\|\phi_{2}\|(b-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma(\alpha_{3}+\alpha_{4}+1)}\\&& +\sum\limits_{i = 1}^{m}\dfrac{|\delta_{i}| \|\phi_{1}\|(\eta_{i}-a)^{\alpha_{1}+\alpha_{2}}}{\Gamma(\alpha_{1}+\alpha_{2}+1)} +\sum\limits_{i = 1}^{m}\dfrac{|\delta_{i}| \Lambda_{1}(\eta_{i}-a)^{\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma(\alpha_{1}+\alpha_{2}+\zeta_{1}+1)} \\&& +\sum\limits_{k = 1}^{n}\dfrac{|\theta_{k}| \|\phi_{1}\|(\xi_{k}-a)^{p_{k}+\alpha_{1}+\alpha_{2}}}{\Gamma(p_{k}+\alpha_{1}+\alpha_{2}+1)} +\sum\limits_{k = 1}^{n}\dfrac{|\theta_{k}| \Lambda_{1}(\xi_{k}-a)^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma(p_{k}+\alpha_{1}+\alpha_{2}+\zeta_{1}+1)} \Big)\Big\}\to 0, \end{eqnarray*}$ $

as $ t_{2}\rightarrow t_{1} $ independently of $ (x, y) \in B_{r}. $ Thus $ |(\mathcal{T}_{1, 1}, \mathcal{T}_{2, 1})(x, y)(t_{2})-(\mathcal{T}_{1, 1}, \mathcal{T}_{2, 1})(x, y)(t_{1})| $ vanishes as $ t_{2}\rightarrow t_{1} $ independently of $ (x, y) \in B_{r}, $ which shows that $ (\mathcal{T}_{1, 1}, \mathcal{T}_{2, 1}) $ is equicontinuous. So we deduce by the Arzelá-Ascoli theorem that $ (\mathcal{T}_{1, 1}, \mathcal{T}_{2, 1}) $ is compact on $ B_{r}. $

It follows from the steps (i)-(iii) that the hypotheses of Krasnosel'ski${\rm{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over i} }}$'s fixed point theorem are satisfied, so its conclusion implies that the problem (1.3) has at least one solution on $ [a, b]. $ This finishes the proof.

Remark 3.1. The conclusion of Theorem 3.2 can also be achieved by assuming $ (H_{1}), (H_{2}), (H_4) $ and the condition: $ (\sigma_{1}+\sigma_{3})K_{1}+(\sigma_{2}+\sigma_{4})K_{2}+(\sigma_{7}+\sigma_{9})L_{1}+(\sigma_{8}+\sigma_{10})L_{2} < 1, $ where $ \sigma_1, \dots, \sigma_4 $ are given in (3.3)–(3.6) and $ \sigma_7, \dots, \sigma_{10} $ are given in (3.9)–(3.12).

In the following result, we prove the existence of solutions for the problem (1.3) by applying the Leray-Schauder alternative ^[29].

Lemma 3.2. (Leray-Schauder alternative ^[29]) Let $ E $ be a Banach space, $ M $ be closed, convex subset of $ E $, $ U $ is an open subset of $ C $ and $ 0 \in U. $ Suppose that $ F:\overline{U}\to C $ is continuous, compact map (that is, $ F(U) $ is a relatively compact subset of $ C $). Then either (i) $ F $ has a fixed point in $ \overline{U} $, or (ii) there are $ u \in \partial U, $ and $ \lambda \in (0, 1) $ with $ u = \lambda F(U) $.

Theorem 3.3. Assume that ($ H_{3}) $ holds. Then there exists at least one solution for the problem (1.3) on $ [a, b] $ provided that

$\begin{equation} (\sigma_{1}+\sigma_{3})u_{i}+(\sigma_{2}+\sigma_{4})v_{i}+(\sigma_{7}+\sigma_{9})\omega_{i}+(\sigma_{8}+\sigma_{10})\tau_{i}+(\sigma_{5}+\sigma_{6}) < 1\; , \; i = 1, 2, \end{equation}$ $

(3.16)

where $ \sigma_1, \dots, \sigma_{10} $ are given in (3.3)–(3.12).

Proof. For all $ (x, y) \in B_{\rho}\subset \mathcal{X} \times \mathcal{X}, $ where $ B_{\rho} $ defined by (3.14), there exist positive constants $ N_{1}, \dots, N_4 $ such that $ |f_{1}(t, x, y)|\leq N_{1}, |f_{2}(t, x, y)|\leq N_{2}, |g_{1}(x, y)|\leq N_{3}, |g_{2}(x, y)|\leq N_{4}, $ Then we show that $ \mathcal{T}: \mathcal{X} \times \mathcal{X} \rightarrow \mathcal{X} \times \mathcal{X} $ is completely continuous. Observe that continuity of $ f_{1}, f_{2}, g_{1}, g_2 $ implies that of the operator $ \mathcal{T}. $ For $ (x, y) \in B_{\rho} $, as in the proof of Theorem 3.1, we have

$\begin{eqnarray*} \|\mathcal{T}_{1}(x, y)\| &\leq& \sigma_{1}N_{1}+\sigma_{2}N_{2}+\sigma_{7}N_{3}+\sigma_{8}N_{4}+\rho \sigma_{5}, \\ \| \mathcal{T}_{2}(x, y)\|&\leq & \sigma_{3}N_{1}+\sigma_{4}N_{2}+\sigma_{9}N_{3}+\sigma_{10}N_{4}+\rho \sigma_{6}. \end{eqnarray*}$ $

From the preceding inequalities, we get

$ \| \mathcal{T}(x, y)\| \leq (\sigma_{1}+ \sigma_{3})N_{1}+(\sigma_{2}+\sigma_{4})N_{2}+(\sigma_{7}+\sigma_{9})N_{3}+(\sigma_{8}+\sigma_{10})N_{4}+\rho (\sigma_{5}+\sigma_{6}), $

which implies that $ \mathcal{T}B_{\rho} $ is uniformly bounded.

Next we show that $ \mathcal{T}B_{\rho} $ is equicontinues. Let $ t_{1}, t_{2} \in [a, b] $ with $ t_{2} > t_{1}. $ Then, for arbitrary $ (x, y) \in B_{\rho}, $ we obtain

$\begin{eqnarray*} &&|\mathcal{T}_{1}(x, y)(t_{2})-\mathcal{T}_{1}(x, y)(t_{1})| \\ &\leq & |\int_{a}^{t_{2}}\dfrac{(t_{2}-s)^{\alpha_{1}+\alpha_{2}-1}}{\Gamma(\alpha_{1}+\alpha_{2})}f_{1}(s, x(s), y(s))ds -\int_{a}^{t_{1}}\dfrac{(t_{1}-s)^{\alpha_{1}+\alpha_{2}-1}}{\Gamma(\alpha_{1}+\alpha_{2})}f_{1}(s, x(s), y(s))ds | \\&& + |\int_{a}^{t_{2}}\dfrac{(t_{2}-s)^{\alpha_{1}+\alpha_{2}+\zeta_1-1}}{\Gamma(\alpha_{1}+\alpha_{2}+\zeta_1)}g_1(x(s), y(s))ds -\int_{a}^{t_{1}}\dfrac{(t_{1}-s)^{\alpha_{1}+\alpha_{2}+\zeta_1-1}}{\Gamma(\alpha_{1}+\alpha_{2}+\zeta_1)}g_1(x(s), y(s))ds | \\&&+ \dfrac{|\lambda_{1}|}{\Gamma(\alpha_{2})} |\int_{a}^{t_{1}}[(t_{2}-s)^{\alpha_{2}-1}-(t_{1}-s)^{\alpha_{2}-1}]x(s)ds+\int_{t_1}^{t_{2}}(t_{2}-s)^{\alpha_{2}-1}x(s)ds |\\&& +\dfrac{ |(t_{2}-a)^{\alpha_{2}+\epsilon_{1}-1}-(t_{1}-a)^{\alpha_{2}+\epsilon_{1}-1}|}{|\Delta| \Gamma(\alpha_{2}+\epsilon_{1})} \{\Omega_{2} ( I_a^{\alpha_{1}+\alpha_{2}+\zeta_1} (|g_1(x(b), y(b)) |)\\&& +I_{a}^{\alpha_{1}+\alpha_{2}}(|f_{1}(b, x(b), y(b))|)+|\lambda_{1}| I_{a}^{\alpha_{2}}|x(b)| +|\lambda_{2} |\sum\limits_{i = 1}^{m}|\mu_{i}| I_{a}^{\alpha_{4}}|y(\eta_{i})|\\&& +\sum\limits_{i = 1}^{m}|\mu_{i}| I_{a}^{\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}(|g_2(x(\eta_{i}), y(\eta_{i})) |) +\sum\limits_{i = 1}^{m}|\mu_{i}| I_{a}^{\alpha_{3}+\alpha_{4}}(|f_{2}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))|)\\&& +\sum\limits_{k = 1}^{n}|\nu_{k}|I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}}(I^{\zeta_2}_{b}(|g_2(x(\xi_{k}), y(\xi_{k}))|) +\sum\limits_{k = 1}^{n}|\nu_{k}|I_{a}^{q_{k}+\alpha_{3}+\alpha_{4}}(| f_{2}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))|)\\&& +|\lambda_{2}| \sum\limits_{k = 1}^{n}|\nu_{k}|I_{a}^{q_{k}+\alpha_{4}}|y(\xi_{k})| ) +\Omega_{4} (|\lambda_{2}| I_{a}^{\alpha_{4}}|y(b)|+I_{a}^{\alpha_{3}+\alpha_{4}}(| f_{2}(b, x(b), y(b))|)\\&& +\sum\limits_{i = 1}^{m}|\delta_{i}|I_a^{\alpha_{1}+\alpha_{2}+\zeta_1} (|g_1(x(\eta_{i}), y(\eta_{i})) |) +\sum\limits_{i = 1}^{m}|\delta_{i}|I_{a}^{\alpha_{1}+\alpha_{2}}(|f_{1}(\eta_{i}, x(\eta_{i}), y(\eta_{i}))|)\\&& +|\lambda_{1}| \sum\limits_{i = 1}^{m}|\delta_{i}|I_{a}^{\alpha_{2}}|x(\eta_{i})| +\sum\limits_{k = 1}^{n}|\theta_{k}|I_a^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_1}(|g_1(x(\xi_{k}), y(\eta_{i})) |)\\&& +\sum\limits_{k = 1}^{n}|\theta_{k}|I_{a}^{p_{k}+\alpha_{1}+\alpha_{2}}(|f_{1}(\xi_{k}, x(\xi_{k}), y(\xi_{k}))|)+|\lambda_{1} |\sum\limits_{k = 1}^{n}|\theta_{k}|I_{a}^{p_{k}+\alpha_{2}}|x(\xi_{k})| ) \}\\ && \leqslant \dfrac{N_{1}}{\Gamma(\alpha_{1}+\alpha_{2}+1)} \{2(t_{2}-t_1)^{\alpha_{1}+\alpha_{2}}+|(t_{2}-a)^{\alpha_{1}+\alpha_{2}}-(t_{1}-a)^{\alpha_{1}+\alpha_{2}}| \}\\&& +\dfrac{N_{3}}{\Gamma(\alpha_{1}+\alpha_{2}+\zeta_1+1)} \{2(t_{2}-t_1)^{\alpha_{1}+\alpha_{2}+\zeta_1}+|(t_{2}-a)^{\alpha_{1}+\alpha_{2}+\zeta_1}-(t_{1}-a)^{\alpha_{1}+\alpha_{2}+\zeta_1}| \}\\&& +\dfrac{|\lambda_{1}|\rho}{\Gamma(\alpha_{2}+1)} \{2(t_{2}-t_1)^{\alpha_{2}}+|(t_{2}-a)^{\alpha_{2}}-(t_{1}-a)^{\alpha_{2}}| \}\\&& +\dfrac{ |(t_{2}-a)^{\alpha_{2}+\epsilon_{1}-1}-(t_{1}-a)^{\alpha_{2}+\epsilon_{1}-1}|}{|\Delta| \Gamma(\alpha_{2}+\epsilon_{1})} \{\Omega_{2} (\dfrac{|\lambda_{1}|\rho(b-a)^{\alpha_{2}}}{\Gamma(\alpha_{2}+1)} +\sum\limits_{i = 1}^{m}\dfrac{|\lambda_{2}\mu_{i}|\rho(\eta_{i}-a)^{\alpha_{4}}}{ \Gamma(\alpha_{4}+1)}\\&& +\sum\limits_{k = 1}^{n}\dfrac{|\lambda_{2} \nu_{k}|\rho(\xi_{k}-a)^{q_{k}+\alpha_{4}}}{\Gamma(q_{k}+\alpha_{4}+1)} + \dfrac{N_{1}(b-a)^{\alpha_{1}+\alpha_{2}}}{\Gamma(\alpha_{1}+\alpha_{2}+1)} +\dfrac{N_{3}(b-a)^{\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma(\alpha_{1}+\alpha_{2}+\zeta_{1}+1)} \\&& +\sum\limits_{i = 1}^{m} \dfrac{|\mu_{i}| N_{2}(\eta_{i}-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma(\alpha_{3}+\alpha_{4}+1)} +\sum\limits_{i = 1}^{m} \dfrac{|\mu_{i}| N_{4}(b-a)^{\zeta_{2}}(\eta_{i}-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma(\alpha_{3}+\alpha_{4}+1)\Gamma(\zeta_{2}+1)}\\&& +\sum\limits_{k = 1}^{n}\dfrac{|\nu_{k}| N_{2}(\xi_{k}-a)^{q_{k}+\alpha_{3}+\alpha_{4}}}{\Gamma(q_{k}+\alpha_{3}+\alpha_{4}+1)} +\sum\limits_{k = 1}^{n}\dfrac{|\nu_{k}| N_{4}(b-a)^{\zeta_{2}}(\xi_{k}-a)^{q_{k} +\alpha_{3}+\alpha_{4}}}{\Gamma(\zeta_{2}+1)\Gamma(q_{k}+\alpha_{3}+\alpha_{4}+1)} )\\&& +|\Omega_{4}| ( \dfrac{N_{2}(b-a)^{\alpha_{3}+\alpha_{4}}}{\Gamma(\alpha_{3}+\alpha_{4}+1)} +\sum\limits_{i = 1}^{m}\dfrac{|\delta_{i}| N_{1}(\eta_{i}-a)^{\alpha_{1}+\alpha_{2}}}{\Gamma(\alpha_{1}+\alpha_{2}+1)}\\&& +\sum\limits_{i = 1}^{m}\dfrac{|\delta_{i}|N_{3}(\eta_{i}-a)^{\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma(\alpha_{1}+\alpha_{2}+\zeta_{1}+1)} +\sum\limits_{k = 1}^{n}\dfrac{|\theta_{k}|N_{1}(\xi_{k}-a)^{p_{k}+\alpha_{1}+\alpha_{2}}}{\Gamma(p_{k}+\alpha_{1}+\alpha_{2}+1)} \\&& +\sum\limits_{k = 1}^{n}\dfrac{|\theta_{k}| N_{3}(\xi_{k}-a)^{p_{k}+\alpha_{1}+\alpha_{2}+\zeta_{1}}}{\Gamma(p_{k}+\alpha_{1}+\alpha_{2}+\zeta_{1}+1)} +\dfrac{|\lambda_{2}|\rho (b-a)^{\alpha_{4}}}{\Gamma(\alpha_{4}+1)}\\&& +\sum\limits_{i = 1}^{m}\dfrac{|\lambda_{1} \delta_{i}|\rho(\eta_{i}-a)^{\alpha_{2}}}{\Gamma(\alpha_{2}+1)} +\sum\limits_{k = 1}^{n}\dfrac{|\lambda_{1} \theta_{k}|\rho(\xi_{k}-a)^{p_{k}+\alpha_{2}}}{\Gamma(p_{k}+\alpha_{2}+1)} ) \}, \end{eqnarray*}$ $

which tends to zero as $ t_{1}\rightarrow t_{2} $ independent of $ (x, y) \in B_{\rho}. $ Similarly, it can be established that $ | \mathcal{T}_{2}(x, y)(t_{2})-\mathcal{T}_{2}(x, y)(t_{1})|\rightarrow 0 $ as $ t_{2}\rightarrow t_{1} $ independently of $ (x, y) \in B_{\rho}. $ Thus the operator $ \mathcal{T} $ is equicontinuous. Hence, by the Arzelá-Ascoli theorem, the operator $ \mathcal{T} $ is completely continuous.

Next we consider the set

$ \Omega = \{(x, y) \in \mathcal{X}\times \mathcal{X}: (x, y) = r \mathcal{T}(x, y), \; \; 0\leq r \leq 1\}, $

and show that it is bounded. Let $ (x, y) \in \Omega $, then $ (x, y) = r \mathcal{T}(x, y) $ implies that $ x(t) = r \mathcal{T}_{1}(x, y)(t), $ and $ y(t) = r \mathcal{T}_{2}(x, y)(t), \forall t \in [a, b]. $ By the condition ($ H_{3} $), we obtain

$\begin{eqnarray*} \|x\|&\leq & \sigma_{1}( u_{0}+u_{1}|x|+u_{2}|y|) +\sigma_{2}( v_{0}+v_{1}|x|+v_{2}|y|)+ \sigma_{7}( \omega_{0}+\omega_{1}|x|+\omega_{2}|y|)\\ &&+\sigma_{8}( \tau_{0}+\tau_{1}|x|+\tau_{2}|y|)+\sigma_{5}(\|x\|+ \|y\| ), \\ \|y\|&\leq & \sigma_{3}( u_{0}+u_{1}|x|+u_{2}|y|) +\sigma_{4}( v_{0}+v_{1}|x|+v_{2}|y|)+ \sigma_{9}( \omega_{0}+\omega_{1}|x|+\omega_{2}|y|)\\ &&+\sigma_{10}( \tau_{0}+\tau_{1}|x|+\tau_{2}|y|) +\sigma_{6}(\|x\|+ \|y\| ). \end{eqnarray*}$ $

Adding the above inequalities, we get

$\begin{eqnarray*} \|x\|+ \|y\| &\leq & (\sigma_{1}+\sigma_{3})( u_{0}+u_{1}|x|+u_{2}|y|) +(\sigma_{2}+\sigma_{4})( v_{0}+v_{1}|x|+v_{2}|y|)\\&&+ (\sigma_{7}+\sigma_{9})( \omega_{0}+\omega_{1}|x|+\omega_{2}|y|)+(\sigma_{8}+\sigma_{10})( \tau_{0}+\tau_{1}|x|+\tau_{2}|y|) \\ &&+(\sigma_{5}+\sigma_{6})(\|x\|+ \|y\| ), \end{eqnarray*}$ $

which leads to

$ \|(x, y)\| \leq \dfrac{(\sigma_{1}+\sigma_{3})u_{0}+(\sigma_{2}+\sigma_{4})v_{0}+(\sigma_{7}+\sigma_{9})\omega_{0}+(\sigma_{8}+\sigma_{10}) \tau_{0}}{\sigma^{*}}, $

where

$\begin{eqnarray*} \sigma^{*}& = & \min\{1-(\sigma_{1}+\sigma_{3})u_{1}-(\sigma_{2}+\sigma_{4})v_{1}-(\sigma_{7}+\sigma_{9})\omega_{1}-(\sigma_{8}+\sigma_{10})\tau_{1}-(\sigma_{5}+\sigma_{7}), \\ &&1-(\sigma_{1}+\sigma_{3})u_{2}-(\sigma_{2}+\sigma_{4})v_{2}-(\sigma_{7}+\sigma_{9})\omega_{2}-(\sigma_{8}+\sigma_{10})\tau_{2}-(\sigma_{5}+\sigma_{7})\} > 0 \end{eqnarray*}$ $

by condition (3.16). Thus the set $ \Omega $ is bounded. Hence it follows by the Leray-Schauder alternative for single-valued maps ^[29] that the problem (1.3) has at least one solution on $ [a, b], $ which completes the proof.

4. Examples

Consider a coupled system of Hilfer fractional differential equations with boundary conditions:

$\begin{equation} \left\{ \begin{array}{ll} ^HD^{1/2, 3/4}(^HD^{1/6, 4/5}+1/90 )x(t) = I_{0^+}^{1/2}g_{1}(x, y)+f_{1}(t, x, y), \, t \in [0, 1], \\[0.4cm] ^HD^{1/2, 3/4}(^HD^{1/2, 1/7}+1/100 )y(t) = I_{1^-}^{1/3}g_{2}(x, y)+f_{2}(t, x, y), \, t \in [0, 1], \\[0.4cm] x(0) = y(0) = 0, \\ x(1) = \dfrac{1}{100}y(1/10)+\dfrac{1}{200}y(1/5)+\dfrac{1}{300}y(3/10)+\dfrac{1}{400}y(2/5)+\dfrac{1}{500}y(1/2) \\[0.4cm] \quad +\dfrac{1}{90}I_{0^+}^{1/2}y(3/5) +\dfrac{1}{70}I_{0^+}^{1/2}y(7/10) +\dfrac{1}{20}I_{0^+}^{1/2}y(4/5), \\[0.4cm] y(1) = \dfrac{1}{35}x(1/10)+\dfrac{1}{100}x(1/5)+\dfrac{1}{21}x(3/10)+\dfrac{1}{70}x(2/5)+\dfrac{1}{500}x(1/2)\\[0.4cm] \quad + \dfrac{1}{100}I_{0^+}^{1/3}x(3/5)+\dfrac{1}{200} I_{0^+}^{1/3}x(7/10) +\dfrac{1}{300}I_{0^+}^{1/3}x(4/5). \end{array} \right. \end{equation}$ $

(4.1)

Here $ \alpha_{1} = 1/2, \alpha_{2} = 1/6, \alpha_{3} = 1/2, \alpha_{4} = 1/2,$ $ \beta_{1} = 3/4, \beta_{2} = 4/5, \beta_{3} = 3/4, \beta_{4} = 1/7, \lambda_{1} = 1/90, \lambda_{2} = 1/100,$ $\epsilon_{1} = 7/8 = \epsilon_{3}, q_{k} = 1/2, p_{k} = 1/3,$ $k = 1, 2, 3, m = 5, n = 3, \eta_{1} = 1/10, \eta_{2} = 1/5, \eta_{3} = 3/10,$ $\eta_{4} = 2/5, \eta_{5} = 1/2, \xi_{1} = 3/5, \xi_{2} = 7/10, \xi_{3} = 4/5,$ $\mu_{1} = 1/100, \mu_{2} = 1/200,$ $\mu_{3} = 1/300, \mu_{4} = 1/400, \mu_{5} = 1/500,$ $ \nu_{1} = 1/90, \nu_{2} = 1/70, \nu_{3} = 1/20,$ $\delta_{1} = 1/35, \delta_{2} = 1/100, \delta_{3} = 1/21, \delta_{4} = 1/70, \delta_{5} = 1/500, \theta_{1} = 1/100,$ $ \theta_{2} = 1/200, \theta_{3} = 1/300, \zeta_{1} = 1/2, \zeta_{2} = 1/3 $.

With the given data, it is found that $ |\Delta| = 1.1419, \sigma_{1} = 2.2278, \sigma_{3} = 0.1836, \sigma_{2} = 0.1096,$ $\sigma_{4} = 2.0132, \sigma_{5} = 0.0009, \sigma_{6} = 0.0025, \sigma_{7} = 1.8560, \sigma_{8} = 0.0475, \sigma_{9} = 0.1328,$ $\sigma_{10} = 1.1252. $

$ ({\bf a)} $ For illustrating Theorem 3.1, we take

$\begin{equation} \left\{ \begin{array}{ll} f_{1}(t, x, y) = \dfrac{2\arctan x+\pi}{14\pi(1+t)}+\dfrac{1}{7(t+\pi)}\sin|y|, \\ f_{2}(t, x, y) = \dfrac{1}{7}\arctan x+\dfrac{3}{(21+t)}\dfrac{|y|}{(1+|y|)}+\dfrac{t^3}{(1+t^2)}, \\ g_{1}(x, y) = \dfrac{1}{12} \Big(\dfrac{|x|}{(1+|x|)}+|y| \Big), \\ g_{2}(x, y) = \dfrac{1}{17} \Big(\sin|x|+\arctan|y| \Big). \end{array} \right. \end{equation}$ $

(4.2)

It can easily be verified that $ f_{1}, f_{2} $ satisfy the condition $ (H_1) $ with $ K_{1} = 1/7\pi, K_{2} = 1/7, $ respectively and $ g_{1}, g_{2} $ satisfy the condition $ (H_2) $ with $ L_{1} = 1/12, L_{2} = 1/17, $ respectively. Furthermore

$ (\sigma_{1}+\sigma_{3})K_{1}+(\sigma_{4}+\sigma_{2})K_{2}+(\sigma_{9}+\sigma_{7})L_{1}+(\sigma_{10}+\sigma_{8})L_{2} +\sigma_{5}+\sigma_{6}\approx 0.65102 < 1. $

Clearly the hypotheses of Theorem 3.1 are satisfied and hence it follows by its conclusion that the system (4.1) with $ f_{1}(t, x, y), f_{2}(t, x, y), g_{1}(x, y) $ and $ g_{2}(x, y) $ given by (4.2) has a unique solution on $ [0, 1]. $ On the other hand, one can deduce that the system (4.1) with (4.2) has at least one solution on $ [0, 1] $ by the application of Remark 3.1 with $ (\sigma_{1}+\sigma_{3})K_{1}+(\sigma_{4}+\sigma_{2})K_{2}+(\sigma_{9}+\sigma_{7})L_{1}+(\sigma_{10}+\sigma_{8})L_{2}\approx 0.6476 < 1. $

$ ({\bf b)} $ As an application of Theorem 3.2, consider

$\begin{equation} \left\{ \begin{array}{ll} f_{1}(t, x, y) = \dfrac{\arctan x}{10(t^2+1)}+ \dfrac{\sin |y|}{17(1+t)} , \\ f_{2}(t, x, y) = \dfrac{2}{\sqrt{t^{2}+2}} + \dfrac{2|x|}{5\pi(8+t)(1+|x|)}, \\ g_{1}(x, y) = \dfrac{|x|}{2(1+|x|)}+ \dfrac{1}{6} \arctan y , \\ g_{2}(x, y) = \dfrac{1}{3} e^{-|x|}+ \dfrac{1}{7} \cos|x|. \end{array} \right. \end{equation}$ $

(4.3)

Using the given values, we find that the assumption ($ H_{4} $) is satisfied since $ |f_{1}(t, x, y)| $ $ \leqslant \dfrac{\pi}{20(1+t^2)}+\dfrac{1}{17(t+1)} = \phi_{1}(t) $ and $ |f_{2}(t, x, y)|\leqslant \dfrac{2}{\sqrt{t^{2}+2}} + \dfrac{2}{5\pi(8+t)} = \phi_{2}(t), $ $ |g_{1}(x, y)|\leqslant (6+\pi)/12 = \Lambda_{1}, $ $ |g_{2}(x, y)|\leqslant 10/21 = \Lambda_{2}. $ Also $ (\sigma_{5}+\sigma_{6})\approx 0.0034 < 1 $ holds true. As all the assumptions of Theorem 3.2 are satisfied, so its conclusion implies that the system (4.1) with the nonlinearities (4.3) has at least one solution on $ [0, 1]. $

$ ({\bf c)} $ In order to demonstrate the application of Theorem 3.3, let us choose

$\begin{equation} \left\{ \begin{array}{ll} f_{1}(t, x, y) = \arctan x+ \dfrac{e^{-t}|x|^2}{17(1+|x|)} +\dfrac{1}{26} y \cos x, \\ f_{2}(t, x, y) = \dfrac{2}{\sqrt{t^{2}+2}}+\dfrac{2}{\pi(8+t)}x \; \arctan y + \dfrac{|x||y|}{5(1+|x|)}, \\ g_{1}(x, y) = \ln 7+ \dfrac{1}{21} x \sin |y| + \dfrac{1}{13} y , \\ g_{2}(x, y) = 3e^{-|x|}+ \dfrac{1}{7} x \cos y+ \dfrac{1}{11} y \arctan x. \end{array} \right. \end{equation}$ $

(4.4)

Obviously ($ H_{3} $) holds true with positive values of $ u_{0}, v_{0}, \omega_{0}, \tau_{0} $ and $ u_{1} = 1/17, u_{2} = 1/26, v_1 = 1/8, v_2 = 1/5, \omega_{1} = 1/21, \omega_{2} = 1/13, \tau_{1} = 1/7, \tau_{2} = \pi/22. $ Also, $ (\sigma_{1}+\sigma_{3})u_{1}+(\sigma_{2}+\sigma_{4})v_{1}+(\sigma_{7}+\sigma_{9})\omega_{1}+(\sigma_{8}+\sigma_{10})\tau_{1}+(\sigma_{5}+\sigma_{6})\approx 0.6728 < 1, $ and $ (\sigma_{1}+\sigma_{3})u_{2}+(\sigma_{2}+\sigma_{4})v_{2}+(\sigma_{7}+\sigma_{9})\omega_{2}+(\sigma_{8}+\sigma_{10})\tau_{2}+(\sigma_{5}+\sigma_{6})\approx 0.8412 < 1. $ As the hypothesis of Theorem 3.3 is verified, therefore we deduce by its conclusion that there exists at least one solution of the system (4.1) with $ f_{1}, f_{2}, g_{1} $ and $ g_2 $ given by (4.4).

5. Conclusions

In the present research work, we investigated the existence and uniqueness of solutions for a new coupled system of multi-term Hilfer fractional differential equations of different orders involving non-integral and autonomous type Riemann-Liouville mixed integral nonlinearities equipped with nonlocal coupled multi-point and Riemann-Liouville integral boundary conditions. Firstly, we proved an auxiliary result concerning the linear variant of the given problem, helping us to transform the problem at hand into a fixed point problem. Then we proved the existence of a unique solution for the given problem by applying Banach's contraction mapping principle and derived the existence results by means of Krasnosel'ski${\rm{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over i} }}$'s fixed point theorem and Leray-Schauder nonlinear alternative. All the obtained results are well illustrated by numerical examples. Our results are new and enrich the literature on nonloacl nonlinear integral boundary value problems for Hilfer fractional differential equations.

Acknowledgments

This project was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, Saudi Arabia under grant no. (KEP-PHD-80-130-42). The authors, therefore, acknowledge with thanks DSR technical and financial support.

Conflict of interest

The authors declare no conflict of interest.

References

[1]	van der Voet GB, Todorov TI, Centeno JA, et al. (2007) Metals and health: a clinical toxicological perspective on tungsten and review of the literature. Mil Med 172: 1002-1005. doi: 10.7205/MILMED.172.9.1002
[2]	Agency for Toxic Substances and Disease Registry (ATSDR) (2005) Toxicological profile for tungsten. Atlanta, GA: U.S. Department of Health and Human Services, Public Health Service.
[3]	United States Fish and Wildlife Service, Nontoxic shot regulations for hunting waterfowl and coots in the U.S.. Available from: www.regulations.gov/#!documentDetail;D=FWS-R9-MB-2011-0077-0009.
[4]	Kraabel BJ, Miller MW, Getzy DM, et al. (1996) Effects of embedded tungsten-bismuth-tin shot and steel shot on mallards (Anas platyrhynchos). J Wildlife Dis 32: 1-8. doi: 10.7589/0090-3558-32.1.1
[5]	Kelly ME, Fitzgerald SD, Aulerich RJ, et al. (1998) Acute effects of lead, steel, tungsten-iron, and tungsten-polymer shot administered to game-farm mallards. J Wildlife Dis 34: 673-687. doi: 10.7589/0090-3558-34.4.673
[6]	Mitchell RR, Fitzgerald SD, Aulerich RJ, et al. (2001) Hematological effects and metal residue concentrations following chronic dosing with tungsten-iron and tungsten-polymer shot in adult game-farm mallards. J Wildlife Dis 37: 459-467. doi: 10.7589/0090-3558-37.3.459
[7]	Mitchell RR, Fitzgerald SD, Aulerich RJ, et al. (2001) Reproductive effects and duckling survivability following chronic dosing with tungsten-iron and tungsten-polymer shot in adult game-farm mallards. J Wildlife Dis 37: 468-474. doi: 10.7589/0090-3558-37.3.468
[8]	Mitchell RR, Fitzgerald SD, Aulerich RJ, et al. (2001) Health effects following chronic dosing with tungsten-iron and tungsten-polymer shot in adult game-farm mallards. J Wildlife Dis 37: 451-458.
[9]	Brewer L, Fairbrother A, Clark J, et al. (2003) Acute toxicity of lead, steel, and an iron-tungsten-nickel shot to mallard ducks (Anas platyrhynchos). J Wildlife Dis 39: 638-648.
[10]	Kerley CR, Easterly CE, Eckerman KF, et al. (1996) Environmental acceptability of high-performance alternatives for depleted uranium penetrators. ORNL/TM-13286. Oak Ridge National Laboratory. Available from: http://www.osti.gov/scitech/biblio/464128-environmental-acceptability-high-performance-alternatives-depleted-uranium-penetrators.
[11]	Kalinich JF, Emond CA, Dalton TK, et al. (2005) Embedded weapons-grade tungsten alloy shrapnel rapidly induces metastatic high-grade rhabdomyosarcomas in F344 rats. Environ Health Perspect 113: 729-734. doi: 10.1289/ehp.7791
[12]	Schuster BE, Roszell LE, Murr LE, et al. (2012) In vivo corrosion, tumor outcome, and microarray gene expression for two types of muscle-implanted tungsten alloys. Toxicol Appl Pharmacol 265: 128-138. doi: 10.1016/j.taap.2012.08.025
[13]	Emond CA, Vergara VB, Lombardini ED, et al. (2015) Induction of rhabdomyosarcoma by embedded military-grade tungsten/nickel/cobalt not by tungsten/nickel/iron in the B6C3F1 mouse. Int J Toxicol 34: 44-54. doi: 10.1177/1091581814565038
[14]	Koutsospyros A, Braida W, Christodoulatos C, et al. (2006) A review of tungsten: From environmental obscurity to scrutiny. J Hazard Mater 136: 1-19. doi: 10.1016/j.jhazmat.2005.11.007
[15]	Gbaruko BC, Igwe JC (2007) Tungsten: Occurrence, chemistry, environmental, and health exposure issues. Global J Environ Res 1: 27-32.
[16]	Ogundipe A, Greenberg B, Braida W, et al. (2006) Morphological characterization and spectroscopic studies of the corrosion behaviour of tungsten heavy alloys. Corr Sci 48: 3281-3297. doi: 10.1016/j.corsci.2005.12.004
[17]	Dermatas D, Braida W, Christodoulatos C, et al. (2004) Solubility, sorption, and soil respiration effects of tungsten and tungsten alloys. Environ Forensics 5: 5-13. doi: 10.1080/15275920490423980
[18]	Bednar AJ, Jones WT, Boyd RE, et al. (2008) Geochemical parameters influencing tungsten mobility in soils. J Environ Qual 37: 229-233. doi: 10.2134/jeq2007.0305
[19]	Strigul N (2010) Does speciation matter for tungsten ecotoxicology? Ecotoxicol Environ Safety 73: 1099-1113.
[20]	Institute of Laboratory Animal Resources (2010) Guide for the Care and Use of Laboratory Animals. 8th edition. Washington, DC: National Academy Press.
[21]	Rao GN (1996) New diet (NTP-2000) for rats in the National Toxicology Program toxicity and carcinogenicity studies. Fund Appl Toxicol 32: 102-108. doi: 10.1006/faat.1996.0112
[22]	Hockley AD, Goldin JH, Wake MJC, et al. (1990) Skull repair in children. Pediatric Neurosurg 16: 271-275.
[23]	Johansson CB, Hansson HA, Albrektsson T (1990) Qualitative interfacial study between bone and tantalum, niobium or commercially pure titanium. Biomaterials 11: 277-280. doi: 10.1016/0142-9612(90)90010-N
[24]	Strecker EP, Hagan B, Liermann D, et al. (1993) Iliac and femoropopliteal vascular occlusive disease treated with flexible tantalum stents. Cardiovasc Intervent Radiol 16: 158-164. doi: 10.1007/BF02641885
[25]	Emond CA, Vergara VB, Lombardini ED, et al. (2015) The role of the component metals in the toxicity of military-grade tungsten alloy. Toxics 3: 499-514. doi: 10.3390/toxics3040499
[26]	Witten ML, Sheppard PR, Witten BL (2012) Tungsten toxicity. Chemico-Biol Interact 196: 87-88. doi: 10.1016/j.cbi.2011.12.002
[27]	Rubin CS, Holmes AK, Belson MG, et al. (2007) Investigating childhood leukemia in Churchill County, Nevada. Environ Health Perspect 115: 151-157.
[28]	Kalinich JF, Vergara VB, Emond CA (2008) Urinary and serum metal levels as indicators of embedded tungsten alloy fragments. Mil Med 173: 754-758. doi: 10.7205/MILMED.173.8.754
[29]	Miller AC, Mog S, McKinney LA, et al. (2001) Neoplastic transformation of human osteoblast cells to the tumorigenic phenotype by heavy metal-tungsten alloy particles: induction of genotoxic effects. Carcinogenesis 22: 115-125. doi: 10.1093/carcin/22.1.115
[30]	Miller AC, Xu J, Prasanna PGS, et al. (2002) Potential late health effects of the heavy metals, depleted uranium and tungsten, used in armor piercing munitions: comparison of neoplastic transformation and genotoxicity using the known carcinogen nickel. Mil Med 167: 120-122.
[31]	Kalinich JF, Kasper CE (2014) Do metals that translocate to the brain exacerbate traumatic brain injury? Med Hypoth 83: 558-562.
[32]	Leggett RW (1997) A model of the distribution and retention of tungsten in the human body. Sci Total Environ 206: 147-165. doi: 10.1016/S0048-9697(97)80006-X
[33]	Emond CA, Kalinich JF (2012) Biokinetics of embedded surrogate radiological dispersal device material. Health Phys 102: 124-136.
[34]	Endoh H, Kaneko T, Nakamura H, et al. (2000) Improved cardiac contractile functions in hypoxia-reoxygenation in rats treated with low concentration Co²⁺. Am J Physiol Heart Circ Physiol 279: H2713-H2719.
[35]	Rakusan K, Cicutti N, Kolar F (2001) Cardiac function, microvascular structure, and capillary hematocrit in hearts of polycythemic rats. Am J Physiol Heart Circ Physiol 281: H2425-H2431.
[36]	Fastje CD, Harper K, Terry C, et al. (2012) Exposure to sodium tungstate and Respiratory Syncytial Virus results in hematological/immunological disease in C57BL/6J mice. Chemico-Biol Interact 196: 89-95. doi: 10.1016/j.cbi.2011.04.008
[37]	Kelly ADR, Lemaire M, Young YK, et al. (2013) In vivo tungsten exposure alters B-cell development and increases DNA damage in murine bone marrow. Toxicol Sci 131: 434-446. doi: 10.1093/toxsci/kfs324
[38]	Jalaguier-Coudray A, Cohen M, Thomassin-Piana J, et al. (2015) Calcifications and tungsten deposits after breast-conserving surgery and intraoperative radiotherapy for breast cancer. Eur J Radiol 84: 2521-2525. doi: 10.1016/j.ejrad.2015.10.004
[39]	Bolt AM, Sabourin V, Molina F, et al. (2015) Tungsten targets the tumor microenvironment to enhance breast cancer metastasis. Toxicol Sci 143: 165-177. doi: 10.1093/toxsci/kfu219
[40]	Harris RM, Williams TD, Hodges NJ, et al. (2011) Reactive oxygen species and oxidative DNA damage mediate the cytotoxicity of tungsten-nickel-cobalt alloys in vitro. Toxicol Appl Pharmacol 250: 19-28. doi: 10.1016/j.taap.2010.09.020
[41]	Bardack S, Dalgard CL, Kalinich JF, et al. (2014) Genotoxic changes to rodent cells exposed in vitro to tungsten, nickel, cobalt, and iron. Int J Environ Res Pub Health 11: 2922-2940. doi: 10.3390/ijerph110302922
[42]	Harris RM, Williams TD, Waring RH, et al. (2015) Molecular basis of carcinogenicity of tungsten alloy particles. Toxicol Appl Pharmacol 283: 223-233. doi: 10.1016/j.taap.2015.01.013
[43]	Laulicht F, Brocato J, Cartularo L, et al. (2015) Tungsten-induced carcinogenesis in human bronchial epithelial cells. Toxicol Appl Pharmacol 288: 33-39. doi: 10.1016/j.taap.2015.07.003
[44]	Guandalini GS, Zhang L, Fornero E, et al. (2011) Tissue distribution of tungsten in mice following oral exposure to sodium tungstate. Chem Res Toxicol 24: 488-493. doi: 10.1021/tx200011k
[45]	McInturf SM, Bekkedal MYV, Wilfong E, et al. (2011) The potential reproductive, neurobehavioral and systemic effects of soluble sodium tungstate exposure in Sprague-Dawley rats. Toxicol Appl Pharmacol 254: 133-137. doi: 10.1016/j.taap.2010.04.021
[46]	Radcliffe PM, Leavens TL, Wagner DJ, et al. (2010) Pharmacokinetics of radiolabeled tungsten (¹⁸⁸W) in male Sprague-Dawley rats following acute sodium tungstate inhalation. Inhalation Toxicol 22: 69-76. doi: 10.3109/08958370902913237
[47]	McDonald JD, Weber WM, Marr R, et al. (2007) Disposition and clearance of tungsten after single-dose oral and intravenous exposure in rodents. J Toxicol Environ Health A 70: 829-836. doi: 10.1080/15287390701211762
[48]	NTP Range-Finding Report: Immunotoxicity of Sodium Tungstate Dihydrate in Female B6C3F1/N Mice (CASRN: 10213-10-2) National Toxicology Program, U.S. Department of Health and Human Services. Available from: ntp.niehs.nih.gov/testing/types/imm/abstract/i03038/index.html.
[49]	U.S. Environmental Protection Agency (2008) Emerging contaminant tungsten. Fact sheet 505-F-070-005. Available from: http://nepis.epa.gov/Exe/ZyPDF.cgi/P1000L3K.PDF?Dockey=P1000L3K.PDF.

This article has been cited by:

1.	Ahmed Alsaedi, Madeaha Alghanmi, Bashir Ahmad, Boshra Alharbi, Uniqueness of solutions for a ψ-Hilfer fractional integral boundary value problem with the p-Laplacian operator, 2023, 56, 2391-4661, 10.1515/dema-2022-0195
2.	Muhammad Bilal Khan, Hakeem A. Othman, Michael Gr. Voskoglou, Lazim Abdullah, Alia M. Alzubaidi, Some Certain Fuzzy Aumann Integral Inequalities for Generalized Convexity via Fuzzy Number Valued Mappings, 2023, 11, 2227-7390, 550, 10.3390/math11030550
3.	Muhammad Bilal Khan, Hakeem A. Othman, Gustavo Santos-García, Muhammad Aslam Noor, Mohamed S. Soliman, Some new concepts in fuzzy calculus for up and down λ-convex fuzzy-number valued mappings and related inequalities, 2023, 8, 2473-6988, 6777, 10.3934/math.2023345
4.	Muhammad Bilal Khan, Hakeem A. Othman, Gustavo Santos-García, Tareq Saeed, Mohamed S. Soliman, On fuzzy fractional integral operators having exponential kernels and related certain inequalities for exponential trigonometric convex fuzzy-number valued mappings, 2023, 169, 09600779, 113274, 10.1016/j.chaos.2023.113274
5.	Bashir Ahmad, Sotiris K. Ntouyas, 2024, Chapter 8, 978-3-031-62512-1, 289, 10.1007/978-3-031-62513-8_8
6.	Haitham Qawaqneh, Hasanen A. Hammad, Hassen Aydi, Exploring new geometric contraction mappings and their applications in fractional metric spaces, 2024, 9, 2473-6988, 521, 10.3934/math.2024028

Reader Comments

Your name:*

Email:*
© 2016 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Environmental Science

1.2 2.7

Metrics

Article views(6370) PDF downloads(1194) Cited by(5)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Tables(7)

AIMS Environmental Science

Tissue distribution patterns of solubilized metals from internalized tungsten alloy in the F344 rat

Related Papers:

Abstract

1. Introduction

2. Preliminaries

3. Existence and uniqueness results

4. Examples

5. Conclusions

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

AIMS Environmental Science

Tissue distribution patterns of solubilized metals from internalized tungsten alloy in the F344 rat

Related Papers:

Abstract

1. Introduction

2. Preliminaries

3. Existence and uniqueness results

4. Examples

5. Conclusions

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog