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]	Dayang Dai, Dabuxilatu Wang . A generalized Liu-type estimator for logistic partial linear regression model with multicollinearity. AIMS Mathematics, 2023, 8(5): 11851-11874. doi: 10.3934/math.2023600
[2]	Muhammad Nauman Akram, Muhammad Amin, Ahmed Elhassanein, Muhammad Aman Ullah . A new modified ridge-type estimator for the beta regression model: simulation and application. AIMS Mathematics, 2022, 7(1): 1035-1057. doi: 10.3934/math.2022062
[3]	Sihem Semmar, Omar Fetitah, Mohammed Kadi Attouch, Salah Khardani, Ibrahim M. Almanjahie . A Bernstein polynomial approach of the robust regression. AIMS Mathematics, 2024, 9(11): 32409-32441. doi: 10.3934/math.20241554
[4]	Yanting Xiao, Wanying Dong . Robust estimation for varying-coefficient partially linear measurement error model with auxiliary instrumental variables. AIMS Mathematics, 2023, 8(8): 18373-18391. doi: 10.3934/math.2023934
[5]	Gaosheng Liu, Yang Bai . Statistical inference in functional semiparametric spatial autoregressive model. AIMS Mathematics, 2021, 6(10): 10890-10906. doi: 10.3934/math.2021633
[6]	Juxia Xiao, Ping Yu, Zhongzhan Zhang . Weighted composite asymmetric Huber estimation for partial functional linear models. AIMS Mathematics, 2022, 7(5): 7657-7684. doi: 10.3934/math.2022430
[7]	Xin Liang, Xingfa Zhang, Yuan Li, Chunliang Deng . Daily nonparametric ARCH(1) model estimation using intraday high frequency data. AIMS Mathematics, 2021, 6(4): 3455-3464. doi: 10.3934/math.2021206
[8]	Emrah Altun, Mustafa Ç. Korkmaz, M. El-Morshedy, M. S. Eliwa . The extended gamma distribution with regression model and applications. AIMS Mathematics, 2021, 6(3): 2418-2439. doi: 10.3934/math.2021147
[9]	Muhammad Amin, Saima Afzal, Muhammad Nauman Akram, Abdisalam Hassan Muse, Ahlam H. Tolba, Tahani A. Abushal . Outlier detection in gamma regression using Pearson residuals: Simulation and an application. AIMS Mathematics, 2022, 7(8): 15331-15347. doi: 10.3934/math.2022840
[10]	Zawar Hussain, Atif Akbar, Mohammed M. A. Almazah, A. Y. Al-Rezami, Fuad S. Al-Duais . Diagnostic power of some graphical methods in geometric regression model addressing cervical cancer data. AIMS Mathematics, 2024, 9(2): 4057-4075. doi: 10.3934/math.2024198

Abstract

1. Introduction

The field of mathematical analysis that deals with the study of arbitrary order integrals and derivatives is known as fractional calculus. Because of its numerous applications across a wide range of fields, this field has increased in importance and recognition over the past few years. According to researchers, this field is the most effective at identifying anomalous kinetics and has numerous uses in a variety of fields. Ordinary differential equations with fractional derivatives can be used to simulate a variety of issues, including statistical, mathematical, engineering, chemical, and biological issues. Several distinct forms of fractional integrals and derivative operators (see e.g., ^[1,2,3,4]), including Riemann-Liouville, Caputo, Riesz, Hilfer, Hadamard, Erdélyi-Kober, Saigo, Marichev-Saigo-Maeda and others, have been thoroughly investigated by researchers. From an application perspective, we suggest the readers to see the work related to the fractional differential equations presented by ^[5,6,7,8]. In ^[9], the authors studied symmetric and antisymmetric solitons in the defocused saturable nonlinearity and the PT-symmetric potential of the fractional nonlinear Schrödinger equation. In ^[10], the fractional exponential function approach is used to study a time-fractional Ablowitz-Ladik model, and bright and dark discrete soliton solutions, discrete exponential solutions, and discrete peculiar wave solutions are discovered. In ^[11], the authors presented the rich vector exact solutions for the coupled discrete conformable fractional nonlinear Schrödinger equations by taking into account the conformable fractional derivative.

On the other hand, special functions like Gamma, Beta, Mittag-Leffler, et al. play a vital part in the foundation of fractional calculus. Moreover, the Mittag-Leffler function is regarded as the fundamental function in fractional calculus. The Prabhakar fractional operator containing a three-parameter version of the aforementioned function in the kernel. The M-L function has been extensively studied to construct solutions of fractional PDEs, such as dynamical characteristic of analytical fractional solitons for the space-time fractional Fokas-Lenells equation, soliton dynamics based on exact solutions of conformable fractional discrete complex cubic Ginzburg-Landau equation and Abundant fractional soliton solutions of a space-time fractional perturbed Gerdjikov-Ivanov equation by a fractional mapping method, see ^[12,13,14]. Strong generalizations of the univariate and bivariate Mittag-Leffler functions, which are known to be important in fractional calculus, are the multivariate Mittag-Leffler functions.

The well-known one-parameter Mittag-Leffler (M-L) function is defined by ^[15,16] as follows

$\begin{equation} \varepsilon_{a}(z_1) = \sum\limits_{l = 0}^{\infty}\frac{z_1^l}{\Gamma(al+1)}\,\,\,\,\,\,\,\,\,(a\in\mathbb{C};\Re(a) > 0,z_1\in\mathbb{C}), \end{equation}$

(1.1)

where $\mathbb{C}$ represents the set of complex numbers and $\Re(a)$ denotes the real part of the complex number.

The generalization of (1.1) with two parameters is defined by ^[17,18] as

$\begin{equation} \varepsilon_{a,b}(z_1) = \sum\limits_{l = 0}^{\infty}\frac{z_1^l}{\Gamma(al+b)}\,\,\,\,\,\,\,\,\,(a,b\in\mathbb{C};\Re(a) > 0,\Re(b) > 0), \end{equation}$

(1.2)

Later on, Agarwal ^[19], Humbert ^[20] and Humbert and Agarwal ^[21] studied the properties and applications of M-L functions. In ^[22], the generalization of (1.1) and (1.2) is defined by

$\begin{equation} \varepsilon_{a,b}^{c}(z_1) = \sum\limits_{l = 0}^{\infty}\frac{(c)_l}{\Gamma(a l+b)}\frac{z_1^l}{l!}\,\,\,\,\,\,\,\,\,(a,b,c\in\mathbb{C};\Re(a) > 0,\Re(b) > 0). \end{equation}$

(1.3)

In ^[23], the following generalization of the M-L function is defined by

$\begin{equation} \varepsilon_{a,b}^{c,q}(z_1) = \sum\limits_{l = 0}^{\infty}\frac{(c)_{lq}}{\Gamma(a l+b)}\frac{z_1^l}{l!}\,\,\,\,\,\,\,\,\,(a,b,c\in\mathbb{C};\Re(a) > 0,\Re(b) > 0, q > 0). \end{equation}$

(1.4)

In ^[24], Rahman et al. proposed the following generalized of M-L function by

$\begin{align} \varepsilon_{a,b,p}^{c,q,d}(z_1) = \sum\limits_{l = 0}^{\infty}\frac{B_p(c +lq; d-c)(d)_{lq}}{B(c, d-c)\Gamma(a l+b)}\frac{z_1^l}{l!}, \end{align}$

(1.5)

where $a, b, c, d\in\mathbb{C}; \Re(c) > 0, \Re(a) > 0, \Re(b) > 0, q > 0$ and $B_p(x, y) = \int_0^1 t^{x-1}(1-t)^{y-1}e^{-t-\frac{p}{t}}dt$ is the extension of beta function (see ^[25]).

