Next Article in Journal
Using Deep Learning to Detect the Need for Forest Thinning: Application to the Lungau Region, Austria
Next Article in Special Issue
Optimal Confidence Regions for Weibull Parameters and Quantiles under Progressive Censoring
Previous Article in Journal
An Optimization Precise Model of Stroke Data to Improve Stroke Prediction
Previous Article in Special Issue
Using Epidemiological Models to Predict the Spread of Information on Twitter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Mathematical Study on a Fractional-Order SEIR Mpox Model: Analysis and Vaccination Influence

1
Department of Mathematics, Al Zaytoonah University, Amman 11733, Jordan
2
Nonlinear Dynamics Research Center (NDRC), Ajman University, Ajman 346, United Arab Emirates
3
Faculty of Computer Studies, Arab Open University, Riyadh 11681, Saudi Arabia
4
School of Mathematical Sciences, Universiti Sains Malaysia, Penang 11800, Malaysia
5
Department of Mathematics and Statistics, American University of Sharjah, Sharjah 26666, United Arab Emirates
*
Author to whom correspondence should be addressed.
Algorithms 2023, 16(9), 418; https://doi.org/10.3390/a16090418
Submission received: 24 July 2023 / Revised: 29 August 2023 / Accepted: 30 August 2023 / Published: 1 September 2023
(This article belongs to the Special Issue Mathematical Models and Their Applications IV)

Abstract

:
This paper establishes a novel fractional-order version of a recently expanded form of the Susceptible-Exposed-Infectious-Recovery (SEIR) Mpox model. This model is investigated by means of demonstrating some significant findings connected with the stability analysis and the vaccination impact, as well. In particular, we analyze the fractional-order Mpox model in terms of its invariant region, boundedness of solution, equilibria, basic reproductive number, and its elasticity. In accordance with an effective vaccine, we study the progression and dynamics of the Mpox disease in compliance with various scenarios of the vaccination ratio through the proposed fractional-order Mpox model. Accordingly, several numerical findings of the proposed model are depicted with the use of two numerical methods; the Fractional Euler Method (FEM) and Modified Fractional Euler Method (MFEM). Such findings demonstrate the influence of the fractional-order values coupled with the vaccination rate on the dynamics of the established disease model.

1. Introduction

Mathematical modeling has turned into an invaluable tool in revealing the behaviors of infectious diseases and formulating effective control strategies. Mpox, a viral disease with similarities to smallpox, has attracted considerable attention in the field of epidemiological models due to its potential impact on human health. As the Mpox disease is a viral zoonotic infection, it might transmit from animals to people. In addition, it might transmit from the environment to a person, as it can transmit from one person to another [1]. In 2022, the Mpox disease began to spread like a pandemic in numerous countries all over the world. It is spreading now in some African countries [2,3].
To gain deeper insights into the transmission dynamics and develop accurate predictive models for Mpox, many researchers have increasingly explored the utilization of fractional calculus. Fractional calculus offers a promising tool for highlighting the intricate characteristics of infectious disease dynamics, as it allows for the incorporation of memory-dependent and non-local effects. Traditional models often suppose an instantaneous mixing by assuming that past interactions have no influence on future outcomes. So, from this point, the notion of fractional calculus might be utilized as an alternative to these models, as it can enable one to formulate a model of long-range dependence with extra freedom of fractional-order values. This would enable a more realistic representation of the spread and evolution of the disease; see [4].
In the context of Mpox, fractional calculus can potentially enhance our understanding of various epidemiological factors such as transmission rates, recovery rates, and the impact of interventions. By incorporating fractional derivatives and integrals into the modeling framework, researchers can address the non-local nature of interactions and the presence of memory effects in disease transmission. These effects can reflect the persistence of immunity, delayed responses to interventions, or even other temporal dependencies; see [5]. Developing a fractional calculus-based model for Mpox involves considering fractional differential equations that can describe the growth of the diseased population over time. The fractional-order derivatives correspond to the memory or long-range dependence, which might be observed in the transmission of disease. Accurate estimation of the fractional-order value requires a combination of empirical evidence, fitting the model to available data, and a deep understanding of the underlying epidemiological processes [6]. The utilization of fractional calculus in modeling Mpox has the potential to yield valuable insights into the disease’s dynamics and control strategies. Incorporating memory effects and non-local interactions in modeling the disease can provide a more accurate representation of its spread and evolution over time. Furthermore, such modeling can assist public health officials in evaluating and optimizing intervention strategies, estimating the impact of vaccination programs, and assessing the effectiveness of different control measures [7].
In this study, we aim to explore the application of fractional calculus in modeling Mpox transmission dynamics. By introducing fractional-order derivatives and considering memory-dependent effects, we attempt to enhance our understanding of the disease’s behavior and provide valuable tools for policymakers and health care professionals in their efforts to mitigate and control Mpox outbreaks. For this purpose, a novel fractional-order version of the SEIR model is established in light of operating the Caputo operator. This model is analyzed with the help of proposing certain novel stability results, including the investigation of its invariant region, equilibria, Basic Reproductive Number (BRN), elasticity, and boundedness of solution, which is observed by using the Adomian decomposition method. Additionally, we also solved this model numerically using the Fractional Euler Method (FEM) and Modified Fractional Euler Method (MFEM) [8] for the aim of revealing the influence of the fractional-order values and the vaccination rate on the dynamics of the model at hand. In summary, this paper is arranged in the following manner. In the next part, some key definitions and basic concepts associated with fractional calculus are recalled briefly, whereas the formulation of the Mpox illness model is proposed in Section 3. In Section 4, we study the stability analysis of the Mpox model. In Section 5, we solve the proposed fractional-order Mpox model using the FEM and MFEM, followed by the conclusion of this work.

2. Preliminary Concepts

