Next Article in Journal
May the Piezoresistivity of GNP-Modified Cement Mortar Be Related to Its Fractal Structure?
Next Article in Special Issue
A New Three-Step Root-Finding Numerical Method and Its Fractal Global Behavior
Previous Article in Journal
Mittag–Leffler Function as an Approximant to the Concentrated Ferrofluid’s Magnetization Curve
Previous Article in Special Issue
On Weighted (k, s)-Riemann-Liouville Fractional Operators and Solution of Fractional Kinetic Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional Dynamics of Typhoid Fever Transmission Models with Mass Vaccination Perspectives

by
Hamadjam Abboubakar
1,*,
Raissa Kom Regonne
2,† and
Kottakkaran Sooppy Nisar
3,*,†
1
Department of Computer Engineering, University Institute of Technology of Ngaoundéré, University of Ngaoundéré, Ngaoundéré P.O. Box 455, Cameroon
2
National School of Agro-Industrial Sciences, University of Ngaoundéré, Ngaoundéré P.O. Box 455, Cameroon
3
Department of Mathematics, College of Arts and Science, Prince Sattam bin Abdulaziz University, Wadi Aldawaser 11991, Saudi Arabia
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Fractal Fract. 2021, 5(4), 149; https://doi.org/10.3390/fractalfract5040149
Submission received: 7 August 2021 / Revised: 20 September 2021 / Accepted: 20 September 2021 / Published: 30 September 2021

Abstract

:
In this work, we formulate and mathematically study integer and fractional models of typhoid fever transmission dynamics. The models include vaccination as a control measure. After recalling some preliminary results for the integer model (determination of the epidemiological threshold denoted by R c , asymptotic stability of the equilibrium point without disease whenever R c < 1 , the existence of an equilibrium point with disease whenever R c > 1 ), we replace the integer derivative with the Caputo derivative. We perform a stability analysis of the disease-free equilibrium and prove the existence and uniqueness of the solution of the fractional model using fixed point theory. We construct the numerical scheme and prove its stability. Simulation results show that when the fractional-order η decreases, the peak of infected humans is delayed. To reduce the proliferation of the disease, mass vaccination combined with environmental sanitation is recommended. We then extend the previous model by replacing the mass action incidences with standard incidences. We compute the corresponding epidemiological threshold denoted by R c and ensure the uniform stability of the disease-free equilibrium, for both new models, when R c < 1 . A new calibration of the new model is conducted with real data of Mbandjock, Cameroon, to estimate R c = 1.4348 . We finally perform several numerical simulations that permit us to conclude that such diseases can possibly be tackled through vaccination combined with environmental sanitation.

1. Introduction

Typhoid fever, caused by a salmonella bacterium (Salmonella typhi), is a tropical disease transmitted by the ingestion of food or/and water contaminated with feces. It is most prevalent in countries located below the equator, in Southeast Asia, and in the Indian subcontinent, where hygiene conditions are poor [1,2]. The principal symptoms of typhoid fever are insomnia, fever, generalized fatigue, headaches, stomach ache, anorexia, constipation or diarrhea, and vomiting. These symptoms can last several weeks. Without effective treatment, typhoid fever can lead to death. According to the World Health Organization (WHO), the number of cases of typhoid fever is estimated to be between 11 and 21 million, with 128,000 to 161,000 deaths annually due to the severity of the disease [1,2]. Vaccination, sanitary measures, and hygiene measures are the best ways to prevent the spread of the disease [2].
Since the work of Sir Ronald Ross on malaria [3], mathematical tools such as differential equations have usually been used to understand and describe the dynamics of infectious diseases [4,5,6,7,8,9,10,11,12]. In [7], the authors proposed a S V I I c R that takes into account some control mechanisms (treatment, education campaigns, and vaccination). Quarantining the infected individuals and their treatment are the main control measures studied in [8]. The author used optimal control methods to conclude that the outbreak can be eliminated or controlled if the control strategies are combined to their highest levels. The same conclusions are given in [10,11]. Recently, Olumuyiwa James et al. [9] formulated and studied an optimal control model for typhoid fever that takes into account both indirect and direct transmission. They compared various proposed strategies using numerical simulations and concluded that the disease burden can be controlled if all the available control measures are combined.
Very recently, many authors have proposed fractional-order models in mathematical epidemiology (ref. [13]), ecology (ref. [14]), plant epidemiology (refs. [15,16]), and psychology (ref. [17]). Indeed, the necessity of the use of fractional derivatives in epidemiology, for example, comes from the fact that these operators have many properties, such as their different types of kernels and the crossover behavior in the model, which can only be described using these operators. Moreover, any real data that have zigzag dynamics (mostly many) that cannot be projected by an integer-order derivative can be solved by a fractional model more clearly.
The principal disadvantage of models with integer derivatives is that they do not permit the definition of memory effects. Replacing integer derivatives with fractional derivatives makes it possible to remedy this problem. Indeed, they offer different ways to forecast data by varying the fractional-order parameter [12,18,19]. Several fractional operators have been defined so far. The most popular are the fractional operators of Caputo, the fractional Caputo–Fabrizio operator, and the fractional operator of Atangana–Baleanu. Each operator explores the dynamics of the studied phenomenon differently, thus helping us to predict more variations in the evolution of the phenomenon. The advantages and disadvantages of each fractional operator and their application domains can be found in [20,21,22].
To the best of our knowledge, there are few mathematical works on typhoid fever using fractional operators [6,12]. In [12], the authors used the Caputo–Fabrizio operator to extend the model proposed in [23]. They provided existence, uniqueness, and stability criteria for the proposed fractional-order typhoid model. More recently, Abboubakar et al. [6] formulated a S I R -B-type compartmental model with both integer and Caputo derivatives. The only control measure was vaccination. They computed the control reproduction number, R c , and performed stability analysis of the disease-free equilibrium point for both models. The present contributions are listed as follows:
  • Using a fractional derivative in place of an integer derivative, as used in our previous model [5], we formulate a new model. To prove the existence and uniqueness of the solutions, we use fixed point theory. The corresponding numerical scheme is obtained through the Adams–Bashforth–Moulton method [24,25]. The stability of this numerical scheme is also proven. Finally, several numerical simulations are carried out from the real values of parameters estimated with real data of Mbandjock, in Cameroon (see [5]).
  • Secondly, we extend the previous models by replacing the mass action incidence law with the standard incidence law. For these new models, we compute the corresponding control reproduction number, R c , and ensure the uniform stability of the equilibrium point without disease. As in [6], model parameters are estimated. With these new parameter values, we finally perform several numerical simulations that permit us to compare the quantitative dynamics of the two types of models.
The rest of the work is organized as follows. We devote Section 2 to preliminary definitions of the fractional derivative in the sense of Caputo and useful results. Formulation of the models with mass incidence law and standard incidence, as well as their mathematical analysis (computation of control reproduction numbers, asymptotic stability analysis of the disease-free equilibrium, existence as well as the uniqueness of solutions, construction of the numerical scheme with its stability), is also described in this section. The calibration of the model with standard incidences and several numerical results is given in Section 3. We end the paper with a discussion and conclusions.

2. Materials and Methods

2.1. Useful Definitions and Results

For over ten years, fractional derivatives have captured the attention of researchers, who use these fractional operators to model physical, chemical, and biological processes. One can cite the dynamics of transmissible diseases [26,27,28]. Before the formulation of the fractional models, it is important to recall their definition, as well as two results that will be used later in the fractional model analysis.
Definition 1.
Let f C 1 l , and we have the following relation:
D τ ν f ( τ ) = d r f ( τ ) d τ r , ν = r N 1 Γ ( r ν ) 0 τ ( τ ι ) r ν 1 f ( r ) ( ι ) d ι , 1 + r < ν < r , r N ,
which represents the Caputo derivative of f.
Lemma 1.
Assume that χ , Q , h , Y > 0 , k h Y with k N , and
y q , m = ( m q ) χ 1 q = 1 , 2 , , m 1 , 0 q m .
Let q = k q = i y q , m | e q | = 0 for k > m 1 .
If
| e m | Q h χ q = 1 m 1 y q , m | e q | + | η 0 | , m = 1 , 2 , , k ,
then
| e k | M | η 0 | , k { 1 , 2 , }
where M R + does not depend on h and k .
Lemma 2.
If 0 < χ < 1 and d N , then there exist positive constants W χ , 1 and W χ , 2 only dependent on χ , such that
( 1 + v ) χ v χ W χ , 1 ( 1 + v ) χ 1 and ( 2 + v ) χ + 1 2 ( 1 + v ) χ + 1 + v χ + 1 W χ , 2 ( 1 + v ) χ 1 .

2.2. Model Dynamics with Mass Action Incidence Law

2.2.1. Model Formulation in ODE Sense and Its Analysis

In a previous work [5], we formulated and analyzed a new mathematical model for the transmission dynamics of typhoid fever with application to the town of Mbandjock, in the central region of Cameroon. The model is divided into seven compartments: susceptible individuals S ( t ) , vaccinated individuals V ( t ) , infected individuals in latent stage E ( t ) , infected individuals without any sign of the disease C ( t ) , symptomatic infected individuals I ( t ) , recovered individuals R ( t ) , and the density of bacteria in the environment B ( t ) . Each individual in each compartment naturally dies at a rate μ h . Susceptible humans are recruited at a constant rate Λ h . The population in compartment S h decreases either by vaccination at a rate ξ , or by infection at an incidence rate ν B ( t ) S ( t ) , where ν is the contact rate. We denote by θ the rate at which vaccinated individuals lose their immunity. The vaccine efficacy is denoted by ϵ . The compartment E of latent individuals, which include infected susceptible individuals and vaccinated individuals, progresses either to the carriers compartment C at a rate q γ 1 or to the compartment I at the rate ( 1 q ) γ 1 . Asymptomatic individuals become symptomatic at the rate ( 1 p ) γ 2 or recover at the rate p γ 2 . δ denotes the disease-induced death rate of symptomatic individuals. Recovered individuals become susceptible at a rate α . With this brief description, the mathematical formulation of the model studied in [5] is presented as follows:
S ˙ ( t ) = Λ h ( k 1 + ν B ( t ) ) S ( t ) + θ V ( t ) + α R ( t ) ,
V ˙ ( t ) = k 2 + ( 1 ϵ ) ν B ( t ) V ( t ) + ξ S ( t ) ,
E ˙ ( t ) = k 3 E ( t ) + ν π V ( t ) + S ( t ) B ( t ) ,
C ˙ ( t ) = q γ 1 E ( t ) k 4 C ( t ) ,
I ˙ ( t ) = q 1 γ 1 E ( t ) + p 1 γ 2 C ( t ) k 5 + σ I ( t ) ,
R ˙ ( t ) = p γ 2 C ( t ) + σ I ( t ) k 6 R ( t ) ,
B ˙ ( t ) = p c C ( t ) + p i I ( t ) μ b B ( t ) ,
where k 1 = ξ + μ h , k 2 = θ + μ h , k 3 = γ 1 + μ h , k 4 = μ h + γ 2 , k 5 = δ + μ h , k 6 = α + μ h , π = 1 + ϵ , q 1 = 1 q , p 1 = p + 1 , k 7 = k 1 k 2 θ ξ = μ h ( k 2 + ξ ) > 0 , k 8 = k 5 + σ .
Model (2) is defined in the following set:
W = S , V , E , C , I , R , B R + 7 : N = V + S + C + E + I + R Λ h μ h ; B ( p i + p c ) Λ h μ h μ b ,
in which a dynamical system is defined, and where N denotes the human population.
Without disease, model (2) has the following equilibrium: Q 0 = S 0 , V 0 , 0 , 0 , 0 , 0 , 0 , where S 0 = Λ h k 2 / ( μ h ( k 2 + ξ ) ) and V 0 = Λ h ξ / ( μ h ( k 2 + ξ ) ) . Using the same approach developed in [29], we obtain the control reproduction number given by
R c = ν Λ h ( k 2 + π ξ ) γ 1 p c q ( σ + k 5 ) + p i ( k 4 ( 1 q ) + γ 2 q ( 1 p ) ) μ b μ h k 3 k 4 ( k 2 + ξ ) ( σ + k 5 ) .
Considering the model without vaccination, R c is equal to the basic reproduction number:
R 0 = ν Λ h γ 1 p c q ( σ + k 5 ) + p i ( k 4 ( 1 q ) + γ 2 q ( 1 p ) ) μ b μ h k 3 k 4 ( σ + k 5 ) .
Thus, it follows that
R c = R 0 ( k 2 + π ξ ) ( k 2 + ξ ) .
Since π ξ = ( 1 ϵ ) ξ ξ , we have ( k 2 + π ξ ) ( k 2 + ξ ) 1 , which means that R c R 0 . This proves that mass vaccination is a useful tool that can be used to effectively tackle this kind of tropical disease.
For typhoid model (2) in the ODE sense, the following results were proven in [5].
Proposition 1
([5]). For model (2), Q 0 is locally and globally asymptotically stable in W if R c < 1 and unstable if R c > 1 .
Proposition 2
([5]). Let us define the following coefficients:
a 2 = R c 4 k 3 2 k 4 2 k 8 2 μ h 2 ( ξ + k 2 ) 2 π × × ( μ h γ 1 α ( k 5 + σ q ) + γ 2 μ h k 8 ( α + γ 1 ) + γ 1 γ 2 α ( 1 p q ) k 5 + [ μ h 2 + ( α + γ 2 + γ 1 ) μ h ] μ h k 8 ) , a 1 = R c 2 k 3 2 k 4 2 k 8 2 μ h 2 ( ξ + k 2 ) 2 γ 1 Λ h k 6 π ( k 4 q 1 + γ 2 p 1 q ) ( R c 2 R b 2 ) , a 0 = γ 1 2 μ h Λ h 2 k 3 k 4 k 6 k 8 ( k 4 q 1 + γ 2 p 1 q ) 2 ( π ξ + k 2 ) 2 ( ξ + k 2 ) ( R c 2 1 ) .
Model (2) with the integer derivative either has (1) only one endemic equilibrium whenever ( a 0 < 0 R c > 1 ) or ( a 1 < 0 and a 0 = 0 or a 1 2 4 a 2 a 0 = 0 ), (2) two endemic equilibrium points if ( a 0 > 0 ( R c < 1 ) , a 1 < 0 ( R c > R b ) and a 1 2 4 a 2 a 0 > 0 ), or (3) no equilibrium otherwise.
Theorem 1
([5]). Model (2) exhibits a supercritical bifurcation at R c = 1 , which implies that whenever R c > 1 , the endemic equilibrium is locally asymptotically stable.
Remark 1.
Proposition 1 combined with Theorem 1 implies that Proposition 2 (iii) will never hold true for model (2). Thus, the condition R c < 1 is sufficient to eradicate the disease.