Gorenflo et al. ^[26] and Haubold et al. ^[27]) studied the various properties of generalized M-L function. In ^[28], a new generalization of M-L function (1.3) is presented by

$\begin{equation} \varepsilon_{a,b,p}^{c}(z_1) = \sum\limits_{l = 0}^{\infty}\frac{(c;p)_l}{\Gamma(a l+b)}\frac{z_1^l}{l!}\,\,\,\,\,\,\,\,\,(p\geq0, a,b,c\in\mathbb{C}; , \Re(a) > 0,\Re(b) > 0,), \end{equation}$

(1.6)

where $(\lambda; p)_l$ is the Pochhammer symbol defined by Srivastava et al. ^[29,30] as

$\begin{equation} (\lambda;p)_{\mu} = \begin{cases} \frac{\Gamma_{p}(\lambda+\mu)}{\Gamma(\lambda)}; &(p > 0, \lambda,\mu\in\mathbb{C} )\\ \\ (\lambda)_{\mu}; & (p = 0,\ {\lambda,\mu}\in \mathbb{C}\setminus\{{0}\}.\ \ \ \ \ \ \ \ \ \ \ \ \end{cases} \end{equation}$

(1.7)

The researchers examined the developments of these extension, (1.6) and (1.7) and studied their related features and applications. In ^[30], Srivastava et al. proposed the following generalized hypergeometric function

$\begin{equation} _sF_t[(\delta_1; p),\cdots, (\delta_s); (\zeta_1),\cdots, (\zeta_t); z_1] = \sum\limits_{l = 0}^{\infty}\frac{(\delta_{1}; p)_l\cdots (\delta_s)_l}{(\zeta_1)_l\cdots (\zeta_t)_l}\ \frac{z_1^l}{l!}\, , \end{equation}$

(1.8)

where $\delta_{j}\in\mathbb{C}$ for j = 1, 2, $\cdots$ , s, $\zeta_{k}\in\mathbb{C}$ for $k = 1, 2, \cdots, t,$ and $\zeta_{k}\neq$ 0, -1, -2, $\cdots$ .

The integral representation of $(\mu; p)_\eta$ is explained by

$\begin{equation} (\mu;p)_{\eta} = \frac{1}{\Gamma(\mu)}\ \int_0^\infty\ s^{\mu+\eta-1}\ e^{-s-\frac{p}{s}} d s, \end{equation}$

(1.9)

where $\Re(\rho) > 0$ and $\Re(\mu+\eta) > 0$ . In particular, the related confluent hypergeometric function $_1F_1$ and the Gauss hypergeometric function $_2F_1$ are given by

$\begin{equation} _2F_1[(\delta_1; p),b;\lambda; z_1] = \sum\limits_{l = 0}^{\infty}\frac{(\delta_{1}; p)_l (b)_l}{(\lambda)_l}\ \frac{z_1^l}{l!}\, , \end{equation}$

(1.10)

and

$\begin{equation} _1F_1[(\delta_1; p);\lambda; z_1] = \Phi[(\delta_1; p);\lambda; z_1] = \sum\limits_{l = 0}^{\infty}\frac{(\delta_{1}; p)_l}{(\lambda)_l}\ \frac{z_1^l}{l!}\, . \end{equation}$

(1.11)

The expansion of the generalised hypergeometric function $_rF_s$ , which was studied by ^[30], has $r$ numerator and $s$ denominator parameters. Researchers recently developed several extensions of special functions, together with their corresponding characteristics and applications. Using extended beta functions as its foundation, Nisar et al. ^[31], Bohner et al. ^[32] and Rahman et al. ^[33] developed an enlargement of fractional derivative operators.

The multivariate M-L function is defined by ^[34] as follows:

$\begin{align} \mathcal{E}_{(a_j), b}^{(c_j)} (z_1, z_2, \dots, z_j)& = \mathcal{E}_{(a_1, a_2,\dots,a_j), b}^{(c_1, c_2 ,\dots, c_j)} (z_1, z_2, \dots z_j) \\ = & \sum\limits_{m_1, m_2, \dots, m_j = 0}^{\infty} \frac{(c_1)_{m_1}(c_2)_{m_2} \dots (c_j)_{m_j} (z_1)^{m_1} \dots (z_j)^{m_j}}{\Gamma(a_{1}m_1+a_2 m_2+\dots a_j m_j + b)m_1! \dots m_j!}, \end{align}$

(1.12)

where $z_i, a_i, b, c_i \in \mathbb{C}$ ; $i = 1, 2, \ldots, j$ , $\Re{(a_i)} > 0$ , $\Re{(b)} > 0$ and $\Re{(c_i)} > 0$ .

In ^{[35,36,37,38,39]}, the authors have studied various properties and applications of different type of generalized M-L functions. For real (complex) valued functions, the Lebesgue measurable space is defined by

$\begin{equation} \textsf{L}(r,s) = \left\{ \begin{array}{c} h: \|h\|_1 = \int_{r}^{s}|h(x)|dx < \infty \\ \end{array} \right\}. \end{equation}$

(1.13)

The left and right sides fractional integral operators of the Riemann-Liouville type are defined by ^[3,4] as follows:

$\begin{equation} (\mathfrak{I}_{r+}^{\lambda}h)(x) = \frac{1}{\Gamma(\lambda)}\int\limits_{r}^{x}\frac{h(\varrho)}{(x-\varrho)^{1-\lambda}}d\varrho,\; \; \; (x > r), \end{equation}$

(1.14)

and

$\begin{equation*} (\mathfrak{I}_{s-}^{\lambda}h)(x) = \frac{1}{\Gamma(\lambda)}\int\limits_{x}^{s}\frac{h(\varrho)}{(\varrho-x)^{1-\lambda}}d\varrho,\; \; \; (x < s), \end{equation*}$

where $h\in \textsf{L}(r, s)$ , $\lambda\in\mathbb{C}$ and $\Re(\lambda) > 0$ .

The left and right sides Riemann-Liouville fractional derivatives for the function $h(x)\in \textsf{L}(r, s)$ , $\lambda\in\mathbb{C}$ , $\Re(\lambda) > 0$ and $n = [\Re(\lambda)]+1$ are defined in ^[3,4] by

$\begin{equation} (\mathfrak{D}_{r+}^{\lambda}h)(x) = (\frac{d}{dx})^n(\mathfrak{I}_{r+}^{n-\lambda}h)(x) \end{equation}$

(1.15)

and

$\begin{equation*} (\mathfrak{D}_{s-}^{\lambda}h)(x) = (-\frac{d}{dx})^n(\mathfrak{I}_{s-}^{n-\lambda}h)(x), \end{equation*}$

respectively. The generalized differential operator $\mathfrak{D}_{r+}^{\lambda, v}$ of order $0 < \lambda < 1$ and type $0 < v < 1$ with respect to $x$ can be found in ^[2,4] as

$\begin{equation} (\mathfrak{D}_{r+}^{\lambda,v}h) = \left( \begin{array}{c} \mathfrak{I}_{r+}^{v(1-\lambda)}\frac{d}{dx}(\mathfrak{I}_{r+}^{(1-v)(1-\lambda)}h) \\ \end{array} \right)(x). \end{equation}$

(1.16)

In particular, if $v = 0,$ then (1.16) will lead to the operator $\mathfrak{D}_{r+}^{\lambda}$ defined in (1.15).

We also take into account the aforementioned well-known results.

Theorem 1.1. In ^[40], the following result for the fractional integral is presented by

$\begin{equation} \mathfrak{I}_{r+}^{\lambda}(\varrho-r)^{\eta-1} = \frac{\Gamma(\eta)}{\Gamma(\lambda+\eta)}(x-r)^{\lambda+\eta-1}, \end{equation}$

(1.17)

where $\lambda$ , $\eta\in\mathbb{C}$ , $\Re(\lambda) > 0$ , $\Re(\eta) > 0$ ,

Theorem 1.2. ^[41] Suppose that the function $h(z)$ has a power series expansion $h(z) = \sum\limits_{k = 0}^{\infty}k_nz^k$ and it is analytic in the disc $|z| < R$ , then we have the following result

$\begin{equation*} \label{a11} \mathfrak{D}_{z}^{\lambda}\{z^{\eta-1}h(z)\} = \frac{\Gamma(\eta)}{\Gamma(\lambda+\eta)}\sum\limits_{n = 0}^{\infty}\frac{a_n(\eta)_n}{(\lambda+\eta)_n}z^n. \end{equation*}$

Lemma 1.1. (Srivastava and Tomovski ^[42]) Suppose that $x > r$ , $\lambda\in(0, 1)$ , $v\in[0, 1]$ , $\Re(\eta) > 0$ and $\Re(\lambda) > 0$ , then we have

$\begin{equation} \mathfrak{D}_{r+}^{\lambda}[(\varrho-r)^{\eta-1}](x) = \frac{\Gamma(\eta)}{\Gamma(\eta-\lambda)}(x-r)^{\eta-\lambda-1}. \end{equation}$

(1.18)

The generalized multivariate M-L function (1.12) is then defined in terms of the modified Pochhammer symbol (1.7) and its different features as well as the accompanying integral operators are examined. This is driven by the aforementioned modifications of special functions.

Motivated by the above results and literature, the paper has the following structure: First, we describe and investigate a novel generalization of the multivariate M-L function using a generalized Pochhammer symbol. Secondly, we offer a few differential and fractional integral formulas for the explored multivariate M-L function. By using the new form of the multivariate M-L function, a new generalization of the fractional integral operator is introduced, and some fundamental characteristics of the operator are discussed.

2. The generalized multivariate of M-L function

We are in a position to present the generalized multivariate M-L function by utilizing the extended Pochhammer symbol in (1.7) as follows:

$\begin{equation} \varepsilon_{(a_j),b;p}^{(c_j)}(z_1,z_2,\cdots,z_j) = \sum\limits_{l_1,\cdots,l_j = 0}^{\infty}\frac{(c_1;p)_{l_1}(c_2)_{l_2}\cdots(c_j)_{l_j}}{\Gamma(a_1 l_1+a_2l_2+\cdots+c_jl_j+b)}\frac{z_1^{l_1} z_2^{l_2}\cdots z_j^{l_j}}{l_1!\cdots l_j!}, \end{equation}$

(2.1)

where $a_i, b, c_i\in\mathbb{C}; \Re(a_i) > 0, \Re(b) > 0, , p\geq0$ for $i = 1, 2, \cdots, j$ . The special case for $a_1 = 1$ and $l_2 = \cdots = l_j = 0$ in (2.1) can be reduced to extended confluent hypergeometric function (1.11) as follows:

$\begin{equation} \varepsilon_{1,b;p}^{c_1}(z_1) = \frac{1}{\Gamma(b)} {_1F_1}[(c_1; p);b; z_1] = \frac{1}{\Gamma(b)} \Phi[(c_1; p);b; z_1]. \end{equation}$

(2.2)

In coming results, we demonstrate some fundamental characteristics and integral representations of the generalized multivariate M-L function.

Theorem 2.1. For the multivariate M-L function defined in (2.1), the following relation holds true:

$\begin{equation} \varepsilon_{(a_j),b;p}^{(c_j)}(z_1,z_2,\cdots,z_j) = b \varepsilon_{(a_j),b+1;p}^{(c_j)}(z_1,z_2,\cdots,z_j) \end{equation}$

(2.3)

$\begin{equation*} +[a_1z_1\frac{d}{dz_1}+\cdots+ a_jz_j\frac{d}{dz_j}] \varepsilon_{(a_j),b+1;p}^{(c_j)}(z_1,\cdots,z_j), \end{equation*}$

where $a_i, b, c_i\in\mathbb{C}; \Re(a_i) > 0, \Re(b) > 0, , p\geq0$ for $i = 1, 2, \cdots, j$ .

Proof. From (2.1), we have

$\begin{align*} &b \varepsilon_{(a_j),b+1,p}^{(c_j)}(z_1,\cdots,z_j)+[a_1z_1\frac{d}{dz_1}+\cdots+ a_jz_j\frac{d}{dz_j}]\varepsilon_{(a_j),b+1;p}^{(c_j)}(z_1,\cdots,z_j)\\ = &b \sum\limits_{l_1,\cdots,l_j = 0}^{\infty}\frac{ (c_1,p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b+1)}\frac{z_1^{l_1}\cdots z_j^{l_j}}{l_1!\cdots l_j!}\\+&[a_1z_1\frac{d}{dz_1}+\cdots+ a_jz_j\frac{d}{dz_j}]\sum\limits_{l_1,\cdots,l_j = 0}^{\infty}\frac{ (c_1,p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b+1)}\frac{z_1^{l_1}\cdots z_j^{l_j}}{l_1!\cdots l_j!}\\ = &b \sum\limits_{l_1,\cdots,l_j = 0}^{\infty}\frac{ (c_1,p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b+1)}\frac{z_1^{l_1}\cdots z_j^{l_j}}{l_1!\cdots l_j!}\\+&[a_1z_1 \frac{d}{dz_1}+\cdots +a_jz_j \frac{d}{dz_j}]\sum\limits_{l_1,\cdots,l_j = 0}^{\infty}\frac{ (c_1,p)_{l_1}(c_2)_{l_2}\cdots(c_j)_{l_j}}{\Gamma(a_1 l_1+\cdots+a_jl_j+b+1)}\frac{z_1^{l_1}\cdots z_jl_j}{l_1!\cdots l_j!}\\ = &b \sum\limits_{l_1,\cdots,l_j = 0}^{\infty}\frac{ (c_1,p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b+1)}\frac{z_1^{l_1}\cdots z_j^{l_j}}{l_1!\cdots l_j!}\\+&\sum\limits_{l_1,\cdots,l_j = 0}^{\infty}\frac{(c_1,p)_{l_1}(c_2)_{l_2}\cdots(c_j)_{l_j} }{\Gamma(a_1 l_1+\cdots+a_jl_j+b+1)}\frac{z_1^{l_1}\cdots z_j^{l_j}}{l_1!\cdots l_j!}(a_1l_1+\cdots+a_jl_j)\\ = & \sum\limits_{l_1,\cdots,l_j = 0}^{\infty}\frac{ (c_1,p)_{l_1}(c_2)_{l_2}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b+1)}\frac{z_1^{l_1}\cdots z_j^{l_j}}{l_1!\cdots l_j!}(a_1l_1+\cdots+a_jl_j+b)\ \ (\mbox{using} \ \Gamma(z_1+1) = z_1 \Gamma(z_1))\\ = &\sum\limits_{l = 0}^{\infty}\frac{(c_1,p)_{l_1}(c_2)_{l_2}\cdots(c_j)_{l_j} }{ \Gamma(a_1l_1+\cdots+a_jl_j+b)}\frac{z_1^{l_1}\cdots z_j^{l_j}}{l_1!\cdots l_j!}\\ = &\varepsilon_{(a_j),b,p}^{(c_j)}(z_1,z_2,\cdots,z_j), \end{align*}$

