A mathematical model for policy of vaccinating recovered people in controlling the spread of COVID-19 outbreak

1.   Introduction

Fractional calculus is an emerging field drawing attention from both theoretical and applied disciplines. In particular, fractional calculus is a powerful tool for explaining problems in ecology, biology, chemistry, physics, mechanics, networks, flow in porous media, electricity, control systems, viscoelasticity, mathematical biology, fitting of experimental data, and so forth. One may see the papers ^[1,2,3,4,5] and the references therein.

Fractional difference calculus or discrete fractional calculus is a very new field for mathematicians. Some real-world phenomena are being studied with the assistance of fractional difference operators. Basic definitions and properties of fractional difference calculus can be found in the book ^[6]. Fractional boundary value problems can be found in the books ^[7,8]. Now, the studies of boundary value problems for fractional difference equations are extended to be more complex. Excellent papers related to discrete fractional boundary value problems can be found in ^{[9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35]} and references cited therein. In particular, there are some recent papers that present the Caputo fractional difference calculus ^{[36,37,38,39,40,41]}. In the literature, there are apparently few research works studying boundary value problems for Caputo fractional difference-sum equations. For example, ^[42] studied a boundary value problem for $p-$Laplacian Caputo fractional difference equations with fractional sum boundary conditions of the forms

In ^[43], investigated a nonlocal fractional sum boundary value problem for a Caputo fractional difference-sum equation of the form

In addition, ^[44] considered a periodic boundary value problem for Caputo fractional difference-sum equations of the form

We aim to fill the gaps related to the boundary value problem of Caputo fractional difference-sum equations. The goal of this paper is to enrich this new research area by using the unknown function of Caputo fractional difference and fractional sum in the problem. So, in this paper, we consider a sequential nonlinear Caputo fractional sum-difference equation with fractional difference boundary value conditions of the form

where $\; \rho_1, \rho_2\in {\mathbb{R}},$ $0 < \alpha, \beta, \gamma, \nu, \mu\leq 1$, $1 < \alpha+\beta\leq 2\;$ are given constants, $\; \mathcal{H}\in C\big({\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta} \times {\mathbb{R}}^3, \mathbb{R}\big)$, and for $\; \varphi \in C\big(\odot\times {\mathbb{R}}^2, [0, \infty)\big)$, $\odot: = \left\lbrace (t, r)\, :\, t, r\in {\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}\rm{\; and }\; r\leq t\right\rbrace$. The operator $\Psi^\mu$ is defined by

The plan of this paper is as follows. In Section 2, we recall some definitions and basic lemmas. Also, we derive the solution of (1.4) by converting the problem to an equivalent equation. In Section 3, we prove existence and uniqueness results of the problem (1.4) using the Banach contraction principle and Schaefer's theorem. Furthermore, we also show the existence of a positive solution to (1.4). An illustrative example is presented in Section 4.

2.   Preliminaries

In the following, there are notations, definitions and lemmas which are used in the main results.

Definition 2.1. ^[10] We define the generalized falling function by $t^{\underline{\alpha}}: = \frac{\Gamma(t+1)}{\Gamma(t+1-\alpha)}$, for any $t$ and $\alpha$ for which the right-hand side is defined. If $t+1-\alpha$ is a pole of the Gamma function and $t+1$ is not a pole, then $t^{\underline{\alpha}} = 0$.

Lemma 2.1. ^[9] Assume the following factorial functions are well defined. If $t\leq r$, then $t^{\underline{\alpha}}\leq r^{\underline{\alpha}}$ for any $\alpha > 0$.

Definition 2.2. ^[10] For $\alpha > 0$ and $f$ defined on $\mathbb{N}_a: = \{a, a+1, \ldots\}$, the $\alpha$-order fractional sum of $f$ is defined by

where $t\in\mathbb{N}_{a+\alpha}$ and $\sigma(s) = s+1$.

Definition 2.3. ^[11] For $\alpha > 0$ and $f$ defined on $\mathbb{N}_a$, the $\alpha$-order Caputo fractional difference of $f$ is defined by

where $t\in \mathbb{N}_{a+N-\alpha}$ and $N \in\mathbb{N}$ is chosen so that $0\leq N-1 < \alpha < N$. If $\alpha = N$, then $\Delta^\alpha_C f(t) = \Delta^N f(t)$.

