On a SEIR-type model of COVID-19 using piecewise and stochastic differential operators undertaking management strategies

Mdi Begum Jeelani; Kamal Shah; Hussam Alrabaiah; Abeer S. Alnahdi; Mdi Begum Jeelani; Kamal Shah; Hussam Alrabaiah; Abeer S. Alnahdi

doi:10.3934/math.20231395

AIMS Mathematics

2023, Volume 8, Issue 11: 27268-27290. doi: 10.3934/math.20231395

Previous Article Next Article

Research article Special Issues

On a SEIR-type model of COVID-19 using piecewise and stochastic differential operators undertaking management strategies

1.
Department of Mathematics and Statistics, College of Science, Imam Mohammad Ibn Saud Islamic University, Riyadh, Saudi Arabia
2.
Department of Mathematics, University of Malakand, Chakdara Dir (Lower), Khyber Pakhtunkhawa, Pakistan
3.
Department of Computer Science and Mathematics, Lebanese American University, Byblos Lebanon
4.
College of Engineering, Al Ain University, Al Ain, UAE
5.
Department of Mathematics, Tafila Technical University, Tafila, Jordan

Received: 03 August 2023 Revised: 23 August 2023 Accepted: 28 August 2023 Published: 26 September 2023
MSC : 03C65, 26A33, 34A08

In this work, an epidemic model of a susceptible, exposed, infected and recovered SEIR-type is established for the distinctive dynamic compartments and epidemic characteristics of COVID-19 as it spreads across a population with a heterogeneous rate. The proposed model is investigated using a novel approach of fractional calculus known as piecewise derivatives. The existence theory is demonstrated through the establishment of sufficient conditions. In addition, result related to Hyers-Ulam stability is also derived for the considered model. A numerical method based on modified Euler procedure is also constructed to simulate the approximate solutions of the proposed model by employing various values of fractional orders. We testified the numerical results by using real available data of Japan. In addition, some results for the SEIR-type model are also presented graphically using the stochastic process, and the obtained results are discussed.

Keywords:

Citation: Mdi Begum Jeelani, Kamal Shah, Hussam Alrabaiah, Abeer S. Alnahdi. On a SEIR-type model of COVID-19 using piecewise and stochastic differential operators undertaking management strategies[J]. AIMS Mathematics, 2023, 8(11): 27268-27290. doi: 10.3934/math.20231395

Related Papers:

[1]	Qi Xiao, Jin Zhong . Characterizations and properties of hyper-dual Moore-Penrose generalized inverse. AIMS Mathematics, 2024, 9(12): 35125-35150. doi: 10.3934/math.20241670
[2]	Waleed Mohamed Abd-Elhameed, Amr Kamel Amin, Nasr Anwer Zeyada . Some new identities of a type of generalized numbers involving four parameters. AIMS Mathematics, 2022, 7(7): 12962-12980. doi: 10.3934/math.2022718
[3]	Faik Babadağ . A new approach to Jacobsthal, Jacobsthal-Lucas numbers and dual vectors. AIMS Mathematics, 2023, 8(8): 18596-18606. doi: 10.3934/math.2023946
[4]	Changsheng Luo, Jiagui Luo . Complete solutions of the simultaneous Pell equations $(a^2+1)y^2-x^2 = y^2-bz^2 = 1$ . AIMS Mathematics, 2021, 6(9): 9919-9938. doi: 10.3934/math.2021577
[5]	Cencen Dou, Jiagui Luo . Complete solutions of the simultaneous Pell's equations $(a^2+2)x^2-y^2 = 2$ and $x^2-bz^2 = 1$ . AIMS Mathematics, 2023, 8(8): 19353-19373. doi: 10.3934/math.2023987
[6]	Faik Babadağ, Ali Atasoy . A new approach to Leonardo number sequences with the dual vector and dual angle representation. AIMS Mathematics, 2024, 9(6): 14062-14074. doi: 10.3934/math.2024684
[7]	Elahe Mehraban, T. Aaron Gulliver, Salah Mahmoud Boulaaras, Kamyar Hosseini, Evren Hincal . New sequences from the generalized Pell $p-$ numbers and mersenne numbers and their application in cryptography. AIMS Mathematics, 2024, 9(5): 13537-13552. doi: 10.3934/math.2024660
[8]	Bin Zhou, Xiujuan Ma, Fuxiang Ma, Shujie Gao . Robustness analysis of random hyper-networks based on the internal structure of hyper-edges. AIMS Mathematics, 2023, 8(2): 4814-4829. doi: 10.3934/math.2023239
[9]	Waleed Mohamed Abd-Elhameed, Anna Napoli . New formulas of convolved Pell polynomials. AIMS Mathematics, 2024, 9(1): 565-593. doi: 10.3934/math.2024030
[10]	Ümit Tokeşer, Tuğba Mert, Yakup Dündar . Some properties and Vajda theorems of split dual Fibonacci and split dual Lucas octonions. AIMS Mathematics, 2022, 7(5): 8645-8653. doi: 10.3934/math.2022483

Abstract

1. Introduction

Dual numbers were first given by Clifford (1845–1879), and some properties of those were studied in the geometrical investigation, and Kotelnikov ^[1] introduced their first applications. Study applied to line geometry and kinematics dual numbers and dual vectors ^[2]. He demonstrated that the directed lines of Euclidean 3-space and the points of the dual unit sphere in $\mathbb{D}^3$ have a one-to-one relationship. Field theory also relies heavily on these numbers ^[3]. The most intriguing applications of dual numbers in field theory are found in a number of Wald publications ^[4]. Dual numbers have contemporary applications in kinematics, dynamics, computer modeling of rigid bodies, mechanism design, and kinematics ^[5,6,7].

Complex numbers have significant advantages in derivative computations. However, the second derivative computations lost these advantages ^[8]. J. A. Fike developed the hyper-dual numbers to solve this issue ^[9]. These numbers may be used to calculate both the first and second derivatives while maintaining the benefits of the first derivative using complex numbers. Furthermore, it is demonstrated that this numerical approach is appropriate for open kinematic chain robot manipulators, sophisticated software, and airspace system analysis and design ^[10].

In the literature, sequences of integers have an important place. The most famous of these sequences have been demonstrated in several areas of mathematics. These sequences have been researched extensively because of their complex characteristics and deep connections to several fields of mathematics. The Fibonacci and Lucas sequences and their related numbers are of essential importance due to their various applications in biology, physics, statistics, and computer science ^[11,12,13]. Many authors were interested in introducing and investigating several generalizations and modifications of Fibonacci and Lucas sequences. The authors investigated two classes that generalize Fibonacci and Lucas sequences, and they utilized them to compute some radicals in reduced forms. Panwar ^[14] defined the generalized $k$ -Fibonacci sequence as

$\begin{align*} F_{k, n} = pkF_{k, n-1}+qF_{k, n-2}, \end{align*}$

with initial conditions $F_{k, 0} = a$ and $F_{k, 1} = b$ . If $a = 0, k = 2, p = q = b = 1$ , the classic Pell sequence and for $a = b = 2, k = 2, p = q = 1$ , Pell-Lucas sequences appear.

The Pell numbers are the numbers of the following integer sequence:

$0, 1, 2, 5, 12, 29, 70, 169, 408, 985, 2378, ...$

The sequence of Pell numbers, which is denoted by ${P}_n$ is defined as the linear reccurence relation

$\begin{align*} P_n = 2P_{n-1}+P_{n-2} , \quad P_0 = 0, P_1 = 1, \ \ n\ge2. \end{align*}$

The integer sequence of Pell-Lucas numbers denoted by ${Q}_n$ is given by

$2, 2, 6, 14, 34, 82, 198, 478, 1154, 2786, 6726, ...$ ,

with the same reccurence relation

$\begin{align*} Q_n = 2Q_{n-1}+Q_{n-2} , \quad Q_0 = Q_1 = 2, \ \ n\ge2. \end{align*}$

The characteristic equation of these numbers is $x^2-2x-1 = 0$ , with roots $\alpha = 1+\sqrt{2}$ and $\beta = 1-\sqrt{2}$ and the Binet's forms of these sequences are given as^{[15,16,17,18]},

$\begin{equation} P_n = \frac{\alpha^n-\beta^n}{\alpha-\beta} \end{equation}$

(1.1)

and

$\begin{equation} Q_n = \alpha^n+\beta^n. \end{equation}$

(1.2)

The set of dual numbers is defined as

$\begin{equation*} \mathbb{D} = \{d = a+\varepsilon a^{*}\mid a, a^{*}\in \mathbb{R}, \varepsilon ^2 = 0, \varepsilon \neq 0\}. \end{equation*}$

The set of hyper-dual numbers is

$\begin{align*} \widetilde{\mathbb{D}} = \left \{\gamma = \gamma_0+\gamma_1\varepsilon +\gamma_2\varepsilon^* +\gamma_3\varepsilon\varepsilon^* \mid \gamma_0, \gamma_1, \gamma_2, \gamma_3\in \mathbb{R}\right\}, \end{align*}$