2.2.2. Fractional-Order Typhoid Model

The following model is obtained when we replace the integer derivative operator in (2) with the non-integer operator in the Caputo sense.
t 0 C D t η S ( t ) = Λ h ( k 1 + ν B ( t ) ) S ( t ) + θ V ( t ) + α R ( t ) ,
t 0 C D t η V ( t ) = k 2 + ( 1 ϵ ) ν B ( t ) V ( t ) + ξ S ( t ) ,
t 0 C D t η E ( t ) = k 3 E ( t ) + ν π V ( t ) + S ( t ) B ( t ) ,
t 0 C D t η C ( t ) = q γ 1 E ( t ) k 4 C ( t ) ,
t 0 C D t η I ( t ) = q 1 γ 1 E ( t ) + p 1 γ 2 C ( t ) k 5 + σ I ( t ) ,
t 0 C D t η R ( t ) = p γ 2 C ( t ) + σ I ( t ) k 6 R ( t ) ,
t 0 C D t η B ( t ) = p c C ( t ) + p i I ( t ) μ b B ( t ) .
with S ( 0 ) 0 , V ( 0 ) 0 , B ( 0 ) 0 , E ( 0 ) 0 , I ( 0 ) 0 , C ( 0 ) 0 , and R ( 0 ) 0 .

Asymptotic Stability of the Disease-Free Equilibrium

Before investigating the stability of the disease-free equilibrium point, we consider the fractional-order system (5) as follows:
t 0 C D t η x ( t ) = F x ( t ) ,
where x ( ζ ) R 7 , F R 7 × R 7 , 0 < η < 1 . The characteristic equation of the matrix F evaluated at any equilibrium (see [30]) is given by
det s ( I ( 1 η ) F ) η F = 0 .
For the fractional model (5), the asymptotic stability of Q 0 is claimed in the following result.
Theorem 2.
The disease-free equilibrium Q 0 is uniformly asymptotically stable if R c < 1 , and unstable otherwise.
Proof. 
The Jacobian matrix of system (5) evaluated at the disease–free equilibrium Q 0 is given by
J ( Q 0 ) = k 1 θ 0 0 0 α ν S 0 ξ k 2 0 0 0 0 ν π V 0 0 0 k 3 0 0 0 ν ( S 0 + π V 0 ) 0 0 q γ 1 k 4 0 0 0 0 0 q 1 γ 1 p 1 γ 2 ( k 5 + σ ) 0 0 0 0 0 p γ 2 σ k 6 0 0 0 0 p c p i 0 μ b = J 1 J 2 J 3 J 4 ,
where J 1 = k 1 θ ξ k 2 , J 2 = 0 0 0 α ν S 0 0 0 0 0 ν π V 0 , J 3 = O R 5 × 2 , and
J 4 = k 3 0 0 0 ν ( S 0 + π V 0 ) q γ 1 k 4 0 0 0 q 1 γ 1 p 1 γ 2 ( k 5 + σ ) 0 0 0 p γ 2 σ k 6 0 0 p c p i 0 μ b .
The characteristic Equation (7) of the typhoid fractional model (5) becomes
det s ( I 2 ( 1 η ) J 1 ) η J 1 = 0 det s ( I 5 ( 1 η ) J 4 ) η J 4 = 0 ,
which is equivalent to
a 2 s 2 + a 1 s + a 0 = 0 , s ( 1 η ) k 6 + 1 + k 6 η s 4 + A 1 A 5 s 3 + A 2 A 5 s 2 + A 3 A 5 s + A 4 A 5 = 0 ,
where a 2 = ( 1 η ) 2 k 7 + ( 1 η ) ( k 1 + k 2 ) + 1 , a 1 = 2 η ( 1 η ) k 7 + η ( k 1 + k 2 ) , a 0 = η 2 k 7 ,
A 5 = η 1 3 k 3 k 4 k 8 μ b η 1 k 9 + q 1 p i + p c q ( 1 R c 2 ) + η 1 3 μ b k 3 k 8 k 8 p c q + p 1 γ 2 p i q + k 3 k 4 q 1 k 4 p i + p 1 γ 2 p i q + k 4 k 8 k 9 + η 1 2 μ b k 8 k 9 + k 4 k 9 + k 3 k 9 + η 1 k 9 μ b + η 1 3 k 3 + η 1 2 k 4 + η 1 2 k 3 + η 1 k 8 + η 1 2 k 3 + η 1 k 4 + η 1 k 3 + 1 k 9 ,
A 1 = 4 k 3 k 4 k 8 k 9 μ b η 1 3 η + 3 η 1 2 k 3 k 8 μ b η q 1 k 4 p i + 3 η 1 2 η k 3 k 4 μ b k 8 p c q ( 1 R c 2 ) + 3 η 1 2 η k 3 k 4 μ b q 1 k 4 p i + p 1 γ 2 p i q + 3 η 1 2 k 3 k 8 μ b η k 8 p c q + p 1 γ 2 p i q + 2 η 1 k 4 k 9 μ b η + ( 2 η 1 k 3 + 1 ) k 9 μ b η + 3 η 1 2 k 4 k 8 k 9 μ b η + 2 η 1 k 8 k 9 μ b η + 3 η 1 2 k 3 + 2 η 1 k 4 + 2 η 1 k 3 + 1 k 8 + 2 η 1 k 3 + 1 k 4 + k 3 k 9 η ,
A 2 = 6 η 1 2 k 3 k 4 k 8 k 9 μ b η 2 + 3 η 1 k 3 k 4 k 8 μ b η 2 p c q + q 1 p i ( 1 R c 2 ) + 3 η 1 η 2 k 3 k 8 μ b k 8 p c q + p 1 γ 2 p i q + k 4 k 8 k 9 μ b + k 3 + k 3 k 4 μ b q 1 k 4 p i + p 1 γ 2 p i q + k 4 k 8 + k 3 k 8 + k 3 k 4 + k 8 μ b + k 4 μ b + k 3 μ b k 9 η 2 ,
A 3 = μ b η 3 k 3 k 4 k 8 4 k 9 η 1 + q 1 p i + p c q ( 1 R c 2 ) + ( μ b + k 3 ) k 4 k 8 k 9 η 3 + k 3 μ b η 3 k 8 q ( k 8 p c + p 1 γ 2 p i ) + k 4 p i ( q 1 k 4 + p 1 γ 2 q ) ,
A 4 = 1 R c 2 η 4 k 3 k 4 k 8 μ b q 1 k 4 p i + k 8 p c q + p 1 γ 2 p i q k 9 ,
and η 1 = 1 η .
Since both coefficients of the first equation of (9) are positive, it follows that the real parts of the solution of (9) are negative. From the second equation of (9), we have that s 1 = k 6 η ( 1 η ) k 6 + 1 is one solution, and the other solutions are the root of T ( s ) : = A 5 s 4 + A 1 s 3 + A 2 s 2 + A 3 s + A 4 . Note that A 4 > 0 R c < 1 , and R c < 1 A 1 > 0 , A 2 > 0 , A 3 > 0 , and A 5 > 0 . Then, it follows that, if R c < 1 , then the disease-free equilibrium Q 0 is asymptotically stable whenever the following Routh–Hurwitz criteria A 1 A 2 A 5 2 A 3 A 5 > 0 and A 1 A 2 A 3 A 5 3 A 1 2 A 4 A 5 3 A 3 2 A 5 2 > 0 are satisfied for polynomial T ( s ) .
For uniform stability, we use the Lyapunov function:
M ( E , C , I , B ) = b 1 E ( t ) + b 2 C ( t ) + b 3 I ( t ) + b 4 B ( t ) ,
where b 1 = 1 , b 2 = k 3 ( k 8 p c + p i p 1 γ 2 ) / [ k 4 p i q 1 γ 1 + q γ 1 ( k 8 p c + p i p 1 γ 2 ) ] , and b 3 = k 3 k 4 p i / [ k 4 p i q 1 γ 1 + q γ 1 ( k 8 p c + p i p 1 γ 2 ) ] and b 4 = k 3 k 4 k 8 / [ k 4 p i q 1 γ 1 + q γ 1 ( k 8 p c + p i p 1 γ 2 ) ] .
The fractional derivative of M is given by
t 0 C D t η M ( E , C , I , B ) t 0 C D t η E b 1 + t 0 C D t η C b 2 + t 0 C D t η I b 3 + t 0 C D t η B b 4 = b 1 ν B S + π V k 3 E + b 2 q γ 1 E k 4 C + b 3 q 1 γ 1 E + p 1 γ 2 C k 8 I + b 4 p c C + p i I μ b B b 1 ν B S 0 + π V 0 k 3 E + b 2 q γ 1 E k 4 C + b 3 q 1 γ 1 E + p 1 γ 2 C k 8 I + b 4 p c C + p i I μ b B = b 1 ν B S 0 + π V 0 b 1 k 3 E + b 2 q γ 1 E b 2 k 4 C + b 3 q 1 γ 1 E + b 3 p 1 γ 2 C b 3 k 8 I + b 4 p c C + b 4 p i I b 4 μ b B = b 1 ν B S 0 + π V 0 b 4 μ b B + b 3 q 1 γ 1 + b 2 q γ 1 b 1 k 3 E + b 4 p c + b 3 p 1 γ 2 b 2 k 4 C + b 4 p i b 3 k 8 I = μ b k 3 k 4 ( k 5 + σ ) p i q 1 γ 1 k 4 + q γ 1 ( p 1 p i γ 2 + p c ( k 5 + σ ) ) R c 2 1 B .
Thus, t 0 C D t η M ( E , C , I , B ) < 0 whenever R c < 1 , and t 0 C D t η M ( E , C , I , B ) = 0 if and only if R c = 1 or B ( ζ ) = 0 . Setting B = 0 in (5), we obtain S = S 0 , V = V 0 , and E = C = I = R = 0 . Thus, lim t S ( t ) , V ( t ) , E ( t ) , C ( t ) , I ( t ) , R ( t ) S 0 , V 0 , 0 , 0 , 0 , 0 , 0 : = Q 0 . Consequently, by [31] (Theorem 2.5), it follows that if R c < 1 , then Q 0 is uniformly asymptotically stable in W . □

Existence and Uniqueness Analysis

Before describing the existence and uniqueness of solutions for the fractional model using the fixed point theorem, we define T as a Banach space of continuous functions defined on an interval P with the norm
X = i = 1 i = 7 X i ,
where X = S , V , E , C , I , R , B , X i = sup { | X i ( t ) | : t P } , and T = M ( P ) × M ( P ) × M ( P ) × M ( P ) × M ( P ) × M ( P ) × M ( P ) .
Let us write system (5) as follows:
C D t η S ( t ) = H 1 ( t , S ) , C D t η V ( t ) = H 2 ( t , V ) , C D t η E ( t ) = H 3 ( t , E ) , C D t η C ( t ) = H 4 ( t , C ) C D t η I ( t ) = H 5 ( t , I ) C D t η R ( t ) = H 6 ( t , R ) C D t η B ( t ) = H 7 ( t , B )
Application of the Caputo fractional integral operator permits us to reduce (11) to the following system:
S ( 0 ) + S ( t ) = 0 t ( t ς ) η 1 H 1 ( ς , S ) d ς 1 Γ ( η ) , V ( 0 ) + V ( t ) = 0 t ( t ς ) η 1 H 2 ( ς , V ) d ς 1 Γ ( η ) , E ( 0 ) + E ( t ) = 0 t ( t ς ) η 1 H 3 ( ς , E ) d ς 1 Γ ( η ) , C ( 0 ) + C ( t ) = 0 t ( t ς ) η 1 H 4 ( ς , C ) d ς 1 Γ ( η ) , I ( 0 ) + I ( t ) = 0 t ( t ς ) η 1 H 5 ( ς , I ) d ς 1 Γ ( η ) , R ( 0 ) + R ( t ) = 0 t ( t ς ) η 1 H 6 ( ς , R ) d ς 1 Γ ( η ) , B ( 0 ) + B ( t ) = 0 t ( t ς ) η 1 H 7 ( ς , R ) d ς 1 Γ ( η ) ,
with 0 < η < 1 .
Now, we will provide the Lipschitz conditions fulfilled by H i , for i = 1 , 2 , , 7 , as well as the contraction conditions. In the following theorem, we only provide the condition for H 1 , the rest being similar.
Theorem 3.
The kernel H 1 satisfies the Lipschitz and contraction conditions provided that 0 ν κ 7 + k 1 < 1 .
Proof. 
For S, we proceed as below:
H 1 ( t , S ) H 1 ( t , S 1 ) = ν B ( S S 1 ) k 1 ( S S 1 ) = ν B + k 1 ( S ( t ) S 1 ( t ) ) ν κ 7 + k 1 S ( t ) S 1 ( t )
where κ 7 is the upper bound of the function B ( t ) . Now, setting W 1 = ν κ 7 + k 1 , we finally obtain
H 1 ( t , S ) H 1 ( t , S 1 ) W 1 S ( t ) S ( t 1 ) ,
which provide the Lipschitz condition. If, additionally, we can have 0 < W 1 = ν κ 7 + k 1 < 1 , then the contraction is also obtained.
As in the case of H 1 , it is easy to obtain the Lipschitz condition for the other kernels. Thus, we have
H 2 ( t , C ) H 2 ( t , C 1 ) W 2 V ( t ) V ( t 1 ) , H 3 ( t , E ) H 3 ( t , E 1 ) W 3 E ( t ) E ( t 1 ) , H 4 ( t , C ) H 4 ( t , C 1 ) W 4 C ( t ) C ( t 1 ) , H 5 ( t , I ) H 5 ( t , I 1 ) W 5 I ( t ) I ( t 1 ) , H 6 ( t , R ) H 6 ( t , R 1 ) W 6 R ( t ) R ( t 1 ) , H 7 ( t , B ) H 7 ( t , B 1 ) W 7 B ( t ) B ( t 1 ) .
 □