In this part, we recall fundamental descriptions and theorems related to fractional calculus, including the Riemann–Liouville integral and derivative, the Caputo derivative, and other relevant concepts [9].
Definition 1.
The Riemann–Liouville integral operator of the function h ( x ) of fractional-order 0 < ρ 1 can be expressed by
I ρ h ( x ) = 1 Γ ( ρ ) 0 x h ( s ) ( x s ) ρ 1 d s , x > 0 .
In the following content, we list below some properties of the Riemann–Liouville integral operator:
I ρ h ( x ) h ( x ) as ρ 0 ,
I ρ ( x a ) ν = Γ ( ν + 1 ) Γ ( ρ + ν + 1 ) ( x a ) ρ + ν , ν 1 , a R ,
I ρ I λ h ( x ) = I λ I ρ h ( x ) , ρ , λ 0 ,
I ρ I λ h ( x ) = I ρ + λ h ( x ) , ρ , λ 0 .
Definition 2.
The Caputo fractional derivative operator of the function h ( x ) of fractional-order 0 < ρ 1 can be expressed by
C D * ρ h ( x ) = 1 Γ ( 1 ρ ) 0 x h ( s ) ( x s ) ρ d s , x > 0 .
Herein, we also list in what follows, some further properties of the Caputo derivative operator:
  • For a constant c, we have C D * ρ c = 0 .
  • For a R , we have
    C D * ρ ( x a ) λ = Γ ( λ + 1 ) Γ ( λ ρ + 1 ) ( x a ) λ ρ , λ > ρ 1 , 0 , o t h e r w i s e .
  • For constants a and b, C D * ρ is a linear operator, i.e.,
    C D * ρ ( a h ( x ) + b g ( x ) ) = a C D * ρ h ( x ) + b C D * ρ g ( x ) .
In the same regard, it is important to remember the following basic property:
I ρ C D * ρ h ( x ) = h ( x ) k = 1 n h k ( 0 + ) x k k ! , x > 0 ,
where n 1 < ρ n such that n N .
Definition 3.
The Caputo operator C D * ρ might be outlined in terms of the Riemann–Liouville integral operator along the following lines:
C D * ρ f = I n ρ D n f ,
where ρ R + and n = ρ .
Definition 4
([9]). The two parameters of the Mittag-Leffler function can be expressed as
E γ , δ ( x ) = j = 0 x k Γ ( γ k + δ ) ,
where γ , δ > 0 and x C .
Theorem 1
([9]). Assume C D * j ρ h ( x ) C ( 0 , T ] , for j = 0 , 1 , 2 , , n + 1 and 0 < ρ 1 . Then, the function h ( x ) might be expanded about x 0 in the following manner:
h ( x ) = j = 0 n ( x x 0 ) j ρ Γ ( j ρ + 1 ) ( C D * j ρ h ) ( x 0 ) + ( x x 0 ) ( n + 1 ) ρ Γ ( ( n + 1 ) ρ + 1 ) ( C D * ( n + 1 ) ρ h ) ( ξ ) ,
x ( 0 , T ] , where 0 < ξ < x .

3. Mpox Model

Mathematical models are extensively engaged in simulating the dynamics and the transmission of infectious diseases. In accordance with an effective vaccine, we study the progression and dynamics of the Mpox disease in compliance with various scenarios of the vaccination rate through a compartmental model of the type SEIR. In this paper, we consider the intervention of vaccination playing a significant role as a parameter, impacting on the model’s dynamics. The equilibrium of the disease-free state will be determined, and then its stability with computing the BRN R 0 * will be examined in Section 4.
The diagram of the compartment model, including the vaccination rate for the aimed Mpox model, is shown in Figure 1, in which the total of the susceptible population is denoted by S, the total of the exposed population is denoted by E, the total of the infected population is denoted by I, and, finally, the total of the recovery population by R.
In the same context, Figure 1 exhibits the parameters ( Λ , σ , ν , μ , α , γ , β ) that are responsible for the dynamics of the considered model. These parameters are outlined in Table 1.
It is noteworthy that the rate of transmission β indicates the transmission’s contact from S to E as a result of interacting with I cases. Additionally, γ 1 indicates the period of incubation, whereas α 1 indicates the period of recovery. In light of the aforementioned considerations, the conventional Mpox model might be described in the following manner [10]:
d S d t = Λ β S ( t ) I ( t ) ( ν + μ ) S ( t ) d E d t = ( μ + γ ) E ( t ) + β S ( t ) I ( t ) d I d t = γ E ( t ) ( μ + σ + α ) I ( t ) d R d t = ν S ( t ) + α I ( t ) μ R ( t ) .
with initial conditions
S ( 0 ) = S 0 , E ( 0 ) = E 0 , I ( 0 ) = I 0 , R ( 0 ) = R 0 .
For more clarification, model (10) is established in accordance with particular conditions, which might be stated as follows:
  • The contact transmission parameter β does not take into account some factors, like climate, age, marital status, or gender.
  • All of the model’s parameters have no negative values.
  • Inter-individual relations within society are distinguished by uniformity and homogeneity. This is because of the inherited suppositions for all of the model’s states in which all persons have the same parameter value of the contact transmission regardless of their situations, such as their health circumstances or age.
  • The size of population N is considered constant at any time, and it meets the following equality:
    S ( t ) + E ( t ) + I ( t ) + R ( t ) = N ( t ) .
  • In the proposed model, the demography is taken into account, in which the natural death μ and natural birth Λ are incorporated.
  • System (10) satisfies the following property:
    d N d t h ( N ) ,
    where h is a function depending only on N [11]. To see this, one can observe
    d N d t = d d t ( S ( t ) + E ( t ) + I ( t ) + R ( t ) ) = Λ μ N σ I .
    Now, due to σ I 0 , we obtain
    d N d t Λ μ N .
    With the use of separation of variables, we can have
    ( Λ μ N ) c e t ,
    or
    N Λ μ + c e t .
    By using N ( 0 ) = N 0 , we obtain
    N ( t ) Λ μ Λ μ N 0 e μ t .
    Now, as t , we obtain
    N ( t ) Λ μ ,
    and, hence, property (13) is held for system (10).