or can be rewritten as

$\begin{align*} \widetilde{\mathbb{D}} = \left \{\gamma = d+\varepsilon^*d^* \mid d, d^*\in \mathbb{D}\right\}, \end{align*}$

where $\varepsilon$ , $\varepsilon^*$ and $\varepsilon\varepsilon^*$ are hyper-dual units that satisfy

$\begin{equation*} (\varepsilon)^2 = ({\varepsilon^*})^2 = 0, \; \varepsilon\neq\varepsilon^*\neq0, \; \varepsilon\varepsilon^* = \varepsilon^*\varepsilon. \end{equation*}$

This set forms commutative and associative algebra over both the dual and real numbers ^[8,9,10].

The square root of a hyper-dual number $\gamma$ can be defined by

$\begin{equation} \sqrt{\gamma} = \sqrt{\gamma_{0}}+\frac{\gamma_{1}}{2\sqrt{\gamma_{0}}}\varepsilon+\frac{\gamma_{2}}{2\sqrt{\gamma_{0}}}\varepsilon^*+(\frac{\gamma_{3}}{2\sqrt{\gamma_{0}}}-\frac{\gamma_{1}\gamma_{2}}{4\gamma_{0}\sqrt{\gamma_{0}}})\varepsilon\varepsilon^*. \end{equation}$

(1.3)

A hyper-dual vector is any vector of the form

$\begin{align*} \vec{\gamma} = \vec{\gamma}_{0}+\vec{\gamma}_{1}\varepsilon+\vec{\gamma}_{2}\varepsilon^*+\vec{\gamma}_{3}\varepsilon\varepsilon^*, \end{align*}$

where $\vec{\gamma}_{0}, \vec{\gamma}_{1}, \vec{\gamma}_{2}, \vec{\gamma}_{3}$ are real vectors, this vector can be rewritten as $\vec{\gamma} = \vec{d}+\varepsilon^*\vec{d}^*$ , where $\vec{d}$ and $\vec{d}^*$ are dual vectors. Let $\vec{\gamma}$ and $\vec{\delta}$ be hyper-dual vectors, then their scalar product is defined as

$\begin{align} \langle \vec{\gamma}, \vec{\delta}\rangle_{HD} = &\langle\vec{\gamma}_{0}, \vec{\delta}_{0}\rangle +\left(\langle \vec{\gamma}_{0}, \vec{\delta}_{1}\rangle+\langle \vec{\gamma}_{1}, \vec{\delta}_{0}\rangle\right)\varepsilon+\left(\langle \vec{\gamma}_{0}, \vec{\delta}_{2}\rangle+\langle \vec{\gamma}_{2}, \vec{\delta}_{0}\rangle\right)\varepsilon^*\\[6pt] &+\left(\langle \vec{\gamma}_{0}, \vec{\delta}_{3}\rangle+\langle \vec{\gamma}_{1}, \vec{\delta}_{2}\rangle+\langle \vec{\gamma}_{2}, \vec{\delta}_{1}\rangle+\langle \vec{\gamma}_{3}, \vec{\delta}_{0}\rangle\right) \varepsilon\varepsilon^*, \end{align}$

(1.4)

which continents inner products of real vectors.

Let $f(x_0+ x_1\varepsilon+ x_2\varepsilon^*+ x_3\varepsilon\varepsilon^*)$ be a hyper-dual function, then

$\begin{align} f(x_0+ x_1\varepsilon+ x_2\varepsilon^*+ x_3\varepsilon\varepsilon^*) = f(x_0)+x_1f(x_0)\varepsilon+x_2f'(x_0)\varepsilon^*+(x_3f'(x_0)+x_1x_2f''(x_0))\varepsilon\varepsilon^*. \end{align}$

(1.5)

Suppose $\vec{\gamma}$ , $\vec{\delta}$ and $\Phi$ be unit hyper-dual vectors and hyper-dual angle respectively then by using (1.5) the scalar product can be written as

$\begin{align} \langle\vec{\gamma}, \vec{\delta}\rangle_{HD} & = \cos{\Phi}\\ & = \cos\phi-\varepsilon^* \phi^{\ast }\sin\phi\\ & = (\cos\psi -\varepsilon \psi^{\ast }\sin\psi)-\varepsilon^* \phi^{\ast }(\sin\psi+\varepsilon \psi^{\ast }\cos\psi), \end{align}$

(1.6)

where $\phi$ and $\psi$ are, respectively, dual and real angles.

The norm of a hyper-dual vector $\vec{\gamma}$ is given by

$\begin{align*} \left\| \vec{\gamma}\right\|_{HD} & = \left\| \vec{\gamma}_0\right\|+\frac{\langle \vec{\gamma}_0, \vec{\gamma}_1\rangle }{\left\| \vec{\gamma}_0\right\|}\varepsilon+\frac{\langle \vec{\gamma}_0, \vec{\gamma}_2\rangle }{\left\| \vec{\gamma}_0\right\|} \varepsilon^* +\left(\frac{\langle \vec{\gamma}_0, \vec{\gamma}_3\rangle }{\left\| \vec{\gamma}_0\right\|}+\frac{\langle \vec{\gamma}_1, \vec{\gamma}_2\rangle }{\left\| \vec{\gamma}_0\right\|}-\frac{\langle \vec{\gamma}_0, \vec{\gamma}_1\rangle \langle \vec{\gamma}_0, \vec{\gamma}_2\rangle}{\left\| \vec{\gamma}_0\right\|^3}\right)\varepsilon\varepsilon^*, \end{align*}$

for $\left\| \vec{\gamma}_0\right\|\neq0$ . If $\left\| \vec{\gamma}\right\|_{HD} = 1$ that is $\left\| \vec{\gamma}_0\right\| = 1$ and $\langle \vec{\gamma}_0, \vec{\gamma}_1\rangle = \langle \vec{\gamma}_0, \vec{\gamma}_2\rangle = \langle \vec{\gamma}_0, \vec{\gamma}_3\rangle = \langle \vec{\gamma}_1, \vec{\gamma}_2\rangle = 0$ , then $\vec{\gamma}$ is a unit hyper-dual vector.

In this paper, we introduce the hyper-dual Pell and the hyper-dual Pell-Lucas numbers, which provide a natural generalization of the classical Pell and Pell-Lucas numbers by using the concept of hyper-dual numbers. We investigate some basic properties of these numbers. We also define a new vector and angle, which are called hyper-dual Pell vector and angle. We give properties of these vectors and angles to exert in the geometry of hyper-dual space.

2. Hyper-dual Pell and Hyper-dual Pell-Lucas numbers

In this section, we define the hyper-dual Pell and hyper-dual Pell-Lucas numbers and then demonstrate their fundamental identities and properties.

Definition 2.1. The $n^{th}$ hyper-dual Pell ${\mathrm{HP}}_n$ and hyper-dual Pell-Lucas ${\mathrm{HQ}}_n$ numbers are defined respectively as

$\begin{align} {\mathrm{HP}}_n = \mathrm{P}_n+\mathrm{P}_{n+1}\varepsilon+\mathrm{P}_{n+2}\varepsilon^*+\mathrm{P}_{n+3}\varepsilon\varepsilon^* \end{align}$

(2.1)

and

$\begin{align} \mathrm{HQ}_n = \mathrm{Q}_n+\varepsilon\mathrm{Q}_{n+1}+\varepsilon^*\mathrm{Q}_{n+2}+\varepsilon\varepsilon^*\mathrm{Q}_{n+3}, \end{align}$

(2.2)

where $\mathrm{P}_n$ and $\mathrm{Q}_n$ are $n^{th}$ Pell and Pell-Lucas numbers.

The few hyper-dual Pell and hyper-dual Pell-Lucas numbers are given as

${\mathrm{HP}}_1 = 1+2\varepsilon+5\varepsilon^*+12\varepsilon\varepsilon^*, {\mathrm{HP}}_2 = 2+5\varepsilon+12\varepsilon^*+29\varepsilon\varepsilon^*, ...$

and

${\mathrm{HQ}}_1 = 2+6\varepsilon+14\varepsilon^*+34\varepsilon\varepsilon^*, {\mathrm{HQ}}_2 = 6+14\varepsilon+34\varepsilon^*+82\varepsilon\varepsilon^*, ...$

Theorem 2.1. The Binet-like formulas of the hyper-dual Pell and hyper-dual Pell-Lucas numbers are given, respectively, by

$\begin{align} {\mathrm{HP}}_n = \frac{\varphi^n\;\underline{\varphi}-\psi^n\;\underline{\psi}}{\varphi-\psi} \end{align}$

(2.3)

and

$\begin{align} {\mathrm{HQ}}_n = \varphi^n\;\underline{\varphi}+\psi^n\;\underline{\psi}, \end{align}$

(2.4)

where