Lemma 2.2. ^[11] Assume that $\alpha > 0$ and $0\leq N-1 < \alpha\leq N.$ Then,

for some $C_i\in\mathbb{R}$, $0\leq i\leq N-1$.

To study the solution of the boundary value problem (1.4), we need the following lemma that deals with a linear variant of the boundary value problem (1.4) and gives a representation of the solution.

Lemma 2.3. Let $\; \Lambda\left(\rho_1-1\right) \neq 0$, $0 < \alpha, \beta, \gamma, \nu, \mu \leq 1$, $1 < \alpha+\beta\leq 2$ and $h\in C\left({\mathbb{N}}_{\alpha+\beta-1, T+\alpha+\beta-1}, \mathbb{R}\right)$ be given. Then, the problem

has the unique solution

where

Proof. Using the fractional sum of order $\alpha \in (0, 1]$ for $(2.1)$ and from Lemma 2.2, we obtain

for $t\in {\mathbb{N}}_{\alpha-1, T+\alpha}$.

Using the fractional sum of order $0 < \beta\leq1$ for $(2.5)$, we obtain

for $t\in{\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}$.

By substituting $t = \alpha+\beta-2, T+\alpha+\beta$ into $(2.6)$ and employing the first condition of $(2.2)$, we obtain

Using the fractional Caputo difference of order $0 < \gamma\leq1$ for $(2.6)$, we obtain

for ${\mathbb{N}}_{\alpha+\beta-\gamma-1, T+\alpha+\beta-\gamma+1}$.

By substituting $t = \alpha+\beta-\gamma-1, T+\alpha+\beta-\gamma+1$ into $(2.8)$ and employing the second condition of $(2.2)$, it implies

The constant $C_2$ can be obtained by substituting $C_1$ into $(2.7)$. Then, we get

where $\Lambda$ is defined by $(2.4)$. Substituting the constants $C_1$ and $C_2$ into $(2.6)$, we obtain $(2.3)$.

3.   Main results

In this section, we wish to establish the existence results for the problem $(1.4)$. We denote $\; {\mathcal{C}} = C({\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}, \mathbb{R})$ as the Banach space of all functions $x$ with the norm defined by

where $\|x\| = \max\limits_{t\in {\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}} |x(t)|$ and $\|\Delta^{\nu}_Cx\| = \max\limits_{t\in {\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}}\big|\Delta^{\nu}_{C}x(t-\nu+1)\big|$.

The following assumptions are assumed:

$\textbf{(A1)}$ ${\mathcal{H}}[t, x, y, z]:{\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta} \times {\mathbb{R}}^3\rightarrow \mathbb{R}$ is a continuous function.

$\textbf{(A2)}$ There exist constants $K_1, K_2 > 0$ such that for each $t\in {\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}$ and all $x_i, y_i, z_i\in {\mathbb{R}}, \; i = 1, 2$, we have

and

where $\; \Psi^\mu (t, 0): = \frac{1}{\Gamma(\mu)}\sum\limits_{s = \alpha+\beta-\mu-2}^{t-\mu}(t-\sigma(s))^{\underline{\mu-1}}\, \varphi\Big(t, s+\mu, 0, 0\Big)$.

$\textbf{(A3)}$ $\varphi:\odot \times {\mathbb{R}}^2\rightarrow {\mathbb{R}}$ is continuious for $(t, s)\in \odot$, and there exists a constant $L > 0$, such that for each $(t, s)\in \odot$ and all $x_i, y_i\in \mathcal{C}, \; i = 1, 2$ we have

Let us define the operator ${\widetilde{\mathcal{H}}}[t, x(t)]$ by

for each $t\in {\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}\;$ and $\; x\in {\mathcal{C}}$.

Note that $\; \Delta^{-\beta} \Delta^{-\alpha}\, {\widetilde{\mathcal{H}}}[t, x(t)]$ and $\Delta^{\nu}_C\, \Delta^{-\beta} \Delta^{-\alpha} \, {\widetilde{\mathcal{H}}}[t, x(t)]\;$ exist when $\; \nu < \alpha+\beta\leq2.$

Lemma 3.1. Assume that (A1)–(A3) hold. Then, the following property holds:

(A4) There exits a positive constant $\Theta$ such that

for each $t\in {\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}\;$ and $\; x_1, x_2\in {\mathcal{C}}$, where

Proof. By $\textbf{(A3)}$, for each $t\in {\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}\;$ and $\; x_1, x_2\in {\mathcal{C}}$, we obtain

and hence

Next, we define the operator ${\mathcal{F}}:\mathcal{C}\longrightarrow\mathcal{C}$ by

where

By Lemma 2.3, we find that any solution of the problem $(1.4)$ is the fixed point of the operator ${\mathcal{F}}$.

Lemma 3.2. Assume that the function ${\mathcal{A}}_{\alpha, \beta, \gamma, \rho_1, \rho_2}(t, s, r, \xi)$ satisfies the following properties:

(A5) ${\mathcal{A}}_{\alpha, \beta, \gamma, \rho_1, \rho_2}(t, s, r, \xi)$ is a continuous function for all $(t, s, r, \xi)\in {\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}\times{\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}\times{\mathbb{N}}_{\alpha-1, T+\alpha+1}\times {\mathbb{N}}_{0, T+2} = :\mathcal{D}$, and there exist two constants $\; \Omega_1, \Omega_2 > 0$, such that

where

and $\; _t\Delta_C^{\nu}\;$ is the Caputo fractional difference with respect to $t$.

Proof. It is obvious that ${\mathcal{A}}_{\alpha, \beta, \gamma, \rho_1, \rho_2}(t, s, r, \xi)$ is a continuous function for all $(t, s, r, \xi)\in {\mathcal{D}}$. Next, we consider

and

Thus, the condition $\textbf{(A5)}$ holds.

In what follows, we consider the existence and uniqueness of a solution to the problem (1.4) using the Banach contraction principle.

Theorem 3.1. Assume that (A1)–(A5) hold. If

where $\Omega_1, \Omega_2$ are defined as (3.5)–(3.6), and

then, the problem $(1.4)$ has a unique solution in ${\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}.$

Proof. Choose a constant $R$ satisfying

We will show that ${\mathcal{F}}(B_R)\subset B_R,$ where $B_R = \{x \in \mathcal{C}: \|x\|_{\mathcal{C}} \leq R\}$. For all $x\in B_R,$ we have

and

Thus,

and hence, ${\mathcal{F}}(B_R)\subset B_R$.

We next show that $\mathcal{F}$ is a contraction. For all $x, y\in \mathcal{C}$ and for each $t\in \mathbb{N}_{\alpha+\beta-2, T+\alpha+\beta}$, we have

and

Thus,

Therefore, ${\mathcal{F}}$ is a contraction. Hence, by using Banach fixed point theorem, we get that ${\mathcal{F}}$ has a fixed point which is a unique solution of the problem (1.4).

We next deduce the existence of a solution to (1.4) by using the following Schaefer's fixed point theorem.

Theorem 3.2. ^[45] (Arzelá-Ascoli Theorem) A set of functions in $C[a, b]$ with the sup norm is relatively compact if and only if it is uniformly bounded and equicontinuous on $[a, b]$.

Theorem 3.3. ^[45] If a set is closed and relatively compact, then it is compact.

Theorem 3.4. ^[46] Let $X$ be a Banach space and $T: X\rightarrow X$ be a continuous and compact mapping. If the set

is bounded, then $T$ has a fixed point.

Theorem 3.5. Suppose that (A1)–(A5) hold. Then, the problem (1.4) has at least one solution on $\mathbb{N}_{\alpha+\beta-2, T+\alpha+\beta}$.

Proof. We shall use Schaefer's fixed point theorem to prove that the operator $F$ defined by $(3.3)$ has a fixed point. It is clear that ${\mathcal{F}}:\mathcal{C}\longrightarrow\mathcal{C}$ is completely continuous. So, it remains to show that the set

Let $u\in E$. Then,

Thus, for each $t\in {\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}$, we have

and

Hence,

This shows that $E$ is bounded. By Schaefer's fixed point theorem, we conclude that the problem (1.4) has at least one solution.

In the sequel, we discuss the positivity of the obtained solution $x\in\mathcal{C}$. To this end, we add adequate assumptions and provide the following theorem.

We note that a positive solution of $(1.4)$ in $\mathcal{C}$ is a function $x(t) > 0$ which has $\Delta^\nu_C\, x(t-\nu+1) > 0$ for all $t\in {\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}$.