Recursively, Equation (12) can be rewritten as follows:
S n ( t ) S ( 0 ) = 1 Γ ( η ) 0 t ( t ς ) η 1 H 1 ( ς , S n 1 ) d ς , V n ( t ) V ( 0 ) = 1 Γ ( η ) 0 t ( t ς ) η 1 H 2 ( ς , V n 1 ) d ς , E n ( t ) E ( 0 ) = 1 Γ ( η ) 0 t ( t ς ) η 1 H 3 ( ς , E n 1 ) d ς , C n ( t ) C ( 0 ) = 1 Γ ( η ) 0 t ( t ς ) η 1 H 4 ( ς , C n 1 ) d ς , I n ( t ) I ( 0 ) = 1 Γ ( η ) 0 t ( t ς ) η 1 H 5 ( ς , I n 1 ) d ς , R n ( t ) R ( 0 ) = 1 Γ ( η ) 0 t ( t ς ) η 1 H 6 ( ς , R n 1 ) d ς , B n ( t ) B ( 0 ) = 1 Γ ( η ) 0 t ( t ς ) η 1 H 7 ( ς , B n 1 ) d ς .
By determing, in a recursive manner, the difference between the successive terms of (11), we obtain
ψ 1 n ( t ) = S n ( t ) S n 1 ( t ) = 1 Γ ( η ) 0 t ( t ς ) η 1 ( H 1 ( ς , S n 1 ) H 1 ( ς , S n 2 ) ) d ς , ψ 2 n ( t ) = V n ( t ) V n 1 ( t ) = 1 Γ ( η ) 0 t ( t ς ) η 1 ( H 2 ( ς , V n 1 ) H 2 ( ς , V n 2 ) ) d ς , ψ 3 n ( t ) = E n ( t ) E n 1 ( t ) = 1 Γ ( η ) 0 t ( t ς ) η 1 ( H 3 ( ς , E n 1 ) H 3 ( ς , E n 2 ) ) d ς , ψ 4 n ( t ) = C n ( t ) C n 1 ( t ) = 1 Γ ( η ) 0 t ( t ς ) η 1 ( H 4 ( ς , C n 1 ) H 4 ( ς , C n 2 ) ) d ς , ψ 5 n ( t ) = I n ( t ) I n 1 ( t ) = 1 Γ ( η ) 0 t ( t ς ) η 1 ( H 5 ( ς , I n 1 ) H 5 ( ς , I n 2 ) ) d ς , ψ 6 n ( t ) = R n ( t ) R n 1 ( t ) = 1 Γ ( η ) 0 t ( t ς ) η 1 ( H 6 ( ς , R n 1 ) H 6 ( ς , R n 2 ) ) d ς , ψ 7 n ( t ) = B n ( t ) B n 1 ( t ) = 1 Γ ( η ) 0 t ( t ς ) η 1 ( H 7 ( ς , B n 1 ) H 7 ( ς , B n 2 ) ) d ς ,
with S 0 ( t ) = S ( 0 ) , V 0 ( t ) = V ( 0 ) , E 0 ( t ) = E ( 0 ) , C 0 ( t ) = C ( 0 ) , I 0 ( t ) = I ( 0 ) , R 0 ( t ) = R ( 0 ) , and B 0 ( t ) = B ( 0 ) .
The norm of ϕ 1 n ( t ) gives
ψ 1 n ( t ) = S n ( t ) S n 1 ( t ) = 1 Γ ( η ) 0 t ( t ς ) η 1 ( H 1 ( ς , S n 1 ) H 1 ( ς , S n 2 ) ) d ς 1 Γ ( η ) 0 t ( t ς ) η 1 ( H 1 ( ς , S n 1 ) H 1 ( ς , S n 2 ) ) d ς .
With the Lipschitz condition (13), we obtain
ψ 1 n ( t ) = S n ( t ) S n 1 ( t ) 1 Γ ( η ) W 1 0 t ( t ς ) η 1 S n 1 S n 2 d ς .
Thus, we have
ψ 1 n ( t ) 1 Γ ( η ) W 1 0 t ( t ς ) η 1 ψ 1 ( n 1 ) ( ς ) d ς .
By proceeding in a similar way, we obtain, for the other ϕ i n ( t ) , i = 2 , , 7 ,
ψ 2 n ( t ) 1 Γ ( η ) W 2 0 t ( t ς ) η 1 ψ 2 ( n 1 ) ( ς ) d ς , ψ 3 n ( t ) 1 Γ ( η ) W 3 0 t ( t ς ) η 1 ψ 3 ( n 1 ) ( ς ) d ς , ψ 4 n ( t ) 1 Γ ( η ) W 4 0 t ( t ς ) η 1 ψ 4 ( n 1 ) ( ς ) d ς , ψ 5 n ( t ) 1 Γ ( η ) W 5 0 t ( t ς ) η 1 ψ 5 ( n 1 ) ( ς ) d ς , ψ 6 n ( t ) 1 Γ ( η ) W 6 0 t ( t ς ) η 1 ψ 6 ( n 1 ) ( ς ) d ς , ψ 7 n ( t ) 1 Γ ( η ) W 7 0 t ( t ς ) η 1 ψ 7 ( n 1 ) ( ς ) d ς .
Each n th term of the state variables of (5) is given by
S n ( t ) = i = 1 n ψ 1 i ( t ) , V n ( t ) = i = 1 n ψ 2 i ( t ) , E n ( t ) = i = 1 n ψ 3 i ( t ) , C n ( t ) = i = 1 n ψ 4 i ( t ) , I n ( t ) = i = 1 n ψ 5 i ( t ) , R n ( t ) = i = 1 n ψ 6 i ( t ) , B n ( t ) = i = 1 n ψ 7 i ( t ) .
The following result guarantees that the solution of the fractional model (5) is unique.
Theorem 4.
If the following inequality holds,
1 Γ ( η ) x η W i < 1 , for i = 1 , 2 , , 7 ,
then the solution of the fractional model (5), for t [ 0 , T ] , is unique.
Proof. 
With Equation (13), inequalities (19) and (20), combined with a recursive technique, we obtain:
ψ 1 n ( t ) S 0 ( t ) W 1 Γ ( η ) x η n , ψ 2 n ( t ) W 0 ( t ) W 2 Γ ( η ) x η n , ψ 3 n ( t ) E 0 ( t ) W 3 Γ ( η ) x η n , ψ 4 n ( t ) C 0 ( t ) W 4 Γ ( η ) x η n , ψ 5 n ( t ) I 0 ( t ) W 5 Γ ( η ) x η n , ψ 6 n ( t ) R 0 ( t ) W 6 Γ ( η ) x η n , ψ 7 n ( t ) B 0 ( t ) W 7 Γ ( η ) x η n .
Therefore, the above sequences satisfy lim n ψ i n ( t ) 0 , j = 1 , 2 , , 7 . The triangle inequality applied to (23) permits us to obtain
S k + n ( t ) S n ( t ) j = n + 1 k + n Z 1 j = Z 1 n + 1 Z 1 k + n + 1 1 Z 1 , V k + n ( t ) V n ( t ) j = n + 1 k + n Z 2 j = Z 2 n + 1 Z 2 k + n + 1 1 Z 2 , E k + n ( t ) E n ( t ) j = n + 1 k + n Z 3 j = Z 3 n + 1 Z 3 k + n + 1 1 Z 3 , C k + n ( t ) C n ( t ) j = n + 1 k + n Z 4 j = Z 4 n + 1 Z 4 k + n + 1 1 Z 4 , I k + n ( t ) I n ( t ) j = n + 1 k + n Z 5 j = Z 5 n + 1 Z 5 k + n + 1 1 Z 5 , R k + n ( t ) R n ( t ) j = n + 1 k + n Z 6 j = Z 6 n + 1 Z 6 k + n + 1 1 Z 6 , B k + n ( t ) B n ( t ) j = n + 1 k + n Z 7 j = Z 7 n + 1 Z 7 n + k + 1 1 Z 7 ,
with 1 Γ ( η ) b η W l < 1 and Z l = 1 Γ ( η ) W l x η n , l = 1 , 2 , , 7 .
Therefore, S n , V n , E n , C n , I n , R n , and B n are uniformly convergent Cauchy sequences (see [32]). With n , it follows that the limit of these sequences represents the unique solution to model (5). □

Numerical Scheme of the Fractional Model and Its Stability Analysis

Several methods have been developed to construct numerical schemes for fractional models. One can cite, among others, the Implicit Quadrature method [24], the Approximate Mittag–Leffler method [33], the Predictor Corrector method [34], and the Adams–Bashforth–Moulton method [25]. The choice of the method depends on several factors, such as the amount of information treated [35] and the accuracy order [36]. The numerical scheme proposed in this work is constructed using the Adams–Bashforth–Moulton method.
Let us consider the following general form of a fractional differential equation [37]:
D t η φ ( t ) = h ( t , φ ( t ) ) , 0 t T , φ ( l ) ( 0 ) = φ 0 l , l = 0 , 1 , 2 , , n 1 , w h e r e n = [ η ] ,
which is equivalent to
φ ( t ) = l = 0 n 1 h 0 l t l l ! + 1 Γ ( η ) 0 t ( t ζ ) η 1 h ( ζ , φ ( ζ ) ) d ζ .
For η [ 0 , 1 ] , 0 t T and setting κ = T / N and t m = m κ , for m = 0 , 1 , 2 , , N Z + , the solution of the fractional model is
S 1 + m = S 0 + κ η Γ ( η + 2 ) Λ h + η R 1 + m p + θ V 1 + m p ( ν B 1 + m p + k 1 ) S 1 + m p + κ η Γ ( η + 2 ) j = 0 m a j , 1 + m Λ h + η R j + θ V j ( ν B j + k 1 ) S j , V 1 + m = V 0 + κ η Γ ( η + 2 ) ξ S 1 + m p ( 1 ϵ ) ν B 1 + m p + k 2 V 1 + m p + κ η Γ ( η + 2 ) j = 0 m a j , 1 + m ξ S j ( 1 ϵ ) ν B j + k 2 V j , E 1 + m = E 0 + κ η Γ ( η + 2 ) ν B 1 + m p S 1 + m p + π V 1 + m p k 3 E 1 + m p + κ η Γ ( η + 2 ) j = 0 m a j , 1 + m ν B j S j + π V j k 3 E j , C 1 + m = C 0 + κ η Γ ( η + 2 ) q γ 1 E 1 + m p k 4 C 1 + m p + κ η Γ ( η + 2 ) j = 0 m a j , 1 + m q γ 1 E j k 4 C j , I 1 + m = I 0 + κ η Γ ( η + 2 ) q 1 γ 1 E 1 + m p + p 1 γ 2 C 1 + m p k 5 + σ I 1 + m p + κ η Γ ( η + 2 ) j = 0 m a j , 1 + m q 1 γ 1 E j + p 1 γ 2 C j k 5 + σ I j , R 1 + m = R 0 + κ η Γ ( η + 2 ) p γ 2 C 1 + m p + σ I 1 + m p k 6 R 1 + m p + κ η Γ ( η + 2 ) j = 0 m a j , 1 + m p γ 2 C j + σ I j k 6 R j , B 1 + m = B 0 + κ η Γ ( η + 2 ) p c C 1 + m p + p i I 1 + m p μ b B 1 + m p + κ η Γ ( η + 2 ) j = 0 m a j , 1 + m p c C j + p i I j μ b B j ,
where
S 1 + m p = S 0 + 1 Γ ( η ) j = 0 m b j , 1 + m Λ h + η R j + θ V j ( ν B j + k 1 ) S j , V 1 + m p = V 0 + 1 Γ ( η ) j = 0 m b j , 1 + m ξ S j ( 1 ϵ ) ν B j + k 2 V j , E 1 + m p = E 0 + 1 Γ ( η ) j = 0 m b j , 1 + m ν B j S j + π V j k 3 E j , C 1 + m p = C 0 + 1 Γ ( η ) j = 0 m b j , 1 + m q γ 1 E j k 4 C j , I 1 + m p = I 0 + 1 Γ ( η ) j = 0 m b j , 1 + m = q 1 γ 1 E j + p 1 γ 2 C j k 5 + σ I j , R 1 + m p = R 0 + 1 Γ ( η ) j = 0 m b j , 1 + m p γ 2 C j + σ I j k 6 R j , B 1 + m p = B 0 + 1 Γ ( η ) j = 0 m b j , 1 + m p c C j + p i I j μ b B j ,
and
a j , 1 + m = m η + 1 ( m + 1 ) ( η + m ) , j = 0 , ( 2 j + m ) η + 1 2 ( m j + 1 ) 1 + η + ( j + m ) 1 + η 1 j m , 1 , j = 1 + m ,
b j , 1 + m = κ η η ( m + 1 j ) η ( m j ) η , 0 j m .
We then claim the following result.
Theorem 5.
Under some conditions, the above numerical scheme (see Equations (27) and (28)) is stable.
Proof. 
Let S 0 , S j ( j = 0 , , 1 + m ) and S 1 + m p ( m = 0 , , N 1 ) be perturbations of S 0 , S j , and S 1 + m p , respectively. By using Equations (19) and (28), the following perturbation equations are obtained:
S 1 + m p = S 0 + 1 Γ ( η ) j = 0 m b j , 1 + m ( G 1 ( t j , S j + S j ) G 1 ( t j , S j ) ) ,
S 1 + m = S 0 + κ η Γ ( η + 2 ) ( G 1 ( t 1 + m , S 1 + m p + S 1 + m P ) G 1 ( t 1 + m , S 1 + m p ) ) + κ η Γ ( η + 2 ) j = 0 m a j , 1 + m ( G 1 ( t j , S j + S j ) G 1 ( t j , S j ) ) .
The Lipschitz condition permits us to obtain
| S 1 + m | ϕ 0 + κ η M Γ ( η + 2 ) | S 1 + m p | + j = 1 m a j , 1 + m | S j | ,
where ϕ 0 = max 0 m N | S 0 | + κ η M a m , 0 Γ ( η + 2 ) | S 0 | . From [32] (Equation (3.18)), it follows that
| S 1 + m p | Θ 0 + M Γ ( η ) j = 1 m b j , 1 + m | S j | ,
where Θ 0 = max 0 m N | S 0 | + M b m , 0 Γ ( η ) | S 0 | . Substituting | S 1 + m p | from Equation (32) into Equation (31) gives
| S 1 + m | ϱ 0 + κ η M Γ ( η + 2 ) M Γ ( η ) j = 1 m b j , 1 + m | S j | + j = 1 m a j , 1 + m | S j | ϱ 0 + κ η M Γ ( η + 2 ) j = 1 m M Γ ( η ) b j , 1 + m + a j , 1 + m | S j | ϱ 0 + κ η M C η , 2 Γ ( η + 2 ) j = 1 m ( m + 1 j ) η 1 | S j | ,
where ϱ 0 = max { ϕ 0 + κ η M a 1 + m , 1 + m Γ ( η + 2 ) Θ 0 } .
Thanks to Lemma 2, we have that C η , 2 > 0 and depends only on η , and κ is assumed to be small enough. A direct application of Lemma 1 implies | S 1 + m | C ϱ 0 . The proof for the other variables is obtained in the same way. This ends the proof. □