$\begin{align} \underline{\varphi} = 1+\varphi\varepsilon+\varphi^{2}\varepsilon^*+\varphi^{3}\varepsilon\varepsilon^*, \; \underline{\psi} = 1+\psi\varepsilon+\psi^{2}\varepsilon^*+\psi^{3}\varepsilon\varepsilon^*. \end{align}$

(2.5)

Proof. From (2.1) and the Binet formula of Pell numbers, we obtain

$\begin{align*} {\mathrm{HP}}_n& = \mathrm{P}_n+\mathrm{P}_{n+1}\varepsilon+\mathrm{P}_{n+2}\varepsilon^*+\mathrm{P}_{n+3}\varepsilon\varepsilon^*\\[6pt] & = \frac{\varphi^n-\psi^n}{\varphi-\psi}+\frac{\varphi^{n+1}-\psi^{n+1}}{\varphi-\psi}\varepsilon+\frac{\varphi^{n+2}-\psi^{n+2}}{\varphi-\psi}\varepsilon^*+\frac{\varphi^{n+3}-\psi^{n+3}}{\varphi-\psi}\varepsilon\varepsilon^*\\[6pt] & = \frac{\varphi^n(1+\varphi\varepsilon+\varphi^{2}\varepsilon^*+\varphi^{3}\varepsilon\varepsilon^*)}{\varphi-\psi}-\frac{\psi^n(1+\psi\varepsilon+\psi^{2}\varepsilon^*+\psi^{3}\varepsilon\varepsilon^*)}{\varphi-\psi} \\[6pt] & = \frac{\varphi^n\; \underline{\varphi}-\psi^n\; \underline{\psi}}{\varphi-\psi}. \end{align*}$

On the other hand, using (2.2) and the Binet formula of Pell-Lucas numbers we obtain

$\begin{align*} {\mathrm{HQ}}_n& = \mathrm{Q}_n+\mathrm{Q}_{n+1}\varepsilon+\mathrm{Q}_{n+2}\varepsilon^*+\mathrm{Q}_{n+3}\varepsilon\varepsilon^*\\[6pt] & = \left( \varphi^n+\psi^n \right)+\left( \varphi^{n+1}+\psi^{n+1} \right)\varepsilon+\left( \varphi^{n+2}+\psi^{n+2} \right)\varepsilon^*+\left( \varphi^{n+3}+\psi^{n+3} \right)\varepsilon\varepsilon^*\\[6pt] & = \varphi^n(1+\varphi\varepsilon+\varphi^{2}\varepsilon^*+\varphi^{3}\varepsilon\varepsilon^*)+\psi^n(1+\psi\varepsilon+\psi^{2}\varepsilon^*+\psi^{3}\varepsilon\varepsilon^*) \\[6pt] & = \varphi^n\; \underline{\varphi}+\psi^n\; \underline{\psi}. \end{align*}$

□ The proof is completed.

Theorem 2.2. (Vajda-like identities) For non-negative integers m, n, and r, we have

$\begin{align*} &{\mathrm{HP}}_{m}{\mathrm{HP}}_{n}-{\mathrm{HP}}_{m-r}{\mathrm{HP}}_{n+r} = (-1)^{n+1}\mathrm{P}_{m-n-r}\mathrm{P}_{r}\left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right), \\[6pt] &{\mathrm{HQ}}_{m}{\mathrm{HQ}}_{n}-{\mathrm{HQ}}_{m-r}{\mathrm{HQ}}_{n+r} = (-1)^{n}\; \mathrm{Q}_{m-n}-(-1)^{n+r}\; \mathrm{Q}_{m-n-2r}\; \left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right). \end{align*}$

Proof. By using the Binet-like formula of hyper-dual Pell numbers, we obtain

$\begin{align*} {\mathrm{HP}}_{m}{\mathrm{HP}}_{n}-{\mathrm{HP}}_{m-r}{\mathrm{HP}}_{n+r}& = \left(\frac{\; \varphi^{m}\underline{\varphi}-\; \psi^{m}\underline{\psi}}{\varphi-\psi}\right)\left(\frac{\; \varphi^{n}\underline{\varphi}-\; \psi^{n}\underline{\psi}}{\varphi-\psi}\right) -\left(\frac{\; \varphi^{m-r}\underline{\varphi}-\; \psi^{m-r}\underline{\psi}}{\varphi-\psi}\right)\left(\frac{\; \varphi^{n+r}\underline{\varphi}-\; \psi^{n+r}\underline{\psi}}{\varphi-\psi}\right)\\[6pt] & = \frac{(\varphi^{r}-\psi^{r})(\varphi^{n}\psi^{m-r}-\psi^{n}\varphi^{m-r})}{(\varphi-\psi)^2}\underline{\varphi}\; \underline{\psi}\\[6pt] & = -\frac{(\varphi^{m-n-r}-\psi^{m-n-r})(\varphi^{r}-\psi^{r})}{(\varphi-\psi)^2}\underline{\varphi}\; \underline{\psi}, \end{align*}$

and by using (1.1), we obtain

$\begin{equation*} {\mathrm{HP}}_{m}{\mathrm{HP}}_{n}-{\mathrm{HP}}_{m-r}{\mathrm{HP}}_{n+r} = (-1)^{n+1}\mathrm{P}_{m-n-r}\mathrm{P}_{r}\left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right). \end{equation*}$

Similarly for hyper-dual Pell-Lucas numbers, we can obtain

$\begin{align*} {\mathrm{HQ}}_{m}{\mathrm{HQ}}_{n}-{\mathrm{HQ}}_{m-r}{\mathrm{HQ}}_{n+r}& = \left(\; \varphi^{m}\underline{\varphi}+\; \psi^{m}\underline{\psi}\right)\left(\; \varphi^{n}\underline{\varphi}+\; \psi^{n}\underline{\psi}\right)\\&\quad-\left(\; \varphi^{m-r}\underline{\varphi}+\; \psi^{m-r}\underline{\psi}\right)\left(\; \varphi^{n+r}\underline{\varphi}+\; \psi^{n+r}\underline{\psi}\right) \\& = \underline{\varphi}\; \underline{\psi}\left(\varphi^{m-n}+\psi^{m-n}-\varphi^{m-n-2r}-\psi^{m-n-2r}\right). \end{align*}$

Using (1.2) and (2.5),

$\begin{align*} {\mathrm{HQ}}_{m}{\mathrm{HQ}}_{n}-{\mathrm{HQ}}_{m-r}{\mathrm{HQ}}_{n+r} = (-1)^{n}\; \mathrm{Q}_{m-n}-(-1)^{n+r}\; \mathrm{Q}_{m-n-2r}\; \left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right). \end{align*}$

Thus, we obtain the desired results. □

Theorem 2.3. (Catalan-like identities) For non negative integers n and r, with $n\ge r$ , we have

$\begin{align*} &{\mathrm{HP}}_{n-r}{\mathrm{HP}}_{n+r}-{{\mathrm{HP}}}_{n}^2 = (-1)^{n-r}\mathrm{P}_{r}^{2}\left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right), \\[6pt] &{\mathrm{HQ}}_{n-r}{\mathrm{HQ}}_{n+r}-{\mathrm{HQ}}_{n}^2 = 8(-1)^{n-r}\; \mathrm{P}_{r}^{2}\; \left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right). \end{align*}$

Proof. From (2.3), we obtain

$\begin{align*} {\mathrm{HP}}_{n-r}{\mathrm{HP}}_{n+r}-{\mathrm{HP}}_{n}^2& = \left(\frac{\; \varphi^{n-r}\underline{\varphi}-\; \psi^{n-r}\underline{\psi}}{\varphi-\psi}\right)\left(\frac{\; \varphi^{n+r}\underline{\varphi}-\; \psi^{n+r}\underline{\psi}}{\varphi-\psi}\right)-\left(\frac{\; \varphi^n\underline{\varphi}-\; \psi^n\underline{\psi}}{\varphi-\psi}\right)^2\\[6pt] & = \frac{\varphi^{n}\psi^{n}}{8}\; \underline{\varphi}\; \underline{\psi}\left(2-\psi^{r}\varphi^{-r}-\psi^{-r}\varphi^{r}\right)\\[6pt] & = (-1)^{n-r}\; \underline{\varphi}\; \underline{\psi}\left(\frac{\varphi^{r}-\psi^{r}}{\varphi-\psi}\right)^{2}, \end{align*}$

and by using (1.1) and (2.5), we will have

$\begin{equation*} {\mathrm{HP}}_{n-r}{\mathrm{HP}}_{n+r}-{\mathrm{HP}}_{n}^2 = (-1)^{n-r}\mathrm{P}_{r}^{2}\left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right). \end{equation*}$

On the other hand, from (2.4) and (2.5) we obtain

