The impact of media on the spatiotemporal pattern dynamics of a reaction-diffusion epidemic model

Xin-You Meng; Tao Zhang; Xin-You Meng; Tao Zhang

doi:10.3934/mbe.2020223

Mathematical Biosciences and Engineering

2020, Volume 17, Issue 4: 4034-4047. doi: 10.3934/mbe.2020223

Previous Article Next Article

Research article Special Issues

The impact of media on the spatiotemporal pattern dynamics of a reaction-diffusion epidemic model

Xin-You Meng ^,,
Tao Zhang

School of Science, Lanzhou University of Technology, Lanzhou, Gansu 730050, China

Received: 03 March 2020 Accepted: 27 May 2020 Published: 03 June 2020

In this paper, a reaction-diffusion SI epidemic model with media impact is considered. The boundedness of system and the existence of the state are given. The local stabilities of the endemic states are analyzed. Sufficient conditions of the occurrence of the Turing pattern are obtained by the center manifold theorem and normal form method. Some numerical simulations are given to check in the theoretical results. We find that the influence of media not only inhibits the spread of infectious diseases, but also effects the spatial steady-state of model.

Keywords:

Citation: Xin-You Meng, Tao Zhang. The impact of media on the spatiotemporal pattern dynamics of a reaction-diffusion epidemic model[J]. Mathematical Biosciences and Engineering, 2020, 17(4): 4034-4047. doi: 10.3934/mbe.2020223

Related Papers:

[1]	M. Latha Maheswari, K. S. Keerthana Shri, Mohammad Sajid . Analysis on existence of system of coupled multifractional nonlinear hybrid differential equations with coupled boundary conditions. AIMS Mathematics, 2024, 9(6): 13642-13658. doi: 10.3934/math.2024666
[2]	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
[3]	Zaid Laadjal, Fahd Jarad . Existence, uniqueness and stability of solutions for generalized proportional fractional hybrid integro-differential equations with Dirichlet boundary conditions. AIMS Mathematics, 2023, 8(1): 1172-1194. doi: 10.3934/math.2023059
[4]	M. Manigandan, Subramanian Muthaiah, T. Nandhagopal, R. Vadivel, B. Unyong, N. Gunasekaran . Existence results for coupled system of nonlinear differential equations and inclusions involving sequential derivatives of fractional order. AIMS Mathematics, 2022, 7(1): 723-755. doi: 10.3934/math.2022045
[5]	Subramanian Muthaiah, Dumitru Baleanu, Nandha Gopal Thangaraj . Existence and Hyers-Ulam type stability results for nonlinear coupled system of Caputo-Hadamard type fractional differential equations. AIMS Mathematics, 2021, 6(1): 168-194. doi: 10.3934/math.2021012
[6]	Manal Elzain Mohamed Abdalla, Hasanen A. Hammad . Solving functional integrodifferential equations with Liouville-Caputo fractional derivatives by fixed point techniques. AIMS Mathematics, 2025, 10(3): 6168-6194. doi: 10.3934/math.2025281
[7]	Cuiying Li, Rui Wu, Ranzhuo Ma . Existence of solutions for Caputo fractional iterative equations under several boundary value conditions. AIMS Mathematics, 2023, 8(1): 317-339. doi: 10.3934/math.2023015
[8]	Abdelkader Amara . Existence results for hybrid fractional differential equations with three-point boundary conditions. AIMS Mathematics, 2020, 5(2): 1074-1088. doi: 10.3934/math.2020075
[9]	Djamila Chergui, Taki Eddine Oussaeif, Merad Ahcene . Existence and uniqueness of solutions for nonlinear fractional differential equations depending on lower-order derivative with non-separated type integral boundary conditions. AIMS Mathematics, 2019, 4(1): 112-133. doi: 10.3934/Math.2019.1.112
[10]	Abdulwasea Alkhazzan, Wadhah Al-Sadi, Varaporn Wattanakejorn, Hasib Khan, Thanin Sitthiwirattham, Sina Etemad, Shahram Rezapour . A new study on the existence and stability to a system of coupled higher-order nonlinear BVP of hybrid FDEs under the $p$ -Laplacian operator. AIMS Mathematics, 2022, 7(8): 14187-14207. doi: 10.3934/math.2022782

Abstract

1. Introduction

Widespread applications of fractional calculus significantly contributed to the popularity of the subject. Fractional order operators are nonlocal in nature and give rise to more realistice and informative mathematical modeling of many real world phenomena, in contrast to their integer-order counterparts, for instance, see ^[13,21,31].

Nonlinear fractional order boundary value problems appear in a variety of fields such as applied mathematics, physical sciences, engineering, control theory, etc. Several aspects of these problems, such as existence, uniqueness and stability, have been explored in recent studies ^{[5,6,7,14,22,24,26,28,32]}.

Coupled nonlinear fractional differential equations find their applications in various applied and technical problems such as disease models ^[8,10,29], ecological models ^[18], synchronization of chaotic systems ^[11,33], nonlocal thermoelasticity ^[30], etc. Hybrid fractional differential equations also received significant attention in the recent years, for example, see ^{[2,3,9,15,16,17,19,20,23]}.

The concept of slits-strips conditions introduced by Ahmad et al. in ^[1] is a new idea and has useful applications in imaging via strip-detectors ^[25] and acoustics ^[27].

In ^[1], the authors investigated the following strips-slit problem:

$^cD^px(t) = f_1(t, x(t)), \ n-1 \lt p\le n, \ t\in [0, 1]\\ x(0) = 0, \ x'(0) = 0, \ x^{\prime\prime}(0) = 0, .... \ x^{(n-2)}(0) = 0, \\ x(\xi) = a_1\int_{0}^{\eta}x(s)ds+a_2\int_{\xi_1}^{1}x(s)ds, \ 0 \lt \eta \lt \xi \lt \xi_1 \lt 1,$

where $^cD^p$ denotes the Caputo fractional derivative of order $p, f_1$ is a given continuous function and $a_1, a_2\in\mathfrak{R}$ .

In 2017, Ahmad et al. ^[4] studied a coupled system of nonlinear fractional differential equations

$^cD^{\alpha}x(t) = f_1(t, x(t), y(t)), \ t\in [0, 1], \ 1 \lt \alpha\le 2, \\ ^cD^{\beta}y(t) = f_2(t, x(t), y(t)), \ t\in [0, 1], \ 1 \lt \beta\le 2,$

supplemented with the coupled and uncoupled boundary conditions of the form:

$x(0) = 0, \ x(a_1) = d_1\int_{0}^{\eta}y(s)ds+d_2\int_{\xi_1}^{1}y(s)ds, \ 0 \lt \eta \lt a_1 \lt \xi_1 \lt 1, \\ y(0) = 0, \ y(a_1) = d_1\int_{0}^{\eta}x(s)ds+d_2\int_{\xi_1}^{1}x(s)ds, \ 0 \lt \eta \lt a_1 \lt \xi_1 \lt 1,$

and

$x(0) = 0, \ x(a_1) = d_1\int_{0}^{\eta}x(s)ds+d_2\int_{\xi_1}^{1}x(s)ds, \ 0 \lt \eta \lt a_1 \lt \xi_1 \lt 1, \\ y(0) = 0, \ y(a_1) = d_1\int_{0}^{\eta}y(s)ds+d_2\int_{\xi_1}^{1}y(s)ds, \ 0 \lt \eta \lt a_1 \lt \xi_1 \lt 1,$

where $^cD^{\alpha}$ and $^cD^{\beta}$ denote the Caputo fractional derivatives of orders $\alpha$ and $\beta$ respectively, $f_1, f_2:[0, 1]\times \mathfrak{R}\times \mathfrak{R}\to \mathfrak{R}$ are given continuous functions and $d_1, d_2$ are real constants.

In this article, motivated by aforementioned works, we introduce and study the following hybrid nonlinear fractional differential equations:

$\begin{array}{l} \left. \begin{array}{cl} ^cD^{\gamma}\left[u(t)-h_1(t, u(t), v(t))\right] = \theta_1(t, u(t), v(t)), \ t\in[0, 1], \ 1 \lt \gamma\le 2, \\ ^cD^{\delta}\left[v(t)-h_2(t, u(t), v(t))\right] = \theta_2(t, u(t), v(t)), \ t\in[0, 1], \ 1 \lt \delta\le 2, \end{array} \right. \end{array}$

(1.1)

equipped with coupled slit-strips-type integral boundary conditions:

$\begin{array}{l} \left. \begin{array}{cl} u(0) = 0, \ u(\eta) = \omega_1\int_{0}^{\xi_1}v(s)ds+\omega_2\int_{\xi_2}^{1}v(s)ds, \ 0 \lt \xi_1 \lt \eta \lt \xi_2 \lt 1, \\ v(0) = 0, \ v(\eta) = \omega_1\int_{0}^{\xi_1}u(s)ds+\omega_2\int_{\xi_2}^{1}u(s)ds, \ 0 \lt \xi_1 \lt \eta \lt \xi_2 \lt 1, \end{array} \right. \end{array}$

(1.2)