In order to propose the aimed fractional-order Mpox model, the Caputo operator C D * ρ is operated on the modified model (10) to obtain
C D * ρ S ( t ) = Λ β S ( t ) I ( t ) ( ν + μ ) S ( t ) C D * ρ E ( t ) = β S ( t ) I ( t ) ( μ + γ ) E ( t ) C D * ρ I ( t ) = γ E ( t ) ( μ + σ + α ) I ( t ) C D * ρ R ( t ) = ν S ( t ) + α I ( t ) μ R ( t ) ,
with initial conditions (11).

4. Stability Analysis

Herein, we aim to examine the primary properties of model (19) in view of various aspects including the invariant region, boundedness of solution, equilibrium point, BRN, and its elasticity.

4.1. The Invariant Region

This subsection intends to study the non-negativity of the solution of model (19). For instance, the solution of model (19) will be validated in accordance with the positive values of the initial conditions. For this purpose, we suppose for all non-negative values of t that the domain of the desired solution ψ has the following form:
ψ = S ( t ) , E ( t ) , I ( t ) , R ( t ) R + 4 : S ( t ) + E ( t ) + I ( t ) + R ( t ) Λ μ ,
where S ( t ) , E ( t ) , I ( t ) , R ( t ) 0 . One could observe that all model parameters have non-negative values (i.e., σ , ν , μ , α , γ , β , Λ , 0 ). Now, in order to show the non-negativity of the solution of model (19), we recall the following two results, which would pave the way to achieve our first target.
Lemma 1
([12,13]). (Generalized Mean Value Theorem). Suppose that Θ ( t ) C [ a 1 , a 2 ] and C D * ρ Θ ( t ) C ( a 1 , a 2 ] , where 0 < ρ 1 , then we have:
Θ ( t ) = Θ ( a 1 ) + 1 Γ ( ρ ) C D * ρ Θ ( ζ ) ( t a 1 ) ρ ,
where 0 ζ t , t ( a 1 , a 2 ] .
Remark 1
([12,13]). Suppose that Θ ( t ) C [ 0 , a 2 ] and C D * ρ Θ ( t ) C ( 0 , a 2 ] , where 0 < ρ 1 . In view of the above Lemma, one can conclude that, if C D * ρ Θ ( t ) 0 , t ( 0 , a 2 ] , then Θ ( t ) will be a non-decreasing vector-valued function, while if C D * ρ Θ ( t ) 0 , t ( 0 , a 2 ] , then such a function will be a non-increasing one.
Theorem 2.
There exists a unique solution Θ ( t ) = [ S ( t ) , E ( t ) , I ( t ) , R ( t ) ] T of system (19), subject to the initial condition Θ ( 0 ) = [ S ( 0 ) , E ( 0 ) , I ( 0 ) , R ( 0 ) ] T R + 4 , and this solution will remain in R + 4 .
Proof. 
Tracking the necessary results deduced by the author in [13] leads us to identify, for definite, a unique solution Θ ( t ) on ( 0 , ) for system (19), subject to the initial condition Θ ( 0 ) . On the other hand, in order to show that the domain R + 4 is indeed the positively invariant region, we have to observe the following assertions:
C D * ρ S ( t ) | S = 0 = Λ β S ( t ) I ( t ) ( ν + μ ) S ( t ) | S = 0 = Λ 0 , C D * ρ E ( t ) | E = 0 = β S ( t ) I ( t ) ( μ + γ ) E ( t ) | E = 0 = β S ( t ) I ( t ) 0 , C D * ρ I ( t ) | I = 0 = γ E ( t ) ( μ + σ + α ) I ( t ) | I = 0 = γ E ( t ) 0 , C D * ρ R ( t ) | R = 0 = ν S ( t ) + α I ( t ) μ R ( t ) | R = 0 = ν S ( t ) + α I ( t ) 0 .
This validates that the solution of system (19) is non-negative over the domain given in (20). □

4.2. Boundedness of Solution

It was reported in [14] that the boundedness of the solution is regarded as a very important trait of this type of population model, with infectious diseases. In the epidemic population models, the solutions must be bounded, as unbounded solutions are physically meaningless. In what follows, we propose to show that the solution given at any specific time for system (19) is bounded by Λ μ , as assumed previously in (20). For this purpose, we should note that
C D * ρ N = C D * ρ S + C D * ρ E + C D * ρ I + C D * ρ R = Λ μ ( S + E + I + R ) σ I ,
where N is given in (12). Accordingly, we can have
C D * ρ N = Λ μ N σ I .
Under the assumed conditions in the domain ψ reported in (20), one could observe that σ I 0 , which yields the following inequality:
C D * ρ N ( t ) Λ μ N ( t ) .
In order to solve the above inequality analytically with the initial conditions
N ( 0 ) = N 0 ,
we use the so-called Adomian decomposition method (or simply ADM), as it is very suitable in this case. To see this, we first operate I ρ on both sides of the inequality (23) to obtain
N ( t ) N 0 + I ρ Λ μ I ρ N ( t ) ,
or
N ( t ) N 0 + Λ Γ ( ρ + 1 ) t ρ μ I ρ N ( t ) ,
Now, in view of the ADM, the general solution of inequality (25) might be supposed as N ( t ) = k = 0 N k ( t ) . This would immediately yield
k = 0 N k ( t ) N 0 + Λ Γ ( ρ + 1 ) t ρ μ I ρ k = 0 N k ( t ) ,
which consequently implies
N 0 ( t ) = N 0 + Λ Γ ( ρ + 1 ) t ρ , N k ( t ) μ I ρ k = 0 N k ( t ) , k 1 .
Thus, based on (26), we can obtain, for instance, N 1 ( t ) as follows:
N 1 ( t ) μ N 0 Γ ( ρ + 1 ) t ρ + Λ Γ ( 2 ρ + 1 ) t 2 ρ .
In the same way, we can obtain
N 2 ( t ) μ 2 N 0 Γ ( 2 ρ + 1 ) t 2 ρ + Λ Γ ( 3 ρ + 1 ) t 3 ρ .
If we continue in this manner, we obtain
N 3 ( t ) μ 3 N 0 Γ ( 3 ρ + 1 ) t 3 ρ + Λ Γ ( 4 ρ + 1 ) t 4 ρ .
Now, due to the solution having the form N ( t ) = k = 0 N k ( t ) , then, with the help of using the Mittag-Leffler function, we can gain
N ( t ) Λ μ N 0 E ρ , 1 μ t ρ .
Consequently, with the use of our assumptions in (20), we can deduce
N ( t ) Λ μ Λ μ N 0 E ρ , 1 μ t ρ .
Therefore, the solution of system (19) is bounded by the term Λ μ as t , and hence this system is well-posed over its domain ψ .