$\begin{align*} {\mathrm{HQ}}_{n-r}{\mathrm{HQ}}_{n+r}-{\mathrm{HQ}}_{n}^2& = \left(\; \varphi^{n-r}\underline{\varphi}+\; \psi^{n-r}\underline{\psi}\right)\left(\; \varphi^{n+r}\underline{\varphi}+\; \psi^{n+r}\underline{\psi}\right)-\left(\; \varphi^n\underline{\varphi}+\; \psi^n\underline{\psi}\right)^2 \\& = \underline{\varphi}\; \underline{\psi}\left(\varphi^{n-r}\psi^{n+r}+\varphi^{n+r}\psi^{n-r}-2\psi^{n}\varphi^{n}\right)\\ & = 8(-1)^{n-r}\; \underline{\varphi}\; \underline{\psi}\left(\frac{\varphi^{r}-\psi^{r}}{\varphi-\psi}\right)^{2}\\ & = 8(-1)^{n-r}\; \mathrm{P}_{r}^{2}\; \left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right). \end{align*}$

□

Corollary 2.1. (Cassini-like identities) For non-negative integer n, we have

$\begin{align*} &\mathrm{HP}_{n-1}\; \mathrm{HP}_{n+1}-\mathrm{HP}_{n}^2 = (-1)^{n-1}\left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right), \\[6pt] &\mathrm{HQ}_{n-1}\; \mathrm{HQ}_{n+1}- \mathrm{HQ}_{n}^2 = 8(-1)^{n-1}\left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right). \end{align*}$

Proof. We can get the result by taking $r = 1$ in Theorem 2.3. □

Theorem 2.4. (d'Ocagne-like identities) For non-negative integers n and m,

$\begin{align*} &\mathrm{HP}_{m+1}\; \mathrm{HP}_{n}-\mathrm{HP}_{m}\; \mathrm{HP}_{n+1} = \left( -1 \right)^{m}\mathrm{P}_{n-m}\left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right), \\[6pt] &\mathrm{HQ}_{m+1}\; \mathrm{HQ}_{n}-\mathrm{HQ}_{m}\; \mathrm{HQ}_{n+1} = 8\left( -1 \right)^{n}\mathrm{P}_{m-n}\left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right). \end{align*}$

Proof. Using (1.1), (2.3), and (2.5), we have

$\begin{align*} \mathrm{HP}_{m+1}\; \mathrm{HP}_{n}-\mathrm{HP}_{m}\; \mathrm{HP}_{n+1}& = \left(\frac{\; \varphi^{m+1}\underline{\varphi}-\; \psi^{m+1}\underline{\psi}}{\varphi-\psi}\right)\left(\frac{\; \varphi^{n}\underline{\varphi}-\; \psi^{n}\underline{\psi}}{\varphi-\psi}\right)\\ &\quad-\left(\frac{\; \varphi^{m}\underline{\varphi}-\; \psi^{m}\underline{\psi}}{\varphi-\psi}\right)\left(\frac{\; \varphi^{n+1}\underline{\varphi}-\; \psi^{n+1}\underline{\psi}}{\varphi-\psi}\right)\\ & = (\varphi-\psi)(\varphi^{n}\psi^{m}-\varphi^{m}\psi^{n}) \underline{\varphi}\; \underline{\psi} \\ & = \left( -1 \right)^{m}\mathrm{P}_{n-m}\left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right). \end{align*}$

Using (1.2), (2.4) and (2.5), we have

$\begin{align*} \mathrm{HQ}_{m+1}\; \mathrm{HQ}_{n}-\mathrm{HQ}_{m}\; \mathrm{HQ}_{n+1} = 8\left( -1 \right)^{n}\mathrm{P}_{m-n}\left( 1+2\varepsilon+6\varepsilon^*+12\varepsilon\varepsilon^* \right). \end{align*}$

□

3. Hyper dual Pell vectors and angle

In this section, we introduce hyper-dual Pell vectors and hyper-dual Pell angle. We will give geometric properties of them.

Definition 3.1. The $n^{th}$ hyper-dual Pell vector is defined as

$\begin{align*} \overrightarrow{{\mathrm{HP}}}_n = \vec{\mathrm{P}}_n+\vec{\mathrm{P}}_{n+1}\varepsilon+\vec{\mathrm{P}}_{n+2}\varepsilon^*+\vec{\mathrm{P}}_{n+3}\varepsilon\varepsilon^*, \end{align*}$

where $\vec{\mathrm{P}}_{n} = (\mathrm{P}_n, \mathrm{P}_{n+1}, \mathrm{P}_{n+2})$ is a real Pell vector. The hyper-dual Pell vector $\overrightarrow{\mathrm{HP}}_n$ can be rewritten in terms of dual Pell vectors $\vec{\mathbb{P}}_n$ and $\vec{\mathbb{P}} ^*_{n}$ as

$\begin{align*} \nonumber \overrightarrow{\mathrm{HP}}_n & = (\vec{\mathrm{P}}_n+\vec{\mathrm{P}}_{n+1}\varepsilon)+(\vec{\mathrm{P}}_{n+2}+\vec{\mathrm{P}}_{n+3}\varepsilon)\varepsilon^* \\ & = \vec{\mathbb{P}}_n+\varepsilon^*\vec{\mathbb{P}} ^*_{n}. \end{align*}$

Theorem 3.1. The scalar product of hyper-dual Pell vectors $\overrightarrow{{\mathrm{HP}}}_n$ and $\overrightarrow{{\mathrm{HP}}}_m$ is

$\begin{align} \langle \overrightarrow{\mathrm{HP}}_n, \overrightarrow{\mathrm{HP}}_m \rangle = &\frac{7\mathrm{Q}_{n+m+2}}{8}-\frac{(-1)^{m}\mathrm{Q}_{n-m}}{8} +(\frac{7\mathrm{Q}_{n+m+3}}{4}-\frac{(-1)^{m}\mathrm{Q}_{n-m}}{4})\varepsilon \\ &+(\frac{7\mathrm{Q}_{n+m+4}}{4}-\frac{3(-1)^{m}\mathrm{Q}_{n-m}}{4})\varepsilon^* +(\frac{7\mathrm{Q}_{n+m+5}}{2}-\frac{3(-1)^{m}\mathrm{Q}_{n-m}}{2})\varepsilon\varepsilon^*. \end{align}$

(3.1)

Proof. By using (1.4), we can write

$\begin{align} \langle \overrightarrow{\mathrm{HP}}_n, \overrightarrow{\mathrm{HP}}_{m}\rangle = &\langle\vec{\mathrm{P}}_{n}, \vec{\mathrm{P}}_{m}\rangle +\left(\langle \vec{\mathrm{P}}_{n}, \vec{\mathrm{P}}_{m+1}\rangle+\langle \vec{\mathrm{P}}_{n+1}, \vec{\mathrm{P}}_{m}\rangle\right)\varepsilon+\left(\langle\vec{\mathrm{P}}_{n}, \vec{\mathrm{P}}_{m+2}\rangle+\langle \vec{\mathrm{P}}_{n+2}, \vec{\mathrm{P}}_{m}\rangle\right)\varepsilon^*\\[6pt] &+\left(\langle \vec{\mathrm{P}}_{n}, \vec{\mathrm{P}}_{m+3}\rangle+\langle \vec{\mathrm{P}}_{n+1}, \vec{\mathrm{P}}_{m+2}\rangle+\langle \vec{\mathrm{P}}_{n+2}, \vec{\mathrm{P}}_{m+1}\rangle+\langle \vec{\mathrm{P}}_{n+3}, \vec{\mathrm{P}}_{m}\rangle\right) \varepsilon\varepsilon^*. \end{align}$

(3.2)

Now we calculate the above inner products for real Pell vectors $\vec{\mathrm{P}}_{n}$ and $\vec{\mathrm{P}}_{m}$ by using Binet's formula of Pell numbers as

$\begin{align*} \langle \vec{\mathrm{P}}_n, \vec{\mathrm{P}}_{m}\rangle& = \mathrm{P}_{n}\mathrm{P}_{m}+\mathrm{P}_{n+1}\mathrm{P}_{m+1}+ \mathrm{P}_{n+2}\mathrm{P}_{m+2} \\ & = \left(\frac{\varphi^{n}-\psi^{n}}{\varphi-\psi}\right)\left(\frac{\varphi^{m}-\psi^{m}}{\varphi-\psi}\right)+\left(\frac{\varphi^{n+1}-\psi^{n+1}}{\varphi-\psi}\right)\left(\frac{\varphi^{m+1}-\psi^{m+1}}{\varphi-\psi}\right) +\left(\frac{\varphi^{n+2}-\psi^{n+2}}{\varphi-\psi}\right)\left(\frac{\varphi^{m+2}-\psi^{m+2}}{\varphi-\psi}\right)\\ & = \frac{\varphi^{n+m}+\psi^{n+m}}{(\varphi-\psi)^2}+\frac{\varphi^{n+m+2}+\psi^{n+m+2}}{(\varphi-\psi)^2} +\frac{\varphi^{n+m+4}+\psi^{n+m+4}}{(\varphi-\psi)^2}-\frac{\left( \varphi^{n}\psi^{m}+\varphi^{m}\psi^{n} \right)\varphi^{-m}\psi^{-m}}{(\varphi-\psi)^2\; \varphi^{-m}\psi^{-m}} \\ & = \frac{1}{8}\left( \mathrm{Q}_{n+m}+\mathrm{Q}_{n+m+2}+\mathrm{Q}_{n+m+4}+(-1)^{m}\mathrm{Q}_{n-m} \right)\\ & = \frac{7\mathrm{Q}_{n+m+2}}{8}-\frac{(-1)^{m}\mathrm{Q}_{n-m}}{8}. \end{align*}$