which is the desired result (2.3).

Theorem 2.2. For the generalized multivariate M-L function defined in (1.12), the following relations hold true:

$\begin{equation} \big(\frac{d}{dz_1}\big)^m\cdots(\frac{d}{dz_j}\big)^m \varepsilon_{(a_j),b;p}^{(c_j)}(z_1,z_2,\cdots,z_j) = (c_1)_m\cdots(c_j)_m \varepsilon_{(a_j),b+(a_j)m;p}^{(c_j)+m}(z_1,\cdots,z_j), \end{equation}$

(2.4)

and

$\begin{equation} \big(\frac{d}{dz_1}\big)^m [z_1^{b-1} \varepsilon_{(a_j),b;p}^{(c_j)}(\varpi_1 z_1^{a_1},\cdots,\varpi_j z_1^{a_j}))] = z_1^{b-m-1} \varepsilon_{(a_j),b-m;p}^{(c_j)}(\varpi_1 z_1^{a_1},\cdots,\varpi_j z_1^{a_j}), \end{equation}$

(2.5)

where $a_i, b, c_i\in\mathbb{C}; \Re(a_i) > 0, \Re(b) > 0, , p\geq0$ for $i = 1, 2, \cdots, j$ , and $\Re(b-m) > 0$ with $m\in\mathbb{N}.$