2.3. Model Dynamics with the Standard Incidence Law

In this section, we extend model (2) by replacing the mass action incidence law with the standard incidence law, and considering direct transmission (human to human). The new typhoid fever transmission dynamics model is thus presented as follows:
S ˙ ( t ) = Λ h + θ V ( t ) + α R ( t ) k 1 S ( t ) β I + C N ( t ) S ( t ) ν B ( t ) B ( t ) + K S ( t ) ,
V ˙ ( t ) = k 2 + π β C + I N ( t ) + π ν B ( t ) K + B ( t ) V ( t ) + ξ S ( t ) ,
E ˙ ( t ) = β C + I N ( t ) + ν B ( t ) B ( t ) + K π V ( t ) + S ( t ) k 3 E ( t ) ,
C ˙ ( t ) = q γ 1 E ( t ) k 4 C ( t ) ,
I ˙ ( t ) = q 1 γ 1 E ( t ) + p 1 γ 2 C ( t ) k 8 I ( t ) ,
R ˙ ( t ) = p γ 2 C ( t ) + σ I ( t ) k 6 R ( t ) ,
B ˙ ( t ) = p c C ( t ) + p i I ( t ) μ b B ( t ) ,
where β is the direct transmission rate, and K represents the half-saturation constant.
The corresponding fractional model is given by
t 0 C D t η S ( t ) = Λ h + θ V ( t ) + α R ( t ) k 1 S ( t ) β I + C N ( t ) S ( t ) ν B ( t ) K + B ( t ) S ( t ) ,
t 0 C D t η V ( t ) = k 2 + π β I + C N ( t ) + π ν B ( t ) K + B ( t ) V ( t ) + ξ S ( t ) ,
t 0 C D t η E ( t ) = β I + C N ( t ) + ν B ( t ) K + B ( t ) π V ( t ) + S ( t ) k 3 E ( t ) ,
t 0 C D t η C ( t ) = q γ 1 E ( t ) k 4 C ( t ) ,
t 0 C D t η I ( t ) = q 1 γ 1 E ( t ) + p 1 γ 2 C ( t ) k 8 I ( t ) ,
t 0 C D t η R ( t ) = p γ 2 C ( t ) + σ I ( t ) k 6 R ( t ) ,
t 0 C D t η B ( t ) = p c C ( t ) + p i I ( t ) μ b B ( t ) .
Without loss of generality, it is evident that the new model (34) (resp. (35)) is also defined in W .
Model (34) has the same disease-free equilibrium as model system (2). Using the same approach developed in [29], we define the next-generation matrix of model system (34) as
N G M = R 1 H 0 p 1 γ 2 β N 0 k 4 k 8 + H 0 β N 0 k 4 H 0 β N 0 k 8 R 4 0 0 0 0 0 0 0 0 R 5 p 1 γ 2 p i k 4 k 8 + p c k 4 p i k 8 0
where R 1 = H 0 β γ 1 N 0 k 3 k 4 p 1 γ 2 q + q 1 k 4 k 8 + q , R 4 = H 0 ν K μ b , R 5 = γ 1 k 3 k 4 p i p 1 γ 2 q + q 1 k 4 k 8 + p c q , and H 0 = S 0 + π V 0 .
Thus, the control reproduction number of model (34), which is the spectral radius of NGM, is given by
R c = R 1 + R 1 2 + 4 R 4 R 5 2 .
The following result is a direct consequence of Theorem 2 in [29] (see Appendix A for the proof).
Proposition 3.
For model (34) (resp. (35)), Q 0 is locally asymptotically stable in W if R c < 1 and unstable if R c > 1 .
Theorem 6.
For the new typhoid model (34) (resp. (35)), the disease-free equilibrium Q 0 is globally asymptotically stable if R c < 1 and unstable otherwise.
Proof. 
See Appendix B. □

3. Results

3.1. Numerical Results of the Fractional Model with Mass Action Incidence Law

We perform several simulations with the parameter values listed in Table 1.
The step size is κ = 10 8 , the initial time is T = 0, and the final time is T > 0 . We begin by showing the impact of fractional order on the dynamics of the disease. To this aim, the fractional-order parameter varies between η = 1 and η = 0.5 .
Figure 1 displays the impact of the Caputo fractional operator on the model dynamics. For different values of the fractional-order parameter η , the infected state profiles are drawn. From Figure 1, it follows that when the fractional order η decreases, the solutions of our fractional model (5) have different behaviors. Indeed, when the fractional order decreases, the number of infected humans in latent, carrier, and symptomatic states increases. This is the same for the compartment B. This phenomenon was also observed in a malaria fractional model studied by [38]. It is important to note that for η = 1 , the solutions of the fractional model converge to the solutions of the integer model.
To evaluate the impact of vaccination on typhoid fever transmission dynamics, we fix the vaccine efficacy at ϵ = 70 % while the vaccine coverage parameter varies between ξ = 0 % and ξ = 90 % ( ξ { 0.90 , 0.50 , 0.20 , 0 } ), with different values of the fractional-order η ( η { 1 , 0.90 , 0.80 , 0.70 , 0.60 , 0.50 } ). The results are displayed in Figure 2, Figure 3, Figure 4 and Figure 5.
Without vaccination coverage ξ = 0 (Figure 2), the peak delayed phenomenon is observable. Indeed, for η = 1 , the peak date corresponds to T = 45 months, with approximately 3000 infected individuals in the latent stage, 400 asymptomatic individuals, 950 symptomatic individuals, and 26,000 free salmonella in the environment. This peak date is delayed when the fractional-order parameter η decreases. Thus, the peak dates are beyond T = 45 months. From Figure 3, Figure 4 and Figure 5 we note that this peak date is forward delayed beyond T = 45 months whenever the fractional-order parameter η decreases. Moreover, the number of infected individuals decreases with the decrease in the fractional-order parameter.
Now, in addition to vaccination, we consider environmental sanitation. To this aim, the bacterial decay rate μ b is modified to μ b : = μ b + ω , where ω { 0 , 0.1 , 0.2 . 0.3 , 0.4 } represents the additional decay rate of free salmonella due to environmental sanitation [6]. The vaccination coverage is fixed at ξ = 32.21 % as reported in Table 1. From Figure 6, Figure 7, Figure 8 and Figure 9 it is evident that mass vaccination combined with environmental sanitation has a positive impact, reducing the disease burden.

3.2. Numerical Results of the Fractional Model with Standard Incidence Law

To simulate the new fractional model (35), we use the parameter values listed in Table 2. Note that the new typhoid model has been calibrated using real data from Mbandjock, Cameroon (see [5,6]). In Figure 10, panel (a) shows the cumulative typhoid cases versus fitted confirmed cases (infectious individuals tested positive), which is equal to ( 1 q ) γ 1 E ( t ) + ( 1 p ) γ 2 C ( t ) , while panel (b) presents the cumulative estimated cases for the next year. The following fractions are used as initial conditions S ( 0 ) = 20,950/32,000, V ( 0 ) = 20/32,000, E ( 0 ) = 200/32,000, C ( 0 ) = 150/32,000, I ( 0 ) = 60/32,000, R ( 0 ) = 1/32,000, and B ( 0 ) = 500 / 10 6 . The relative change is r = 1.83 × 10 7 and the function tolerance is equal to 10 6 .
First, we observe the general dynamics of the new fractional model. The results are displayed in Figure 11. As for the case of the fractional model with mass incidence law (5), Figure 11 reveals that when the fractional order η decreases, the solutions of our fractional model (35) have different behaviors. The number of typhoid cases decreases and the peak is delayed when the fractional order decreases.
Vaccination coverage impact is studied numerically. From Figure 12, Figure 13, Figure 14 and Figure 15, it follows that the more ξ increases, the fewer individuals are infected. This shows that mass vaccination plays a important role in reducing the spread of the disease.

4. Discussion and Conclusions

In this work, we extended our previous S V E I R -B compartmental model [5] by replacing the integer derivative with fractional derivatives, to evaluate the memory effect on the transmission dynamics of typhoid fever. We began by recalling some previous results on the integer model (the control reproduction number R c , existence and stability of equilibrium points). In order to describe the non-local character as well as long-term memory effects in the typhoid fever transmission dynamics, we replaced the integer derivative with the fractional derivative in the Caputo sense and studied the asymptotic stability of the disease-free equilibrium. Using fixed point theory, we proved the existence as well as the uniqueness of the solutions of the fractional model. We used the Adams–Bashforth method to construct the numerical scheme of the proposed fractional model. We then established the stability of this proposed numerical scheme. We simulated our fractional model using the Adams–Bashforth–Moulton scheme implemented by [39]. Using parameter values for Mbandjock, a city in the central region of Cameroon, we simulates the model by varying the fractional-order parameter, the vaccination coverage, and the bacterial decay rate. Apart from the fact that the solutions of the fractional model converged to the solutions of the integer model when the fractional-order approached one ( η = 1 ), the simulation results showed that the expected date of the disease peak was forward delayed when the fractional-order parameter decreased. In addition, combining vaccination with environmental sanitation can permit a considerable reduction in the disease’s spread.
We then extended the previous models by replacing the mass action incidence law with the standard incidence. The analysis of the models showed that the disease-free equilibrium is also globally asymptotically stable whenever the corresponding reproduction number R c is less than one. Due to the complexity of the newly proposed models, we could not prove the existence and uniqueness of the endemic equilibrium. However, numerical simulations showed that it is possible that the new typhoid fever models permit a unique endemic equilibrium that is globally stable whenever R c > 1 , and no equilibrium otherwise. We also found that, from a quantitative point of view, the disease burden was overestimated with the models the with mass incidence law compared to the one with the standard incidence law. Indeed, for the models with mass incidence, the control reproduction number was estimated at 2.4750, while the one with the standard incidence was estimated at 1.4348. This was in accordance with our previous work in which we considered the standard incidence law. In [6], the control reproduction number was estimated at 1.3722. As for the models with mass action incidences, we observed a delay in the disease peaks whenever the fractional-order derivative decreased.
It was observed that mass vaccination can overcome this disease. In fact, if the means are put in place to finance and implement vaccination campaigns in rural areas, it is possible to eradicate typhoid fever. Moreover, these vaccination campaigns must be accompanied by awareness campaigns among the population in order to combat this type of disease, as well as instructing citizens on ways to protect their environment against the proliferation of salmonella.
Our main contribution in this paper consisted in the formulation, using both integer and fractional derivatives, of new transmission dynamics typhoid fever models that incorporate the standard incidence rates and mass vaccination. The values of the control reproductive number differ from the model with mass action incidence and those with the standard incidences. Indeed, for the model with mass action incidences, R c = 2.4750 , while, for those with standard incidences, R c * = 1.4348 . This proves that mass action incidence overestimates the disease burden.

Author Contributions

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

Funding

The APC was funded by Kottakkaran Sooppy Nisar. No other funds have been received for this work.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to calibrate our model were taken from Mbandjock district hospital.

Acknowledgments

The authors thank the LESIA laboratory manager of the National School of Agro-industrial Sciences of the University of Ngaoundéré for their hospitality during the performance of the numerical simulations. The authors thank the Handling Editor and the anonymous reviewers for their comments and suggestions, which enabled us to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proof of Proposition 3