□ Similarly,

$\begin{align*} \langle \vec{\mathrm{P}}_n, \vec{\mathrm{P}}_{m+1} \rangle& = \frac{7\mathrm{Q}_{n+m+3}}{8}+\frac{(-1)^{m}\mathrm{Q}_{n-m-1}}{8}, \\ \langle \vec{\mathrm{P}}_{n+1}, \vec{\mathrm{P}}_{m} \rangle& = \frac{7\mathrm{Q}_{n+m+3}}{8}-\frac{(-1)^{m}\mathrm{Q}_{n-m+1}}{8}, \\ \langle \vec{\mathrm{P}}_{n}, \vec{\mathrm{P}}_{m+2} \rangle& = \frac{7\mathrm{Q}_{n+m+4}}{8}-\frac{(-1)^{m}\mathrm{Q}_{n-m-2}}{8}, \\ \langle \vec{\mathrm{P}}_{n+2}, \vec{\mathrm{P}}_{m} \rangle& = \frac{7\mathrm{Q}_{n+m+4}}{8}-\frac{(-1)^{m}\mathrm{Q}_{n-m+2}}{8}, \\ \langle \vec{\mathrm{P}}_{n}, \vec{\mathrm{P}}_{m+3} \rangle& = \frac{7\mathrm{Q}_{n+m+5}}{8}+\frac{(-1)^{m}\mathrm{Q}_{n-m-3}}{8}, \\ \langle \vec{\mathrm{P}}_{n+1}, \vec{\mathrm{P}}_{m+2} \rangle& = \frac{7\mathrm{Q}_{n+m+5}}{8}-\frac{(-1)^{m}\mathrm{Q}_{n-m-1}}{8}, \\ \langle \vec{\mathrm{P}}_{n+2}, \vec{\mathrm{P}}_{m+1} \rangle& = \frac{7\mathrm{Q}_{n+m+5}}{8}+\frac{(-1)^{m}\mathrm{Q}_{n-m+1}}{8}, \\ \langle \vec{\mathrm{P}}_{n+3}, \vec{\mathrm{P}}_{m} \rangle& = \frac{7\mathrm{Q}_{n+m+5}}{8}-\frac{(-1)^{m}\mathrm{Q}_{n-m+3}}{8}. \end{align*}$

By substituting these equalities in (3.2), we obtain the result.

Example 3.1. Let $\overrightarrow{{\mathrm{HP}}}_1 = (1, 2, 5)+(2, 5, 12)\varepsilon+(5, 12, 29)\varepsilon^*+(12, 29, 70)\varepsilon\varepsilon^*$ and $\overrightarrow{{\mathrm{HP}}}_0 = (0, 1, 2)+(1, 2, 5)\varepsilon+(2, 5, 12)\varepsilon^*+(5, 12, 29)\varepsilon\varepsilon^*$ be the hyper-dual Pell vectors. The scalar product of $\overrightarrow{{\mathrm{HP}}}_1$ and $\overrightarrow{{\mathrm{HP}}}_0$ are

$\begin{align*} \langle \overrightarrow{\mathrm{HP}}_1, \overrightarrow{\mathrm{HP}}_0 \rangle & = \frac{7\mathrm{Q}_{3}-Q_1}{8} +\frac{7\mathrm{Q}_{4}-Q_1}{4}\varepsilon +\frac{7\mathrm{Q}_{5}-3Q_1}{4}\varepsilon^* +\frac{7\mathrm{Q}_{6}-3Q_1}{2}\varepsilon\varepsilon^*\\ & = 12+59\varepsilon+142\varepsilon^*+690\varepsilon\varepsilon^*. \end{align*}$

By the other hand

$\begin{align*} \langle \overrightarrow{\mathrm{HP}}_1, \overrightarrow{\mathrm{HP}}_{0}\rangle \nonumber & = \langle\vec{\mathrm{P}}_{1}, \vec{\mathrm{P}}_{0}\rangle +\left(\langle \vec{\mathrm{P}}_{1}, \vec{\mathrm{P}}_{1}\rangle+\langle \vec{\mathrm{P}}_{2}, \vec{\mathrm{P}}_{0}\rangle\right)\varepsilon+\left(\langle\vec{\mathrm{P}}_{1}, \vec{\mathrm{P}}_{2}\rangle+\langle \vec{\mathrm{P}}_{3}, \vec{\mathrm{P}}_{0}\rangle\right)\varepsilon^*\\[6pt] &\quad+\left(\langle \vec{\mathrm{P}}_{1}, \vec{\mathrm{P}}_{3}\rangle+\langle \vec{\mathrm{P}}_{2}, \vec{\mathrm{P}}_{2}\rangle+\langle \vec{\mathrm{P}}_{3}, \vec{\mathrm{P}}_{1}\rangle+\langle \vec{\mathrm{P}}_{4}, \vec{\mathrm{P}}_{0}\rangle\right) \varepsilon\varepsilon^*\\ & = 12+(30+29)\varepsilon+(72+70)\varepsilon^*+(174+173+174+169)\varepsilon\varepsilon^*\\ & = 12+59\varepsilon+142\varepsilon^*+690\varepsilon\varepsilon^*. \end{align*}$

The results are the same as we expected.

Corollary 3.1. The norm of $\overrightarrow{{\mathrm{HP}}}_n$ is

$\begin{align} \| \overrightarrow{\mathrm{HP}}_n \|^2 = \langle \overrightarrow{\mathrm{HP}}_n, \overrightarrow{\mathrm{HP}}_n \rangle = &\frac{7\mathrm{Q}_{2n+2}}{8}-\frac{(-1)^{n}}{4} +(\frac{7\mathrm{Q}_{2n+3}}{4}-\frac{(-1)^{n}}{2})\varepsilon \\&+(\frac{7\mathrm{Q}_{2n+4}}{4}-\frac{3(-1)^{n}}{2})\varepsilon^* +(\frac{7\mathrm{Q}_{2n+5}}{2}-3(-1)^{n}\mathrm)\varepsilon\varepsilon^*. \end{align}$

(3.3)

Proof. The proof is clear from taking $m = n$ in (3.1). □

Example 3.2. Find the norm of $\overrightarrow{{\mathrm{HP}}}_1 = (1, 2, 5)+(2, 5, 12)\varepsilon+(5, 12, 29)\varepsilon^*+(12, 29, 70)\varepsilon\varepsilon^*$ .

If we take $n = 1$ in (3.3) and use (1.3), then we will get

$\begin{align*} \| \overrightarrow{\mathrm{HP}}_1 \| & = \sqrt{\frac{7\mathrm{Q}_{4}}{8}+\frac{1}{4} +(\frac{7\mathrm{Q}_{5}}{4}+\frac{1}{2})\varepsilon +(\frac{7\mathrm{Q}_{6}}{4}+\frac{3}{2})\varepsilon^* +(\frac{7\mathrm{Q}_{7}}{2}+3\mathrm)\varepsilon\varepsilon^*}\\ & = \sqrt{30+144\varepsilon+348\varepsilon^*+1676\varepsilon\varepsilon^*}\\ & = \sqrt{30}+\frac{72}{\sqrt{30}}\varepsilon+\frac{174}{\sqrt{30}}\varepsilon^*+\frac{734}{5\sqrt{30}}\varepsilon\varepsilon^*. \end{align*}$

From (1.6) and (3.1), the following cases can be given for the scalar product of hyper-dual Pell vectors $\overrightarrow{\mathrm{HP}}_n$ and $\overrightarrow{\mathrm{HP}}_m$ .

Case 3.1. Assume that $\cos\phi = 0$ and $\phi^{\ast }\neq 0$ , then $\psi = \frac{\pi}{2}$ , $\psi^* = 0$ , therefore

$\begin{align*} \langle \overrightarrow{\mathrm{HP}}_n, \overrightarrow{\mathrm{HP}}_m \rangle = -\varepsilon^* \phi^*& = (\frac{7\mathrm{Q}_{m+n+4}}{4}-\frac{3(-1)^{m}\mathrm{Q}_{n-m}}{4})\varepsilon^*+(\frac{7\mathrm{Q}_{m+n+5}}{2}-\frac{3(-1)^{m}\mathrm{Q}_{n-m}}{2}\mathrm)\varepsilon\varepsilon^*, \end{align*}$

then, we get