Proof. Differentiating (1.12) $m$ times with respect to $z_1, z_2, \cdots, z_j$ respectively, we get

$\begin{align*} & (\frac{d}{dz_1}\big)^m\cdots(\frac{d}{dz_j}\big)^m \varepsilon_{(a_j),b;p}^{(c_j)}(z_1,\cdots,z_j) = (\frac{d}{dz_1}\big)^m\cdots(\frac{d}{dz_j}\big)^m \sum\limits_{l_1 = l_2 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p)_{l_1}(c_2)_{l_2}\cdots(c_j)_{l_j}}{\Gamma(a_1 l_1+\cdots+a_jl_j+b)}\frac{z_1^{l_1}\cdots z_j^{l_j}}{l_1!\cdots l_j!}\\ = &\sum\limits_{l_1 = \cdots = l_j = m}^{\infty}\frac{(c_1;p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b)}\frac{l_1!\cdots l_j!\ z_1^{l_1-m}\cdots z_j^{l_j-m}}{(l_1-m)!\cdots(l_j-m)!\ l1!\cdots l_j!}\\ & = \sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p)_{l_1+m}\cdots(c_j)_{l_j+m}}{\Gamma(a_1(l_1+m)+\cdots a_j(l_j+m)+b)}\frac{z_1^{l_1}\cdots z_j^{l_j}}{l_1!\cdots l_j!}\ \ (\mbox{Replacing}\ l_i \ \mbox{by}\ l_i+m) \\ & = \sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1)_m\cdots(c_j)_m\ (c_1+m;p)_{l_1}\cdots(c_j+m)_{l_j}}{\Gamma(a_1l_1+\cdots a_jl_j+b+(a_1+\cdots+a_j)m)}\frac{z_1^{l_1}\cdots z_j^{l_j}}{l_1!\cdots l_j!}. \end{align*}$

Now using $(\lambda; \sigma)_{\mu+p} = (\lambda)_{\mu} (\lambda+\mu; \sigma)_p$ and $(\lambda)_{\mu+p} = (\lambda)_{\mu} (\lambda+\mu)_p$ , we get

$\begin{align*} & (\frac{d}{dz_1}\big)^m\cdots(\frac{d}{dz_j}\big)^m \varepsilon_{(a_j),b;p}^{(c_j)}(z_1,\cdots,c_j)\\ = &(c_1)_m \cdots(c_j)_m \sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1+m;p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots a_jl_j+b+(a_1+\cdots +a_j)m)}\frac{z_1^{l}\cdots z_j^{l_j}}{l_1!\cdots l_j!}\\ & = (c_1)_m\cdots (c_j)_{m} \ \varepsilon_{(a_j),b+(a_j)m;p}^{(c_j)+m}(z_1,z_2,\cdots, z_j), \end{align*}$

which is the desired result (2.4). Similarly, to prove (2.5), we have

$\begin{align*} &\big(\frac{d}{dz_1}\big)^m [z_1^{b-1} \varepsilon_{(a_j),b;p}^{(c_j)}(\varpi_1 z_1^{a_1},\cdots \varpi_j z_j^{a_j})]\\ = &\big(\frac{d}{dz_1}\big)^m z_1^{b-1}\sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b)}\frac{(\varpi_1 z_1^{a_1})^{l_1}\cdots(\varpi_jz_1^{a_j})^{l_j}}{l_1!\cdots l_j!}\\ = &\big(\frac{d}{dz_1}\big)^m \sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+ a_jl_j+b)} \frac{z_1^{b-1+a_1l_1+\cdots+a_jl_j}}{l_1!\cdots l_j!} \varpi_1^{l_1}\cdots\varpi_j^{l_j}\\ = &\sum\limits_{l_1 = \cdots l_j = 0}^{\infty}\frac{(c_1;p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b)} \frac{\varpi_1^{l_1}\cdots\varpi_j^{l_j}}{l_1!\cdots l_j!} \frac{(a_1l_1+\cdots+a_jl_j+b-1)!}{(a_1l_1+\cdots+a_jl_j+b-m-1)!}\ z_1^{a_1l_1+\cdots+a_jl_j+b-m-1}. \end{align*}$

Differentiating $m$ times and using the relation $l(l-1)! = l!$ , we get

$\begin{align*} &\big(\frac{d}{dz_1}\big)^m [z_1^{b-1} \varepsilon_{(a_j),b;p}^{(c_j)}(\varpi_1 z_1^{a_1},\cdots \varpi_j z_j^{a_j})]\\ = &\sum\limits_{l_1 = \cdots l_j = 0}^{\infty}\frac{(c_1;p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b)} \frac{\Gamma(a_1l_1+\cdots+a_jl_j+b)}{\Gamma(a_1l_1+\cdots+a_jl_j+b-m)}\, \varpi_1^{l_1}\cdots\varpi_j^{l_j} \, \frac{z_1^{a_1l_1+\cdots+a_jl_j+b-1-m}}{l_1!\cdots l_j!} \\ = & z_1^{b-m-1}\sum\limits_{l_1 = \cdots l_j = 0}^{\infty}\frac{(c_1;p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b-m)}\frac{(\varpi_1 z_1^{a_1})^{l_1}\cdots(\varpi_1 z_1^{a_j})^{l_j}}{l_1!\cdots l_j!}\\ = &z_1^{b-m-1}\varepsilon_{(a_j),b-m;p}^{(c_j)}(\varpi_1 z_1^{a_1},\cdots,\varpi_j z_1^{a_j}). \end{align*}$

The proof is completed.

Corollary 2.1. The generalized multivariate M-L function has the following integral representations:

$\begin{equation*} \int_{0}^{z_1}t^{b-1}\varepsilon_{(a_j),b;p}^{(c_j)}(\varpi_1 t^{a_1},\cdots,\varpi_j t^{a_j}) d t = z_1^{b}\varepsilon_{(a_j),b+1;p}^{(c_j)}(\varpi_1 z_1^{a_1},\cdots,\varpi_jz_1^{a_j}), \end{equation*}$

where $a_i, b, c_i, \varpi_i\in\mathbb{C}; \Re(a_i) > 0, \Re(b) > 0, p\geq0$ for $i = 1, 2, \cdots, j$ .

3. Fractional integral and differential formulas of generalized M-L function

In this section, we present some fractional integration and differentiation formulas of generalized M-L function given in (2.1).

Theorem 3.1. Suppose $x > r\, (r\in\mathbb{R_+} = [0, \infty))$ , $a_i$ , $b$ , $c_i$ , $\varpi\in\mathbb{C}$ , $\Re(a_i) > 0$ and $\Re(c_i) > 0$ , $\Re(b) > 0$ and $\Re(\lambda) > 0$ , then the following relations hold true:

$\begin{align} &\mathfrak{I}_{r+}^{\lambda}\left[(\varrho-r)^{b-1}\varepsilon_{(a_j),b;p}^{(c_j)}(\varpi_1(\varrho-r)^{a_1},\cdots,\varpi_j(\varrho-r)^{a_j})\right](x) \\ = &(x-r)^{\lambda+b-1} \varepsilon_{(a_i), b+\lambda;p}^{(c_i)}(\varpi_1(x-r)^{a_1},\cdots,\varpi_j(x-r)^{a_j}), \end{align}$

(3.1)

$\begin{align} &\mathfrak{D}_{r+}^{\lambda}\left[(\varrho-r)^{b-1}\varepsilon_{(a_i),b;p}^{( c_i)}(\varpi_1(\varrho-r)^{a_1},\cdots,\varpi_j(\varrho-r)^{a_j})\right](x) \\ = &(x-r)^{b-\lambda-1} \varepsilon_{(a_i),b-\lambda;p}^{(c_i)}(\varpi_1(x-r)^{a_1},\cdots,\varpi_j(x-r)^{a_j}) \end{align}$

(3.2)

and

$\begin{align} &\mathfrak{D}_{r+}^{\lambda,v}\left[(\varrho-r)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(x-r)^{a_1},\cdots,\varpi_j(x-r)^{a_j})\right](x) \\ = &(x-r)^{b-\lambda-1} \varepsilon_{(a_i),b-\lambda;p}^{(c_i)}(\varpi_1(x-r)^{a_1},\cdots,\varpi_j(x-r)^{a_j}). \end{align}$

(3.3)

Proof. Consider

$\begin{align*} &\mathfrak{I}_{r+}^{\lambda}\left[(\varrho-r)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(x-r)^{a_1},\cdots,\varpi_j(x-r)^{a_j})\right](x)\\ = &\frac{1}{\Gamma(\lambda)}\int_{r}^{x}\frac{(x-r)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(\varrho-r)^{a_1},\cdots,\varpi_1(\varrho-r)^{a_j})} {(x-\varrho)^{1-\lambda}}d\varrho\\ = &\frac{1}{\Gamma(\lambda)}\sum\limits_{n = 0}^{\infty}\frac{(c_1;p,v)_{l_1}\cdots(c_j)_{l_n}\varpi^{l_1}\cdots\varpi^{l_j}}{\Gamma(a_1 l_1+\cdots+a_jl_j+b)l_1!\cdots l_j!}\int_{r}^{x}(\varrho-r)^{b+a_1 l_1+\cdots+a_jl_j-1}(x-\varrho)^{\lambda-1}d\varrho\\ = &\sum\limits_{n = 0}^{\infty}\frac{(c_1;p,v)_{l_1}\cdots(c_j)_{l_n}\varpi^{l_1}\cdots\varpi^{l_j}}{\Gamma(a_1 l_1+\cdots+a_jl_j+b)l_1!\cdots l_j!} \left( \mathfrak{I}_{r+}^{\lambda}[(\varrho-r)^{b+a_1 l_1+\cdots+a_jl_j-1}] \right). \end{align*}$

By the use of (1.17), we have

$\begin{align*} &\mathfrak{I}_{r+}^{\lambda}\left[(\varrho-r)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(x-r)^{a_1},\cdots,\varpi_j(x-r)^{a_j})\right](x)\\ = &\sum\limits_{n = 0}^{\infty}\frac{(c_1;p,v)_{l_1}\cdots(c_j)_{l_n}\varpi^{l_1}\cdots\varpi^{l_j}}{\Gamma(a_1 l_1+\cdots+a_jl_j+b)l_1!\cdots l_j!} (x-r)^{b+\lambda+a_1l_1+\cdots+a_jl_j-1}. \frac{\Gamma(a_1l_1+\cdots+a_jl_j+b)}{\Gamma(a_1l_1+\cdots+a_jl_j+b+\lambda)}\\ = &(x-r)^{b+\lambda-1}\sum\limits_{n = 0}^{\infty}\frac{(c_1;p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b+\lambda)} \frac{[\varpi_1^{l_1}(x-r)^{a_1l_1}\cdots\varpi_j^{l_j}(x-r)^{a_jl_j}]}{l_1!\cdots l_j!}\\ = &(x-r)^{b+\lambda-1}\varepsilon_{(a_i),b+\lambda;p}^{(c_i)}(\varpi_1(x-r)^{a_1},\cdots,\varpi_j(x-r)^{a_j}), \end{align*}$

which gives the proof of (3.1).

Next, we have

$\begin{align*} &\mathfrak{D}_{r+}^{\lambda}\left[(\varrho-r)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(\varrho-r)^{a_1},\cdots,\varpi_j(\varrho-r)^{a_j})\right]\\ = &(\frac{d}{dx})^n \left\{ \mathfrak{I}_{r+}^{n-\lambda}[(\varrho-r)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(\varrho-r)^{a_1},\cdots,\varpi_j(\varrho-r)^{a_j})] \right\}, \end{align*}$

which on using (3.1) takes the following form:

$\begin{align*} &\mathfrak{D}_{r+}^{\lambda}\left[(\varrho-r)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(\varrho-r)^{a_1},\cdots,\varpi_j(\varrho-r)^{a_j})\right] \\ = &(\frac{d}{dx})^n \left\{ (x-r)^{b-\lambda+n-1}\varepsilon_{(a_i),b-\lambda+n;p}^{(c_i)}(\varpi_1(x-r)^{a_1},\cdots,\varpi_j(x-r)^{a_j}) \right\}. \end{align*}$

Applying (2.5), we get

$\begin{align*} &\mathfrak{D}_{r+}^{\lambda}\left[(\varrho-r)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(x-r)^{a_1},\cdots,\varpi_j(x-r)^{a_j})\right](x) \\ = &\left\{ (x-r)^{\eta-\lambda-1}\varepsilon_{(a_i),b-\lambda;p}^{(c_i)}(\varpi_1(x-r)^{a_1},\cdots, \varpi_j(x-r)^{a_j}) \right\}, \end{align*}$

which gives the proof of (3.2).

To obtain (3.3), we have

$\begin{align*} &\left( \mathfrak{D}_{r+}^{\lambda,v}\left[(\varrho-r)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(\varrho-r)^{a_1},\cdots,\varpi_j(\varrho-r)^{a_j})\right] \right)(x)\\ = &\left( \mathfrak{D}_{r+}^{\lambda,v}\left[\sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p,v)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j +b)}\frac{\varpi^{l_1}\cdots\varpi^{l_j}}{l_1!\cdots l_j!}(\varrho-r)^{a_1 l_1+\cdots+a_jl_j+b-1}\right] \right)(x)\\ = &\sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p,v)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j +b)}\frac{\varpi^{l_1}\cdots\varpi^{l_j}}{l_1!\cdots l_j!}\\ \times&\left( \mathfrak{D}_{r+}^{\lambda,v}[(\varrho-r)^{a_1 l_1+\cdots+a_jl_j+b-1}] \right)(x). \end{align*}$

By applying (1.18), we get