Theorem 3.6. Suppose that (A1)–(A5) are fulfilled in ${\mathbb{R}^+}$, where $\; \mathcal{H}\in C\big({\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta} \times {\mathbb{R}}^+\times {\mathbb{R}}^+\times {\mathbb{R}}^+, \mathbb{R}^{+}\big)\;$ and $\; \varphi \in C\big(\odot\times {\mathbb{R}^{+}}\times {\mathbb{R}^{+}}, {\mathbb{R}}^{+}\big)$. If condition $(3.7)$ is satisfied, for $\; \alpha, \beta, \gamma, \nu, \mu \in (0, 1),$ and in addition

then a solution in $\mathcal{C}$ of the problem $(1.4)$ is positive.

Proof. By Theorem $3.1$ and the fact that, for $\; \rho_1 > 1 \; \; \rm{and}\; \; \rho_2 > \frac{\Gamma(T+\gamma+4)}{\Gamma(2-\gamma)\Gamma(T+3)}$, the condition $(3.7)$ is a particular case, the problem $(1.4)$ admits a unique solution in $\mathcal{C}$.

Moreover, since $\; \alpha, \beta, \gamma, \nu, \mu \in (0, 1)$, we obtain for each $(t, s, r, \xi)\in {\mathcal{D}}$,

and

It results that the unique solution $x(t)$ of problem $(1.4)$ which satisfies with (3.3) is positive for each $t\in {\mathbb{N}}_{\alpha+\beta-2, T+\alpha+\beta}$.

4.   An example

In this section, we present an example to illustrate our results.

Example. Consider the following fractional difference boundary value problem:

By letting $\; \alpha = \frac{2}{3}, \; \beta = \frac{5}{6}, \; \gamma = \frac{1}{3}, \; \nu = \frac{1}{2}, \; \mu = \frac{1}{4}, \; T = 4, \; \rho_1 = 2, \; \rho_2 = 20$, ${\mathcal{H}}[t, x, y, z] = \frac{1}{(t+5)^5}\left[\frac{x}{1+|x|} + \frac{y}{1+|y|} + z \right] {\rm{\; and }}\; \varphi[t+\frac{1}{2}, s+\frac{1}{4}, x, y] = \frac{e^{-s}}{(t+5)^2}\left[\frac{x+1}{3+|x|} + \frac{y+1}{3+|y|} \right],$ we can show that

Observe that (A1)–(A5) hold for all $x_i, y_i, z_i\in \mathbb{R}, \; i = 1, 2$, and for each $t\in \mathbb{N}_{\frac{-1}{2}, \frac{11}{2}}$, we obtain

So, $K_1 = \left(\frac{2}{11}\right)^5\approx 0.000199, \;$ and $\; K_2 = \max\limits_{t\in \mathbb{N}_{\frac{-1}{2}, \frac{11}{2}}}{\widetilde{\mathcal{H}}[t, 0]}\approx0.0000394.$

Next, for all $x_i, y_i\in \mathbb{R}, \; i = 1, 2,$ and each $(t, s)\in \mathbb{N}_{\frac{-1}{2}, \frac{11}{2}}\times \mathbb{N}_{\frac{-1}{2}, \frac{11}{2}}$, we obtain

So, $K_2 = e^{-\frac{1}{2}}\left(\frac{2}{11}\right)^5\approx 0.000121.$

Finally, we can show that

Hence, by Theorem 3.1, the problem (4.1) has a unique solution.

Moreover, by Theorem 3.5, the problem (4.1) has at least one solution on $\mathbb{N}_{\frac{-1}{2}, \frac{11}{2}}$.

Furthermore, ${\mathcal{H}}, \varphi \in {\mathbb{R}}^+,$ and

Therefore, the solution of the problem (4.1) is positive on $\mathbb{N}_{\frac{-1}{2}, \frac{11}{2}}$ by Theorem 3.6.

5.   Conclusions

In the present research, we considered a sequential nonlinear Caputo fractional sum-difference equation with fractional difference boundary conditions. Notice that the unknown function of this problem is in the form of Caputo fractional difference and fractional sum with different orders, which expands the research scope of the problems in ^[42,43,44]. Existence results are established by a Banach contraction principle and Schaefer's fixed point theorem. {The results of the paper are new and enrich the subject of boundary value problems for Caputo fractional difference-sum equations. In future work, we may extend this work by considering new boundary value problems.

Acknowledgments

This research was funded by the National Science, Research and Innovation Fund (NSRF) and Suan Dusit University with Contract no. 64-FF-06.

Conflict of interest

The authors declare no conflict of interest.