where $^cD^{\gamma}$ , $^cD^{\delta}$ denote the Caputo fractional derivatives of orders $\gamma$ and $\delta$ respectively, $\theta_i, h_i:[0, 1]\times \mathfrak{R}\times\mathfrak{R}\to\mathfrak{R}$ are given continuous functions with $h_i(0, u(0), v(0)) = 0, \; i = 1, 2$ and $\omega_1, \omega_2$ are real constants.

We arrange the rest of the paper as follows. In section 2, we present some definitions and obtain an auxiliary result, while section 3 contains the main results for the problems (1.1) and (1.2). Section 4 is devoted to the illustrative examples for the derived results.

2. An auxiliary lemma

Let us first recall some related definitions ^[21].

Definition 2.1. For a locally integrable real-valued function $g_1:{[a, \infty)}\longrightarrow \mathbb{R}$ , we define the Riemann-Liouville fractional integral of order $\sigma > 0$ as

$I^{\sigma}g_1(t) = \frac{1}{\Gamma(\sigma)}\int_{0}^{t}\frac{g_1(\tau)}{(t-\tau)^{1-\sigma}}d \tau, \ \sigma \gt 0,$

where $\Gamma$ is the Euler's gamma function.

Definition 2.2. The Caputo derivative of order $\sigma$ for an $n$ -times continuously differentiable function $g_1:[0, \infty)\to \mathfrak{R}$ is defined by

$^cD^{\sigma}g_1(t) = \frac{1}{\Gamma(n-\sigma)}\int_{0}^{t}(t-\tau)^{n-\sigma-1}g_1^{(n)}(\tau)d \tau, \ n-1 \lt \sigma \lt n, \ n = [\sigma]+1$

where $[\sigma]$ is the integer part of a real number.

Lemma 2.1. For $\chi_i, \Phi_i \in C([0, 1], \mathfrak{R})$ with $\chi_i(0) = 0, i = 1, 2,$ the following linear system of equations:

$\begin{array}{l} \left. \begin{array}{cl} ^cD^{\gamma}\left[u(t)-\chi_1(t)\right] = \Phi_1(t), \ t\in[0, 1], \ 1 \lt \gamma\le 2, \\ ^cD^{\delta}\left[v(t)-\chi_2(t)\right] = \Phi_2(t), \ t\in[0, 1], \ 1 \lt \delta\le 2, \end{array} \right. \end{array}$

(2.1)

equipped with coupled slit-strips-type integral boundary conditions (1.2), is equivalent to the integral equations:

$\begin{array}{l} u(t) = \frac{t}{\eta^2-\Delta^2}\Bigg[\eta\Bigg\{\omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\Phi_2(\tau)d\tau+\chi_2(s)\Bigg)ds \\ + \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\Phi_2(\tau)d\tau+\chi_2(s)\Bigg)ds -\int_{0}^{\eta}\frac{(\eta-s)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(s)ds-\chi_1(\eta)\Bigg\} \\ + \Delta\Bigg\{\omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(\tau)d\tau+\chi_1(s)\Bigg)ds \\ + \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(\tau)d\tau+\chi_1(s)\Bigg)ds -\int_{0}^{\eta}\frac{(\eta-s)^{\delta-1}}{\Gamma(\delta)}\Phi_2(s)ds-\chi_2(\eta)\Bigg\} \Bigg] \\ + \int_{0}^{t}\frac{(t-s)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(s)ds+\chi_1(t), \end{array}$

(2.2)

$\begin{array}{l} v(t) = \frac{t}{\eta^2-\Delta^2}\Bigg[\Delta\Bigg\{\omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\Phi_2(\tau)d\tau+\chi_1(s)\Bigg)ds \\ + \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\Phi_2(\tau)d\tau+\chi_1(s)\Bigg)ds -\int_{0}^{\eta}\frac{(\eta-\tau)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(s)ds-\chi_2(\eta)\Bigg\} \\ + \eta\Bigg\{\omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(\tau)d\tau+\chi_2(s)\Bigg)ds \\ + \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(\tau)d\tau+\chi_2(s)\Bigg)ds -\int_{0}^{\eta}\frac{(\eta-s)^{\delta-1}}{\Gamma(\delta)}\Phi_2(s)ds-\chi_2(\eta)\Bigg\} \Bigg] \\ + \int_{0}^{t}\frac{(t-s)^{\delta-1}}{\Gamma(\delta)}\Phi_2(s)ds+\chi_2(t), \end{array}$

(2.3)

where it is assumed that

$\begin{array}{l} \Delta = \frac{1}{2}\left(\omega_1\xi_1^2+\omega_2(1-\xi_2^2)\right)\neq 0. \end{array}$

(2.4)

Proof. Solving the fractional differential equations in (2.1), we get

$\begin{array}{l} u(t) = c_0+c_1t+\int_{0}^{t}\frac{(t-s)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(s)ds+\chi_1(t) \end{array}$

(2.5)

and

$\begin{array}{l} v(t) = c_2+c_3 t+\int_{0}^{t}\frac{(t-s)^{\delta-1}}{\Gamma(\delta)}\Phi_2(s)ds+\chi_2(t), \end{array}$

(2.6)

where $c_0, c_1, c_2, c_3$ $\in \mathfrak{R}$ are arbitrary constants.

Using the conditions $u(0) = 0$ and $v(0) = 0$ in (2.5) and (2.6), we find that $c_0 = 0$ and $c_2 = 0$ . Thus (2.5) and (2.6) become

$\begin{array}{l} u(t) = c_1 t+\int_{0}^{t}\frac{(t-s)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(s)ds+\chi_1(t), \end{array}$

(2.7)

$\begin{array}{l} v(t) = c_3 t+\int_{0}^{t}\frac{(t-s)^{\delta-1}}{\Gamma(\delta)}\Phi_2(s)ds+\chi_2(t), \end{array}$

(2.8)

Making use of the coupled slit-strips-type integral boundary conditions given by (1.2) in (2.7) and (2.8) together with the notation (2.4), we obtain a system of equations:

$\begin{array}{l} \omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\Phi_2(\tau)d\tau+\chi_2(s)\Bigg)ds+ \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\Phi_2(\tau)d\tau+\chi_2(s)\Bigg)ds \\ - \int_{0}^{\eta}\frac{(\eta-s)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(s)ds-\chi_1(\eta)) = c_1\eta-\Delta c_3, \end{array}$

(2.9)

$\begin{array}{l} \omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(\tau)d\tau+\chi_1(s)\Bigg)ds + \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(\tau)d\tau+\chi_1(s)\Bigg)ds \\ - \int_{0}^{\eta}\frac{(\eta-s)^{\delta-1}}{\Gamma(\delta)}\Phi_2(s)ds-\chi_2(\eta) = c_3\eta-\Delta c_1. \end{array}$

(2.10)

Solving the systems (2.9)–(2.10) for $c_1$ and $c_3$ , we find that

$c_1 = \frac{t}{\eta^2-\Delta^2}\Bigg[\eta\Bigg\{\omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\Phi_2(\tau)d\tau+\chi_2(s)\Bigg)ds \nonumber\\ + \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\Phi_2(\tau)d\tau+\chi_2(s)\Bigg)ds -\int_{0}^{\eta}\frac{(\eta-s)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(s)ds-\chi_1(\eta)\Bigg\}\nonumber \\ + \Delta\Bigg\{\omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(\tau)d\tau+\chi_1(s)\Bigg)ds + \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(\tau)d\tau+\chi_1(s)\Bigg)ds \nonumber \\ -\int_{0}^{\eta}\frac{(\eta-s)^{\delta-1}}{\Gamma(\delta)}\Phi_2(s)ds-\chi_2(\eta)\Bigg\} \Bigg]$

and

$c_3 = \frac{t}{\eta^2-\Delta^2}\Bigg[\Delta\Bigg\{\omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\Phi_2(\tau)d\tau+\chi_1(s)\Bigg)ds \nonumber\\ + \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\Phi_2(\tau)d\tau+\chi_1(s)\Bigg)ds -\int_{0}^{\eta}\frac{(\eta-s)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(s)ds-\chi_2(\eta)\Bigg\} \nonumber \\ + \eta\Bigg\{\omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(\tau)d\tau+\chi_2(s)\Bigg)ds + \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\Phi_1(\tau)d\tau+\chi_2(s)\Bigg)ds \nonumber\\ -\int_{0}^{\eta}\frac{(\eta-s)^{\delta-1}}{\Gamma(\delta)}\Phi_2(s)ds-\chi_1(\eta)\Bigg\} \Bigg]$

Inserting the values of $c_1$ and $c_3$ in (2.7) and (2.8) leads to the integral equations (2.2) and (2.3). By direct computation, one can obtain the converse of the lemma. The proof is finished.

3. Main results

Let $\mathcal{W} = \{\tilde{w}(t) : \tilde{w}(t)\in C([0, 1])\}$ be a Banach space equipped with the norm $\left\|\tilde{w}\right\| = \max\{|\tilde{w}(t)|, t\in [0, 1]\}$ , Then the product space $(\mathcal{W}\times\mathcal{W}, \left\|(u, v)\right\|)$ endowed with the norm $\left\|(u, v)\right\| = \left\|u\right\|+\left\|v\right\|$ , $(u, v)\in \mathcal{W}\times\mathcal{W}$ is also a Banach space.