4.3. The Disease-Free Equilibrium (DFE)

The DFE is the case in which society is clear of the Mpox illness. Such a point is P 0 = ( S 0 , E 0 , I 0 , R 0 ) , where S 0 = S ( 0 ) is the initial susceptible state, E 0 = E ( 0 ) is the initial exposed state, I 0 = I ( 0 ) is the initial infected state, and R 0 = R ( 0 ) is the initial recovered state. Thus, finding a solution for model (19), while considering E 0 = 0 and I 0 = 0 , leads us to the following DFE:
P 0 = ( S 0 , E 0 , I 0 , R 0 ) = Λ μ + ν , 0 , 0 , ν Λ μ ( μ + ν ) ,
where S 0 = Λ ( μ + ν ) and R 0 = ν Λ μ ( μ + ν ) .

4.4. The BRN R 0 *

The BRN is regarded a very significant index used to indicate the harshness of the infectious illness. It can be employed to describe the occurrence of the pandemic and the potential of outbreak. This index is regarded as an epidemiological measure to determine the hazard of the infectious transmission. With the aim of computing the BRN, we use the next-generation matrix method (NGMM), which was suggested in [15] and afterward developed in [16]. Such a method has been broadly employed for evaluating the BRN in a number of epidemiological models. Additionally, for the purpose of finding the BRN, we employ the NGMM by taking into account the exposed and infected classes. Afterward, we break down d Y d t = Q G , in which Q = d E d t and G = d I d t . This, consequently, gives
d Y d t = d E d t d I d t , Q = β S ( t ) I ( t ) 0 , G = ( μ + γ ) E γ E + ( μ + σ + α ) I .
Now, by computing the Jacobian matrices of Q and G at P 0 declared in (29), we obtain
G = γ + μ 0 γ α + μ + σ , Q = 0 β Λ μ + ν 0 0 .
Thus, we can evaluate R 0 * by finding the spectral radius Φ ( · ) of the (NGMM) F = Q G 1 . In other words, we have
R 0 * = Φ ( F ) = Φ ( Q G 1 ) = β γ Λ ( γ + μ ) ( μ + ν ) ( α + μ + σ ) .

4.5. Local Sensitivity and Elasticity Analysis of R 0 *

In the field of epidemiological models, the parameters’ values have a vital role in computing R 0 * . Under the assumption that all parameters are retained at specific values (estimated) including the parameter ϱ , we call this measure “local sensitivity” [17]. In such a measure, the formula used to evaluate the effect of changing one parameter ϱ on R 0 * is written as
S ϱ R 0 * = R 0 * ϱ .
However, although this method does not consider a simultaneous changing of all parameters, it is influenced highly by the values of ϱ and R 0 * . Hence, a better indicator will be the elasticity which is defined by the change in the quantity of R 0 * , as a percentage, due to the value change of the parameter ϱ , as a percentage. The formula of elasticity is written as
Υ ϱ R 0 * = R 0 * ϱ ϱ R 0 * % Δ R 0 * % Δ ϱ ,
where a positive elasticity indicates an increment in R 0 * in terms of the parameter ϱ , while a negative elasticity indicates a decrement in R 0 * , in terms of the parameter ϱ .
From Equation (32) of R 0 * and using the elasticity Equation (34), we investigate the elasticity for all parameters of the model. From this point of view, we can calculate the local elasticity of R 0 * in relation to the vaccination rate ( ν ) as
Υ ν R 0 * = R 0 * ν ν R 0 * = β γ Λ ( γ + μ ) ( μ + ν ) 2 ( α + μ + σ ) ν R 0 * = ν μ + ν < 0 .
Since the value of ν indicates the vaccination rate, then the increment of this rate leads to the decrement of R 0 * and, thus, negatively affects the spreading of the disease. To see this, we plot Figure 2 based on the data reported in Table 2. This figure clearly shows the influence of the vaccination rate v on R 0 * and, hence, confirms our aforesaid claim.
On the other hand, the local elasticity of R 0 * , in relation to the contact transmission rate ( β ), can be computed as
Υ β R 0 * = R 0 * β β R 0 * = γ Λ ( γ + μ ) ( μ + ν ) ( α + μ + σ ) β R 0 * = 1 > 0 .
Thus, the increment of the transmission rate β by 1 % leads to a rise in R 0 * . This is the case of spreading infectious diseases due to the increment in personal contact and communications between individuals and due to the pathogen specifications. However, we plot Figure 3 on the basis of the same values reported in Table 2, but here we take ν = 0.15 and allow β to be varied from 0.045 to 0.18. Such a figure also shows the influence of the transmission rate β on R 0 * , which completely coincides with our algebraic calculations mentioned previously.
In the same manner, the local elasticity of R 0 * , in relation to the transmission rate ( α ) from the I state to the R state, can be computed as
Υ α R 0 * = R 0 * α α R 0 * = β γ Λ ( γ + μ ) ( μ + ν ) ( α + μ + σ ) 2 α R 0 * = α α + μ + σ < 0 .
The inversevalue of α , which is α 1 , is the recovery period. So, the increment in α (decrement in recovery period) leads to a decrement in R 0 * . This reflects the movement rate rise from the I state to the R state, and, consequently, the individuals become immunized. This would negatively affect the spreading of disease. Now, by taking ν = 0.15 and β = 0.05 , and by keeping the remaining parameters as they are in Table 2, except for changing α within the interval 1 32 , 1 15 , we plot Figure 4, which shows the influence of the transmission rate α on R 0 * as well. This also confirms our assertion discussed previously.
It should be mentioned here, that based on the data reported in Table 2, we can easily compute the BRN as R 0 * = 0.000611 , and due to this value being less than 1, then one can deduce that the Mpox illness will not promptly spread.