$\begin{align*} &\left( \mathfrak{D}_{r+}^{\lambda,v}[(\varrho-r)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(\varrho-r)^{a_1},\cdots,\varpi_j(\varrho-r)^{a_j})] \right)(x)\\ = &\sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p,v)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j +b)}\frac{\varpi^{l_1}\cdots\varpi^{l_j}}{l_1!\cdots l_j!}\\ \times&\frac{\Gamma(a_1l_1+\cdots+ a_jl_j+b)}{\Gamma(a_1l_1+\cdots+a_jl_j +b-\lambda)}(x-r)^{a_1l_1+ \cdots+a_jl_j+b-\lambda-1}\\ = &(x-r)^{b-\lambda-1}\sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p,v)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j +b-\lambda)}\frac{\varpi^{l_1}(x-r)^{a_1}\cdots\varpi^{l_j}(x-r)^{a_j}}{l_1!\cdots l_j!}\\ = &(x-r)^{b-\lambda-1}\varepsilon_{(a_i),b-\lambda;p}^{(c_i)}(\varpi^{l_1}(x-r)^{a_1},\cdots,\varpi^{l_j}(x-r)^{a_j}), \end{align*}$

which completes the required proof.

Remark 3.1. Applying Theorem 3.1 for $p = 0$ , then we obtain the result presented in ^[34].

4. Generalization of fractional integral and its properties

In this section, we define a fractional integral involving newly defined multivariate M-L function and discuss its properties.

Definition 4.1. Let $b, a_i, c_i, \varpi_i\in\mathbb{C}$ , $\Re(c_i) > 0$ , $\Re(a_i) > 0$ and $\Re(b) > 0$ and $h\in L(r, s)$ . Then the generalized left and right sided fractional integrals are defined by

$\begin{align} &\left(\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}h\right)(x) = \int_{r}^{x}(x-\varrho)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)} (\varpi_1(x-\varrho)^{a_1},\cdots,\varpi_j(x-\varrho)^{a_j})h(\varrho)d\varrho, (x > r) \end{align}$

(4.1)

and

$\begin{align} &\left(\mathfrak{R}_{s-;(a_i),b;p}^{(\varpi_i);(c_i)}h\right)(x) = \int_{x}^{s}(\varrho-x)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)} (\varpi_1(\varrho-x)^{a_1},\cdots,\varpi_j(\varrho-x)^{a_j})h(\varrho)d\varrho, (x < s), \end{align}$

(4.2)

respectively.

Remark 4.1. If we consider $p = 0$ , then the operators defined in (4.1) and (4.2) will take the form defined earlier by ^[34]. Similarly, if we consider $p = 0$ and $j = 1$ , then the operators defined in (4.1) and (4.2) will take the form defined by ^[22]. If we take $j = 1$ , then the work done in this paper will lead to the work presented by ^[28]. Also, if we consider one of $\varpi_i = 0$ , for $i = 1, 2, \cdots, j$ , then the operators defined in (4.1) and (4.2) will take the form of the classical operators.

Next, we prove the following properties of integral operator defined in (4.1).

Theorem 4.1. Suppose that $b, a_i, \lambda, c_i, \varpi_i\in\mathbb{C}$ , $\Re(a_i) > 0$ , $\Re(b) > 0$ , $\Re(\lambda) > 0$ , $p\geq0$ and $\Re(c_i) > 0$ for $i = 1, 2, \cdots, j$ , then the following result holds true:

$\begin{align*} \left(\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}[(\varrho-r)^{\lambda-1}] \right) (x) = (x-r)^{\lambda+b-1}\Gamma(\lambda)\varepsilon_{(a_i),b+\lambda}^{(c_i);p}(\varpi_1(x-r)^{a_1},\cdots,\varpi_j(x-r)^{a_j}). \end{align*}$

Proof. By the use of definition (4.1), we have

$\begin{equation*} \left( \mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}h \right)(x) = \int_{r}^{x}(x-\varrho)^{b-1}\varepsilon_{(a_i),b}^{(c_i)}(\varpi_1(x-\varrho)^{a_1},\cdots,\varpi_j(x-\varrho)^{a_j})h(\varrho)d\varrho. \end{equation*}$

Therefore, we get

$\begin{align*} &\left( \mathfrak{R}_{r+;(a_i), b;p}^{(\varpi_i);(c_i)}[(\varrho-r)^{\lambda-1}] \right) (x) = \int_{r}^{x}(x-\varrho)^{b-1}(\varrho-r)^{\lambda-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(x-\varrho)^{a_1},\cdots,\varpi_j(x-\varrho)^{a_j})d\varrho\\ = &\sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b))}\frac{\varpi_1^{l_1}\cdots\varpi_j^{l_j}} {l_1!\cdots l_j!} \left( \int_{r}^{x}(\varrho-r)^{\lambda-1}(x-\varrho)^{\lambda+a_1l_1+ \cdots+a_jl_j-1}d\varrho \right)\\ = &\sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b))}\frac{\varpi_1^{l_1}\cdots\varpi_j^{l_j}} {l_1!\cdots l_j!} \mathfrak{I}_{r+}^{a_1l_1+\cdots+a_jl_j+b}[(\varrho-r)^{\lambda-1}]\\ & = (x-r)^{b+\lambda-1} \sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p)_{l_1}\cdots(c_j)_{l_j}}{\Gamma(a_1l_1+\cdots+a_jl_j+b))} \frac{[\varpi_1(x-r)^{a_1l_1}\cdots\varpi_j(x-r)^{a_jl_j}]}{l_1!\cdots l_j!}\\ \times& \frac{\Gamma(\lambda)\Gamma(a_1l_1+\cdots+a_jl_j +b)}{\Gamma(a_1l_1+\cdots+a_jl_j+b+\lambda)}\\ = &(x-r)^{b+\lambda-1}\Gamma(\lambda)\varepsilon_{(a_i),b+\lambda;p}^{(c_i)}(\varpi_1(x-r)^{a_1},\cdots,\varpi_j(x-r)^{a_j}), \end{align*}$

which gives the desired proof.

Theorem 4.2. Suppose that $c_i, a_i, b, \varpi_i \in\mathbb{C}$ , $\Re(a_i) > 0$ , $\Re(b) > 0,$ $p\geq0$ for $i = 1, 2, \cdots, j$ , then the following result holds true:

$\begin{equation*} \|\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}\Phi\|_1\leq K\|\Phi\|_1. \end{equation*}$

Where

$\begin{align*} &K: = (s-r)^{\mathfrak{Re}(b)}\sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{|(c_1;p)_{l_1}\cdots(c_j)_{l_j}|}{\Gamma(a_1l_1+\cdots+a_jl_j +b)(\Re(b)+\Re(a_1)l_1+\cdots+\Re(a_j)l_j)}\\ \times&\frac{|\varpi_1^{l_1}(s-r)^{a_1l_1}\cdots\varpi_j^{l_j}(s-r)^{a_jl_j}|}{l_1!\cdots l_j!}. \end{align*}$

Proof. By the use of (1.13) and (4.1) and by interchanging integration and summation order, we have

$\begin{align*} \|\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i) }\Phi\|_1& = \int\limits_{r}^{s}|\int_{r}^{x}(x-\varrho)^{b-1} \varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(x-\varrho)^{a_1},\cdots,\varpi_j(x-\varrho)^{a_j}) \Phi(\varrho)d\varrho|dx\\ \leq&\int_{r}^{s}\left[ \int_{\varrho}^{x} (x-\varrho)^{\Re(b)-1}|\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(x-\varrho)^{a_1},\cdots,\varpi_j(x-\varrho)^{a_j})|dx \right]|\Phi(\varrho)|d\varrho\\ = &\int_{r}^{s}\left[ \int_{0}^{x-\varrho} u^{\Re(b)-1}|\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1 u^{a_1},\cdots,\varpi_j u^{a_j})|du \right]|\Phi(\varrho)|d\varrho, \end{align*}$

by setting $u = x-\varrho.$ After simplification, we obtain