We need the following assumptions to derive the main results.

(A1) Let $\theta_1, \theta_2:[0, 1]\times\mathfrak{R}^2\to\mathfrak{R}$ be continuous and bounded functions and there exists constants $m_i, n_i$ such that, for all $t\in [0, 1]$ and $x_i, y_i\in \mathfrak{R}, i = 1, 2,$

$\left|\theta_1(t, x_1, x_2)-\theta_1(t, y_1, y_2)\right|\le m_1|x_1-y_1|+m_2|x_2-y_2|, \\ \left|\theta_2(t, x_1, x_2)-\theta_2(t, y_1, y_2)\right|\le n_1|x_1-y_1|+n_2|x_2-y_2|.$

(A2) For continuous and bounded functions $h_i,$ i = 1, 2, there exist real constants $\mu_i, \beta_i, \sigma_i > 0$ such that, for all $x_i, y_i \in \mathfrak {R},$ $|h_i(t, x, y)|\le \mu_i$ for all $(t, x, y) \in [0, 1]\times \mathfrak{R}\times\mathfrak{R}$ and

(A3) $\sup\limits_{t\in [0, 1]}\theta_1(t, 0, 0) = \mathcal{N}_1 < \infty \ \ \mbox{and} \, \sup\limits_{t\in [0, 1]}\theta_2(t, 0, 0) = \mathcal{N}_2 < \infty.$

$(A4)$ For the sake of computational convenience, we set

$\mathcal{M}_1 = \frac{1}{\Gamma(\gamma+1)} + \frac{1}{|\eta^2-\Delta^2|}\left[\frac{\eta^{\gamma+1}}{\Gamma(\gamma+1)}+|\Delta||\omega_1|\frac{\xi_1^{\gamma+1}}{\Gamma(\gamma+2)}+\Delta|\omega_2|\frac{1-\xi_2^{\gamma+1}}{\Gamma(\gamma+2) }\right], \\ \mathcal{M}_2 = \frac{1}{|\eta^2-\Delta^2|}\left[\eta|\omega_1|\frac{\xi_1^{\delta+1}}{\Gamma(\delta+2)}+\eta|\omega_2|\frac{1-\xi_2^{\delta+1}}{\Gamma(\delta+2)}+|\Delta|\frac{\eta^{\delta}}{\Gamma(\delta+1)}\right], \\ \mathcal{M}_3 = \frac{1}{|\eta^2-\Delta^2|}\left[\eta|\omega_1|\frac{\xi_1^{\gamma+1}}{\Gamma(\gamma+2)}+\eta|\omega_2|\frac{1-\xi_2^{\gamma+1}}{\Gamma(\gamma+2)}+|\Delta|\frac{\eta^{\gamma}}{\Gamma(\gamma+1)}\right], \\ \mathcal{M}_4 = \frac{1}{\Gamma(\delta+1)} + \frac{1}{|\eta^2-\Delta^2|}\left[\frac{\eta^{\delta+1}}{\Gamma(\delta+1)}+|\Delta||\omega_1|\frac{\xi_1^{\delta+1}}{\Gamma(\delta+2)}+\Delta|\omega_2|\frac{1-\xi_2^{\delta+1}}{\Gamma(\delta+2)}\right], \\ \mathcal{N}_3 = \frac{\eta}{|\eta^2-\Delta^2|}\left[|\omega_1|\xi_1\mu_2+|\omega_2|\mu_2(1-\xi_2)+\mu_1\right]+|\Delta|\left[|\omega_1|\mu_1\xi_1+\mu_1|\omega_2|(1-\xi_2)+\mu_2\right]+\mu_1, \\ \mathcal{N}_4 = \frac{1}{|\eta^2-\Delta^2|}\left[|\Delta||\omega_1|\mu_1\xi_1+|\omega_2|(1-\xi_2)\mu_1+\mu_2\right]+|\eta|\left[|\omega_1|\mu_2\xi_1+|\omega_2|(1-\xi_2)\mu_2+\mu_1\right]+\mu_2, \\ \mathcal{N}_5 = \frac{1}{|\eta^2-\Delta^2|}\left[\eta|\omega_1|\xi_1+|\omega_2|\eta(1-\xi_2)+|\Delta|\right], \\ \mathcal{N}_6 = \frac{1}{|\eta^2-\Delta^2|}\left[|\Delta|\xi_1|\omega_1|+|\Delta||\omega_2|(1-\xi_2)+\eta\right]+1, \\ \mathcal{N}_7 = \frac{1}{|\eta^2-\Delta^2|}\left[|\Delta|+|\eta||\omega_1|\xi_1+|\eta||\omega_2|(1-\xi_2)+|\omega_2|+|\Delta|\right],$

and

$\begin{array}{l} \mathcal{M}_k = \min\Big\{1-(\mathcal{M}_1+\mathcal{M}_3)k_1-(\mathcal{M}_2+\mathcal{M}_4)\lambda_1, \\ 1-(\mathcal{M}_1+\mathcal{M}_3)k_2-(\mathcal{M}_2+\mathcal{M}_4)\lambda_2\Big\}, k_i, \lambda_i\ge 0, \ i = 1, 2. \end{array}$

(3.1)

(A5) $(\mathcal{M}_1+\mathcal{M}_3)(m_1+m_2)+(\mathcal{M}_2+\mathcal{M}_4)(n_1+n_2)+(\mathcal{N}_5+\mathcal{N}_6)(\sigma_1+\sigma_2)+(\mathcal{N}_7+\mathcal{N}_8)(\beta_1+\beta_2) < 1.$

In view of Lemma 1, we define an operator $\mathcal{T}: \mathcal{W}\times\mathcal{W}\to \mathcal{W}\times\mathcal{W}$ associated with the problems (1.1) and (1.2) as follows:

$\begin{array}{l} \mathcal{T}(u, v)(t) = \left( \begin{array}{cl} \mathcal{T}_1(u, v)(t) \\ \mathcal{T}_2(u, v)(t) \end{array} \right), \end{array}$

(3.2)

where

$\mathcal{T}_1(u, v)(t) = \frac{t}{\eta^2-\Delta^2}\Bigg[\eta\Bigg\{\omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\theta_2(\tau, u(\tau), v(\tau))d\tau+h_2(s, u(s), v(s)\Bigg)ds \nonumber\\ + \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\theta_2(\tau, u(\tau), v(\tau))d\tau+h_2(s, u(s), v(s)\Bigg)ds \nonumber \\ -\int_{0}^{\eta}\frac{(\eta-s)^{\gamma-1}}{\Gamma(\gamma)}\theta_1(s, u(s), v(s))ds-h_1(\eta, u(\eta), v(\eta))\Bigg\} \nonumber \\ + \Delta\Bigg\{\omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\theta_1(\tau, u(\tau), v(\tau))d\tau+h_1(s, u(s), v(s)\Bigg)ds\nonumber \\ + \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\theta_1(\tau, u(\tau), v(\tau))d\tau+h_1(s, u(s), v(s)\Bigg)ds \nonumber\\ -\int_{0}^{\eta}\frac{(\eta-s)^{\delta-1}}{\Gamma(\delta)}\theta_2(s, u(s), v(s))ds-h_2(\eta, u(\eta), v(\eta))\Bigg\} \Bigg]\nonumber \\ + \int_{0}^{t}\frac{(t-s)^{\gamma-1}}{\Gamma(\gamma)}\theta_1(s, u(s), v(s))ds+h_1(t, u(t), v(t))$

and

$\mathcal{T}_2(u, v)(t) = \frac{t}{\eta^2-\Delta^2}\Bigg[\Delta\Bigg\{\omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\theta_2(\tau, u(\tau), v(\tau))d\tau+h_1(s, u(s), v(s)\Bigg)ds \nonumber\\ + \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\theta_2(\tau, u(\tau), v(\tau))d\tau+h_1(s, u(s), v(s)\Bigg)ds \nonumber \\ -\int_{0}^{\eta}\frac{(\eta-\tau)^{\gamma-1}}{\Gamma(\gamma)}\theta_1(s, u(s), v(s))ds-h_2(\eta, u(\eta), v(\eta))\Bigg\} \nonumber \\ + \eta\Bigg\{\omega_1\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\theta_1(\tau, u(\tau), v(\tau))d\tau+h_2(s, u(s), v(s)\Bigg)ds\nonumber \\ + \omega_2\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}\theta_1(\tau, u(\tau), v(\tau))d\tau+h_2(s, u(s), v(s)\Bigg)ds \nonumber\\ -\int_{0}^{\eta}\frac{(\eta-s)^{\delta-1}}{\Gamma(\delta)}\theta_2(s, u(s), v(s))ds-h_1(\eta, u(\eta), v(\eta))\Bigg\} \Bigg]\nonumber \\ + \int_{0}^{t}\frac{(t-s)^{\gamma-1}}{\Gamma(\gamma)}\theta_2(s, u(s), v(s))ds+h_2(t, u(t), v(t)).$