5. Numerical Findings

In the following content, we intend to employ two numerical methods (the FEM and MFEM) to obtain approximate solutions of the fractional-order Mpox system (19). In fact, these methods represent two fractional versions of the conventional Euler method. One can refer to [8] for the purpose of obtaining a comprehensive view about such methods. In this context, Theorem 1 can lay the foundation to recall the FEM and MFEM, which are considered for the following fractional problem:
C D * ρ z ( x ) = g ( t , z ( t ) ) , z ( 0 ) = z 0 , t > 0 , 0 < ρ 1 .
To address problem (38), we assume [ 0 , a ] as the interval over which we wish to obtain the desired numerical solution. In general, we are not capable to find z ( t ) , which represents the analytical solution to the problem at hand. As an alternative, we can establish a set of points ( t i , z ( t i ) ) , and then employ them in finding the approximate solution [8,18]. For simplification purposes, we divide [ 0 , a ] into k sub-intervals [ t i , t i + 1 ] with h = a / k through t i = i h , for i = 0 , 1 , 2 , , k . Now, we assume that z ( t ) , C D * ρ z ( t ) , and C D * 2 ρ z ( t ) are continuous on ( 0 , a ] . Therefore, by means of Theorem 1, one might expand z ( t ) about t = t i as
z ( t ) = z ( t i ) + ( t t i ) ρ Γ ( ρ + 1 ) C D * ρ z ( t i ) + ( t t i ) 2 ρ Γ ( 2 ρ + 1 ) C D * 2 ρ z ( ξ ) ,
for some ξ ( t i , t i + 1 ) . Thus, if we substitute t i + 1 instead of t in (39), we obtain
z ( t i + 1 ) = z ( t i ) + h ρ Γ ( ρ + 1 ) C D * ρ z ( t i ) + ( t i + 1 t i ) 2 ρ Γ ( 2 ρ + 1 ) C D * 2 ρ z ( ξ ) ,
for some ξ ( t i , t i + 1 ) . Now, if one chooses h = t i + 1 t i too small, then the last term of (40) can be eliminated to obtain
z ( t i + 1 ) = z ( t i ) + h ρ Γ ( ρ + 1 ) C D * ρ z ( t i ) .
Actually, Equation (41) represents the primary formula of the FEM. In the same regard, if one substitutes
C D * ρ z ( t i + 1 ) = g t i + h ρ Γ ( ρ + 1 ) , z ( t i ) + h ρ Γ ( ρ + 1 ) g ( t i , z ( t i ) )
instead of (41), the result of z ( t i + 1 ) will be consequently yielded. In other words, we obtain
z ( t i + 1 ) = z ( t i ) + h ρ Γ ( ρ + 1 ) × g t i + h ρ Γ ( ρ + 1 ) , z ( t i ) + h ρ Γ ( ρ + 1 ) g ( t i , z ( t i ) ) ,
for i = 0 , 1 , 2 , , k 1 . Equation (42) represents the main formula of the MFEM that would be compared next with the FEM’s formula. In what follows, we attempt to implement only the MFEM on the fractional-order Mpox system (19), as the FEM is similar to the MFEM. To do so, we reconsider such a system again as
C D * ρ S ( t ) = g 1 ( t , S ( t ) , E ( t ) , I ( t ) , R ( t ) ) C D * ρ E ( t ) = g 2 ( t , S ( t ) , E ( t ) , I ( t ) , R ( t ) ) C D * ρ I ( t ) = g 3 ( t , S ( t ) , E ( t ) , I ( t ) , R ( t ) ) C D * ρ R ( t ) = g 4 ( t , S ( t ) , E ( t ) , I ( t ) , R ( t ) ) ,
where
g 1 ( t , S ( t ) , E ( t ) , I ( t ) , R ( t ) ) = Λ β S ( t ) I ( t ) ( ν + μ ) S ( t ) g 2 ( t , S ( t ) , E ( t ) , I ( t ) , R ( t ) ) = β S ( t ) I ( t ) ( μ + γ ) E ( t ) g 3 ( t , S ( t ) , E ( t ) , I ( t ) , R ( t ) ) = γ E ( t ) ( μ + σ + α ) I ( t ) g 4 ( t , S ( t ) , E ( t ) , I ( t ) , R ( t ) ) = ν S ( t ) + α I ( t ) μ R ( t ) .
More precisely, with the aim of generating ( t k , S ( t k ) ) in relation to the compartment S, one should suppose that S ( t ) , C D * ρ S ( t ) , and C D * 2 ρ S ( t ) are all continuous on ( 0 , T ] . From this perspective, if we assume that
g 1 ( t , S ( t ) , E ( t ) , I ( t ) , R ( t ) ) = Λ β S ( t ) I ( t ) ( ν + μ ) S ( t ) ,
so that
C D * ρ S ( t ) = g 1 ( t , S ( t ) , E ( t ) , I ( t ) , R ( t ) ) ,
then, with the use of (42), one might obtain
S ( t i + 1 ) = S ( t i ) + h ρ Γ ( ρ + 1 ) × g 1 t i + h ρ Γ ( ρ + 1 ) , S ( t i ) + h ρ Γ ( ρ + 1 ) g 1 ( t i , S ( t i ) ) ,
for i = 0 , 1 , 2 , , k 1 .
In a similar manner, the aforesaid approach might be applied for the remaining classes to obtain their approximate solutions. In the long run, we can infer the following approximations of model (19):
S ( t i + 1 ) = S ( t i ) + h ρ Γ ( ρ + 1 ) × g 1 t i + h ρ Γ ( ρ + 1 ) , S ( t i ) + h ρ Γ ( ρ + 1 ) g 1 ( t i , S ( t i ) ) , E ( t i + 1 ) = E ( t i ) + h ρ Γ ( ρ + 1 ) × g 2 t i + h ρ Γ ( ρ + 1 ) , E ( t i ) + h ρ Γ ( ρ + 1 ) g 2 ( t i , E ( t i ) ) , I ( t i + 1 ) = I ( t i ) + h ρ Γ ( ρ + 1 ) × g 3 t i + h ρ Γ ( ρ + 1 ) , I ( t i ) + h ρ Γ ( ρ + 1 ) g 3 ( t i , I ( t i ) ) , R ( t i + 1 ) = R ( t i ) + h ρ Γ ( ρ + 1 ) × g 4 t i + h ρ Γ ( ρ + 1 ) , R ( t i ) + h ρ Γ ( ρ + 1 ) g 4 ( t i , R ( t i ) ) ,
where g 1 , g 2 , g 3 , g 4 are already outlined in (44), i = 0 , 1 , 2 , , k 1 .
In the following content, we propose to depict certain numerical findings that demonstrate the dynamics of the fractional-order Mpox system (19). In this connection, we consider Table 3, which takes very close data to certain available data captured out of the Indian community, on the basis of reference [19,20].
On the basis of the previous discussion, we compare, in Figure 5, Figure 6, Figure 7 and Figure 8, between several MFEMs’ solutions for the I and E of system (19), in accordance with several fractional-order values of ρ . The primary goal of these simulations is to notice the impact of the vaccination rate on the exposed and infected states.
For further description and in accordance with ρ = 0.75 , we simulate, once again, the infected and exposed sizes of system (19), as shown in Figure 9 and Figure 10, by using the MFEM with and without consideration of the vaccination rate ν . Accordingly, one might obviously notice that, if the vaccination rate is raised, then the amounts of the exposed and infected states are reduced, validating the impact of such a rate on the model at hand.
For the same connection, it was demonstrated in [8] that the MFEM is a modified numerical scheme for the FEM. Thus, with the aim of validating the numerical results generated by using the MFEM, we make some comparisons between its numerical results and the FEM’s numerical results (see Figure 11, Figure 12, Figure 13 and Figure 14).
In the same regard, for the purpose of showing the validity of the domain ψ reported in Section 4.1 and the validity of the boundedness of solution of system (19) discussed in Section 4.2, one can easily check, based on the data reported in Table 3, that the sum of all classes satisfies the inequality declared in (20). In other words, one can check
S ( t ) + E ( t ) + I ( t ) + R ( t ) 1.79121 .
This bound is further confirmed numerically by plotting the sum of S ( t ) + E ( t ) + I ( t ) + R ( t ) according to different fractional-order values, as shown in Figure 15.