Proof. 
The Jacobian matrix of (34) (resp. (35)) evaluated at the disease-free equilibrium Q 0 is given by
J ( Q 0 ) = k 1 θ 0 S 0 β N 0 S 0 β N 0 α S 0 ν K ξ k 2 0 π V 0 β N 0 π V 0 β N 0 0 π V 0 ν K 0 0 k 3 H 0 β N 0 H 0 β N 0 0 H 0 ν K 0 0 γ 1 q k 4 0 0 0 0 0 q 1 γ 1 p 1 γ 2 k 8 0 0 0 0 0 p γ 2 σ k 6 0 0 0 0 p c p i 0 μ b = J 1 J 2 J 3 J 4 ,
where J 1 = k 1 θ ξ k 2 , J 4 = k 3 H 0 β N 0 H 0 β N 0 0 H 0 ν K γ 1 q k 4 0 0 0 q 1 γ 1 p 1 γ 2 k 8 0 0 0 p γ 2 σ k 6 0 0 p c p i 0 μ b ,
J 2 = 0 S 0 β N 0 S 0 β N 0 α S 0 ν K 0 π V 0 β N 0 π V 0 β N 0 0 π V 0 ν K , and J 3 = 0 R 5 × 2 . The eigenvalues of J ( Q 0 ) are those of J 1 and J 4 . It is evident that the eigenvalues of J 1 have negative real parts. Indeed, the characteristic polynomial of J 1 is T ( x ) = det ( J 1 x I 2 ) = x 2 + ( k 1 + k 2 ) x + k 7 . Since all its coefficients are positive, it follows that all its roots have negative real parts. A trivial eigenvalue of J 4 is x = k 6 . The others are the roots of the following polynomial: I ( x ) = x 4 + a 1 x 3 + a 2 x 2 + a 3 x + a 4 , with a 1 = μ b + k 8 + k 4 + k 3 , a 4 = k 3 k 4 k 8 μ b 1 R c R c R 1 + 1 ,
a 2 = 1 k 8 q + p 1 γ 2 q + q 1 k 4 k 8 2 μ b q + k 4 k 8 μ b q + k 3 k 8 μ b q + p 1 γ 2 k 8 μ b q + p 1 γ 2 k 4 μ b q + p 1 γ 2 k 3 μ b q + k 4 k 8 2 q + k 3 k 8 2 q + k 3 k 4 k 8 q 1 R 1 + p 1 γ 2 k 4 k 8 q + p 1 γ 2 k 3 k 8 q + p 1 γ 2 k 3 k 4 q + q 1 k 4 k 8 μ b + q 1 k 4 2 μ b + q 1 k 3 k 4 μ b + q 1 k 4 2 k 8 + q 1 k 3 k 4 k 8 1 R 1 + q 1 k 3 k 4 2 ,
and
a 3 = 1 ( ( p 1 γ 2 k 8 + p 1 2 γ 2 2 ) p i + ( k 8 2 + p 1 γ 2 k 8 ) p c ) q 2 + ( ( q 1 k 4 k 8 + 2 p 1 q 1 γ 2 k 4 ) p i + q 1 k 4 k 8 p c ) q + q 1 2 k 4 2 p i × × ( ( ( ( p 1 γ 2 k 4 + p 1 γ 2 k 3 ) k 8 2 + ( ( ( 1 R 1 ) p 1 γ 2 k 3 + p 1 2 γ 2 2 ) k 4 + p 1 2 γ 2 2 k 3 ) k 8 + p 1 2 γ 2 2 k 3 k 4 ) μ b + ( 1 R 1 ) p 1 γ 2 k 3 k 4 k 8 2 + ( 1 R 1 ) p 1 2 γ 2 2 k 3 k 4 k 8 ) p i + K 1 p c ) q 2 + ( K 2 p i + ( q 1 k 3 k 4 2 k 8 μ b 1 R c 2 + R 1 R c + q 1 k 4 2 k 8 2 μ b + ( 1 R 1 ) q 1 k 3 k 4 k 8 2 μ b + ( 1 R 1 ) q 1 k 3 k 4 2 k 8 2 ) p c ) q + ( q 1 2 k 3 k 4 2 k 8 μ b 1 R c ( 1 R 1 + R c ) + q 1 2 k 4 3 k 8 μ b + q 1 2 k 3 k 4 3 μ b + ( 1 R 1 ) q 1 2 k 3 k 4 3 k 8 ) p i .
K 1 = k 8 1 R c 1 R 1 + R c + p 1 γ 2 1 R c 2 + R 1 R c k 3 k 4 k 8 μ b + ( k 4 + k 3 ) ( k 8 + p 1 γ 2 ) k 8 2 μ b + ( 1 R 1 ) k 8 + p 1 γ 2 k 3 k 4 k 8 2
K 2 = ( k 8 + p 1 γ 2 ) q 1 k 3 k 4 k 8 μ b 1 R c 2 + R 1 R c + ( 1 R 1 ) p 1 γ 2 q 1 k 3 k 4 k 8 μ b + ( ( 1 R 1 ) q 1 k 3 + 2 p 1 γ 2 q 1 ) k 4 2 k 8 μ b 2 + 2 p 1 γ 2 q 1 k 3 k 4 2 μ b + q 1 k 4 2 k 8 2 μ b + q 1 k 3 k 4 2 k 8 ( 1 R 1 ) k 8 + 2 p 1 γ 2
It is clear that a 1 is always positive, and a i , i 2 , 3 , 4 are positive if R c < 1 . Indeed, it is important to note that
R c < 1 R 1 < 1 ,
which implies that K 1 > 0 and K 2 > 0 .
Thus, all coefficients of the polynomial I ( x ) are always positive whenever R c < 1 . It follows that, if R c < 1 , then the disease-free equilibrium is locally asymptotically stable if and only if the following conditions hold (because of the length of the expressions, we omit them here):
a 1 a 2 a 3 > 0 and a 1 a 2 a 3 a 1 2 a 4 a 3 2 > 0 .
This ends the proof. □
It remains now to prove the corresponding result for the new fractional model (35). To this aim, let us define the following equation:
det r I ( 1 η ) J ( Q 0 ) η J ( Q 0 ) = 0 ,
which is the characteristic equation of
J ( Q 0 ) = k 1 ϑ 0 S 0 β N 0 S 0 β N 0 0 S 0 ν K ξ k 2 0 V 0 β π N 0 V 0 β π N 0 0 V 0 ν π N 0 0 0 k 3 H 0 β N 0 H 0 β N 0 0 H 0 ν K 0 0 γ 1 q k 4 0 0 0 0 0 q 1 γ 1 p 1 γ 2 k 8 0 0 0 0 0 γ 2 p σ k 6 0 0 0 0 p c p i 0 μ b .
From [4,30], it follows that Q 0 is asymptotically stable, for the new fractional model, if all solutions of (A3) have negative real parts.
Setting D : = r I ( 1 η ) J ( Q 0 ) η Q 0 = D 1 0 R 5 × 2 D 4 , with
D 1 : = η 1 k 1 + 1 r + k 1 η η 1 r ϑ η ϑ η 1 r ξ η ξ η 1 k 2 + 1 r + k 2 η and
D 4 : = η 1 k 3 + 1 r + k 3 η H 0 η 1 β r N 0 H 0 β η N 0 H 0 η 1 β r N 0 H 0 β η N 0 0 H 0 η 1 ν r K H 0 η ν K η 1 γ 1 q r γ 1 η q η 1 k 4 + 1 r + k 4 η 0 0 0 η 1 q 1 γ 1 r q 1 γ 1 η η 1 p 1 γ 2 r p 1 γ 2 η η 1 k 8 + 1 r + k 8 η 0 0 0 η 1 γ 2 p r γ 2 η p η 1 r σ η σ η 1 k 6 + 1 r + k 6 η 0 0 η 1 p c r p c η η 1 p i r η p i 0 η 1 μ b + 1 r + μ b η ,
it follows that the solutions of (A3) are the solutions of det ( D 1 ) = 0 and det ( D 4 ) = 0 . From the Proof of Theorem 2, it follows that the solutions of det ( D 1 ) = 0 have negative real parts. It thus remains to show that the same is true for det ( D 4 ) = 0 . Note that r = k 8 η η 1 k 8 + 1 < 0 is a solution of det ( D 4 ) = 0 . The others are the solutions of det ( D 4 ) = 0 , where
D 4 : = η 1 k 3 + 1 r + k 3 η H 0 η 1 β r N 0 H 0 β η N 0 H 0 η 1 β r N 0 H 0 β η N 0 H 0 η 1 ν r K H 0 η ν K η 1 γ 1 q r γ 1 η q η 1 k 4 + 1 r + k 4 η 0 0 η 1 q 1 γ 1 r q 1 γ 1 η η 1 p 1 γ 2 r p 1 γ 2 η η 1 k 8 + 1 r + k 8 η 0 0 η 1 p c r p c η η 1 p i r η p i η 1 μ b + 1 r + μ b η .
After some straightforward algebraic computations, we obtain that
det ( D 4 ) = 0 r 4 + A 2 A 1 r 3 + A 3 A 1 r 2 + A 4 A 1 r + A 5 A 1 = 0 ,
where
A 1 = ( ( ( ( p 1 γ 2 k 3 k 4 k 8 2 p 1 2 γ 2 2 k 3 k 4 k 8 ) μ b R c 2 + ( R 1 p 1 γ 2 k 3 k 4 k 8 2 + R 1 p 1 2 γ 2 2 k 3 k 4 k 8 ) μ b R c + ( ( 1 R 1 ) p 1 γ 2 k 3 k 4 k 8 2 + ( 1 R 1 ) p 1 2 γ 2 2 k 3 k 4 k 8 ) μ b ) η 4 + ( ( ( p 1 γ 2 k 4 + p 1 γ 2 k 3 ) k 8 2 + ( ( ( 1 R 1 ) p 1 γ 2 k 3 + p 1 2 γ 2 2 ) k 4 + p 1 2 γ 2 2 k 3 ) k 8 + p 1 2 γ 2 2 k 3 k 4 ) μ b + ( 1 R 1 ) p 1 γ 2 k 3 k 4 k 8 2 + ( 1 R 1 ) p 1 2 γ 2 2 k 3 k 4 k 8 ) η 3 + ( ( p 1 γ 2 k 8 2 + ( p 1 γ 2 k 4 + p 1 γ 2 k 3 + p 1 2 γ 2 2 ) k 8 + p 1 2 γ 2 2 k 4 + p 1 2 γ 2 2 k 3 ) μ b + ( p 1 γ 2 k 4 + p 1 γ 2 k 3 ) k 8 2 + ( ( ( 1 R 1 ) p 1 γ 2 k 3 + p 1 2 γ 2 2 ) k 4 + p 1 2 γ 2 2 k 3 ) k 8 + p 1 2 γ 2 2 k 3 k 4 ) η 2 + ( ( p 1 γ 2 k 8 + p 1 2 γ 2 2 ) μ b + p 1 γ 2 k 8 2 + ( p 1 γ 2 k 4 + p 1 γ 2 k 3 + p 1 2 γ 2 2 ) k 8 + p 1 2 γ 2 2 k 4 + p 1 2 γ 2 2 k 3 ) η + p 1 γ 2 k 8 + p 1 2 γ 2 2 ) p i + ( ( k 3 k 4 k 8 3 p 1 γ 2 k 3 k 4 k 8 2 ) μ b R c 2 + ( R 1 k 3 k 4 k 8 3 + R 1 p 1 γ 2 k 3 k 4 k 8 2 ) μ b R c + ( ( 1 R 1 ) k 3 k 4 k 8 3 + ( 1 R 1 ) p 1 γ 2 k 3 k 4 k 8 2 ) μ b ) p c η 4 + ( ( k 3 k 4 k 8 2 p 1 γ 2 k 3 k 4 k 8 ) μ b R c 2 + ( R 1 k 3 k 4 k 8 2 + R 1 p 1 γ 2 k 3 k 4 k 8 ) μ b R c + ( ( k 4 + k 3 ) k 8 3 + ( ( ( 1 R 1 ) k 3 + p 1 γ 2 ) k 4 + p 1 γ 2 k 3 ) k 8 2 + p 1 γ 2 k 3 k 4 k 8 ) μ b + ( 1 R 1 ) k 3 k 4 k 8 3 + ( 1 R 1 ) p 1 γ 2 k 3 k 4 k 8 2 ) p c η 3 + ( ( k 8 3 + ( k 4 + k 3 + p 1 γ 2 ) k 8 2 + ( p 1 γ 2 k 4 + p 1 γ 2 k 3 ) k 8 ) μ b + ( k 4 + k 3 ) k 8 3 + ( ( ( 1 R 1 ) k 3 + p 1 γ 2 ) k 4 + p 1 γ 2 k 3 ) k 8 2 + p 1 γ 2 k 3 k 4 k 8 ) p c η 2 + ( ( k 8 2 + p 1 γ 2 k 8 ) μ b + k 8 3 + ( k 4 + k 3 + p 1 γ 2 ) k 8 2 + ( p 1 γ 2 k 4 + p 1 γ 2 k 3 ) k 8 ) p c η + ( k 8 2 + p 1 γ 2 k 8 ) p c ) q 2 + ( ( ( ( q 1 k 3 k 4 2 k 8 2 2 p 1 q 1 γ 2 k 3 k 4 2 k 8 ) μ b R c 2 + ( R 1 q 1 k 3 k 4 2 k 8 2 + 2 R 1 p 1 q 1 γ 2 k 3 k 4 2 k 8 ) μ b R c + ( ( 1 R 1 ) q 1 k 3 k 4 2 k 8 2 + ( 2 2 R 1 ) p 1 q 1 γ 2 k 3 k 4 2 k 8 ) μ b ) η 4 + ( ( q 1 k 3 k 4 k 8 2 p 1 q 1 γ 2 k 3 k 4 k 8 ) μ b R c 2 + ( R 1 q 1 k 3 k 4 k 8 2 + R 1 p 1 q 1 γ 2 k 3 k 4 k 8 ) μ b R c + ( ( q 1 k 4 2 + q 1 k 3 k 4 ) k 8 2 + ( ( ( 1 R 1 ) q 1 k 3 + 2 p 1 q 1 γ 2 ) k 4 2 + ( 2 R 1 ) p 1 q 1 γ 2 k 3 k 4 ) k 8 + 2 p 1 q 1 γ 2 k 3 k 4 2 ) μ b + ( 1 R 1 ) q 1 k 3 k 4 2 k 8 2 + ( 2 2 R 1 ) p 1 q 1 γ 2 k 3 k 4 2 k 8 ) η 3 + ( ( q 1 k 4 k 8 2 + ( q 1 k 4 2 + ( q 1 k 3 + 2 p 1 q 1 γ 2 ) k 4 ) k 8 + 2 p 1 q 1 γ 2 k 4 2 + 2 p 1 q 1 γ 2 k 3 k 4 ) μ b + ( q 1 k 4 2 + q 1 k 3 k 4 ) k 8 2 + ( ( ( 1 R 1 ) q 1 k 3 + 2 p 1 q 1 γ 2 ) k 4 2 + ( 2 R 1 ) p 1 q 1 γ 2 k 3 k 4 ) k 8 + 2 p 1 q 1 γ 2 k 3 k 4 2 ) η 2 + ( ( q 1 k 4 k 8 + 2 p 1 q 1 γ 2 k 4 ) μ b + q 1 k 4 k 8 2 + ( q 1 k 4 2 + ( q 1 k 3 + 2 p 1 q 1 γ 2 ) k 4 ) k 8 + 2 p 1 q 1 γ 2 k 4 2 + 2 p 1 q 1 γ 2 k 3 k 4 ) η + q 1 k 4 k 8 + 2 p 1 q 1 γ 2 k 4 ) p i + ( q 1 k 3 k 4 2 k 8 2 μ b R c 2 + R 1 q 1 k 3 k 4 2 k 8 2 μ b R c + ( 1 R 1 ) q 1 k 3 k 4 2 k 8 2 μ b ) p c η 4 + ( q 1 k 3 k 4 2 k 8 μ b R c 2 + R 1 q 1 k 3 k 4 2 k 8 μ b R c + ( ( q 1 k 4 2 + ( 1 R 1 ) q 1 k 3 k 4 ) k 8 2 + q 1 k 3 k 4 2 k 8 ) μ b + ( 1 R 1 ) q 1 k 3 k 4 2 k 8 2 ) p c η 3 + ( ( q 1 k 4 k 8 2 + ( q 1 k 4 2 + q 1 k 3 k 4 ) k 8 ) μ b + ( q 1 k 4 2 + ( 1 R 1 ) q 1 k 3 k 4 ) k 8 2 + q 1 k 3 k 4 2 k 8 ) p c η 2 + ( q 1 k 4 k 8 μ b + q 1 k 4 k 8 2 + ( q 1 k 4 2 + q 1 k 3 k 4 ) k 8 ) p c η + q 1 k 4 k 8 p c ) q + ( ( q 1 2 k 3 k 4 3 k 8 μ b R c 2 + R 1 q 1 2 k 3 k 4 3 k 8 μ b R c + ( 1 R 1 ) q 1 2 k 3 k 4 3 k 8 μ b ) η 4 + ( q 1 2 k 3 k 4 2 k 8 μ b R c 2 + R 1 q 1 2 k 3 k 4 2 k 8 μ b R c + ( ( q 1 2 k 4 3 + ( 1 R 1 ) q 1 2 k 3 k 4 2 ) k 8 + q 1 2 k 3 k 4 3 ) μ b + ( 1 R 1 ) q 1 2 k 3 k 4 3 k 8 ) η 3 + ( ( q 1 2 k 4 2 k 8 + q 1 2 k 4 3 + q 1 2 k 3 k 4 2 ) μ b + ( q 1 2 k 4 3 + ( 1 R 1 ) q 1 2 k 3 k 4 2 ) k 8 + q 1 2 k 3 k 4 3 ) η 2 + ( q 1 2 k 4 2 μ b + q 1 2 k 4 2 k 8 + q 1 2 k 4 3 + q 1 2 k 3 k 4 2 ) η + q 1 2 k 4 2 ) p i ,
A 4 = ( ( 4 ( k 8 + p 1 γ 2 ) p 1 γ 2 k 3 k 4 k 8 μ b 1 R c ( 1 + R c R 1 ) η 4 + ( ( ( p 1 γ 2 k 4 + p 1 γ 2 k 3 ) k 8 2 + ( ( ( 1 R 1 ) p 1 γ 2 k 3 + p 1 2 γ 2 2 ) k 4 + p 1 2 γ 2 2 k 3 ) k 8 + p 1 2 γ 2 2 k 3 k 4 ) μ b + ( 1 R 1 ) p 1 γ 2 k 3 k 4 k 8 2 + ( 1 R 1 ) p 1 2 γ 2 2 k 3 k 4 k 8 ) η 3 ) p i + 4 ( k 8 + p 1 γ 2 ) ( 1 R c ) ( 1 + R c R 1 ) k 3 k 4 k 8 2 μ b p c η 4 + ( ( k 8 + p 1 γ 2 ) k 3 k 4 k 8 μ b R c R 1 R c + ( ( k 4 + k 3 ) k 8 3 + ( ( ( 1 R 1 ) k 3 + p 1 γ 2 ) k 4 + p 1 γ 2 k 3 ) k 8 2 + p 1 γ 2 k 3 k 4 k 8 ) μ b + ( 1 R 1 ) k 3 k 4 k 8 3 + ( 1 R 1 ) p 1 γ 2 k 3 k 4 k 8 2 ) p c η 3 ) q 2 + ( ( ( ( 4 q 1 k 3 k 4 2 k 8 2 8 p 1 q 1 γ 2 k 3 k 4 2 k 8 ) μ b R c 2 + ( 4 R 1 q 1 k 3 k 4 2 k 8 2 + 8 R 1 p 1 q 1 γ 2 k 3 k 4 2 k 8 ) μ b R c + ( ( 4 4 R 1 ) q 1 k 3 k 4 2 k 8 2 + ( 8 8 R 1 ) p 1 q 1 γ 2 k 3 k 4 2 k 8 ) μ b ) η 4 + ( ( q 1 k 3 k 4 k 8 2 p 1 q 1 γ 2 k 3 k 4 k 8 ) μ b R c 2 + ( R 1 q 1 k 3 k 4 k 8 2 + R 1 p 1 q 1 γ 2 k 3 k 4 k 8 ) μ b R c + ( ( q 1 k 4 2 + q 1 k 3 k 4 ) k 8 2 + ( ( ( 1 R 1 ) q 1 k 3 + 2 p 1 q 1 γ 2 ) k 4 2 + ( 2 R 1 ) p 1 q 1 γ 2 k 3 k 4 ) k 8 + 2 p 1 q 1 γ 2 k 3 k 4 2 ) μ b + ( 1 R 1 ) q 1 k 3 k 4 2 k 8 2 + ( 2 2 R 1 ) p 1 q 1 γ 2 k 3 k 4 2 k 8 ) η 3 ) p i + ( 4 q 1 k 3 k 4 2 k 8 2 μ b R c 2 + 4 R 1 q 1 k 3 k 4 2 k 8 2 μ b R c + ( 4 4 R 1 ) q 1 k 3 k 4 2 k 8 2 μ b ) p c η 4 + ( q 1 k 3 k 4 2 k 8 μ b R c 2 + R 1 q 1 k 3 k 4 2 k 8 μ b R c + ( ( q 1 k 4 2 + ( 1 R 1 ) q 1 k 3 k 4 ) k 8 2 + q 1 k 3 k 4 2 k 8 ) μ b + ( 1 R 1 ) q 1 k 3 k 4 2 k 8 2 ) p c η 3 ) q + ( ( 4 q 1 2 k 3 k 4 3 k 8 μ b R c 2 + 4 R 1 q 1 2 k 3 k 4 3 k 8 μ b R c + ( 4 4 R 1 ) q 1 2 k 3 k 4 3 k 8 μ b ) η 4 + ( q 1 2 k 3 k 4 2 k 8 μ b R c 2 + R 1 q 1 2 k 3 k 4 2 k 8 μ b R c + ( ( q 1 2 k 4 3 + ( 1 R 1 ) q 1 2 k 3 k 4 2 ) k 8 + q 1 2 k 3 k 4 3 ) μ b + ( 1 R 1 ) q 1 2 k 3 k 4 3 k 8 ) η 3 ) p i ,
A 5 = k 3 k 4 k 8 μ b R c R 1 + 1 η 4 k 8 q + p 1 γ 2 q + q 1 k 4 p 1 γ 2 p i q + k 8 p c q + q 1 k 4 p i 1 R c ,
A 2 = ( ( ( ( 4 p 1 γ 2 k 3 k 4 k 8 2 4 p 1 2 γ 2 2 k 3 k 4 k 8 ) μ b R c 2 + ( 4 R 1 p 1 γ 2 k 3 k 4 k 8 2 + 4 R 1 p 1 2 γ 2 2 k 3 k 4 k 8 ) μ b R c + ( ( 4 4 R 1 ) p 1 γ 2 k 3 k 4 k 8 2 + ( 4 4 R 1 ) p 1 2 γ 2 2 k 3 k 4 k 8 ) μ b ) η 4 + ( ( ( 3 p 1 γ 2 k 4 + 3 p 1 γ 2 k 3 ) k 8 2 + ( ( ( 3 3 R 1 ) p 1 γ 2 k 3 + 3 p 1 2 γ 2 2 ) k 4 + 3 p 1 2 γ 2 2 k 3 ) k 8 + 3 p 1 2 γ 2 2 k 3 k 4 ) μ b + ( 3 3 R 1 ) p 1 γ 2 k 3 k 4 k 8 2 + ( 3 3 R 1 ) p 1 2 γ 2 2 k 3 k 4 k 8 ) η 3 + ( ( 2 p 1 γ 2 k 8 2 + ( 2 p 1 γ 2 k 4 + 2 p 1 γ 2 k 3 + 2 p 1 2 γ 2 2 ) k 8 + 2 p 1 2 γ 2 2 k 4 + 2 p 1 2 γ 2 2 k 3 ) μ b + ( 2 p 1 γ 2 k 4 + 2 p 1 γ 2 k 3 ) k 8 2 + ( ( ( 2 2 R 1 ) p 1 γ 2 k 3 + 2 p 1 2 γ 2 2 ) k 4 + 2 p 1 2 γ 2 2 k 3 ) k 8 + 2 p 1 2 γ 2 2 k 3 k 4 ) η 2 + ( ( p 1 γ 2 k 8 + p 1 2 γ 2 2 ) μ b + p 1 γ 2 k 8 2 + ( p 1 γ 2 k 4 + p 1 γ 2 k 3 + p 1 2 γ 2 2 ) k 8 + p 1 2 γ 2 2 k 4 + p 1 2 γ 2 2 k 3 ) η ) p i + ( ( 4 k 3 k 4 k 8 3 4 p 1 γ 2 k 3 k 4 k 8 2 ) μ b R c 2 + ( 4 R 1 k 3 k 4 k 8 3 + 4 R 1 p 1 γ 2 k 3 k 4 k 8 2 ) μ b R c + 4 ( 1 R 1 ) ( k 8 + p 1 γ 2 ) k 3 k 4 k 8 2 μ b ) p c η 4 + ( ( 3 k 3 k 4 k 8 2 3 p 1 γ 2 k 3 k 4 k 8 ) μ b R c 2 + ( 3 R 1 k 3 k 4 k 8 2 + 3 R 1 p 1 γ 2 k 3 k 4 k 8 ) μ b R c + ( 3 ( k 4 + k 3 ) k 8 3 + ( ( 3 ( 1 R 1 ) k 3 + 3 p 1 γ 2 ) k 4 + 3 p 1 γ 2 k 3 ) k 8 2 + 3 p 1 γ 2 k 3 k 4 k 8 ) μ b + 3 ( 1 R 1 ) k 3 k 4 k 8 3 + 3 ( 1 R 1 ) p 1 γ 2 k 3 k 4 k 8 2 ) p c η 3 + ( ( 2 k 8 3 + ( 2 k 4 + 2 k 3 + 2 p 1 γ 2 ) k 8 2 + ( 2 p 1 γ 2 k 4 + 2 p 1 γ 2 k 3 ) k 8 ) μ b + ( 2 k 4 + 2 k 3 ) k 8 3 + ( ( ( 2 2 R 1 ) k 3 + 2 p 1 γ 2 ) k 4 + 2 p 1 γ 2 k 3 ) k 8 2 + 2 p 1 γ 2 k 3 k 4 k 8 ) p c η 2 + ( ( k 8 2 + p 1 γ 2 k 8 ) μ b + k 8 3 + ( k 4 + k 3 + p 1 γ 2 ) k 8 2 + ( p 1 γ 2 k 4 + p 1 γ 2 k 3 ) k 8 ) p c η ) q 2 + ( ( ( ( 4 q 1 k 3 k 4 2 k 8 2 8 p 1 q 1 γ 2 k 3 k 4 2 k 8 ) μ b R c 2 + ( 4 R 1 q 1 k 3 k 4 2 k 8 2 + 8 R 1 p 1 q 1 γ 2 k 3 k 4 2 k 8 ) μ b R c + ( ( 4 4 R 1 ) q 1 k 3 k 4 2 k 8 2 + ( 8 8 R 1 ) p 1 q 1 γ 2 k 3 k 4 2 k 8 ) μ b ) η 4 + ( ( 3 q 1 k 3 k 4 k 8 2 3 p 1 q 1 γ 2 k 3 k 4 k 8 ) μ b R c 2 + ( 3 R 1 q 1 k 3 k 4 k 8 2 + 3 R 1 p 1 q 1 γ 2 k 3 k 4 k 8 ) μ b R c + ( ( 3 q 1 k 4 2 + 3 q 1 k 3 k 4 ) k 8 2 + ( ( ( 3 3 R 1 ) q 1 k 3 + 6 p 1 q 1 γ 2 ) k 4 2 + ( 6 3 R 1 ) p 1 q 1 γ 2 k 3 k 4 ) k 8 + 6 p 1 q 1 γ 2 k 3 k 4 2 ) μ b + ( 3 3 R 1 ) q 1 k 3 k 4 2 k 8 2 + ( 6 6 R 1 ) p 1 q 1 γ 2 k 3 k 4 2 k 8 ) η 3 + ( ( 2 q 1 k 4 k 8 2 + ( 2 q 1 k 4 2 + ( 2 q 1 k 3 + 4 p 1 q 1 γ 2 ) k 4 ) k 8 + 4 p 1 q 1 γ 2 k 4 2 + 4 p 1 q 1 γ 2 k 3 k 4 ) μ b + ( 2 q 1 k 4 2 + 2 q 1 k 3 k 4 ) k 8 2 + ( ( ( 2 2 R 1 ) q 1 k 3 + 4 p 1 q 1 γ 2 ) k 4 2 + ( 4 2 R 1 ) p 1 q 1 γ 2 k 3 k 4 ) k 8 + 4 p 1 q 1 γ 2 k 3 k 4 2 ) η 2 + ( ( q 1 k 4 k 8 + 2 p 1 q 1 γ 2 k 4 ) μ b + q 1 k 4 k 8 2 + ( q 1 k 4 2 + ( q 1 k 3 + 2 p 1 q 1 γ 2 ) k 4 ) k 8 + 2 p 1 q 1 γ 2 k 4 2 + 2 p 1 q 1 γ 2 k 3 k 4 ) η ) p i + 4 ( 1 R 1 + R 1 R c R c 2 ) q 1 k 3 k 4 2 k 8 2 μ b p c η 4 + ( 3 q 1 k 3 k 4 2 k 8 μ b R c 2 + 3 R 1 q 1 k 3 k 4 2 k 8 μ b R c + ( ( 3 q 1 k 4 2 + ( 3 3 R 1 ) q 1 k 3 k 4 ) k 8 2 + 3 q 1 k 3 k 4 2 k 8 ) μ b + ( 3 3 R 1 ) q 1 k 3 k 4 2 k 8 2 ) p c η 3 + ( ( 2 q 1 k 4 k 8 2 + ( 2 q 1 k 4 2 + 2 q 1 k 3 k 4 ) k 8 ) μ b + ( 2 q 1 k 4 2 + ( 2 2 R 1 ) q 1 k 3 k 4 ) k 8 2 + 2 q 1 k 3 k 4 2 k 8 ) p c η 2 + ( q 1 k 4 k 8 μ b + q 1 k 4 k 8 2 + ( q 1 k 4 2 + q 1 k 3 k 4 ) k 8 ) p c η ) q + ( ( 4 q 1 2 k 3 k 4 3 k 8 μ b R c 2 + 4 R 1 q 1 2 k 3 k 4 3 k 8 μ b R c + ( 4 4 R 1 ) q 1 2 k 3 k 4 3 k 8 μ b ) η 4 + ( 3 q 1 2 k 3 k 4 2 k 8 μ b R c 2 + 3 R 1 q 1 2 k 3 k 4 2 k 8 μ b R c + ( ( 3 q 1 2 k 4 3 + ( 3 3 R 1 ) q 1 2 k 3 k 4 2 ) k 8 + 3 q 1 2 k 3 k 4 3 ) μ b + ( 3 3 R 1 ) q 1 2 k 3 k 4 3 k 8 ) η 3 + ( ( 2 q 1 2 k 4 2 k 8 + 2 q 1 2 k 4 3 + 2 q 1 2 k 3 k 4 2 ) μ b + ( 2 q 1 2 k 4 3 + ( 2 2 R 1 ) q 1 2 k 3 k 4 2 ) k 8 + 2 q 1 2 k 3 k 4 3 ) η 2 + ( q 1 2 k 4 2 μ b + q 1 2 k 4 2 k 8 + q 1 2 k 4 3 + q 1 2 k 3 k 4 2 ) η ) p i ,
A 3 = ( ( ( ( 6 p 1 γ 2 k 3 k 4 k 8 2 6 p 1 2 γ 2 2 k 3 k 4 k 8 ) μ b R c 2 + ( 6 R 1 p 1 γ 2 k 3 k 4 k 8 2 + 6 R 1 p 1 2 γ 2 2 k 3 k 4 k 8 ) μ b R c + ( ( 6 6 R 1 ) p 1 γ 2 k 3 k 4 k 8 2 + ( 6 6 R 1 ) p 1 2 γ 2 2 k 3 k 4 k 8 ) μ b ) η 4 + ( ( ( 3 p 1 γ 2 k 4 + 3 p 1 γ 2 k 3 ) k 8 2 + ( ( ( 3 3 R 1 ) p 1 γ 2 k 3 + 3 p 1 2 γ 2 2 ) k 4 + 3 p 1 2 γ 2 2 k 3 ) k 8 + 3 p 1 2 γ 2 2 k 3 k 4 ) μ b + ( 3 3 R 1 ) p 1 γ 2 k 3 k 4 k 8 2 + ( 3 3 R 1 ) p 1 2 γ 2 2 k 3 k 4 k 8 ) η 3 + ( ( p 1 γ 2 k 8 2 + ( p 1 γ 2 k 4 + p 1 γ 2 k 3 + p 1 2 γ 2 2 ) k 8 + p 1 2 γ 2 2 k 4 + p 1 2 γ 2 2 k 3 ) μ b + ( p 1 γ 2 k 4 + p 1 γ 2 k 3 ) k 8 2 + ( ( ( 1 R 1 ) p 1 γ 2 k 3 + p 1 2 γ 2 2 ) k 4 + p 1 2 γ 2 2 k 3 ) k 8 + p 1 2 γ 2 2 k 3 k 4 ) η 2 ) p i + ( ( 6 k 3 k 4 k 8 3 6 p 1 γ 2 k 3 k 4 k 8 2 ) μ b R c 2 + ( 6 R 1 k 3 k 4 k 8 3 + 6 R 1 p 1 γ 2 k 3 k 4 k 8 2 ) μ b R c + ( 6 6 R 1 ) ( k 3 k 4 k 8 3 + p 1 γ 2 k 3 k 4 k 8 2 ) μ b ) p c η 4 + ( ( 3 k 3 k 4 k 8 2 3 p 1 γ 2 k 3 k 4 k 8 ) μ b R c 2 + ( 3 R 1 k 3 k 4 k 8 2 + 3 R 1 p 1 γ 2 k 3 k 4 k 8 ) μ b R c + ( ( 3 k 4 + 3 k 3 ) k 8 3 + ( ( ( 3 3 R 1 ) k 3 + 3 p 1 γ 2 ) k 4 + 3 p 1 γ 2 k 3 ) k 8 2 + 3 p 1 γ 2 k 3 k 4 k 8 ) μ b + 3 ( 1 R 1 ) k 3 k 4 k 8 2 k 8 + p 1 γ 2 ) p c η 3 + ( ( k 8 3 + ( k 4 + k 3 + p 1 γ 2 ) k 8 2 + p 1 γ 2 ( k 4 + k 3 ) k 8 ) μ b + ( k 4 + k 3 ) k 8 3 + ( ( ( 1 R 1 ) k 3 + p 1 γ 2 ) k 4 + p 1 γ 2 k 3 ) k 8 2 + p 1 γ 2 k 3 k 4 k 8 ) p c η 2 ) q 2 + ( ( ( ( 6 q 1 k 3 k 4 2 k 8 2 12 p 1 q 1 γ 2 k 3 k 4 2 k 8 ) μ b R c 2 + ( 6 R 1 q 1 k 3 k 4 2 k 8 2 + 12 R 1 p 1 q 1 γ 2 k 3 k 4 2 k 8 ) μ b R c + ( ( 6 6 R 1 ) q 1 k 3 k 4 2 k 8 2 + ( 12 12 R 1 ) p 1 q 1 γ 2 k 3 k 4 2 k 8 ) μ b ) η 4 + ( ( 3 q 1 k 3 k 4 k 8 2 3 p 1 q 1 γ 2 k 3 k 4 k 8 ) μ b R c 2 + ( 3 R 1 q 1 k 3 k 4 k 8 2 + 3 R 1 p 1 q 1 γ 2 k 3 k 4 k 8 ) μ b R c + ( ( 3 q 1 k 4 2 + 3 q 1 k 3 k 4 ) k 8 2 + ( ( ( 3 3 R 1 ) q 1 k 3 + 6 p 1 q 1 γ 2 ) k 4 2 + ( 6 3 R 1 ) p 1 q 1 γ 2 k 3 k 4 ) k 8 + 6 p 1 q 1 γ 2 k 3 k 4 2 ) μ b + ( 3 3 R 1 ) q 1 k 3 k 4 2 k 8 2 + ( 6 6 R 1 ) p 1 q 1 γ 2 k 3 k 4 2 k 8 ) η 3 + ( ( q 1 k 4 k 8 2 + ( q 1 k 4 2 + ( q 1 k 3 + 2 p 1 q 1 γ 2 ) k 4 ) k 8 + 2 p 1 q 1 γ 2 k 4 2 + 2 p 1 q 1 γ 2 k 3 k 4 ) μ b + ( q 1 k 4 2 + q 1 k 3 k 4 ) k 8 2 + ( ( ( 1 R 1 ) q 1 k 3 + 2 p 1 q 1 γ 2 ) k 4 2 + ( 2 R 1 ) p 1 q 1 γ 2 k 3 k 4 ) k 8 + 2 p 1 q 1 γ 2 k 3 k 4 2 ) η 2 ) p i + ( 6 q 1 k 3 k 4 2 k 8 2 μ b R c 2 + 6 R 1 q 1 k 3 k 4 2 k 8 2 μ b R c + ( 6 6 R 1 ) q 1 k 3 k 4 2 k 8 2 μ b ) p c η 4 + ( 3 q 1 k 3 k 4 2 k 8 μ b R c 2 + 3 R 1 q 1 k 3 k 4 2 k 8 μ b R c + ( ( 3 q 1 k 4 2 + ( 3 3 R 1 ) q 1 k 3 k 4 ) k 8 2 + 3 q 1 k 3 k 4 2 k 8 ) μ b + ( 3 3 R 1 ) q 1 k 3 k 4 2 k 8 2 ) p c η 3 + ( ( q 1 k 4 k 8 2 + ( q 1 k 4 2 + q 1 k 3 k 4 ) k 8 ) μ b + ( q 1 k 4 2 + ( 1 R 1 ) q 1 k 3 k 4 ) k 8 2 + q 1 k 3 k 4 2 k 8 ) p c η 2 ) q + ( ( 6 q 1 2 k 3 k 4 3 k 8 μ b R c 2 + 6 R 1 q 1 2 k 3 k 4 3 k 8 μ b R c + ( 6 6 R 1 ) q 1 2 k 3 k 4 3 k 8 μ b ) η 4 + ( 3 q 1 2 k 3 k 4 2 k 8 μ b R c 2 + 3 R 1 q 1 2 k 3 k 4 2 k 8 μ b R c + ( ( 3 q 1 2 k 4 3 + ( 3 3 R 1 ) q 1 2 k 3 k 4 2 ) k 8 + 3 q 1 2 k 3 k 4 3 ) μ b + ( 3 3 R 1 ) q 1 2 k 3 k 4 3 k 8 ) η 3 + ( ( q 1 2 k 4 2 k 8 + q 1 2 k 4 3 + q 1 2 k 3 k 4 2 ) μ b + ( q 1 2 k 4 3 + ( 1 R 1 ) q 1 2 k 3 k 4 2 ) k 8 + q 1 2 k 3 k 4 3 ) η 2 ) p i .
It is possible to show that all the above coefficients are positive whenever R c < 1 . Then, it follows that, if R c < 1 , then the disease-free equilibrium Q 0 is asymptotically stable whenever the following Routh–Hurwitz criteria, A 1 A 2 A 5 2 A 3 A 5 > 0 and A 1 A 2 A 3 A 5 3 A 1 2 A 4 A 5 3 A 3 2 A 5 2 > 0 , are satisfied for polynomial det ( D 4 ) . (Given the heaviness of these coefficients, we do not present the Routh–Hurwitz conditions here.)