Theorem 3.1. Assume that conditions (A1) to (A5) are satisfied. Then there exists a unique solution for the problems (1.1) and (1.2) on $[0, 1]$ .

Proof. In the first step, we establish that $\mathcal{T}\bar{B}_r\subset\bar{B}_r,$ where $\bar{B}_r = \{(u, v)\in\mathcal{W}\times\mathcal{W}:\left\|(u, v)\right\|\le r\}$ is a closed ball with

$r\ge \dfrac{(\mathcal{M}_1+\mathcal{M}_3)\mathcal{N}_1+(\mathcal{M}_2+\mathcal{M}_4)\mathcal{N}_2+\mathcal{N}_3\mu}{1-\left[(\mathcal{M}_1+\mathcal{M}_3)(m_1+m_2)+(\mathcal{M}_2+\mathcal{M}_4)(n_1+n_2) +\mathcal{N}_3(\beta_1+\beta_2)\right]},$

and the operator $\mathcal{T}: \mathcal{W}\times\mathcal{W}\to \mathcal{W}\times\mathcal{W}$ is defined by (3.2). For $(u, v)\in \bar{B}_r$ and $t \in [0, 1]$ , it follows by (A1) that

$\left|\theta_1(t, u(t), v(t))\right| \le \left|\theta_1(t, u(t), v(t))-\theta_1(t, 0, 0)\right| \le m_1||u||+m_2||v||.$

Similarly one can find that $\left|\theta_2(t, u(t), v(t))\right|\le n_1||u||+n_2||v||.$ Then we have

$|\mathcal{T}_1(u, v)(t)| \\ \le\max\limits_{t\in[0, 1]}\Bigg[\frac{t}{|\eta^2-\Delta^2|}\Bigg[\eta\Bigg\{|\omega_1|\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)} |\theta_2(\tau, u(\tau), v(\tau))|d\tau+|h_2(s, u(s), v(s)|\Bigg)ds \nonumber\\ + |\omega_2|\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}|\theta_2(\tau, u(\tau), v(\tau))|d\tau+|h_2(s, u(s), v(s)|\Bigg)ds \nonumber \\ +\int_{0}^{\eta}\frac{(\eta-s)^{\gamma-1}}{\Gamma(\gamma)} |\theta_1(s, u(s), v(s))|ds+|h_1(\eta, u(\eta), v(\eta))|\Bigg\} \nonumber \\ + |\Delta|\Bigg\{|\omega_1|\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}|\theta_1(\tau, u(\tau), v(\tau))|d\tau +|h_1(s, u(s), v(s)|\Bigg)ds\nonumber \\ + |\omega_2|\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)}|\theta_1(\tau, u(\tau), v(\tau))|d\tau +|h_1(s, u(s), v(s)|\Bigg)ds \nonumber\\ +\int_{0}^{\eta}\frac{(\eta-s)^{\delta-1}}{\Gamma(\delta)}|\theta_2(s, u(s), v(s))|ds+|h_2(\eta, u(\eta), v(\eta))|\Bigg\} \Bigg]\nonumber \\ + \int_{0}^{t}\frac{(t-s)^{\gamma-1}}{\Gamma(\gamma)}|\theta_1(s, u(s), v(s))|ds+|h_1(t, u(t), v(t))|\Bigg]\nonumber \\ \le \frac{1}{|\eta^2-\Delta^2|}\Bigg[\eta\Bigg\{|\omega_1|\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)} \left(n_1||u||+n_2||v||+\mathcal{N}_2\right)d\tau+\mu_2\Bigg)ds \nonumber\\ + |\omega_2|\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}| \left(n_1||u||+n_2||v||+\mathcal{N}_2\right)d\tau+\mu_2 \Bigg)ds \nonumber \\ +\int_{0}^{\eta}\frac{(\eta-s)^{\gamma-1}}{\Gamma(\gamma)} \left(m_1||u||+m_2||v||+\mathcal{N}_1\right)d\tau+\mu_1 \Bigg\} \nonumber \\ + |\Delta|\Bigg\{|\omega_1|\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)} \left(m_1||u||+m_2||v||+\mathcal{N}_1\right)d\tau+\mu_1 \Bigg)ds\nonumber \\ + |\omega_2|\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)} \left(m_1||u||+m_2||v||+\mathcal{N}_1\right)d\tau+\mu_1 \Bigg)ds \nonumber\\ +\int_{0}^{\eta}\frac{(\eta-s)^{\delta-1}}{\Gamma(\delta)} \left(n_1||u||+n_2||v||+\mathcal{N}_2\right)ds+\mu_2 \Bigg\}\Bigg]\nonumber \\ + \int_{0}^{t}\frac{(t-s)^{\gamma-1}}{\Gamma(\gamma)} \left(m_1||u||+m_2||v||+\mathcal{N}_1\right)ds + \mu_1 \\ le \frac{1}{|\eta^2-\Delta^2|}\Bigg[\eta|\omega_1|\frac{\xi_1^{\delta+1}}{\Gamma(\delta+2)}+\eta|\omega_2|\frac{1-\xi_2^{\delta+1}}{\Gamma(\delta+2)}+|\Delta|\frac{\eta^{\delta}}{\Gamma(\delta+1)}\Bigg] \left(n_1||u||+n_2||v||+\mathcal{N}_2\right) \nonumber\\ + \Bigg[\frac{1}{|\eta^2-\Delta^2|}\Bigg(\frac{\eta^{\gamma+1}}{\Gamma(\gamma+1)}+|\Delta| |\omega_1|\frac{\xi_1^{\gamma+1}}{\Gamma(\gamma+2)}+|\Delta||\omega_2|\frac{1-\xi_2^{\gamma+1}}{\Gamma(\gamma+2)}\Bigg)+\frac{1}{\Gamma(\gamma+1)}\Bigg] \left(m_1||u||+m_2||v||+\mathcal{N}_1\right) \nonumber\\ + \frac{\eta}{|\eta^2-\Delta^2|}\left(|\omega_1|\mu_2\xi_1+|\omega_2|\mu_2(1-\xi_2)+\mu_1\right)+|\Delta|\left(|\omega_1|\mu_1\xi_1+|\omega_2|\mu_1(1-\xi_2)+\mu_2\right)+\mu_1 \nonumber \\ \le \left(\mathcal{M}_2n_1+\mathcal{M}_1m_1+ \mathcal{M}_2n_2+\mathcal{M}_1m_2\right)r+ \mathcal{M}_2\mathcal{N}_2+\mathcal{M}_1\mathcal{N}_1+\mathcal{N}_3.$

Analogously, one can find that

$|\mathcal{T}_2(u, v)(t)| \le \left(\mathcal{M}_4n_1+\mathcal{M}_3m_1+\mathcal{M}_4n_2+\mathcal{M}_3m_2 \right)r+ \mathcal{M}_4\mathcal{N}_2+\mathcal{M}_3\mathcal{N}_1+\mathcal{N}_4.$

From the foregoing estimates for $\mathcal{T}_1$ and $\mathcal{T}_2$ , we obtain $||\mathcal{T}(u, v)(t)||\le r$ .

Next, for $(u_1, v_1), (u_2, v_2)\in \mathcal{W}\times\mathcal{W}$ and $t\in [0, 1]$ , we get

$\left|\mathcal{T}_1(u_2, v_2)(t)-\mathcal{T}_1(u_1, v_1)(t)\right| \nonumber \\ \le \frac{1}{|\eta^2-\Delta^2|}\Bigg[\eta\Bigg\{|\omega_1|\Bigg(\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)} \left|\theta_2(\tau, u_2(\tau), v_2(\tau))-\theta_2(\tau, u_1(\tau), v_1(\tau))\right|\Bigg)d\tau \nonumber \\ \quad +\left|h_2(s, u_2(s), v_2(s))-h_2(s, u_1(s), v_1(s))\right|ds\Bigg) \nonumber\\ + |\omega_2|\Bigg(\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\delta-1}}{\Gamma(\delta)}\left|\theta_2(\tau, u_2(\tau), v_2(\tau))-\theta_2(\tau, u_1(\tau), v_1(\tau))\right|\Bigg)d\tau \nonumber \\ \quad +\left|h_2(s, u_2(s), v_2(s))-h_2(s, u_1(s), v_1(s))\right|ds\Bigg) \nonumber\\ +\int_{0}^{\eta}\frac{(\eta-s)^{\gamma-1}}{\Gamma(\gamma)} \left|\theta_1(\tau, u_2(\tau), v_2(\tau))-\theta_1(\tau, u_1(\tau), v_1(\tau))\right|ds \nonumber \\ \quad +\left|h_1(\eta, u_2(\eta), v_2(\eta))-h_1(\eta, u_1(\eta), v_1(\eta))\right|\Bigg\} \nonumber\\$

$+ |\Delta|\Bigg\{|\omega_1|\Bigg(\int_{0}^{\xi_1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)} \left|\theta_1(\tau, u_2(\tau), v_2(\tau))-\theta_1(\tau, u_1(\tau), v_1(\tau))\right|\Bigg)d\tau \nonumber \\ \quad +\left|h_1(s, u_2(s), v_2(s))-h_1(s, u_1(s), v_1(s))\right|ds\Bigg) \nonumber\\ + |\omega_2|\Bigg(\int_{\xi_2}^{1} \Bigg(\int_{0}^{s}\frac{(s-\tau)^{\gamma-1}}{\Gamma(\gamma)} \left|\theta_1(\tau, u_2(\tau), v_2(\tau))-\theta_1(\tau, u_1(\tau), v_1(\tau))\right|\Bigg)d\tau \nonumber \\ \quad +\left|h_1(s, u_2(s), v_2(s))-h_1(s, u_1(s), v_1(s))\right|ds\Bigg) \nonumber\\ +\int_{0}^{\eta}\frac{(\eta-s)^{\delta-1}}{\Gamma(\delta)} \left|\theta_1(\tau, u_2(\tau), v_2(\tau))-\theta_1(\tau, u_2(\tau), v_2(\tau)\right|ds \nonumber \\ \quad +\left|h_2(\eta, u_2(\eta), v_2(\eta))-h_2(\eta, u_1(\eta), v_1(\eta))\right|ds\Bigg\}\Bigg]\nonumber\\ + \int_{0}^{t}\frac{(t-s)^{\gamma-1}}{\Gamma(\gamma)} \Bigg(\left|\theta_1(s, u_2(s), v_2(s))-\theta_1(s, u_1(s), v_1(s))\right|\Bigg)ds \nonumber \\ \quad +\left|h_1(t, u_2(t), v_2(t))-h_1(t, u_1(t), v_1(t))\right| \nonumber\\ \le \frac{1}{|\eta^2-\Delta^2|} \Bigg[\eta|\omega_1|\frac{\xi_1^{\delta+1}}{\Gamma(\delta+2)}+\eta|\omega_2|\frac{1-\eta^{\delta+1}}{\Gamma(\delta+2)}+|\Delta|\frac{\eta^{\delta}}{\Gamma(\delta+1)}\Bigg]\times \left(n_1||u_2-u_1||+n_2||v_2-v_1||\right) \nonumber\\ + \Bigg(\frac{1}{|\eta^2-\Delta^2|}\Bigg[\frac{\eta^{\gamma+1}}{\Gamma(\gamma+1)}+|\Delta| |\omega_1|\frac{\xi_1^{\gamma+1}}{\Gamma(\gamma+2)}+|\Delta||\omega_2|\frac{1-\xi_2^{\gamma+1}}{\Gamma(\gamma+2)}\Bigg]\nonumber \\ +\frac{1}{\Gamma(\gamma+1)}\Bigg)\times\left(m_1||u_2-u_1||+m_2||v_2-v_1||\right) \nonumber\\ + \frac{1}{|\eta^2-\Delta^2|}\Bigg[(\eta|\omega_1|\xi_1+\eta|\omega_2|(1-\xi_2)+|\Delta|)(\sigma_1||u_2-u_1||+\sigma_2||v_2-v_1||)\\ +(\eta+|\Delta||\omega_1|\xi_1+|\Delta||\omega_2|(1-\xi_2)+1)(\beta_1||u_2-u_1||+\beta_2||v_2-v_1||)\Bigg] \nonumber\\ le \mathcal{M}_2\left(n_1||u_2-u_1||+n_2||v_2-v_1||\right) + \mathcal{M}_1\left(m_1||u_2-u_1||+m_2||v_2-v_1||\right) \nonumber \\ \ + \mathcal{N}_5\left(\sigma_1||u_2-u_1||+\sigma_2||v_2-v_1||\right)+\mathcal{N}_6\left(\beta_1||u_2-u_1||+\beta_2||v_2-v_1||\right) \nonumber \\ = \left(\mathcal{M}_2n_1+\mathcal{M}_1m_1+\mathcal{N}_5\sigma_1+\mathcal{N}_6\beta_1\right)||u_2-u_1|| + \left(\mathcal{M}_2n_2+\mathcal{M}_1m_2+\mathcal{N}_5\sigma_2\right)||v_2-v_1||)$

which implies that

$\begin{array}{l} \left\|\mathcal{T}_1(u_2, v_2)(t)-\mathcal{T}_1(u_1, v_1)(t)\right\| \\ le \left(\mathcal{M}_2n_1+\mathcal{M}_1m_1+\mathcal{N}_5\sigma_1+\mathcal{N}_6\beta_1+\mathcal{M}_2n_2+\mathcal{M}_1m_2+\mathcal{N}_5\sigma_2+\mathcal{N}_6\beta_2\right)\left(||u_2-u_1||+||v_2-v_1||\right). \end{array}$

(3.3)

Likewise, we have

$\begin{array}{l} \left\|\mathcal{T}_2(u_2, v_2)(t)-\mathcal{T}_2(u_1, v_1)(t)\right\| \\ le \left(\mathcal{M}_4n_1+\mathcal{M}_3m_1+\mathcal{N}_6\sigma_1+\mathcal{N}_7\beta_1+\mathcal{M}_4n_2+\mathcal{M}_3m_2+\mathcal{N}_6\sigma_2+\mathcal{N}_7\beta_2\right)\left(||u_2-u_1||+||v_2-v_1||\right). \end{array}$

(3.4)

From (3.3) and (3.4), we deduce that

$\left\|\mathcal{T}(u_2, v_2)(t)-\mathcal{T}(u_1, v_1)(t)\right\| \nonumber \\ \le \left[(\mathcal{M}_1+\mathcal{M}_3)(m_1+m_2)+(\mathcal{M}_2+\mathcal{M}_4)(n_1+n_2)+(\mathcal{N}_7+\mathcal{N}_8)(\beta_1+\beta_2)+(\mathcal{N}_5+\mathcal{N}_6)(\sigma_1+\sigma_2)\right]\\ \times\left(||u_2-u_1||+||v_2-v_1||\right),$

which shows that $\mathcal{T}$ is a contraction by the assumption (A5) and hence it has a unique fixed point by Banach fixed point theorem. This leads to the conclusion that there exists a unique solution for the problems (1.1) and (1.2) on $[0, 1]$ . The proof is complete.

Now, we discuss the existence of solutions for the problems (1.1) and (1.2) by means of Leray-Schauder alternative (^[12], p. 4).

Theorem 3.2. Assume that there exists real constants $\tilde{k}_0 > 0$ , $\tilde{\lambda}_0 > 0$ and $\tilde{k}_i, \tilde{\lambda}_i\ge 0, \ i = 1, 2$ such that, for any $u_i\in \mathfrak{R}, \ i = 1, 2$

$\left|\theta_1(t, u_1, u_2)\right| \le \tilde{k}_0+\tilde{k}_1|u_1|+\tilde{k}_2|u_2|, \, \, \left|\theta_2(t, u_1, u_2)\right| \le \tilde{\lambda}_0+\tilde{\lambda}_1|u_1|+\tilde{\lambda}_2|u_2|.$

In addition,

$(\mathcal{M}_1+\mathcal{M}_3)\tilde{k}_1+(\mathcal{M}_2+\mathcal{M}_4)\tilde{\lambda}_1 \lt 1, \, \, (\mathcal{M}_1+\mathcal{M}_3)\tilde{k}_2+(\mathcal{M}_2+\mathcal{M}_4)\tilde{\lambda}_2 \lt 1,$

where $\mathcal{M}_i, \ i = 1, 2, 3, 4$ are given in (A4). Then the problems (1.1) and (1.2) have at least one solution on $[0, 1]$ .

Proof. The proof consists of two steps. First we show that the operator $\mathcal{T}:\mathcal{W}\times\mathcal{W}\to\mathcal{W}\times\mathcal{W}$ defined by (3.2) is completely continuous. Observe that continuity of the operator $\mathcal{T}$ follows from that of $\theta_1$ and $\theta_2$ . Consider a bounded set $\Omega \subset \mathcal{W}\times\mathcal{W}$ so that we can find positive constants $l_1$ and $l_2$ such that $\left|\theta_1(t, u(t), v(t))\right|\le l_1$ and $\left|\theta_2(t, u(t), v(t))\right|\le l_2$ for every $(u, v)\in \Omega$ . Hence, for any $(u, v)\in \Omega$ , we find that

$\left|\mathcal{T}_1(u, v)(t)\right| \le \frac{1}{|\eta^2-\Delta^2|} \Bigg[\eta|\omega_1|\frac{\xi_1^{\delta+1}}{\Gamma(\delta+2)}+\eta|\omega_2|\frac{1-\xi_2^{\delta+1}}{\Gamma(\delta+2)}+|\Delta|\frac{\eta^{\delta}}{\Gamma(\delta+1)}\Bigg]l_2 \nonumber\\ + \Bigg\{\frac{1}{|\eta^2-\Delta^2|}\Bigg[\frac{\eta^{\gamma}}{\Gamma(\gamma+1)}+|\Delta| |\omega_1|\frac{\xi_1^{\gamma+1}}{\Gamma(\gamma+2)}+|\Delta||\omega_2|\frac{1-\xi_2^{\gamma+1}}{\Gamma(\gamma+2)}+\frac{1}{\Gamma(\gamma+1)}\Bigg]\Bigg\}l_1 \nonumber\\ +\frac{1}{|\eta^2-\Delta^2|}\Big\{\eta\left[|\omega_1|\xi_1\mu_2+|\omega_2|\mu_2(1-\xi_2)+\mu_1\right]+\eta\mu_1+|\Delta|\left[\mu_1\xi_1|\omega_1|+\mu_1|\omega_2|(1-\xi_2)+\mu_2\right]\Big\}+\mu_1 \\ = \mathcal{M}_2l_2+\mathcal{M}_1l_1+\mathcal{N}_3.$

Thus we deduce that $\left\|\mathcal{T}_1(u, v)\right\|\le \mathcal{M}_2l_2+\mathcal{M}_1l_1+\mathcal{N}_3.$ In a similar fashion, it can be found that $\left\|\mathcal{T}_2(u, v)\right\|\le \mathcal{M}_4l_2+\mathcal{M}_3l_1+\mathcal{N}_4.$ Hence, it follows from the foregoing inequalities that $\mathcal{T}_1$ and $\mathcal{T}_2$ are uniformly bounded and hence the operator $\mathcal{T}$ is uniformly bounded. In order to show that $\mathcal{T}$ is equicontinuous, we take $0 < r_1 < r_2 < 1$ . Then, for any $(u, v)\in \Omega$ , we obtain

$\left|\mathcal{T}_1(u(r_2), v(r_2))-\mathcal{T}_1(u(r_1), v(r_1))\right| \nonumber \\ \le \frac{l_1}{\Gamma(\gamma)}\int_{0}^{r_1}\left[(r_2-s)^{\gamma-1}-(r_1-s)^{\gamma-1}\right]ds+\frac{l_1}{\Gamma(\gamma)}\int_{r_1}^{r_2}(r_2-s)^{\gamma-1}ds \nonumber \\ + \frac{r_2-r_1}{|\eta^2-\Delta^2|} \Bigg\{\Bigg[\eta|\omega_1|\frac{\xi_1^{\delta+1}}{\Gamma(\delta+2)}+\eta|\omega_2|\frac{1-\xi_2^{\delta+1}}{\Gamma(\delta+2)}+|\Delta|\frac{\eta^{\delta}}{\Gamma(\delta+1)}\Bigg]l_2 \nonumber\\ + \Bigg[\frac{\eta^{\gamma+1}}{\Gamma(\gamma+1)}+|\Delta| |\omega_1|\frac{\xi_1^{\gamma+1}}{\Gamma(\gamma+2)}+|\Delta||\omega_2|\frac{1-\xi_2^{\gamma+1}}{\Gamma(\gamma+2)}+\frac{1}{\Gamma(\gamma+2)}\Bigg]l_1+N_3\Bigg\},$

$\left|\mathcal{T}_2(u(r_2), v(r_2))-\mathcal{T}_2(u(r_1), v(r_1))\right| \nonumber \\ \le \frac{l_2}{\Gamma(\delta)}\int_{0}^{r_1}\left[(r_2-s)^{\delta-1}-(r_1-s)^{\delta-1}\right]ds+\frac{l_2}{\Gamma(\delta)}\int_{r_1}^{r_2}(r_2-s)^{\delta-1}ds \nonumber \\ \ \ + \frac{r_2-r_1}{|\eta^2-\Delta^2|} \Bigg\{\Bigg[\eta|\omega_1|\frac{\xi_1^{\delta+1}}{\Gamma(\delta+2)}+\eta|\omega_2|\frac{1-\xi_2^{\delta+1}}{\Gamma(\delta+2)}+\frac{\eta\dot{.}\eta^{\delta+1}}{\Gamma(\delta+1)}\Bigg]l_2 \nonumber\\ + \Bigg[|\Delta|\frac{\eta^{\gamma}}{\Gamma(\gamma+1)}+\eta |\omega_1|\frac{\xi_1^{\gamma+1}}{\Gamma(\gamma+2)}+\eta|\omega_2|\frac{1-\xi_2^{\gamma+1}}{\Gamma(\gamma+2)}\Bigg]l_1+\mathcal{N}_4\Bigg\},$

which imply that the operator $\mathcal{T}(u, v)$ is equicontinuous. In view of the foregoing arguments, we deduce that operator $\mathcal{T}(u, v)$ is completely continuous.

Next, we consider a set $\mathcal{P} = \{(u, v)\in \mathcal{W}\times\mathcal{W} : (u, v) = \lambda\mathcal{T}(u, v), \ 0\le \lambda\le 1\}$ and show that it is bounded. Let us take $(u, v)\in \mathcal{P}$ and $t\in [0, 1]$ . Then it follows from $u(t) = \lambda \mathcal{T}_1(u, v)(t)$ and $v(t) = \lambda \mathcal{T}_2(u, v)(t)$ , together with the given assumptions that

$\|u\|\le \mathcal{M}_1\left(\tilde{k}_0+\tilde{k}_1||u||+\tilde{k}_2||v||\right) +\mathcal{M}_2\left(\tilde{\lambda}_0+\tilde{\lambda}_1||u||+\tilde{\lambda}_2||v||\right)+\mathcal{N}_3, \\ |v\| \le \mathcal{M}_3\left(\tilde{k}_0+\tilde{k}_1||u||+\tilde{k}_2||v||\right)+\mathcal{M}_4\left(\tilde{\lambda}_0+\tilde{\lambda}_1||u||+\tilde{\lambda}_2||v||\right)+\mathcal{N}_4,$

which lead to

$\|u\|+\|v\| \le \left[(\mathcal{M}_1+\mathcal{M}_3)\tilde{k}_0+(\mathcal{M}_2+\mathcal{M}_4)\tilde{\lambda}_0 +\mathcal{N}_3+\mathcal{N}_4\right] \nonumber \\ + \left[(\mathcal{M}_1+\mathcal{M}_3)\tilde{k}_1+(\mathcal{M}_2+\mathcal{M}_4)\tilde{\lambda}_1\right]\|u\|+ \left[(\mathcal{M}_1+\mathcal{M}_3)\tilde{k}_2+(\mathcal{M}_2+\mathcal{M}_4)\tilde{\lambda}_2\right]\|v\|.$

Thus

$\|(u, v)\| \le \frac{(\mathcal{M}_1+\mathcal{M}_3)\tilde{k}_0+(\mathcal{M}_2+\mathcal{M}_4)\tilde{\lambda}_0+\mathcal{N}_3+\mathcal{N}_4}{\mathcal{M}_k},$

where $\mathcal{M}_{k}$ is defined by (3.1). Consequently the set $\mathcal {P}$ is bounded. Hence, it follows by Leray-Schauder alternative (^[12], p. 4) that the operator $\mathcal{T}$ has at least one fixed point. Therefore, the problems (1.1) and (1.2) have at least one solution on $[0, 1].$ This finishes the proof.

4. Examples

Example 4.1. Consider a coupled boundary value problem of fractional differential equations with slit-strips-type conditions given by

$\begin{array}{l} ^cD^{3/2}\left(u(t)-\frac{\sin t|u(t)|}{2(2+|u(t)|)}\right) = \frac{1}{56}u(t)+\frac{2}{7}\frac{v(t)} {1+v(t)}+\frac{5}{7}, \\ ^cD^{5/4}\left(v(t)-\frac{\sin t|v(t)|} {2(2+|v(t)|)}\right) = \frac{1}{39}\frac{|\cos u(t)|} {1+|\cos u(t)|}+\frac{1}{28}\sin v(t)+\frac{3}{7}, \end{array}$

(4.1)

$\begin{array}{l} u(0) = 0, \ u\left(\frac{1}{2}\right) = \int_{0}^{1/5}v(s)ds+\int_{4/5}^{1}v(s)ds, \\ v(0) = 0, \ v\left(\frac{1}{2}\right) = \int_{0}^{1/5}u(s)ds+\int_{4/5}^{1}u(s)ds. \end{array}$

(4.2)

Here $\gamma = \frac{3}{2}$ , $\delta = \frac{5}{4}$ , $\omega_1 = 1$ , $\omega_2 = 1$ , $\eta = \frac{1}{2}$ , $\xi_1 = \frac{1}{5}$ , $\xi_2 = \frac{4}{5}$ . From the given data, we find that $\Delta = -0.11$ , $m_1 = \frac{1}{56}$ , $m_2 = \frac{2}{71}$ , $n_1 = \frac{1}{39}$ , $n_2 = \frac{1}{28}$ , $\mathcal{M}_1\simeq 1.44716$ , $\mathcal{M}_2\simeq 0.51905$ , $\mathcal{M}_3\simeq 0.4046$ , $\mathcal{M}_4\simeq 2.51887$ , $\mathcal{N}_5\simeq 2.94238$ , $\mathcal{N}_6\simeq 7.3223$ , $\mathcal{N}_7\simeq 5.6164$ , $\mathcal{N}_8\simeq 5.2206$ , and $(\mathcal{M}_1+\mathcal{M}_3)(m_1+m_2)+(\mathcal{M}_2+\mathcal{M}_4)(n_1+n_2)+(\mathcal{N}_5+\mathcal{N}_6)(\sigma_1+\sigma_2)+(\mathcal{N}_7+\mathcal{N}_8)(\beta_1+\beta_2)\simeq 0.8030305 < 1$ .

Clearly all the conditions of Theorem 3.1 are satisfied. In consequence, the conclusion of Theorem 3.1 applies to the problems (4.1)–(4.2).

Example 4.2. We consider the problems (4.1)–(4.2) with

$\begin{array}{l} \theta_1(t, u(t), v(t))\ = \frac{1}{2}+\frac{2}{39}tanu(t)+\frac{2}{41}v(t), \, \, \theta_2(t, u(t), v(t))\ = \frac{2}{5}+\frac{1}{9} sinu(t)+\frac{1}{17}v(t). \end{array}$

(4.3)

Observe that

$\left|\theta_1(t, u, v)\right| \le \tilde{k}_0+\tilde{k}_1|u|+\tilde{k}_2|v|, \, \, \left|\theta_2(t, u, v)\right| \le \tilde{\lambda}_0+\tilde{\lambda}_1|u|+\tilde{\lambda}_2|v|$

with $\tilde{k}_0 = \frac{1}{2}$ , $\tilde{k}_1 = \frac{2}{39}$ , $\tilde{k}_2 = \frac{2}{41}$ $\tilde{\lambda}_0 = \frac{2}{5}$ , $\tilde{\lambda}_1 = \frac{1}{9}$ , $\tilde{\lambda}_2 = \frac{1}{17}$ . Furthermore,

$(\mathcal{M}_1+\mathcal{M}_3)\tilde{k}_1+(\mathcal{M}_2+\mathcal{M}_4)\tilde{\lambda}_1 \simeq 0.432507777 \lt 1, (\mathcal{M}_1+\mathcal{M}_3)\tilde{k}_2+(\mathcal{M}_2+\mathcal{M} _4)\tilde{\lambda}_2 \simeq 0.269030756 \lt 1.$

Thus all the conditions of Theorem 3.2 hold true and hence there exists at least one solution for the problems (4.1)–(4.2) with $\theta_1(t, u, v)$ and $\theta_2(t, u, v)$ given by (4.3).

Acknowledgment

The authors thank the reviewers for their useful remarks on our paper.

Conflict of interest

All authors declare no conflicts of interest in this paper.

References

[1]	A. L. Lloyd, V. A. A. Jansen, Spatiotemporal dynamics of epidemics: Synchrony in metapopulation models, Math. Biosci., 188 (2004), 1-16. doi: 10.1016/j.mbs.2003.09.003
[2]	M. Su, C. Hui, The effect of predation on the prevalence and aggregation of pathogens in prey, Biosystems, 105 (2011), 300-306. doi: 10.1016/j.biosystems.2011.05.012
[3]	C. Bret, Modeling and inference for infectious disease dynamics: A likelihood-based approach, Stat. Sci., 33 (2018), 57-69.
[4]	G. Q. Sun, J. H. Xie, S. H. Huang, Z. Jin, M.-T. Li, L. Liu, Transmission dynamics of cholera: Mathematical modeling and control strategies, Commun. Nonlin. Sci. Numer. Si., 45 (2017), 235-244. doi: 10.1016/j.cnsns.2016.10.007
[5]	Y. Kazuo, Threshold dynamics of reaction-diffusion partial differential equations model of Ebola virus disease, Int. J. Biomath., 11 (2018), 1850108. doi: 10.1142/S1793524518501085
[6]	J. R. Petrella, W. R. Hao, A. Rao, P. M. Doraiswamy, Computational causal modeling of the dynamic biomarker cascade in alzheimer's disease, Comput. Math. Methods Med., 2019 (2019), 6216530.
[7]	Q. Yang, H. F. Huo, Dynamics of an edge-based SEIR model for sexually transmitted diseases, Math. Biosci. Eng., 17 (2020), 669-699. doi: 10.3934/mbe.2020035
[8]	D. M. Xiao, S. G. Ruan, Global analysis of an epidemic model with nonmonotone incidence rate, Math. Biosci., 208 (2007), 419-429. doi: 10.1016/j.mbs.2006.09.025
[9]	M. Haque, J. Zhen, E. Venturino, An ecoepidemiological predator-prey model with standard disease incidence, Math. Method Appl. Sci., 32 (2010), 875-898.
[10]	L. I. Bo, S. L. Yuan, W. G. Zhang, Analysis on an epidemic model with a ratio-dependent nonlinear incidence rate, Int. J. Biomath., 4 (2011), 227-239. doi: 10.1142/S1793524511001374
[11]	Y. L. Cai, W. M. Wang, Spatiotemporal dynamics of a reaction-diffusion epidemic model with nonlinear incidence rate, J. Stat. Mech. Theory. E, 2011(2011), P02025.
[12]	X. B. Zhang, S. Q. Chang, H. Huo, Dynamic behavior of a stochastic SIR epidemic model with vertical transmission, Electron. J. Differ. Equat., 2019(2019), 1-20. doi: 10.1186/s13662-018-1939-6
[13]	D. Wodarz, D. N. Levy, N. L. Komarova, Multiple infection of cells changes the dynamics of basic viral evolutionary processes, Evolution, 3 (2019), 104-115.
[14]	D. X. Jia, T. H. Zhang, S. L. Yuan, Pattern dynamics of a diffusive toxin producing phytoplanktonzooplankton model with three-dimensional patch, Int. J. Bifurcat. Chaos, 29 (2019), 1930011. doi: 10.1142/S0218127419300118
[15]	M. Bendahmane, M. Langlais, A reaction-diffusion system with cross-diffusion modeling the spread of an epidemic disease, J. Evol. Equ., 10 (2010), 883-904. doi: 10.1007/s00028-010-0074-y
[16]	W. Z. Huang, M. A. Han, K. Y. Liu, Dynamics of an SIS reaction-diffusion epidemic model for disease transmission, Math. Biosci. Eng., 7 (2017), 51-66.
[17]	S. N. Zhao, S. L. Yuan, H. Wang, Threshold behavior in a stochastic algal growth model with stoichiometric constraints and seasonal variation, J. Differ. Equat., 268 (2020), 5113-5139. doi: 10.1016/j.jde.2019.11.004
[18]	A. D'Onofrio, Stability properties of pulse vaccination strategy in SEIR epidemic model, Math. Biosci., 179(2002), 57-72.
[19]	M. DelaSen, S. Alonso-Quesada, R. Nistal, A. Ibeas, On the existence of equilibrium points, boundedness, oscillating behavior and positivity of a sveirs epidemic model under constant and impulsive vaccination, Advances Differ. Equat., 2011 (2011), 339-370.
[20]	A. L. Shane, R. K. Mody, J. A. Crump, L. K. Pickering, 2017 infectious diseases society of America clinical practice guidelines for the diagnosis and management of infectious diarrhea, Clin. Infect. Dis., 65(2017), 1963-1973. doi: 10.1093/cid/cix959
[21]	C. Kuttler, Modeling to inform infectious disease control, Biometrics, 74 (2018), 382-383. doi: 10.1111/biom.12852
[22]	S. K. Nandi, S. Jana, M. Mandal, K. K. Tapan, Complex dynamics and optimal treatment of an epidemic model with two infectious diseases, Biosystems, 5 (2019), 5-29.
[23]	C. Y. Zheng, C. Y. Xia, Q. T. Guo, M. Dehmer, Interplay between SIR-based disease spreading and awareness diffusion on multiplex networks, J. Parallel Distr. Com., 115 (2018), 20-28. doi: 10.1016/j.jpdc.2018.01.001
[24]	L. J. S. Allen, B. M. Bolker, Y. Lou, A. L. Nevai, Asymptotic profiles of the steady states for an SIS epidemic patch model, SIAM. J. Appl. Math., 67 (2007), 1283-1309. doi: 10.1137/060672522
[25]	R. H. Cui, Y. Lou, A spatial SIS model in advective heterogeneous environments, J. Differ. Equat., 261 (2016), 3305-3343. doi: 10.1016/j.jde.2016.05.025
[26]	G. Q. Sun, J. Zhang, L. P. Song, Z. Jin, B. L. Li, Pattern formation of a spatial predator-prey system, Appl. Math. Comput., 218 (2012), 11151-11162.
[27]	G. Q. Sun, L. Li, Z. K. Zhang, Spatial dynamics of a vegetation model in an arid flat environment, Nonlin. Dynam., 73 (2013), 2207-2219. doi: 10.1007/s11071-013-0935-3
[28]	D. Kalajdzievska, M. Y. Li, Modeling the effects of carriers on transmission dynamics of infectious diseases, Math. Biosci. Eng., 8 (2017), 711-722.
[29]	F. Berezovskaya, G. Karev, B. J. Song, C. Castillo-Chvez, A simple epidemic model with surprising dynamics, Math. Biosci. Eng., 2 (2005), 133-152.
[30]	X. W. Yu, S. L. Yuan, Asymptotic properties of a stochastic chemostat model with two distributed delays and nonlinear perturbation, Discrete Cont. Dyn. Syst. B, 25 (2020), 2373-2390.
[31]	X. Y. Meng, H. F. Huo, X. B. Zhang, Stability and global Hopf bifurcation in a Leslie-Gower predator-prey model with stage structure for prey, J. Appl. Math. Comput., 60 (2019), 1-25. doi: 10.1007/s12190-018-1201-0
[32]	X. Y. Meng, J. Li, Stability and Hopf bifurcation analysis of a delayed phytoplankton-zooplankton model with Allee effect and linear harvesting, Math. Biosci. Eng., 17 (2020), 1973-2002. doi: 10.3934/mbe.2020105
[33]	X. Y. Meng, Y. Q. Wu, Dynamical analysis of a fuzzy phytoplankton-zooplankton model with refuge, fishery protection and harvesting, J. Appl. Math. Comput., 63 (2020), 361-389.
[34]	W. M. Wang, Y. L. Cai, M. J. Wu, K. F. Wang, Z. Q. Li, Complex dynamics of a reaction-diffusion epidemic model, Nonlin. Anal. Real Word Appl., 13 (2012), 2240-2258. doi: 10.1016/j.nonrwa.2012.01.018
[35]	T. Wang, Dynamics of an epidemic model with spatial diffusion, Phys. A, 409 (2014), 119-129. doi: 10.1016/j.physa.2014.04.028
[36]	W. M. Wang, H. Y. Liu, Y. L. Cai, Z. Q. Li, Turing pattern selection in a reaction-diffusion epidemic model, Chinese Phys. B, 20 (2011), 074702. doi: 10.1016/j.physa.2014.04.028
[37]	M. Fras, M. Gosak, Spatiotemporal patterns provoked by environmental variability in a predatorprey model, Biosystems, 114 (2013), 172-177. doi: 10.1016/j.biosystems.2013.09.004
[38]	Z. P. Ma, H. F. Huo, H. Xiang, Spatiotemporal patterns induced by delay and cross fractional diffusion in a predator-prey model describing intraguild predation, Math. Methods Appl. Sci., 43 (2020), 5179-5196. doi: 10.1002/mma.6259
[39]	A. Onofrio, M. Banerjee, P. Manfredi, Spatial behavioural responses to the spread of an infectious disease can suppress Turing and Turing-Hopf patterning of the disease, Phys. A, 545 (2020), 123773. doi: 10.1016/j.physa.2019.123773
[40]	Y. F. Li, J. A. Cui, The effect of constant and pulse vaccination on SIS epidemic models incorporating media coverage, Commun. Nonlin. Sci. Numer. Simul., 14 (2008), 2353-2365.
[41]	M. J. Ma, S. Y. Liu, J. Li, Does media coverage influence the spread of drug addiction, Commun. Nonlinear Sci. Numer. Simul., 50(2017), 169-179. doi: 10.1016/j.cnsns.2017.03.002
[42]	M. A. Khan, S. Islam, G. Zaman, Media coverage campaign in hepatitis b transmission model, Appl. Math. Comput., 331 (2018), 378-393.
[43]	R. S. Liu, J. H. Wu, H. P. Zhu, Media/psychological impact on multiple outbreaks of emerging infectious diseases, Comput. Math. Methods Med., 8 (2007), 153-164. doi: 10.1080/17486700701425870
[44]	J. A. Cui, X. Tao, H. P. Zhu, An SIS infection model incorporating media coverage, Rocky Mountauin J. Math., 38 (2008), 1323-1334.
[45]	Y. P. Liu, J. A. Cui, The impact of media coverage on the dynamics of infections disease, Int. J. Biomath, 1 (2008), 65-74. doi: 10.1142/S1793524508000023
[46]	Y. F. Li, J. A. Cui, The effect of constant and pulse vaccination on SIS epidemic models incorporating media coverage, Commun. Nonlin. Sci. Numer. Simul., 14 (2009), 2353-2365. doi: 10.1016/j.cnsns.2008.06.024
[47]	J. A. Cui, Y. H. Sun, H. P. Zhu, The impact of media on the control of infectious diseases, J. Dyn. Differ. Equ., 20 (2008), 31-53. doi: 10.1007/s10884-007-9075-0
[48]	Y. N. Xiao, S. Y. Tang, J. H. Wu, Media impact switching surface during an infectious disease outbreak, Sci. Rep., 5 (2015), 7838. doi: 10.1038/srep07838
[49]	Q. L. Yan, S. Y. Tang, S. Gabriele, J. H. Wu, Media coverage and hospital notifications: Correlation analysis and optimal media impact duration to manage a pandemic, J. Theor. Biol., 390 (2016), 1-13. doi: 10.1016/j.jtbi.2015.11.002
[50]	J. M. Tchuenche, N. Dube, B. P. Claver, J. S. Robert, T. B. Chris, The impact of media coverage on the transmission dynamics of human influenza, BMC Public Health, 11(2011), S5.
[51]	C. J. Sun, W. Yang, J. Arino, K. Kamran, Effect of media-induced social distancing on disease transmission in a two patch setting, Math. Biosci., 230 (2011), 87-95. doi: 10.1016/j.mbs.2011.01.005
[52]	P. V. D. Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29-48. doi: 10.1364/JOSA.48.000784
[53]	X. Y. Meng, J. G. Wang, Analysis of a delayed diffusive model with Beddington-Deangelis functional response, Int. J. Biomath., 12 (2019), 1950047. doi: 10.1142/S1793524519500475
[54]	Y. Zhao, J. Li, X. Ma, Stochastic periodic solution of a susceptible-infective epidemic model in a polluted environment under environmental fluctuation, Comput. Math. Methods Med., 2018 (2018), 1-15.

This article has been cited by:

1.	Wafa Shammakh, Hadeel Z. Alzumi, Zahra Albarqi, On multi-term proportional fractional differential equations and inclusions, 2020, 2020, 1687-1847, 10.1186/s13662-020-03104-y
2.	Karthikeyan Buvaneswari, Panjaiyan Karthikeyan, Dumitru Baleanu, On a system of fractional coupled hybrid Hadamard differential equations with terminal conditions, 2020, 2020, 1687-1847, 10.1186/s13662-020-02790-y
3.	Zhiwei Lv, Ishfaq Ahmad, Jiafa Xu, Akbar Zada, Analysis of a Hybrid Coupled System of ψ-Caputo Fractional Derivatives with Generalized Slit-Strips-Type Integral Boundary Conditions and Impulses, 2022, 6, 2504-3110, 618, 10.3390/fractalfract6100618
4.	Said Mesloub, Eman Alhazzani, Hassan Eltayeb Gadain, A Two-Dimensional Nonlocal Fractional Parabolic Initial Boundary Value Problem, 2024, 13, 2075-1680, 646, 10.3390/axioms13090646
5.	Haroon Niaz Ali Khan, Akbar Zada, Ioan-Lucian Popa, Sana Ben Moussa, The Impulsive Coupled Langevin ψ-Caputo Fractional Problem with Slit-Strip-Generalized-Type Boundary Conditions, 2023, 7, 2504-3110, 837, 10.3390/fractalfract7120837
6.	Pengyan Yu, Guoxi Ni, Chengmin Hou, Existence and uniqueness for a mixed fractional differential system with slit-strips conditions, 2024, 2024, 1687-2770, 10.1186/s13661-024-01942-3
7.	Karthikeyan Buvaneswari, Panjaiyan Karthikeyan, Kulandhivel Karthikeyan, Ozgur Ege, Existence and uniqueness results on coupled Caputo-Hadamard fractional differential equations in a bounded domain, 2024, 38, 0354-5180, 1489, 10.2298/FIL2404489B
8.	Haroon Niaz Ali Khan, Akbar Zada, Ishfaq Khan, Analysis of a Coupled System of $\psi$ -Caputo Fractional Derivatives with Multipoint–Multistrip Integral Type Boundary Conditions, 2024, 23, 1575-5460, 10.1007/s12346-024-00987-0

Reader Comments

Your name:*

Email:*
© 2020 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

Mathematical Biosciences and Engineering

4.4

Metrics

Article views(4501) PDF downloads(276) Cited by(5)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Figures(4)

Mathematical Biosciences and Engineering

The impact of media on the spatiotemporal pattern dynamics of a reaction-diffusion epidemic model

Related Papers:

Abstract

1. Introduction

2. An auxiliary lemma

3. Main results

4. Examples

Acknowledgment

Conflict of interest

References

This article has been cited by:

Reader Comments

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

Metrics

Figures and Tables

Other Articles By Authors

Catalog

Mathematical Biosciences and Engineering

The impact of media on the spatiotemporal pattern dynamics of a reaction-diffusion epidemic model

Related Papers:

Abstract

1. Introduction

2. An auxiliary lemma

3. Main results

4. Examples

Acknowledgment

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