6. Conclusions

In this study, a novel fractional-order version of Mpox disease has been established in light of the Caputo operator. Accordingly, several innovative findings connecting to the stability analysis have been addressed. In the same connection, the proposed Mpox model has been numerically solved with the use of the FEM and MFEM. As a consequence of concentrating on the numerical findings, we can obviously notice that the behavior of the proposed fractional-order Mpox model is affected by performing any variation in the value of the fractional-order. This enables us to gain more extra degrees of freedom for the established model. In addition, we observe that the behavior of the MFEM’s solutions is completely coincided with the FEM’s ones, confirming the validity of the computational methods used in this work. Based on these solutions, it can be clearly concluded that, if the rate of vaccination is raised, then the exposed and infectious sizes are reduced gradually, confirming the influence of this rate. More precisely, if the vaccination aspect is regarded, then the infectious and exposed cases will vanish in the long run. This will consequently cause Mpox illness to be under control eventually.

Author Contributions

Conceptualization, I.M.B. and I.H.J.; methodology, A.A.A. and M.A.; software, S.B.A.-S.; validation, K.M.; formal analysis, I.M.B.; investigation, A.A.A. and M.A.; resources, S.B.A.-S.; data curation, I.H.J.; writing—original draft preparation, K.M.; writing—review and editing, I.M.B.; visualization, A.A.A. and M.A.; supervision, I.H.J.; project administration, S.B.A.-S.; funding acquisition, K.M. All authors have read and agreed to the published version of the manuscript.

Funding