Appendix B. Proof of Theorem 6

Let us rewrite system (34) as
d E d t d C d t d I d t d B d t = k 3 H 0 β N 0 H 0 β N 0 H 0 ν K γ 1 q k 4 0 0 q 1 γ 1 p 1 γ 2 k 8 0 0 p c p i μ b E C I B N ( S , V , E , C , I , R , B ) ,
where N ( S , V , E , C , I , R , B ) = β ( C + I ) H 0 N 0 H N + ν B H 0 K H K + B 0 0 0 .
In W , H = S + π V < H 0 = S 0 + π V 0 for t > 0 ; thus, N ( S , V , E , C , I , R , B ) O R 4 . Note that Proposition 3 ensures that the following matrix
J ( Q 0 ) = k 3 H 0 β N 0 H 0 β N 0 H 0 ν K γ 1 q k 4 0 0 q 1 γ 1 p 1 γ 2 k 8 0 0 p c p i μ b
has all its eigenvalues with negative real parts. It follows that from the comparison theorem [40], E , C , I , B ( 0 , 0 , 0 , 0 ) and S , V , R ( S 0 , V 0 , 0 ) as t + . Thus, S , V , E , C , I , R , B Q 0 = ( S 0 , V 0 , 0 , 0 , 0 , 0 , 0 ) as t + . We finally conclude that the disease–free equilibrium is globally asymptotically stable in W if R c < 1 .