$\begin{equation*} \phi^* = \left( -1 \right)^{m} \left( \frac{3}{2}+ \varepsilon\right)-\frac{7}{4}\left( \mathrm{Q}_{m+n+4}+2\; \varepsilon\; \mathrm{Q}_{m+n+5}\right) \end{equation*}$

and corresponding dual lines $\mathrm{}{d}_1$ and ${d}_2$ are perpendicular such that they do not intersect each other; see Figure 1.

Figure 1. Geometric representation of hyper-dual angle between the directed dual lines

${d}_1$ and

${d}_2$ .

DownLoad: Full-Size Img PowerPoint

Case 3.2. Assume that $\phi^{\ast } = 0$ and $\phi\neq0$ , then we obtain

$\begin{align*} \langle \overrightarrow{\mathrm{HP}}_n, \overrightarrow{\mathrm{HP}}_m \rangle = \cos\phi& = \left( \frac{7\mathrm{Q}_{m+n+2}}{8}-\frac{(-1)^{m}\mathrm{Q}_{n-m}}{8} \right)+\left( \frac{7\mathrm{Q}_{m+n+3}}{4}-\frac{(-1)^{m}\mathrm{Q}_{n-m}}{4} \right)\varepsilon, \end{align*}$

therefore

$\begin{equation*} \phi = \arccos\left(\left( \frac{7\mathrm{Q}_{m+n+2}}{8}-\frac{(-1)^{m}{Q}_{n-m}}{8} \right)+\left( \frac{7\mathrm{Q}_{m+n+3}}{4}-\frac{(-1)^{m}{Q}_{n-m}}{4} \right)\varepsilon\right), \end{equation*}$

and corresponding dual lines ${d}_1$ and ${d}_2$ intersect each other; see Figure 2.

Figure 2. Intersection of dual lines

${d}_1$ and

${d}_2$ .

DownLoad: Full-Size Img PowerPoint

Case 3.3. Assume that $\cos\phi = 0$ and $\phi^* = 0$ , then $\psi = \frac{\pi}{2}$ and $\psi^{\ast } = 0$ , therefore

$\begin{align*} \langle \overrightarrow{\mathrm{HP}}_n, \overrightarrow{\mathrm{HP}}_m \rangle & = 0, \end{align*}$

and dual lines ${d}_1$ and ${d}_2$ intersect each other at a right angle; see Figure 3.

Figure 3. Perpendicular intersection of dual lines

${d}_1$ and

${d}_2$ .

DownLoad: Full-Size Img PowerPoint

Case 3.4. Assume that $\phi = 0$ and $\phi^* = 0$ , then

$\begin{align*} \langle \overrightarrow{\mathrm{HP}}_n, \overrightarrow{\mathrm{HP}}_m \rangle & = 1, \end{align*}$

in this case corresponding dual lines ${d}_1$ and ${d}_2$ are parallel; see Figure 4.

Figure 4. Parallel of dual lines

${d}_1$ and

${d}_2$ .

DownLoad: Full-Size Img PowerPoint

4. Conclusions

In the present study, we introduce two families of hyper-dual numbers with components containing Pell and the Pell-Lucas numbers. First, we define hyper-dual Pell and Pell-Lucas numbers. Afterwards, by means of the Binet's formulas of Pell and Pell-Lucas numbers, we investigate identities such as the Binet-like formulas, Vajda-like, Catalan-like, Cassini-like, and d'Ocagne-like identities. After that, we define hyper-dual Pell vector and angle with some properties and geometric applications related to them. In the future it would be valuable to replicate a similar exploration and development of our findings on hyper-dual numbers with Pell and Pell-Lucas numbers. These results can trigger further research on the subjects of the hyper-dual numbers, vector, and angle to carry out in the geometry of dual and hyper-dual space.

Author contributions

Faik Babadağ and Ali Atasoy: Conceptualization, writing-original draft, writing-review, editing. All authors have read and approved the final version of the manuscript for publication.

Conflict of interest

The authors declare that they have no conflict of interest.

References

[1]	Z. Ali, F. Rabiei, K. Shah, T. Khodadadi, Fractal-fractional order dynamical behavior of an HIV/AIDS epidemic mathematical model, Eur. Phys. J. Plus, 136 (2021), 36. https://doi.org/10.1140/epjp/s13360-020-00994-5 doi: 10.1140/epjp/s13360-020-00994-5
[2]	M. M. Amirian, Y. Jamali, The concepts and applications of fractional order differential calculus in modeling of viscoelastic systems: a primer, Crit. Rev. Biomed. Eng., 47 (2019), 249–276. https://doi.org/10.1615/CritRevBiomedEng.2018028368 doi: 10.1615/CritRevBiomedEng.2018028368
[3]	A. J. Arenas, G. Gonzalez-Parra, B. M. Chen-Charpentier, Construction of nonstandard finite difference schemes for the SI and SIR epidemic models of fractional order, Math. Comput. Simulat., 121 (2016), 48–63. https://doi.org/10.1016/j.matcom.2015.09.001 doi: 10.1016/j.matcom.2015.09.001
[4]	M. Arfan, K. Shah, A. Ullah, Fractal-fractional mathematical model of four species comprising of prey-predation, Phys. Scripta, 96 (2021), 124053. DOI 10.1088/1402-4896/ac2f37 doi: 10.1088/1402-4896/ac2f37
[5]	J. K. K. Asamoah, M. A. Owusu, Z. Jin, F. T. Oduro, A. Abidemi, E. O. Gyasi, Global stability and cost-effectiveness analysis of COVID-19 considering the impact of the environment: using data from Ghana, Chaos Soliton. Fract., 140 (2020), 110103. https://doi.org/10.1016/j.chaos.2020.110103 doi: 10.1016/j.chaos.2020.110103
[6]	J. K. K. Asamoah, E. Okyere, A. Abidemi, S. E. Moore, G. Q. Sun, Z. Jin, et al., Optimal control and comprehensive cost-effectiveness analysis for COVID-19, Results Phys., 33 (2022), 105177. https://doi.org/10.1016/j.rinp.2022.105177 doi: 10.1016/j.rinp.2022.105177
[7]	J. K. K. Asamoah, Z. Jin, G. Q. Sun, B. Seidu, E. Yankson, A. Abidemi, et al., Sensitivity assessment and optimal economic evaluation of a new COVID-19 compartmental epidemic model with control interventions, Chaos Soliton. Fract., 146 (2021), 110885. https://doi.org/10.1016/j.chaos.2021.110885 doi: 10.1016/j.chaos.2021.110885
[8]	J. K. K. Asamoah, C. S. Bornaa, B. Seidu, Z. Jin, Mathematical analysis of the effects of controls on transmission dynamics of SARS-CoV-2, Alex. Eng. J., 59 (2020), 5069–5078. https://doi.org/10.1016/j.aej.2020.09.033 doi: 10.1016/j.aej.2020.09.033
[9]	J. K. K. Asamoah, Fatmawati, A fractional mathematical model of heartwater transmission dynamics considering nymph and adult amblyomma ticks, Chaos Soliton. Fract., 174 (2023), 113905. https://doi.org/10.1016/j.chaos.2023.113905 doi: 10.1016/j.chaos.2023.113905
[10]	J. K. K. Asamoah, E. Okyere, E. Yankson, A. A. Opoku, A. Adom-Konadu, E. Acheampong, et al., Non-fractional and fractional mathematical analysis and simulations for Q fever, Chaos Soliton. Fract., 156 (2022), 111821. https://doi.org/10.1016/j.chaos.2022.111821 doi: 10.1016/j.chaos.2022.111821
[11]	J. K. K. Asamoah, Fractal-fractional model and numerical scheme based on Newton polynomial for Q fever disease under Atangana-Baleanu derivative, Results Phys., 34 (2022), 105189. https://doi.org/10.1016/j.rinp.2022.105189 doi: 10.1016/j.rinp.2022.105189
[12]	A. Atangana, S. I. Araz, Nonlinear equations with global differential and integral operators:existence, uniqueness with application to epidemiology, Results Phys., 20 (2021), 103593. https://doi.org/10.1016/j.rinp.2020.103593 doi: 10.1016/j.rinp.2020.103593
[13]	A. Atangana, D. Baleanu, New fractional derivatives with non-local and non-singular kernel theory and application to heat transfer model, Therm. Sci., 20 (2016), 763–769. https://doi.org/10.48550/arXiv.1602.03408 doi: 10.48550/arXiv.1602.03408
[14]	M. A. Abdulwasaa, M. S. Abdo, K. Shah, T. A. Nofal, S. K. Panchal, S. V. Kawale, et al., Fractal-fractional mathematical modeling and forecasting of new cases and deaths of COVID-19 epidemic outbreaks in India, Results Phys., 20 (2021), 103702. https://doi.org/10.1016/j.rinp.2020.103702 doi: 10.1016/j.rinp.2020.103702
[15]	R. P. Agarwal, S. Arshad, D. Regan, V. Lupulescu, A Schauder fixed point theorem in semilinear spaces and applications, Fixed Point Theory Appl., 2013 (2013), 306. https://doi.org/10.1186/1687-1812-2013-306 doi: 10.1186/1687-1812-2013-306
[16]	A. Atangana, S. I. Araz, New concept in calculus: piecewise differential and integral operators, Chaos Soliton. Fract., 145 (2021), 110638. https://doi.org/10.1016/j.chaos.2020.110638 doi: 10.1016/j.chaos.2020.110638
[17]	A. J. Arenas, G. González-Parra, B. M. Chen-Charpentier, Construction of nonstandard finite difference schemes for the SI and SIR epidemic models of fractional order, Math. Comput. Simulat., 121 (2016), 48–63. https://doi.org/10.1016/j.matcom.2015.09.001 doi: 10.1016/j.matcom.2015.09.001
[18]	M. S. Arshad, D. Baleanu, M. B. Riaz, M. Abbas, A novel 2-stage fractional Runge-Kutta method for a time fractional logistic growth model, Discrete Dyn. Nat. Soc., 2020 (2020), 1020472. https://doi.org/10.1155/2020/1020472 doi: 10.1155/2020/1020472
[19]	S. Boccaletti, W. Ditto, G. Mindlin, A. Atangana, Modeling and forecasting of epidemic spreading: The case of Covid-19 and beyond, Chaos Soliton. Fract., 135 (2020), 109794. https://doi.org/10.1016/j.chaos.2020.109794 doi: 10.1016/j.chaos.2020.109794
[20]	S. Banihashemi, H. Jafari, A. Babaei, A stable collocation approach to solve a neutral delay stochastic differential equation of fractional order, J. Comput. Appl. Math., 403 (2022), 113845. https://doi.org/10.1016/j.cam.2021.113845 doi: 10.1016/j.cam.2021.113845
[21]	I. I. Bogoch, A. Watts, A. Thomas-Bachli, C. Huber, M. U. Kraemer, K. Khan, Pneumonia of unknown aetiology in Wuhan, China: potential for international spread via commercial air travel, J. Travel Med., 27 (2020), taaa008. https://doi.org/10.1093/jtm/taaa008 doi: 10.1093/jtm/taaa008
[22]	C. Celauro, C. Fecarotti, A. Pirrotta, A. C. Collop, Experimental validation of a fractional model for creep/recovery testing of asphalt mixtures, Constr. Build. Mater., 36 (2012), 458–466. https://doi.org/10.1016/j.conbuildmat.2012.04.028 doi: 10.1016/j.conbuildmat.2012.04.028
[23]	Y. Chen, F. Liu, Q. Yu, T. Li, Review of fractional epidemic models, Appl. Math. Model., 97 (2021), 281–307. https://doi.org/10.1016/j.apm.2021.03.044 doi: 10.1016/j.apm.2021.03.044
[24]	A. Carpinteri, F. Mainardi, Fractals and fractional calculus in continuum mechanics, Vienna: Springer, 1997. https://doi.org/10.1007/978-3-7091-2664-6
[25]	C.T. Deressa, G. F. Duressa, Analysis of Atangana-Baleanu fractional-order SEAIR epidemic model with optimal control, Adv. Differ. Equ., 2021 (2021), 174. https://doi.org/10.1186/s13662-021-03334-8 doi: 10.1186/s13662-021-03334-8
[26]	Z. Dai, Y. Peng, H. A. Mansy, R. H. Sandler, T. J. Royston, A model of lung parenchyma stress relaxation using fractional viscoelasticity, Med. Eng. Phys., 37 (2015), 752–758. https://doi.org/10.1016/j.medengphy.2015.05.003 doi: 10.1016/j.medengphy.2015.05.003
[27]	Z. J. Fu, Z. C. Tang, H. T. Zhao, P. W. Li, T. Rabczuk, Numerical solutions of the coupled unsteady nonlinear convection-diffusion equations based on generalized finite difference method, Eur. Phys. J. Plus, 134 (2019), 272. https://doi.org/10.1140/epjp/i2019-12786-7 doi: 10.1140/epjp/i2019-12786-7
[28]	A. B. Gumel, S. Ruan, T. Day, J. Watmough, F. Brauer, P. van den Driessche, et al., Modelling strategies for controlling SARS out breaks, Proc. R. Soc. Lond. B., 271 (2004), 2223–2232. https://doi.org/10.1098/rspb.2004.2800 doi: 10.1098/rspb.2004.2800
[29]	E. F. D. Goufo, Application of the Caputo-Fabrizio fractional derivative without singular kernel to Korteweg-de Vries-Burgers equation, Math. Model. Anal., 21 (2016), 188–198. https://doi.org/10.3846/13926292.2016.1145607 doi: 10.3846/13926292.2016.1145607
[30]	D. S. Hui, E. I. Azhar, T. A. Madani, F. Ntoumi, R. Kock, O. Dar, et al., The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health–The latest 2019 novel coronavirus outbreak in Wuhan, China, B. Math. Biol., 91 (2020), 264–66. https://doi.org/10.1016/j.ijid.2020.01.009 doi: 10.1016/j.ijid.2020.01.009
[31]	S. Hussain, E. N. Madi, H. Khan, H. Gulzar, S. Etemad, S. Rezapour, et al., On the stochastic modeling of COVID-19 under the environmental white noise, J. Funct. Space., 2022 (2022), 4320865. doilinkhttps://doi.org/10.1155/2022/4320865 doi: 10.1155/2022/4320865
[32]	A. A. Hamou, E. Azroul, Z. Hammouch, A. L. Alaoui, On dynamics of fractional incommensurate model of Covid-19 with nonlinear saturated incidence rate, MedRxiv, 2021 (2021), 07, https://doi.org/10.1101/2021.07.18.21260711 doi: 10.1101/2021.07.18.21260711
[33]	M. T. Hoang, O. F. Egbelowo, Dynamics of a fractional-order hepatitis B epidemic model and its solutions by nonstandard numerical schemes, In: Mathematical modelling and analysis of infectious diseases, Cham: Springer, 2020,127–153. https://doi.org/10.1007/978-3-030-49896-2_5
[34]	G. Jumarie, Stochastic differential equations with fractional Brownian motion input, Int. J. Syst. Sci., 24 (1993), 1113–1131. https://doi.org/10.1080/00207729308949547 doi: 10.1080/00207729308949547
[35]	S. Kumar, A. Kumar, B. Samet, H. Dutta, A study on fractional host-parasitoid population dynamical model to describe insect species, Numer. Meth. Part. D. E., 37 (2021), 1673–1692. https://doi.org/10.1002/num.22603 doi: 10.1002/num.22603
[36]	M. M. Khalsaraei, An improvement on the positivity results for 2-stage explicit Runge-Kutta methods, J. Comput. Appl. Math., 235 (2010), 137–143. https://doi.org/10.1016/j.cam.2010.05.020 doi: 10.1016/j.cam.2010.05.020
[37]	R. Kahn, I. Holmdahl, S. Reddy, J. Jernigan, M. J. Mina, R. B. Slayton, Mathematical modeling to inform vaccination strategies and testing approaches for voronavirus disease 2019 (COVID-19) in nursing homes, Clin. Infect. Dis., 74 (2022), 597–603. https://doi.org/10.1093/cid/ciab517 doi: 10.1093/cid/ciab517
[38]	M. A. Khan, A. Atangana, Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative, Alex. Eng. J., 59 (2020), 2379–2389. https://doi.org/10.1016/j.aej.2020.02.033 doi: 10.1016/j.aej.2020.02.033
[39]	M. A. Khan, A. Atangana, E. Alzahrani, E. Fatmawati, The dynamics of COVID-19 with quarantined and isolation, Adv. Differ. Equ., 2020 (2020), 425. https://doi.org/10.1186/s13662-020-02882-9 doi: 10.1186/s13662-020-02882-9
[40]	A. M. Lopes, J. T. Machado, Fractional order models of leaves, J. Vib. Control, 20 (2014), 998–1008. https://doi.org/10.1177/1077546312473323 doi: 10.1177/1077546312473323
[41]	R. Li, S. Zhong, C. Swartz, An improvement of the Arzela-Ascoli theorem, Topol. Appl., 159 (2012), 2058–2061. http://doi.org/10.1016/j.topol.2012.01.014 doi: 10.1016/j.topol.2012.01.014
[42]	R. Lewandowski, Z. Pawlak, Dynamic analysis of frames with viscoelastic dampers modelled by rheological models with fractional derivatives, J. Sound Vib., 330 (2011), 923–936. https://doi.org/10.1016/j.jsv.2010.09.017 doi: 10.1016/j.jsv.2010.09.017
[43]	B. Li, H. Liang, L. Shi, Q. He, Complex dynamics of Kopel model with nonsymmetric response between oligopolists, Chaos Soliton. Fract., 156 (2022), 111860. https://doi.org/10.1016/j.chaos.2022.111860 doi: 10.1016/j.chaos.2022.111860
[44]	F. Liu, K. Burrage, Novel techniques in parameter estimation for fractional dynamical models arising from biological systems, Comput. Math. Appl., 62 (2011), 822–833. https://doi.org/10.1016/j.camwa.2011.03.002 doi: 10.1016/j.camwa.2011.03.002
[45]	B. Li, H. Liang, Q. He, Multiple and generic bifurcation analysis of a discrete Hindmarsh-Rose model, Chaos Soliton. Fract., 146 (2021), 110856. https://doi.org/10.1016/j.chaos.2021.110856 doi: 10.1016/j.chaos.2021.110856
[46]	J. Mondal, S. Khajanchi, Mathematical modeling and optimal intervention strategies of the COVID-19 outbreak, Nonlinear Dyn., 109 (2022), 177–202. https://doi.org/10.1007/s11071-022-07235-7 doi: 10.1007/s11071-022-07235-7
[47]	J. T. Machado, V. Kiryakova, F. Mainardi, Recent history of fractional calculus, Commun. Nonlinear Sci., 16 (2011), 1140–1153. https://doi.org/10.1016/j.cnsns.2010.05.027 doi: 10.1016/j.cnsns.2010.05.027
[48]	F. C. Meral, T. J. Royston, R. Magin, Fractional calculus in viscoelasticity: an experimental study, Commun. Nonlinear Sci., 15 (2010), 939–945. https://doi.org/10.1016/j.cnsns.2009.05.004 doi: 10.1016/j.cnsns.2009.05.004
[49]	R. L. Magin, Fractional calculus in bioengineering, Crit. Rev. Biomed. Eng., 32 (2004), 1–104. https://doi.org/10.1615/critrevbiomedeng.v32.i1.10 doi: 10.1615/critrevbiomedeng.v32.i1.10
[50]	R. L. Magin, Fractional Calculus in bioengineering, Redding: Begell House, 2006.
[51]	F. Mainardi, An historical perspective on fractional calculus in linear viscoelasticity, Fract. Calc. Appl. Anal., 15 (2012), 712–717. https://doi.org/10.2478/s13540-012-0048-6 doi: 10.2478/s13540-012-0048-6
[52]	I. Nesteruk, Statistics based predictions of coronavirus 2019-nCoV spreading in mainland China, MedRxiv, 4 (2020), 1988–1989. https://doi.org/10.1101/2020.02.12.20021931 doi: 10.1101/2020.02.12.20021931
[53]	O. A. Omar, R. A. Elbarkouky, H. M. Ahmed, Fractional stochastic modelling of COVID-19 under wide spread of vaccinations: Egyptian case study, Alex. Eng. J., 61 (2022), 8595–8609. https://doi.org/10.1016/j.aej.2022.02.002. doi: 10.1016/j.aej.2022.02.002
[54]	J. C. Pedjeu, G. S. Ladde, Stochastic fractional differential equations: Modeling, method and analysis, Chaos Soliton. Fract., 45 (2012), 279–293. https://doi.org/10.1016/j.chaos.2011.12.009 doi: 10.1016/j.chaos.2011.12.009
[55]	A. Y. Rossikhin, M. V. Shitikova, Applications of fractional calculus to dynamic problems of linear and nonlinear hereditary mechanics of solids, Appl. Mech. Rev., 50 (1997), 15–67. https://doi.org/10.1115/1.3101682 doi: 10.1115/1.3101682
[56]	A. Radulescu, C. Williams, K. Cavanagh, Management strategies in a SEIR-type model of COVID 19 community spread, Sci. Rep., 10 (2020), 21256. https://doi.org/10.1038/s41598-020-77628-4 doi: 10.1038/s41598-020-77628-4
[57]	Y. B. Sang, Critical Kirchhoff-Choquard system involving the fractional p-Laplacian operator and singular nonlinearities, Topol. Method. Nonl. An., 58 (2021), 233–274. https://doi.org/10.12775/TMNA.2020.070 doi: 10.12775/TMNA.2020.070
[58]	K. Shah, R. Din, W. Deebani, P. Kumam, Z. Shah, On nonlinear classical and fractional order dynamical system addressing COVID-19, Results Phys., 24 (2021), 104069. https://doi.org/10.1016/j.rinp.2021.104069 doi: 10.1016/j.rinp.2021.104069
[59]	M. Shimizu, W. Zhang, Fractional calculus approach to dynamic problems of viscoelastic materials, JSME International Journal Series C Mechanical Systems, Machine Elements and Manufacturing, 42 (1999), 825–837. https://doi.org/10.1299/jsmec.42.825 doi: 10.1299/jsmec.42.825
[60]	Y. B. Sang, S. H. Liang, Fractional Kirchhoff-Choquard equation involving Schrodinger term and upper critical exponent, J. Geom. Anal., 32 (2022), 5. https://doi.org/10.1007/s12220-021-00747-5 doi: 10.1007/s12220-021-00747-5
[61]	L. Stella, A. P. Martínez, D. Bauso, P. Colaneri, The role of asymptomatic infections in the COVID-19 epidemic via complex networks and stability analysis, SIAM J. Control Optim., 60 (2022), S119–S144. https://doi.org/10.1137/20M1373335 doi: 10.1137/20M1373335
[62]	D. Valério, J. T. Machado, V. Kiryakova, Some pioneers of the applications of fractional calculus, Fract. Calc. Appl. Anal., 17 (2014), 552–578. https://doi.org/10.2478/s13540-014-0185-1 doi: 10.2478/s13540-014-0185-1
[63]	Y. Wang, Z. Wei, J. Cao, Epidemic dynamics of influenza-like diseases spreading in complex networks, Nonlinear Dyn., 101 (2020), 1801–1820. https://doi.org/10.1007/s11071-020-05867-1 doi: 10.1007/s11071-020-05867-1
[64]	J. T. Wu, K. Leung, G. M. Leung, Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study, Lancet, 395 (2020), 689–697. https://doi.org/10.1016/S0140-6736(20)30260-9 doi: 10.1016/S0140-6736(20)30260-9
[65]	B. Wang, L. Li, Y. Wang, An efficient nonstandard finite difference scheme for chaotic fractional-order Chen system, IEEE Access, 8 (2020), 98410–98421. https://doi.org/10.1109/ACCESS.2020.2996271 doi: 10.1109/ACCESS.2020.2996271
[66]	G. C. Wu, M. Luo, L. L. Huang, S. Banerjee, Short memory fractional differential equations for new memristor and neural network design, Nonlinear Dyn., 100 (2020), 3611–3623. https://doi.org/10.1007/s11071-020-05572-z doi: 10.1007/s11071-020-05572-z
[67]	A. Zeb, A. Atangana, Z. A. Khan, S. Djillali, A robust study of a piecewise fractional order COVID-19 mathematical model, Alex. Eng. J., 61 (2022), 5649–5665. https://doi.org/10.1016/j.aej.2021.11.039 doi: 10.1016/j.aej.2021.11.039
[68]	S. Zhao, Q. Lin, J. Ran, S. S. Musa, G. Yang, W. Wang, et al., Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak, Int. J. Infect. Dis., 92 (2020), 214–217. https://doi.org/10.1016/j.ijid.2020.01.050 doi: 10.1016/j.ijid.2020.01.050
[69]	S. Zhao, S. S. Musa, Q. Lin, J. Ran, G. Yang, W. Wang, et al., Estimating the unreported number of novel coronavirus (2019-nCoV) cases in China in the first half of January 2020: A data-driven modelling analysis of the early outbreak, J. Clin. Med., 9 (2020), 388. http://doi.org/10.3390/jcm9020388 doi: 10.3390/jcm9020388
[70]	P. Zhou, X. L. Yang, X. G. Wang, B. Hu, L. Zhang, W. Zhang, et al., A pneumonia outbreak associated with a new coronavirus of probable bat origin, Nature, 579 (2020), 270–273. https://doi.org/10.1038/s41586-020-2012-7 doi: 10.1038/s41586-020-2012-7
[71]	Y. Zhang, H. Sun, H. H. Stowell, M. Zayernouri, S. E. Hansen, A review of applications of fractional calculus in Earth system dynamics, Chaos Soliton. Fract., 102 (2017), 29–46. https://doi.org/10.1016/j.chaos.2017.03.051 doi: 10.1016/j.chaos.2017.03.051
[72]	Naming the coronavirus disease (COVID-19) and the virus that causes it, Available from: World Health Organization (WHO), 2019. https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance/naming-the-coronavirus-disease-(covid-2019)-and-the-virus-that-causes-it.
[73]	World Health Organization, Japan situation, 2020. Available from: https://covid19.who.int/region/wpro/country/jp.
[74]	Japan COVID coronavirus statistics, 2023. Available from: https://www.worldometers.info/coronavirus/country/japan/.
[75]	Japan COVID cases, 2021. Available from: https://www.nytimes.com/interactive/2021/world/japan-covid-cases.html.

Reader Comments

Your name:*

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