The authors extend their appreciation to Arab Open University for funding this research through AOU research fund No. AOURG-2023-008.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the Arab Open University and Al-Zaytoonah University for providing the necessary scientific research supplies to implement the research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Suvvari, T.K.; Ghosh, A.; Lopinti, A.; Islam, M.A.; Bhattacharya, P. Hematological manifestations of Monkeypox virus (MPOX) and impact of human MPOX disease on blood donation—What we need to know? New Microbes New Infect. 2023, 52, 101108. [Google Scholar] [CrossRef] [PubMed]
  2. Adetifa, I.; Muyembe, J.-J.; Bausch, D.G.; Heymann, D.L. Mpox neglect and the smallpox niche: A problem for Africa, a problem for the world. Lancet 2023, 401, 1822–1824. [Google Scholar] [CrossRef]
  3. Sam-Agudu, N.A.; Martyn-Dickens, C.; Ewa, A.U. A global update of mpox (monkeypox) in children. Curr. Opin. Pediatr. 2023, 35, 193–200. [Google Scholar] [CrossRef] [PubMed]
  4. Jahanshahi, H.; Munoz-Pacheco, J.M.; Bekiros, S.; Alotaibi, N.D. A fractional-order SIRD model with time-dependent memory indexes for encompassing the multi-fractional characteristics of the COVID-19. Chaos Solitons Fractals 2021, 143, 110632. [Google Scholar] [CrossRef] [PubMed]
  5. Stuart, C.I.J.M.; Takahashi, Y.; Umezawa, H. On the stability and non-local properties of memory. J. Theor. Biol. 1978, 71, 605–618. [Google Scholar] [CrossRef]
  6. Peter, O.J.; Oguntolu, F.A.; Ojo, M.M.; Oyeniyi, A.O.; Jan, R.; Khan, I. Fractional order mathematical model of monkeypox transmission dynamics. Phys. Scr. 2022, 97, 084005. [Google Scholar] [CrossRef]
  7. Rexma Sherine, V.; Chellamani, P.; Ismail, R.; Avinash, N.; Xavier, G.B.A. Estimating the Spread of Generalized Compartmental Model of Monkeypox Virus Using a Fuzzy Fractional Laplace Transform Method. Symmetry 2022, 14, 2545. [Google Scholar] [CrossRef]
  8. Batiha, I.M.; Bataihah, A.; Al-Nana, A.A.; Alshorm, S.; Jebril, I.H.; Zraiqat, A. A numerical scheme for dealing with fractional initial value problem. Int. J. Innov. Comput. Inf. Control 2023, 19, 763–774. [Google Scholar] [CrossRef]
  9. Kilbas, A.A. Theory and Application of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  10. Almuzini, M.; Batiha, I.M.; Momani, S. A study of fractional-order monkeypox mathematical model with its stability analysis. In Proceedings of the 2023 International Conference on Fractional Differentiation and Its Applications (ICFDA), Ajman, United Arab Emirates, 14–16 March 2023; pp. 1–6. [Google Scholar] [CrossRef]
  11. Garba, S.M.; Gumel, A.B.; Lubuma, J.M.-S. Dynamically-consistent non-standard finite difference method for an epidemic model. Math. Comput. Model. 2011, 53, 131–150. [Google Scholar] [CrossRef]
  12. Lin, W. Global existence theory and chaos control of fractional differential equations. J. Math. Anal. Appl. 2007, 332, 709–726. [Google Scholar] [CrossRef]
  13. Naik, P.A.; Zu, J.; Owolabi, K.M. Global dynamics of a fractional order model for the transmission of HIV epidemic with optimal control. Chaos Solitons Fractals 2020, 138, 109826. [Google Scholar] [CrossRef] [PubMed]
  14. Iqbal, Z.; Ahmed, N.; Baleanu, D.; Adel, W.; Rafiq, M.; Aziz-ur Rehman, M.; Alshomrani, A.S. Positivity and boundedness preserving numerical algorithm for the solution of fractional nonlinear epidemic model of HIV/AIDS transmission. Chaos Solitons Fractals 2020, 134, 109706. [Google Scholar] [CrossRef]
  15. Diekmann, O.; Heesterbeek, J.A.P.; Metz, J.A. On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J. Math. Biol. 1990, 28, 365–382. [Google Scholar] [CrossRef] [PubMed]
  16. Van den Driessche, P.; Watmough, J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2002, 180, 29–48. [Google Scholar] [CrossRef] [PubMed]
  17. Morio, J. Global and local sensitivity analysis methods for a physical system. Eur. J. Phys. 2011, 32, 1577. [Google Scholar] [CrossRef]
  18. Zraiqat, A.; Paikrayb, S.K.; Duttac, H. A Certain Class of Deferred Weighted Statistical B-Summability Involving (p,q)-Integers and Analogous Approximation Theorems. Filomat 2019, 33, 1425–1444. [Google Scholar]
  19. Statisticstimes.com. Countries by GDP Growth. 2020. Available online: https://statisticstimes.com (accessed on 15 March 2023).
  20. van den Driessche, P. Reproduction numbers of infectious disease models. Infect. Dis. Model. 2017, 2, 288–303. [Google Scholar] [PubMed]