References

  1. World Health Organization. Typhoid vaccines: WHO position paper. Wkly. Epidemiol. Rec. 2008, 83, 49–59. [Google Scholar]
  2. World Health Organization. Typhoid vaccines: WHO position paper–March 2018–Vaccins antityphoïdiques: Note de synthèse de l’OMS–mars 2018. Wkly. Epidemiol. Rec. 2018, 93, 153–172. [Google Scholar]
  3. Ross, R. Some quantitative studies in epidemiology. Nature 1911, 87, 466–467. [Google Scholar] [CrossRef]
  4. Abboubakar, H.; Kumar, P.; Erturk, V.S.; Kumar, A. A mathematical study of a tuberculosis model with fractional derivatives. Int. J. Model. Simul. Sci. Comput. 2021, 12, 2150037. [Google Scholar] [CrossRef]
  5. Abboubakar, H.; Racke, R. Mathematical modeling, forecasting, and optimal control of typhoid fever transmission dynamics. Chaos Solitons Fractals 2021, 149, 111074. [Google Scholar] [CrossRef]
  6. Abboubakar, H.; Kombou, L.K.; Koko, A.D.; Fouda, H.P.E.; Kumar, A. Projections and fractional dynamics of the typhoid fever: A case study of Mbandjock in the Centre Region of Cameroon. Chaos Solitons Fractals 2021, 150, 111129. [Google Scholar] [CrossRef]
  7. Edward, S.; Nyerere, N. Modelling typhoid fever with education, vaccination and treatment. Eng. Math. 2016, 1, 44–52. [Google Scholar]
  8. Mushayabasa, S. Modeling the impact of optimal screening on typhoid dynamics. Int. J. Dyn. Control 2016, 4, 330–338. [Google Scholar] [CrossRef]
  9. Peter, O.J.; Ibrahim, M.O.; Edogbanya, H.O.; Oguntolu, F.A.; Oshinubi, K.; Ibrahim, A.A.; Ayoola, T.A.; Lawal, J.O. Direct and Indirect Transmission of Typhoid Fever Model with Optimal Control. Results Phys. 2021, 27, 104463. [Google Scholar] [CrossRef]
  10. Tilahun, G.T.; Makinde, O.D.; Malonza, D. Modelling and optimal control of typhoid fever disease with cost-effective strategies. Comput. Math. Methods Med. 2017, 2017, 2324518. [Google Scholar] [CrossRef] [Green Version]
  11. Tilahun, G.T.; Makinde, O.D.; Malonza, D. Co-dynamics of pneumonia and typhoid fever diseases with cost effective optimal control analysis. Appl. Math. Comput. 2018, 316, 438–459. [Google Scholar] [CrossRef]
  12. Shaikh, A.S.; Nisar, K.S. Transmission dynamics of fractional order Typhoid fever model using Caputo–Fabrizio operator. Chaos Solitons Fractals 2019, 128, 355–365. [Google Scholar] [CrossRef]
  13. Erturk, V.S.; Kumar, P. Solution of a COVID-19 model via new generalized Caputo-type fractional derivatives. Chaos Solitons Fractals 2020, 139, 110280. [Google Scholar] [CrossRef]
  14. Kumar, P.; Erturk, V.S. Environmental persistence influences infection dynamics for a butterfly pathogen via new generalised Caputo type fractional derivative. Chaos Solitons Fractals 2021, 144, 110672. [Google Scholar] [CrossRef]
  15. Kumar, P.; Erturk, V.S.; Almusawa, H. Mathematical structure of mosaic disease using microbial biostimulants via Caputo and Atangana–Baleanu derivatives. Results Phys. 2021, 24, 104186. [Google Scholar] [CrossRef]
  16. Kumar, P.; Suat Ertürk, V.; Nisar, K.S. Fractional dynamics of huanglongbing transmission within a citrus tree. Math. Methods Appl. Sci. 2021, 44, 11404–11424. [Google Scholar] [CrossRef]
  17. Kumar, P.; Erturk, V.S.; Murillo-Arcila, M. A complex fractional mathematical modeling for the love story of Layla and Majnun. Chaos Solitons Fractals 2021, 150, 111091. [Google Scholar] [CrossRef]
  18. Loverro, A. Fractional Calculus: History, Definitions and Applications for the Engineer; Rapport Technique; Department of Aerospace and Mechanical Engineering, Univeristy of Notre Dame: Notre Dame, IN, USA, 2004; pp. 1–28. [Google Scholar]
  19. Nabi, K.N.; Abboubakar, H.; Kumar, P. Forecasting of COVID-19 pandemic: From integer derivatives to fractional derivatives. Chaos Solitons Fractals 2020, 141, 110283. [Google Scholar] [CrossRef] [PubMed]
  20. Angstmann, C.N.; Jacobs, B.A.; Henry, B.I.; Xu, Z. Intrinsic Discontinuities in Solutions of Evolution Equations Involving Fractional Caputo–Fabrizio and Atangana–Baleanu Operators. Mathematics 2020, 8, 2023. [Google Scholar] [CrossRef]
  21. Tarasov, V.E. On chain rule for fractional derivatives. Commun. Nonlinear Sci. Numer. Simul. 2016, 30, 1–4. [Google Scholar] [CrossRef]
  22. Tarasov, V.E. No violation of the Leibniz rule. No fractional derivative. Commun. Nonlinear Sci. Numer. Simul. 2013, 18, 2945–2948. [Google Scholar] [CrossRef] [Green Version]
  23. Peter, O.; Ibrahim, M.; Oguntolu, F.; Akinduko, O.; Akinyemi, S. Direct and indirect transmission dynamics of typhoid fever model by differential transform method. ATBU J. Sci. Technol. Educ. (JOSTE) 2018, 6, 167–177. [Google Scholar]
  24. Diethelm, K. An algorithm for the numerical solution of differential equations of fractional order. Electron. Trans. Numer. Anal 1997, 5, 1–6. [Google Scholar]
  25. Diethelm, K.; Ford, N.J. Analysis of fractional differential equations. J. Math. Anal. Appl. 2002, 265, 229–248. [Google Scholar] [CrossRef] [Green Version]
  26. Boccaletti, S.; Ditto, W.; Mindlin, G.; Atangana, A. Modeling and forecasting of epidemic spreading: The case of Covid-19 and beyond. Chaos Solitons Fractals 2020, 135, 109794. [Google Scholar] [CrossRef]
  27. Khan, M.A.; Atangana, A.; Alzahrani, E. The dynamics of COVID-19 with quarantined and isolation. Adv. Differ. Equ. 2020, 2020, 425. [Google Scholar] [CrossRef]
  28. Schmidt, A.; Gaul, L. On the numerical evaluation of fractional derivatives in multi-degree-of-freedom systems. Signal Process. 2006, 86, 2592–2601. [Google Scholar] [CrossRef]
  29. van den Driessche, P.; Watmough, J. Reproduction numbers and the sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2002, 180, 29–48. [Google Scholar] [CrossRef]
  30. Li, H.; Cheng, J.; Li, H.B.; Zhong, S.M. Stability analysis of a fractional-order linear system described by the Caputo–Fabrizio derivative. Mathematics 2019, 7, 200. [Google Scholar] [CrossRef] [Green Version]
  31. Wojtak, W.; Silva, C.J.; Torres, D.F. Uniform asymptotic stability of a fractional tuberculosis model. Math. Model. Nat. Phenom. 2018, 13, 9. [Google Scholar] [CrossRef]
  32. Li, C.; Zeng, F. The finite difference methods for fractional ordinary differential equations. Numer. Funct. Anal. Optim. 2013, 34, 149–179. [Google Scholar] [CrossRef]
  33. Diethelm, K.; Luchko, Y. Numerical solution of linear multi-term initial value problems of fractional order. J. Comput. Anal. Appl. 2004, 6, 243–263. [Google Scholar]
  34. Diethelm, K.; Ford, N.J.; Freed, A.D. A predictor-corrector approach for the numerical solution of fractional differential equations. Nonlinear Dyn. 2002, 29, 3–22. [Google Scholar] [CrossRef]
  35. Pieskä, J.; Laitinen, E.; Lapin, A. Predictor-corrector methods for solving continuous casting problem. In Domain Decomposition Methods in Science and Engineering; Springer: Berlin/Heidelberg, Germany, 2005; pp. 677–684. [Google Scholar]
  36. Butcher, J.C. Numerical methods for ordinary differential equations in the 20th century. J. Comput. Appl. Math. 2000, 125, 1–29. [Google Scholar] [CrossRef] [Green Version]
  37. Li, C.; Tao, C. On the fractional Adams method. Comput. Math. Appl. 2009, 58, 1573–1588. [Google Scholar] [CrossRef] [Green Version]
  38. Abboubakar, H.; Kumar, P.; Rangaig, N.A.; Kumar, S. A Malaria Model with Caputo-Fabrizio and Atangana-Baleanu Derivatives. Int. J. Model. Simul. Sci. Comput. 2021, 12, 2150013. [Google Scholar] [CrossRef]
  39. Garrappa, R. Predictor-Corrector PECE Method for Fractional Differential Equations. MATLAB, Central File Exchange [File ID: 32918]. Available online: https://www.mathworks.com/matlabcentral/fileexchange/32918-predictor-corrector-pece-method-for-fractional-differential-equations (accessed on 23 September 2021).
  40. Lakshmikantham, V.; Leela, S.; Martynyuk, A.A. Stability Analysis of Non Linear Systems; Springer: Berlin/Heidelberg, Germany, 1989. [Google Scholar]
Figure 1. Simulation results showing the fractional dynamics on the infected state variable profiles for different values of the fractional-order parameter η .
Figure 1. Simulation results showing the fractional dynamics on the infected state variable profiles for different values of the fractional-order parameter η .
Fractalfract 05 00149 g001
Figure 2. Simulation results showing the infected state variable profiles without vaccination ( ξ = 0 ).
Figure 2. Simulation results showing the infected state variable profiles without vaccination ( ξ = 0 ).
Fractalfract 05 00149 g002
Figure 3. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 20 % with different values of the fractional order.
Figure 3. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 20 % with different values of the fractional order.
Fractalfract 05 00149 g003
Figure 4. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 50 % with different values of the fractional order.
Figure 4. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 50 % with different values of the fractional order.
Fractalfract 05 00149 g004
Figure 5. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 90 % with different values of the fractional order.
Figure 5. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 90 % with different values of the fractional order.
Fractalfract 05 00149 g005aFractalfract 05 00149 g005b
Figure 6. Simulation results showing the infected state variable profiles when vaccination is combined with environmental sanitation, for fractional order η = 1 .
Figure 6. Simulation results showing the infected state variable profiles when vaccination is combined with environmental sanitation, for fractional order η = 1 .
Fractalfract 05 00149 g006
Figure 7. Simulation results showing the infected state variable profiles when vaccination is combined with environmental sanitation, for fractional order η = 0.90 .
Figure 7. Simulation results showing the infected state variable profiles when vaccination is combined with environmental sanitation, for fractional order η = 0.90 .
Fractalfract 05 00149 g007
Figure 8. Simulation results showing the infected state variable profiles when vaccination is combined with environmental sanitation, for fractional order η = 0.80 .
Figure 8. Simulation results showing the infected state variable profiles when vaccination is combined with environmental sanitation, for fractional order η = 0.80 .
Fractalfract 05 00149 g008
Figure 9. Simulation results showing the infected state variable profiles when vaccination is combined with environmental sanitation, for fractional order η = 0.70 .
Figure 9. Simulation results showing the infected state variable profiles when vaccination is combined with environmental sanitation, for fractional order η = 0.70 .
Fractalfract 05 00149 g009aFractalfract 05 00149 g009b
Figure 10. Parameter estimation and forecasting of cumulative new cases of typhoid fever in Mbandjock, Cameroon, from 1 July 2019 to 31 August 2020, for the new model (35). Red bullets denote real data (see [5,6]). The fitted model is represented with the blue line, and new forecasted cases are represented by the dotted line.
Figure 10. Parameter estimation and forecasting of cumulative new cases of typhoid fever in Mbandjock, Cameroon, from 1 July 2019 to 31 August 2020, for the new model (35). Red bullets denote real data (see [5,6]). The fitted model is represented with the blue line, and new forecasted cases are represented by the dotted line.
Fractalfract 05 00149 g010aFractalfract 05 00149 g010b
Figure 11. Simulation results showing the fractional dynamics on the infected state variable profiles for different values of the fractional-order parameter η .
Figure 11. Simulation results showing the fractional dynamics on the infected state variable profiles for different values of the fractional-order parameter η .
Fractalfract 05 00149 g011
Figure 12. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 0 % with different values of the fractional order.
Figure 12. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 0 % with different values of the fractional order.
Fractalfract 05 00149 g012
Figure 13. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 20 % with different values of the fractional order.
Figure 13. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 20 % with different values of the fractional order.
Fractalfract 05 00149 g013
Figure 14. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 50 % with different values of the fractional order.
Figure 14. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 50 % with different values of the fractional order.
Fractalfract 05 00149 g014aFractalfract 05 00149 g014b
Figure 15. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 90 % with different values of the fractional order.
Figure 15. Simulation results showing the infected state variable profiles when the vaccination coverage ξ = 90 % with different values of the fractional order.
Fractalfract 05 00149 g015
Table 1. Parameter values of (2) taken from [5].
Table 1. Parameter values of (2) taken from [5].
ParameterValueParameterValueParameterValue
Λ h 1 μ b 0.149990 ν 3.2618 × 10 6
γ 1 0.2145 ϵ 0.9495 δ 0.1499
γ 2 0.1498 ξ 0.3221 σ 0.49992
α 0.0834 θ 0.0833 R c 2.4750
Table 2. Estimated parameter values of the new typhoid model (35).
Table 2. Estimated parameter values of the new typhoid model (35).
ParameterValuesSourceParameterValuesSource
Λ h 3Fitted μ b 0.0015Fitted
γ 1 0.1512Fitted ϵ 0.9497Fitted
γ 2 0.3039Fitted ξ 0.1538Fitted
δ 0.1382Fitted β 0.60Fitted
ν 0.00050FittedK995.7957Fitted
σ 0.4992Fitted R c 1.4348Estimated
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abboubakar, H.; Kom Regonne, R.; Sooppy Nisar, K. Fractional Dynamics of Typhoid Fever Transmission Models with Mass Vaccination Perspectives. Fractal Fract. 2021, 5, 149. https://doi.org/10.3390/fractalfract5040149

AMA Style

Abboubakar H, Kom Regonne R, Sooppy Nisar K. Fractional Dynamics of Typhoid Fever Transmission Models with Mass Vaccination Perspectives. Fractal and Fractional. 2021; 5(4):149. https://doi.org/10.3390/fractalfract5040149

Chicago/Turabian Style

Abboubakar, Hamadjam, Raissa Kom Regonne, and Kottakkaran Sooppy Nisar. 2021. "Fractional Dynamics of Typhoid Fever Transmission Models with Mass Vaccination Perspectives" Fractal and Fractional 5, no. 4: 149. https://doi.org/10.3390/fractalfract5040149

APA Style

Abboubakar, H., Kom Regonne, R., & Sooppy Nisar, K. (2021). Fractional Dynamics of Typhoid Fever Transmission Models with Mass Vaccination Perspectives. Fractal and Fractional, 5(4), 149. https://doi.org/10.3390/fractalfract5040149

Article Metrics

Back to TopTop