$\begin{align*} &\|\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}\Phi\|_1\\ \leq&\int_{r}^{s} \Bigr[ \sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{|(c_1;p)_{l_1}\cdots(c_j)_{l_j}|}{\Gamma(a_1l_1+\cdots+a_jl_j+b)}\frac{|\varpi_1^{a_1}\cdots\varpi_j^{l_j}|} {l_1!\cdots l_j!}\\ \times& \left( \frac{(u)^{\Re(b)+\Re(a_1)l_1+\cdots+\Re(a_j)l_j}}{(\Re(b)+\Re(a_1)l_1+\cdots+\Re(a_j)l_j)} \right)_{0}^{s-r} \Bigr] |\Phi(\varrho)|d\varrho. \end{align*}$

It follows that

$\begin{align*} &\|\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}\Phi\|_1\leq \Bigr\{ (s-r)^{\Re(b)} \sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{|(c_1;p)_{l_1}\cdots(c_j)_{l_j}|}{\Gamma(a_1l_1+\cdots+a_jl_j +b)(\Re(b)+\Re(a_1)l_1+\cdots+\Re(a_j)l_j)}\\ \times&\frac{|\varpi_1^{l_1}(s-r)^{a_1l_1}\cdots\varpi_j^{l_j}(s-r)^{a_jl_j}|}{l_1!\cdots l_j!} \Bigr\} \int\limits_{r}^{s}|\Phi(\varrho)|d\varrho\\ = &K||\Phi||_1, \end{align*}$

where

$\begin{align*} &K = (s-r)^{\mathfrak{Re}(b)}\sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{|(c_1;p)_{l_1}\cdots(c_j)_{l_j}|}{\Gamma(a_1l_1+\cdots+a_jl_j +b)(\Re(b)+\Re(a_1)l_1+\cdots+\Re(a_j)l_j)}\\ \times&\frac{|\varpi_1^{l_1}(s-r)^{a_1l_1}\cdots\varpi_j^{l_j}(s-r)^{a_jl_j}|}{l_1!\cdots l_j!}, \end{align*}$

which gives the desired result.

Corollary 4.1. If we take $a_i, b, c_i, \varpi_i \in \mathbb{C}$ , $\Re(a_i) > 0$ , $\Re(b) > 0$ , $\Re(c_i) > 0$ with $i = 1, 2, \cdots, j$ , then the following result holds true:

$\begin{align*} \left(\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}1\right)(x) = (x-r)^{b}\varepsilon_{(a_i),b+1;p}^{(c_i)}(\varpi_1(x-r)^{a_1},\cdots,\varpi_j(x-r)^{a_j}). \end{align*}$

Proof. Consider

$\begin{align*} &\left(\mathfrak{R}_{r+;(a_i),b}^{(\varpi_i);(c_i)}1\right)(x) = \int_{r}^x(x-\varrho)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(x-\varrho)^{a_1},\cdots,\varpi_j(x-r)^{a_j})d\varrho\\ = &\int_{r}^x(x-\varrho)^{b-1}\left(\sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p)_{l_1}\cdots(c_j)_{l_j}\varpi_1^{l_1}(x-\varrho)^{a_1l_1}\cdots\varpi_j^{l_j}(x-\varrho)^{a_jl_j}} {\Gamma(a_1l_1+\cdots+a_jl_j+b)l_1!\cdots l_j!}\right)d\varrho. \end{align*}$

It follows that

$\begin{align*} \left(\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}1\right)(x) & = \sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p)_{l_1}\cdots(c_j)_{l_j}\varpi_1^{l_1}\cdots\varpi_j^{l_j}} {\Gamma(a_1l_1+\cdots+a_jl_j+b)l_1!\cdots l_j!} \int_{r}^x(x-\varrho)^{b+a_1l_1+\cdots+a_jl_j-1}d\varrho\\ = &(x-r)^{b}\sum\limits_{l_1 = \cdots = l_j = 0}^{\infty}\frac{(c_1;p)_{l_1}\cdots(c_j)_{l_j}\varpi_1^{l_1}(x-r)^{a_1l_1}\cdots\varpi_j^{l_j}(x-r)^{a_jl_j}} {\Gamma(a_1l_1+\cdots+a_jl_j+b)(a_1l_1+\cdots+a_jl_j+b)l_1!\cdots l_j!}\\ = &(x-r)^{b}\,\varepsilon_{(a_i),b+1;p}^{(c_i)}(\varpi_1(x-r)^{a_1},\cdots,\varpi_j(x-r)^{a_j}), \end{align*}$

which gives the desired proof.

Theorem 4.3. The generalized fractional operator can be represented in term of Riemann–Liouville fractional integrals for $c_i$ , $a_i$ , $b$ , $\varpi_i\in\mathbb{C}$ with $\Re(a_i) > 0$ , $\Re(b) > 0$ , $\Re(c_i) > 0$ for $i = 1, 2, \cdots, j$ , $p\geq0$ and $x > r$ as follows:

$\begin{align*} \left(\mathfrak{R}_{r+;(a_i),b}^{(\varpi_i);(c_i)}h\right)(x) = \sum\limits_{l_1 = \cdots = l_j = 0}^\infty \frac{\Gamma(c_1+l_1;p)(c_2)_{l_2}\cdots(c_j)_{l_j}\varpi_1^{a_1}\cdots\varpi_j^{a_j}} {\Gamma(c_1)l_1!\cdots l_j!}\mathfrak{I}_{r+}^{a_1l_1+\cdots+a_jl_j +b}h(x). \end{align*}$

Proof. By utilizing (2.1) in (4.1) and then interchanging the order of summation and integration, we have

$\begin{align*} \left(\mathfrak{R}_{r+;(a_i),b}^{(\varpi_i);(c_i)}h\right)(x)& = \int_{r}^{x}(x-\varrho)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)} (\varpi_1(x-\varrho)^{a_1},\cdots,\varpi_j(x-\varrho)^{a_j})h(\varrho)d\varrho\\ = &\int_{r}^{x}(x-\varrho)^{b-1}\sum\limits_{l_1 = \cdots = l_j = 0}^\infty\frac{\Gamma(c_1+l_1;p)(c_2)_{l_2}\cdots(c_j)_{l_j}\varpi_1^{l_1} (x-\varrho)^{a_1l_1}\cdots\varpi_j^{l_j} (x-\varrho)^{a_jl_j}} {\Gamma(c_1)\Gamma(a_1l_1+\cdots+a_jl_j+b)l_1!\cdots l_j!}h(\varrho)d\varrho\\ = &\sum\limits_{l_1 = \cdots = l_j = 0}^\infty\frac{\Gamma(c_1+l_1;p)(c_2)_{l_2}\cdots(c_j)_{l_j}\varpi_1^{a_1l_1}\cdots\varpi_j^{a_jl_j} }{\Gamma(c_1)l_1!\cdots l_j!} \frac{1}{\Gamma(a_1l_1+\cdots+a_jl_j +b)}\\ \times&\int_{r}^{x}(x-\varrho)^{a_1l_1+\cdots+a_jl_j+b-1}h(\varrho)d\varrho\\ = &\sum\limits_{l_1 = \cdots = l_j = 0}^\infty\frac{\Gamma(c_1+l_1;p)(c_2)_{l_2}\cdots(c_j)_{l_j}\varpi_1^{a_1l_1}\cdots\varpi_j^{a_jl_j} }{\Gamma(c_1)l_1!\cdots l_j!}\mathfrak{I}_{r+}^{a_1l_1 +\cdots+a_jl_j+b}h(x), \end{align*}$

which gives the desired proof.

Theorem 4.4. For $\lambda$ , $c_i$ , $a_i$ , $b$ , $\varpi_i\in\mathbb{C}$ with $\Re(a_i) > 0$ , $\Re(b) > 0$ , $\Re(c_i) > 0$ , $\Re(\lambda) > 0$ , for $i = 1, 2, \cdots, j$ , $p\geq0$ and $x > r$ , then the following result holds true:

$\begin{equation} (\mathfrak{I}_{r+}^{\lambda}[\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}h])(x) = (\mathfrak{R}_{r+;(a_i),b+\lambda}^{(\varpi_i);(c_i)}h)(x) = (\mathfrak{R}_{r+;(a_i),b}^{(\varpi_i);(c_i)}[\mathfrak{I}_{r+}^{\lambda}h])(x), \end{equation}$

(4.3)

where $h\in L(r, s)$ .

Proof. By employing (1.14) and (4.1), we have

$\begin{align*} &(\mathfrak{I}_{r+}^{\lambda}[\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}h])(x)\\ = &\frac{1}{\Gamma(\lambda)}\int_{r}^{x}\frac{[(\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}h)(\varrho)]}{(x-\varrho) ^{1-\lambda}}d\varrho\\ = &\frac{1}{\Gamma(\lambda)}\int_{r}^{x}(x-\varrho)^{\lambda-1} \left[ \int_{r}^{\varrho} (\varrho-u)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(\varrho-u)^{a_1},\cdots,\varpi_j(\varrho-u)^{a_j})h(u)du \right]d\varrho. \end{align*}$

It follows that

$\begin{align*} &(\mathfrak{I}_{r+}^{\lambda}[\mathfrak{R}_{r+;(a_i),b}^{(\varpi_i);(c_i)}h])(x) \\ = &\int_{r}^{x} \left[ \frac{1}{\Gamma(\lambda)}\int_{u}^{x} (x-\varrho)^{\lambda-1}(\varrho-u)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i) }(\varpi_1(\varrho-u)^{a_1},\cdots,\varpi_j(\varrho-u)^{a_j})d\varrho \right]h(u)du. \end{align*}$

By considering $\varrho-u = \theta$ , we get

$\begin{align*} \begin{aligned} &(\mathfrak{I}_{r+}^{\lambda}[\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}h])(x)\\ = &\int_{r}^{x} \left[ \frac{1}{\Gamma(\lambda)}\int_{0}^{x-u} (x-u-\theta)^{\lambda-1}\theta^{b-1}\varepsilon_{(a_i),b;p}^{(c_i) } (\varpi_1\theta^{a_1},\cdots,\varpi_j\theta^{a_j})d\theta \right]h(u)du\\ = &\int_{r}^{x} \left[ \frac{1}{\Gamma(\lambda)}\int_{0}^{x-u} \frac{\theta^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1\theta^{a_1},\cdots,\varpi_j\theta^{a_j})}{(x-u-\theta)^{1-\lambda}}d\theta \right]h(u)du.\end{aligned} \end{align*}$

Hence, from (1.14) and applying (3.1), we obtain

$\begin{align*} &(\mathfrak{I}_{r+}^{\lambda}[\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}h])(x)\\ = &\int_{r}^{x}\left[ \theta^{\lambda+b-1}\varepsilon_{(a_i),b+\lambda;p}^{(c_i)}(\varpi_1\theta^{a_1},\cdots,\varpi_j\theta^{a_j}) \right]h(u)du\\ = &\int_{r}^{x} (x-u)^{\lambda+b-1}\varepsilon_{(a_i),b+\lambda}^{(c_i)}(\varpi_1(x-u)^{a_1},\cdots,\varpi_j(x-u)^{a_j})h(u)du. \end{align*}$

Thus, we have

$\begin{equation} (\mathfrak{I}_{r+}^{\lambda}[\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}h])(x) = \left( \mathfrak{R}_{r+;(a_i),b+\lambda}^{(\varpi_i);(c_i)}h \right)(x). \end{equation}$

(4.4)

Next, consider the right hand side of (4.3) and employing (4.1) to derive the second part, we have

$\begin{align*} &(\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}[\mathfrak{I}_{r+}^{\lambda}h])(x)\\ = &\int_{r}^{x}(x-\varrho)^{b-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(x-\varrho)^{a_1},\cdots,\varpi_j(x-\varrho)^{a_j}) [\mathfrak{I}_{r+}^{\lambda}h] (\varrho)d\varrho\\ = &\int_{r}^{x}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(x-\varrho)^{a_1},\cdots,\varpi_j(x-\varrho)^{a_j})\left( \frac{1}{\Gamma(\lambda)}\int\limits_{r}^{\varrho}\frac{h(u)}{(\varrho-u)^{1-\lambda}}du \right)d\varrho. \end{align*}$

It follows that

$\begin{align*} &(\mathfrak{R}_{r+;(a_i),b}^{(\varpi_i);(c_i)}[\mathfrak{I}_{r+}^{\lambda}h])(x)\\ = &\int\limits_{r}^{x}\frac{1}{\Gamma(\lambda)} \left[ \int_{u}^{x}(x-\varrho)^{b-1}(\varrho-u)^{\lambda-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1(x-\varrho)^{a_1},\cdots,\varpi_j(x-\varrho)^{a_j})d\varrho \right] h(u)du. \end{align*}$

By setting $x-\varrho = \theta$ , we get

$\begin{align*} &(\mathfrak{R}_{r+;(a_i),b}^{(\varpi_i);(c_i)}[\mathfrak{I}_{r+}^{\lambda}h])(x)\\ = &\int_{r}^{x}\frac{1}{\Gamma(\lambda)} \left[ \int_{x-u}^{0}\theta^{b-1}(x-\theta-u)^{\lambda-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1\theta^{a_1},\cdots,\varpi_j\theta^{a_j}) (-d\theta) \right] h(u)du\\ = &\int\limits_{r}^{x}\frac{1}{\Gamma(\lambda)} \left[ \int_{0}^{x-u}\theta^{b-1}(x-\theta-u)^{\lambda-1}\varepsilon_{(a_i),b;p}^{(c_i)}(\varpi_1\theta^{a_1},\cdots,\varpi_j\theta^{a_j})d\theta \right] h(u)du. \end{align*}$

Further, by using (1.14) and applying (3.1), we obtain

$\begin{equation} (\mathfrak{R}_{r+;(a_i),b;p}^{(\varpi_i);(c_i)}[\mathfrak{I}_{r+}^{\lambda}h])(x) = \left( \mathfrak{R}_{r+;(a_i),b+\lambda}^{(\varpi_i);(c_i)}h \right)(x). \end{equation}$

(4.5)

Thus, (4.4) and (4.5) gives the desired proof.

5. Conclusions

Nowadays, the theories are developed very rapidly. The scientists are introducing more advanced and generalized forms of the classical ones. In this present study, we introduced a generalized form of the multivariate M-L function (2.1) by employing the generalized Pochhammer symbol and its properties. By using this more extended form of M-L, we introduced a fractional integral operator and studied some of the basic properties of this operator. The special cases of the main results are if we take $p = 0$ , then the operators defined in (4.1) and (4.2) will reduce to the work done by ^[34]. Similarly, if we take $j = 1$ and $p = 0$ , then the operators defined in (4.1) and (4.2) will lead to the work done by ^[22]. If we take $j = 1$ , then the work done in this paper will lead to the work presented by ^[28]. Moreover, if we consider one of $\varpi_i = 0$ , for $i = 1, 2, \cdots, j$ , then the operators defined in (4.1) and (4.2) will reduce to the classical R-L operators. We believe that our proposed operator will be more applicable in the fields of fractional integral inequalities and operator theory.

Acknowledgments

The author T. Abdeljawad would like to thank Prince Sultan University for supporting through TAS research lab. Manar A. Alqudah: Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R14), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflict of interest

The authors declare 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)