Figure 1. The compartment model, including vaccination rate.
Figure 1. The compartment model, including vaccination rate.
Algorithms 16 00418 g001
Figure 2. Influence of the vaccination rate v on R 0 * .
Figure 2. Influence of the vaccination rate v on R 0 * .
Algorithms 16 00418 g002
Figure 3. Influence of the transmission rate β on R 0 * .
Figure 3. Influence of the transmission rate β on R 0 * .
Algorithms 16 00418 g003
Figure 4. Influence of the transmission rate α on R 0 * .
Figure 4. Influence of the transmission rate α on R 0 * .
Algorithms 16 00418 g004
Figure 5. Exposed sizes E × 10 2 of system (19) for ρ = 1 , 0.9 , 0.8 , 0.7 , 0.6 through MFEM with ν = 0 .
Figure 5. Exposed sizes E × 10 2 of system (19) for ρ = 1 , 0.9 , 0.8 , 0.7 , 0.6 through MFEM with ν = 0 .
Algorithms 16 00418 g005
Figure 6. Exposed sizes E × 10 2 of system (19) for ρ = 1 , 0.9 , 0.8 , 0.7 , 0.6 through MFEM with ν = 0.5 .
Figure 6. Exposed sizes E × 10 2 of system (19) for ρ = 1 , 0.9 , 0.8 , 0.7 , 0.6 through MFEM with ν = 0.5 .
Algorithms 16 00418 g006
Figure 7. Infected sizes I × 10 2 of system (19) for ρ = 1 , 0.9 , 0.8 , 0.7 , 0.6 through MFEM with ν = 0 .
Figure 7. Infected sizes I × 10 2 of system (19) for ρ = 1 , 0.9 , 0.8 , 0.7 , 0.6 through MFEM with ν = 0 .
Algorithms 16 00418 g007
Figure 8. Infected sizes I × 10 2 of system (19) for ρ = 1 , 0.9 , 0.8 , 0.7 , 0.6 through MFEM with ν = 0.5 .
Figure 8. Infected sizes I × 10 2 of system (19) for ρ = 1 , 0.9 , 0.8 , 0.7 , 0.6 through MFEM with ν = 0.5 .
Algorithms 16 00418 g008
Figure 9. Exposed sizes E × 10 2 of system (19) for ρ = 0.75 with and without considering ν .
Figure 9. Exposed sizes E × 10 2 of system (19) for ρ = 0.75 with and without considering ν .
Algorithms 16 00418 g009
Figure 10. Infected sizes I × 10 2 of system (19) for ρ = 0.75 with and without considering ν .
Figure 10. Infected sizes I × 10 2 of system (19) for ρ = 0.75 with and without considering ν .
Algorithms 16 00418 g010
Figure 11. Susceptible sizes S × 10 2 of system (19) for ρ = 0.6 , 0.7 , 0.8 , 0.9 , 1 using FEM and MFEM.
Figure 11. Susceptible sizes S × 10 2 of system (19) for ρ = 0.6 , 0.7 , 0.8 , 0.9 , 1 using FEM and MFEM.
Algorithms 16 00418 g011
Figure 12. Exposed sizes E × 10 2 of system (19) for ρ = 0.6 , 0.7 , 0.8 , 0.9 , 1 using FEM and MFEM.
Figure 12. Exposed sizes E × 10 2 of system (19) for ρ = 0.6 , 0.7 , 0.8 , 0.9 , 1 using FEM and MFEM.
Algorithms 16 00418 g012
Figure 13. Infected sizes I × 10 2 of system (19) for ρ = 0.6 , 0.7 , 0.8 , 0.9 , 1 using FEM and MFEM.
Figure 13. Infected sizes I × 10 2 of system (19) for ρ = 0.6 , 0.7 , 0.8 , 0.9 , 1 using FEM and MFEM.
Algorithms 16 00418 g013
Figure 14. Recovered sizes R × 10 2 of system (19) for ρ = 0.6 , 0.7 , 0.8 , 0.9 , 1 using FEM and MFEM.
Figure 14. Recovered sizes R × 10 2 of system (19) for ρ = 0.6 , 0.7 , 0.8 , 0.9 , 1 using FEM and MFEM.
Algorithms 16 00418 g014
Figure 15. The sizes of S ( t ) + E ( t ) + I ( t ) + R ( t ) for ρ = 0.6 , 0.8 , 1 using MFEM.
Figure 15. The sizes of S ( t ) + E ( t ) + I ( t ) + R ( t ) for ρ = 0.6 , 0.8 , 1 using MFEM.
Algorithms 16 00418 g015
Table 1. Model parameters.
Table 1. Model parameters.
ParameterDenotationValue
Λ Rate of birthDepending on country
γ Rate of transmission (E to I) 1 17 , 1 7
β Rate of transmission (S to E) [ 0.045 , 0.18 ]
α Rate of transmission (I to R) 1 32 , 1 15
σ Rate of fatality [ 0.03 , 0.06 ]
ν Rate of vaccinationVariable [ 0 , 1 ]
μ Rate of natural deathDependent on country
Table 2. The parameters of the fractional-order Mpox model.
Table 2. The parameters of the fractional-order Mpox model.
ParameterValue
Λ 163/10,000
β 0.5
γ 0.05
α 1/17
μ 91/10,000
ν Changes within [ 0 , 1 ]
σ 0.03
Table 3. Parameters of system (19) [19].
Table 3. Parameters of system (19) [19].
ParameterValue
S ( 0 ) 23
E ( 0 ) 6
I ( 0 ) 4
R ( 0 ) 0
Λ 0.01630
β 0.05000
γ 0.05882
α 0.03125
μ 0.00910
ν [ 0 , 1 ]
σ 0.0300
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Batiha, I.M.; Abubaker, A.A.; Jebril, I.H.; Al-Shaikh, S.B.; Matarneh, K.; Almuzini, M. A Mathematical Study on a Fractional-Order SEIR Mpox Model: Analysis and Vaccination Influence. Algorithms 2023, 16, 418. https://doi.org/10.3390/a16090418

AMA Style

Batiha IM, Abubaker AA, Jebril IH, Al-Shaikh SB, Matarneh K, Almuzini M. A Mathematical Study on a Fractional-Order SEIR Mpox Model: Analysis and Vaccination Influence. Algorithms. 2023; 16(9):418. https://doi.org/10.3390/a16090418

Chicago/Turabian Style

Batiha, Iqbal M., Ahmad A. Abubaker, Iqbal H. Jebril, Suha B. Al-Shaikh, Khaled Matarneh, and Manal Almuzini. 2023. "A Mathematical Study on a Fractional-Order SEIR Mpox Model: Analysis and Vaccination Influence" Algorithms 16, no. 9: 418. https://doi.org/10.3390/a16090418

APA Style

Batiha, I. M., Abubaker, A. A., Jebril, I. H., Al-Shaikh, S. B., Matarneh, K., & Almuzini, M. (2023). A Mathematical Study on a Fractional-Order SEIR Mpox Model: Analysis and Vaccination Influence. Algorithms, 16(9), 418. https://doi.org/10.3390/a16090418

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop