# Application of a New Generalized Fractional Derivative and Rank of Control Measures on Cholera Transmission Dynamics

Authors: Kumama Regassa Cheneke; Koya Purnachandra Rao; Geremew Kenassa Edessa
Journal: International Journal of Mathematics and Mathematical Sciences
 (2021)

---

## Abstract

In this study, the mathematical model of the cholera epidemic is formulated and analyzed to show the impact ofVibrio cholerae in reserved freshwater. Moreover, the results obtained from applying the new fractional derivative method show that, as the order of the fractional derivative increases, cholera-preventing behaviors also increase. Also, the finding of our study shows that the dynamics of Vibrio cholerae can be controlled if continuous treatment is applied in reserved freshwater used for drinking purposes so that the intrinsic growth rate of Vibrio cholerae in water is less than the natural death of Vibrio cholerae. We have applied the stability theory of differential equations and proved that the disease-free equilibrium is asymptotically stable if R0<1, and the intrinsic growth rate of the Vibrio cholerae bacterium population is less than its natural death rate. The center manifold theory is applied to show the existence of forward bifurcation at the point R0=1 and the local stability of endemic equilibrium if R0>1. Furthermore, the performed numerical simulation results show that, as the rank of control measures applied increases from no control, weak control, and strong control measures, the recovered individuals are 55.02, 67.47, and 674.7, respectively. Numerical simulations are plotted using MATLAB software package.


---

## Body

## 1. Introduction

Cholera is a waterborne acute gastrointestinal infection characterized by diarrhea and vomiting, which kills within an hour if not treated [1]. Cholera is caused by drinking water or eating food contaminated by a bacterium called Vibrio cholerae (V. cholerae) [2, 3]. Most of the diseases are controlled by public health organizations including cholera which is characterized with severe vomiting and diarrhea [4]. The incubation period of cholera is less than five days [5]. However, cholera outbreak occurs twice in endemic areas, and it becomes the major issue in developing countries. Cholera infection can be transmitted from human to human or human to environment through effective contacts with bacterium V. cholerae [6].Different mathematical models have been carried out to describe the transmission dynamics of cholera infection [7, 8]. Codeco developed the first basic mathematical model for cholera infection. Furthermore, a new fractional derivative was developed by Hattaf [9] and applied to analyze the memory effect on the HIV-infected population. However, the new fractional derivative has not been applied to study the cases in the cholera infection.In this study, we modified the model developed in [5], we focus on the logistic growth model of Vibrio cholerae concentration in reserved freshwater used for drinking purposes, and the contribution of infected individuals to the environment is assumed to be properly managed. That is, all waste disposals from cholera-infected patients are properly managed, and the contributions of humans to the environment are restricted. Hence, we developed an optimal control mathematical model that takes into consideration the loss of immunity after cholera infection recovery. Our model considers the reimmunization system for those individuals that lose immunity. Furthermore, we have applied a new fractional derivative and analyzed the relationship between the order of the derivative and the extinction status of cholera infection.


## 2. Mathematical Model Developments

In this study, we have developed the deterministic mathematical model of cholera with optimal control strategies. We have divided the total population as (i) susceptibleS: these are individuals free of the infection but can be infected with the infection if exposed to the infectious people or environment contaminated with the cholera-causing bacterium; (ii) infected I: these are individuals that get infected with cholera and can transmit infection to others with possible contacts that cause infection; (iii) recovered R: these are individuals recovered from the cholera infection with treatment. The recovered individuals who lose immunity become susceptible; (iv) concentration of the Vibrio cholerae bacterium Bt: this covers an environment contaminated with Vibrio cholerae that cause the cholera epidemic.Furthermore, the following assumptions are incorporated into the formulation of the model:(i)
The total population is considered to be nonconstant(ii)
All humans die due to cholera at the disease-induced death rate ofφ(iii)
All waste products containVibrio cholerae which are assumed to be managed, and there is no contribution of the infected person to reserved freshwater(iv)
All human populations are assumed to die at the natural death rateμ(v)
All new infected human individuals are recruited to the susceptible at the recruitment rateλ(vi)
Treatment is performed in reserved freshwater(vii)
Cholera-infected individuals recover at a recovery rate ofη(viii)
Cholera infection-recovered individuals lose immunity at the rateγ(ix)
Cholera infection is transmitted from human to human at the rateβ(x)
V. cholerae is ingested from the environment at the ingestion rate βc(xi)
The saturated incidence rate is used in the transmission of cholera infection from the environment to humans(xii)
The concentration ofV. cholerae in water that causes 50% chance of cholera infection is denoted by κ(xiii)
The carrying capacity ofV. cholerae is K(xiv)
V. cholerae die at the rate of ν(xv)
Cholera infection transmission preventive control effort isu1(xvi)
Immunity loss control effort isu2(xvii)
V. cholerae clearing control effort is u3In general, the stated assumptions can be described by the flow diagram.The behavior of the population described in Figure1 can be described by the mathematical model as(1)dSdt=λ−1−u1βcSBk+B+βSI+1−u2γR–μS,dIdt=1−u1βcSBk+B+βSI−η+φ+μI,dRdt=ηI−1−u2γ+μR,dBdt=1−u3rB1−BK−νB,with nonnegative initial conditions.Figure 1
The flow diagram of the SIR-B model of the cholera epidemic with control measures.


## 3. Mathematical Analysis of Only Cholera Model without Optimal Control

Formulated model (1) of cholera without control measures is reduced to the following form:(2)dSdt=λ−βcSBk+B−βSI+γR–μS,dIdt=βcSBk+B+βSI−η+φ+μI,dRdt=ηI−γ+μR,dBdt=rB1−BK−νB,with nonnegative initial conditions.


### 3.1. Invariant Regions

Theorem 1.
Solutions of model (1) are invariant in the region Ω⊂ℜ+4 such that(3)Ω=ΩH×ΩB⊂ℜ+3×ℜ+=S,I,R,B:S,I,R∈ℜ+3,B∈ℜ+,Nh=S+I+R≤λμ,B≤K.Proof.
As in [10], to prove this theorem, we add the first three equations of model (1) so that(4)dNhdt=λ−μNh−φI≤λ−μNh.
Applying a definite integral over time interval [0,t], the preceding inequality gives(5)Nh≤λμ−λμ−N0e−μt.
Here, as timet gets larger and larger, Nh approaches to least upper bound λ/μ. Also, Nh≤λ/μ, for all time t. Similarly, from the last equation of model (2), we obtain B≤K for all time t. Therefore, the invariant region Ω⊂ℜ+4 is given by(6)Ω=S,I,R,B:Nh=S+I+R≤λμ,B≤K.


### 3.2. Positivity Property of Solutions

Theorem 2.
For given positive initial conditions, all possible solutions of model (1) are positive in the invariant Ω⊂ℜ+4 such that(7)Ω=S,I,R,B:S0>0,I0≥0,R0≥0,B0≥0.Proof.
To prove the positivity property, we apply the methods employed in [11]. Consider the first equation of model (1):(8)dSdt=λ−βcSBk+B−βSI+γR–μS.
Rearranging the preceding equation, we get(9)dSdt+βcBk+B+βI+μS=λ+γR.
Also, we have(10)ddtStexpμt+∫0tβcBζk+Bζ+βIζdζ=λ+γRexpμt+∫0tβcBζk+Bζ+βIζdζ.
Applying integration on the preceding equation over time interval [0,t], we get(11)Stexpμt+∫0tβcBζk+Bζ+βIζdζ−S0=∫0tλ+γRζexpμt+∫0tβcBζk+Bζ+βIζdζdt.
Hence,(12)St=S0exp−μt+∫0tβcBζk+Bζ+βIζdζ+exp−μt+∫0tβcBζk+Bζ+βIζdζ∫0tλ+γRζexpμt+∫0tβcBζk+Bζ+βIζdζdt>0.
Since the expressions(13)λ+γRζ>0,exp−μt+∫0tβcBζk+Bζ+βIζdζ+exp−μt+∫0tβcBζk+Bζ+βIζdζ,∫0tλ+γRζexpμt+∫0tβcBζk+Bζ+βIζdζdtare positive, the state variable St>0 for all time t. Similarly, all other state variables are positive. Therefore, the set that contains all solution variables of model (1) is a positively invariant set.


### 3.3. Disease-Free Equilibrium (DFE)

The disease-free equilibriumE0 of model (1) is given by(14)E0=λμ,0,0,0.


### 3.4. Endemic Equilibrium (EE)

The computed endemic equilibriumE1=S∗,I∗,R∗,B∗ of model (2) is given by(15)E1=S∗,I∗,ηγ+μI∗,B∗,where I∗=λ−μS∗/η+φ+μ,B∗=K1−ν/r=B, and S∗=S is the solution of the equation(16)μ2R0λS2−βcBK+B+μR0+1S+λ=0.


### 3.5. Reproduction Number

In this section, the reproduction numberR0 is computed by employing the next-generation matrix FV−1 stated in research works [1, 5, 12]. Thus, from model (2), the reproduction number can be computed as(17)R0=ρFV−1,where F=βλ/μβcλ/μk0randV=η+φ+μ00ν.Hence,(18)FV−1=βλbμβckλμν0rν.The eigenvalues of the preceding next-generation matrix are(19)λ1=βλ/bμ,λ2=r/ν.Therefore,(20)R0=maxβλbμ,rν=βλbμ,where b=η+φ+μ.


### 3.6. Local Stability of Disease-Free Equilibrium

Theorem 3.
The disease-free equilibriumE0 of model (2) is locally asymptotically stable for the reproduction number less than one and intrinsic growth rate r of bacteria less than the natural death rate. On the contrary, E0 is unstable for the reproduction number greater than one.Proof.
To approve the local stability ofE0, we compute the eigenvalues of the Jacobian matrix obtained from model (2) and evaluated at E0. Hence, it follows that(21)JE0=−μ−βλμγ−βcλμk0η+φ+μR0−10βcλμk00−γ+μ0000r−ν.
The computed eigenvalues of the preceding matrix are given by(22)λ1=−μ,λ2=−μ,λ3=η+φ+μR0−1,λ4=r−ν.
From the stability principle on the eigenvalues computed from the Jacobian matrix evaluated atE0, we conclude that the disease-free equilibrium E0 is locally asymptotically stable if R0<1 and r<ν and unstable if either R0>1 or r>ν.


### 3.7. Global Stability of Disease-Free Equilibrium

Model (2) can be rewritten in the following form:(23)dXdt=HX,Y,dYdt=GX,Y,GX,0=0.The disease-free equilibriumE01 of the preceding system is(24)E01=X0,0,where X0 is the disease-free equilibrium of the disease-free system.According to Biswas et al. [13], to guarantee global asymptotic stability, we verify the following conditions H1 and H2 to be satisfied:H1: fordX/dt=HX,0,X0 is the globally asymptotically stable equilibriumH2:GX,Y=PY−G^X,Y,G^X,Y≥0 for X,Y∈ΩHere,P=DYGX,0 satisfy the condition of the Metzler matrix, and Ω is a region of feasible solutions. Now, we state the following theorem.Theorem 4.
The disease-free equilibrium of the cholera dynamic model is globally asymptotically stable ifR0<1 and r<ν if conditions H1 and H2 are satisfied and unstable whenever R0>1.Proof.
From model (2), we have(25)HX,0=λ−μS=HS,0.
SolvingHX,0=0, we obtain S=λ/μ. Hence, X0=λ/μ,0.
Here,X0 is the globally stable equilibrium of equation(26)dXdt=HX,0.
From the infected compartments of model (2), we obtain(27)GX,Y=−η+φ+μ00r−νIB−−βcSBk+B−βSIrB2K.
At the disease-free equilibrium, the preceding equation reduces to the following form:(28)GX,Y=GX,Y=−η+φ+μ00r−νIB−00=PY−G^X,Y,where P=−η+φ+μ00r−ν and G^X,Y=00≥0,∀X,Y.
Now, it follows that the formulated model satisfied the hypothesis conditions required asGX,Y=PY−G^X,Y, where G^X,Y≥0,∀X,Y. Therefore, E0 is globally asymptotically stable if R0<1.


### 3.8. Bifurcation Analysis

Theorem 5.
The endemic equilibriumE1 of model (2) exhibits backward bifurcation at R0=1.Proof.
Letβ∗ be a bifurcation parameter corresponding to parameter β and obtained by setting R0=1 and solving for β∗. Therefore, the bifurcation parameter β∗ is given by(29)β∗=μη+φ+μλ.
To determine types of bifurcation, we construct Jacobian matrixJE0,β∗ from model (2) and compute the eigenvalues of JE0,β∗. Thus, we have(30)JE0,β∗=−μ−η+φ+μγ−βcλμk000βcλμk00−γ+μ0000r−ν.
The computed eigenvalues of the preceding matrix are given by(31)λ1=−μ,λ2=0,λ3=−μ,λ4=r−ν.
Here, all eigenvalues are negative exceptλ2 which is simply 0 with condition r<ν. Therefore, according to Akanni et al. [14], model (2) exhibits bifurcation at bifurcation parameter β∗. To determine the type of bifurcation, we compute left eigenvector v=v1,v2,v3,v4 and right eigenvector w=w1,w2,w3,w4 that satisfy the conditions v⋅JE0,β∗=0,JE0,β∗⋅w=0,andv⋅w=1. Hence, the computed left eigenvalues from equation v⋅JE0,β∗=0 are given by(32)v1=0,v2=v2,v3=0,v4=−μkβcλr−ν.
And the computed right eigenvalues from equationJE0,β∗⋅w=0 are given by(33)w1=−η+φ+μw2μ,w2=w2,w3=0,w4=0.
Also, from the conditionv⋅w=1, we have(34)v2=1w2.
Let us choosev2=w2=1; we compute the following second-order partial derivatives of functions evaluated at the disease-free equilibrium and bifurcation parameter.
Now, setx=x1,x2,x3,x4=S,I,R,B and(35)fx=f1f2f3f4x=λ−βcx1x4k+x4−βx1x2+γx3–μx1βcx1x4k+x4+βx1x2−η+φ+μx2ηx2−γ+μx3rx41−x4K−νx4,∂2f1∂x12=0,∂2f1∂x22=0,∂2f1∂x32=0,∂2f1∂x42=2βckx1k+x43,∂2f2∂x12=0,∂2f2∂x22=0,∂2f2∂x32=0,∂2f2∂x42=−2βckx1k+x43,∂2f3∂x12=0,∂2f3∂x22=0,∂2f3∂x32=0,∂2f3∂x42=0,∂2f4∂x12=0,∂2f4∂x22=0,∂2f4∂x32=0,∂2f4∂x42=−2rk,∂2f1∂x2∂x1=−β,∂2f1∂x3∂x1=0,∂2f1∂x4∂x1=−βckk+x42,∂2f2∂x2∂x1=β,∂2f2∂x3∂x1=0,∂2f2∂x4∂x1=βckk+x42,∂2f3∂x2∂x1=0,∂2f3∂x3∂x1=0,∂2f3∂x4∂x1=0,∂2f4∂x2∂x1=0,∂2f4∂x3∂x1=0,∂2f4∂x4∂x1=0,∂2f2∂x1∂βc∗=kx2+x2x4+x4k+x4,∂2f2∂x2∂β∗=x1,∂2f2∂x3∂βc∗=0,∂2f2∂x4∂βc∗=kx1k+x42.
Next, the bifurcation coefficientsa and b can be computed as follows:(36)a=∑k,i,j=14vkwiwj∂2fk∂xi∂xjE0,β∗=−2βcλμk2−2η+φ+μβcμ<0,b=∑k,i=14vkwi∂2fk∂xi∂β∗E0,β∗=v2w2x1=λμ>0.
Hence, according to Akanni et al. [14], the model shows supercritical (forward) bifurcation at R0=1.


### 3.9. Stability of Endemic EquilibriumE1

Theorem 6.
Endemic equilibrium is locally asymptotically stable ifR0>1.Proof.
It follows from the center manifold theory, a model that shows forward bifurcation has endemic equilibrium that exhibits locally asymptotic stability ifR0>1 [14].Theorem 7.
Endemic equilibrium is globally asymptotically stable ifR0>1.Proof.
It follows from the center manifold theory, a model that shows forward bifurcation has endemic equilibrium that exhibits globally asymptotic stability ifR0>1 [14].


## 3.1. Invariant Regions

Theorem 1.
Solutions of model (1) are invariant in the region Ω⊂ℜ+4 such that(3)Ω=ΩH×ΩB⊂ℜ+3×ℜ+=S,I,R,B:S,I,R∈ℜ+3,B∈ℜ+,Nh=S+I+R≤λμ,B≤K.Proof.
As in [10], to prove this theorem, we add the first three equations of model (1) so that(4)dNhdt=λ−μNh−φI≤λ−μNh.
Applying a definite integral over time interval [0,t], the preceding inequality gives(5)Nh≤λμ−λμ−N0e−μt.
Here, as timet gets larger and larger, Nh approaches to least upper bound λ/μ. Also, Nh≤λ/μ, for all time t. Similarly, from the last equation of model (2), we obtain B≤K for all time t. Therefore, the invariant region Ω⊂ℜ+4 is given by(6)Ω=S,I,R,B:Nh=S+I+R≤λμ,B≤K.


## 3.2. Positivity Property of Solutions

Theorem 2.
For given positive initial conditions, all possible solutions of model (1) are positive in the invariant Ω⊂ℜ+4 such that(7)Ω=S,I,R,B:S0>0,I0≥0,R0≥0,B0≥0.Proof.
To prove the positivity property, we apply the methods employed in [11]. Consider the first equation of model (1):(8)dSdt=λ−βcSBk+B−βSI+γR–μS.
Rearranging the preceding equation, we get(9)dSdt+βcBk+B+βI+μS=λ+γR.
Also, we have(10)ddtStexpμt+∫0tβcBζk+Bζ+βIζdζ=λ+γRexpμt+∫0tβcBζk+Bζ+βIζdζ.
Applying integration on the preceding equation over time interval [0,t], we get(11)Stexpμt+∫0tβcBζk+Bζ+βIζdζ−S0=∫0tλ+γRζexpμt+∫0tβcBζk+Bζ+βIζdζdt.
Hence,(12)St=S0exp−μt+∫0tβcBζk+Bζ+βIζdζ+exp−μt+∫0tβcBζk+Bζ+βIζdζ∫0tλ+γRζexpμt+∫0tβcBζk+Bζ+βIζdζdt>0.
Since the expressions(13)λ+γRζ>0,exp−μt+∫0tβcBζk+Bζ+βIζdζ+exp−μt+∫0tβcBζk+Bζ+βIζdζ,∫0tλ+γRζexpμt+∫0tβcBζk+Bζ+βIζdζdtare positive, the state variable St>0 for all time t. Similarly, all other state variables are positive. Therefore, the set that contains all solution variables of model (1) is a positively invariant set.


## 3.3. Disease-Free Equilibrium (DFE)

The disease-free equilibriumE0 of model (1) is given by(14)E0=λμ,0,0,0.


## 3.4. Endemic Equilibrium (EE)

The computed endemic equilibriumE1=S∗,I∗,R∗,B∗ of model (2) is given by(15)E1=S∗,I∗,ηγ+μI∗,B∗,where I∗=λ−μS∗/η+φ+μ,B∗=K1−ν/r=B, and S∗=S is the solution of the equation(16)μ2R0λS2−βcBK+B+μR0+1S+λ=0.


## 3.5. Reproduction Number

In this section, the reproduction numberR0 is computed by employing the next-generation matrix FV−1 stated in research works [1, 5, 12]. Thus, from model (2), the reproduction number can be computed as(17)R0=ρFV−1,where F=βλ/μβcλ/μk0randV=η+φ+μ00ν.Hence,(18)FV−1=βλbμβckλμν0rν.The eigenvalues of the preceding next-generation matrix are(19)λ1=βλ/bμ,λ2=r/ν.Therefore,(20)R0=maxβλbμ,rν=βλbμ,where b=η+φ+μ.


## 3.6. Local Stability of Disease-Free Equilibrium

Theorem 3.
The disease-free equilibriumE0 of model (2) is locally asymptotically stable for the reproduction number less than one and intrinsic growth rate r of bacteria less than the natural death rate. On the contrary, E0 is unstable for the reproduction number greater than one.Proof.
To approve the local stability ofE0, we compute the eigenvalues of the Jacobian matrix obtained from model (2) and evaluated at E0. Hence, it follows that(21)JE0=−μ−βλμγ−βcλμk0η+φ+μR0−10βcλμk00−γ+μ0000r−ν.
The computed eigenvalues of the preceding matrix are given by(22)λ1=−μ,λ2=−μ,λ3=η+φ+μR0−1,λ4=r−ν.
From the stability principle on the eigenvalues computed from the Jacobian matrix evaluated atE0, we conclude that the disease-free equilibrium E0 is locally asymptotically stable if R0<1 and r<ν and unstable if either R0>1 or r>ν.


## 3.7. Global Stability of Disease-Free Equilibrium

Model (2) can be rewritten in the following form:(23)dXdt=HX,Y,dYdt=GX,Y,GX,0=0.The disease-free equilibriumE01 of the preceding system is(24)E01=X0,0,where X0 is the disease-free equilibrium of the disease-free system.According to Biswas et al. [13], to guarantee global asymptotic stability, we verify the following conditions H1 and H2 to be satisfied:H1: fordX/dt=HX,0,X0 is the globally asymptotically stable equilibriumH2:GX,Y=PY−G^X,Y,G^X,Y≥0 for X,Y∈ΩHere,P=DYGX,0 satisfy the condition of the Metzler matrix, and Ω is a region of feasible solutions. Now, we state the following theorem.Theorem 4.
The disease-free equilibrium of the cholera dynamic model is globally asymptotically stable ifR0<1 and r<ν if conditions H1 and H2 are satisfied and unstable whenever R0>1.Proof.
From model (2), we have(25)HX,0=λ−μS=HS,0.
SolvingHX,0=0, we obtain S=λ/μ. Hence, X0=λ/μ,0.
Here,X0 is the globally stable equilibrium of equation(26)dXdt=HX,0.
From the infected compartments of model (2), we obtain(27)GX,Y=−η+φ+μ00r−νIB−−βcSBk+B−βSIrB2K.
At the disease-free equilibrium, the preceding equation reduces to the following form:(28)GX,Y=GX,Y=−η+φ+μ00r−νIB−00=PY−G^X,Y,where P=−η+φ+μ00r−ν and G^X,Y=00≥0,∀X,Y.
Now, it follows that the formulated model satisfied the hypothesis conditions required asGX,Y=PY−G^X,Y, where G^X,Y≥0,∀X,Y. Therefore, E0 is globally asymptotically stable if R0<1.


## 3.8. Bifurcation Analysis

Theorem 5.
The endemic equilibriumE1 of model (2) exhibits backward bifurcation at R0=1.Proof.
Letβ∗ be a bifurcation parameter corresponding to parameter β and obtained by setting R0=1 and solving for β∗. Therefore, the bifurcation parameter β∗ is given by(29)β∗=μη+φ+μλ.
To determine types of bifurcation, we construct Jacobian matrixJE0,β∗ from model (2) and compute the eigenvalues of JE0,β∗. Thus, we have(30)JE0,β∗=−μ−η+φ+μγ−βcλμk000βcλμk00−γ+μ0000r−ν.
The computed eigenvalues of the preceding matrix are given by(31)λ1=−μ,λ2=0,λ3=−μ,λ4=r−ν.
Here, all eigenvalues are negative exceptλ2 which is simply 0 with condition r<ν. Therefore, according to Akanni et al. [14], model (2) exhibits bifurcation at bifurcation parameter β∗. To determine the type of bifurcation, we compute left eigenvector v=v1,v2,v3,v4 and right eigenvector w=w1,w2,w3,w4 that satisfy the conditions v⋅JE0,β∗=0,JE0,β∗⋅w=0,andv⋅w=1. Hence, the computed left eigenvalues from equation v⋅JE0,β∗=0 are given by(32)v1=0,v2=v2,v3=0,v4=−μkβcλr−ν.
And the computed right eigenvalues from equationJE0,β∗⋅w=0 are given by(33)w1=−η+φ+μw2μ,w2=w2,w3=0,w4=0.
Also, from the conditionv⋅w=1, we have(34)v2=1w2.
Let us choosev2=w2=1; we compute the following second-order partial derivatives of functions evaluated at the disease-free equilibrium and bifurcation parameter.
Now, setx=x1,x2,x3,x4=S,I,R,B and(35)fx=f1f2f3f4x=λ−βcx1x4k+x4−βx1x2+γx3–μx1βcx1x4k+x4+βx1x2−η+φ+μx2ηx2−γ+μx3rx41−x4K−νx4,∂2f1∂x12=0,∂2f1∂x22=0,∂2f1∂x32=0,∂2f1∂x42=2βckx1k+x43,∂2f2∂x12=0,∂2f2∂x22=0,∂2f2∂x32=0,∂2f2∂x42=−2βckx1k+x43,∂2f3∂x12=0,∂2f3∂x22=0,∂2f3∂x32=0,∂2f3∂x42=0,∂2f4∂x12=0,∂2f4∂x22=0,∂2f4∂x32=0,∂2f4∂x42=−2rk,∂2f1∂x2∂x1=−β,∂2f1∂x3∂x1=0,∂2f1∂x4∂x1=−βckk+x42,∂2f2∂x2∂x1=β,∂2f2∂x3∂x1=0,∂2f2∂x4∂x1=βckk+x42,∂2f3∂x2∂x1=0,∂2f3∂x3∂x1=0,∂2f3∂x4∂x1=0,∂2f4∂x2∂x1=0,∂2f4∂x3∂x1=0,∂2f4∂x4∂x1=0,∂2f2∂x1∂βc∗=kx2+x2x4+x4k+x4,∂2f2∂x2∂β∗=x1,∂2f2∂x3∂βc∗=0,∂2f2∂x4∂βc∗=kx1k+x42.
Next, the bifurcation coefficientsa and b can be computed as follows:(36)a=∑k,i,j=14vkwiwj∂2fk∂xi∂xjE0,β∗=−2βcλμk2−2η+φ+μβcμ<0,b=∑k,i=14vkwi∂2fk∂xi∂β∗E0,β∗=v2w2x1=λμ>0.
Hence, according to Akanni et al. [14], the model shows supercritical (forward) bifurcation at R0=1.


## 3.9. Stability of Endemic EquilibriumE1

Theorem 6.
Endemic equilibrium is locally asymptotically stable ifR0>1.Proof.
It follows from the center manifold theory, a model that shows forward bifurcation has endemic equilibrium that exhibits locally asymptotic stability ifR0>1 [14].Theorem 7.
Endemic equilibrium is globally asymptotically stable ifR0>1.Proof.
It follows from the center manifold theory, a model that shows forward bifurcation has endemic equilibrium that exhibits globally asymptotic stability ifR0>1 [14].


## 4. Application of a New Generalized Fractional Derivative with the Nonsingular Kernel

In this section, using the works done in [9], we study the dynamics of the human population with a new generalized definition of the fractional derivative with a nonnegative singular kernel in the presence of cholera infection. The rate of change of human population, from model (2), can be described as(37)dSdt=λ−βcSBk+B−βSI+γR–μS,dIdt=βcSBk+B+βSI−η+φ+μI,dRdt=ηI−γ+μR.LetNht be the size of human populations at time t and defined as(38)Nht=St+It+Rt.Differentiating both sides of the preceding equation with respect to timet and using (37), we get(39)dNhdt=λ–μNh−φI≤λ–μNh.Now, considering the upper bound, the preceding inequality reduces to the form(40)dNhdt=λ–μNh.The solution of the preceding equation over time interval [0,t] is given by(41)Nh=λμ−λμ−Nh0e−μt.Definition 1.
(new generalized fractional derivative with nonsingular kernel; see [9]).
Letα∈0,1,β,γ>0, and f∈H1a,b. In the sense of Caputo, the α-order new generalized fractional derivative of function ft with respect to the weight function wt is defined as(42)DCa,t,wα,β,γft=Nα1−α1wt∫atEβ−μαt−xγddxwfxdx,where w∈C1a,b,w,w′>0 on a,b, Nα is a normalized function obeying N0=N1=1,μα=α/1−α, and Eβt=∑k=0+∞tk/Γβk+1 is the Mittag-Leffler function of parameter β.
For simplicity, we symbolizeDCa,t,wα,β,γft by Da,t,wα,β,γft. Assuming that when a cholera infection spreads within a community, the individuals get knowledge about the infection, we replace the classical derivative by Da,t,wα,β,γ. Then, (41) becomes(43)Da,t,wα,β,γNht=λ–μNht.
Applying the Laplace transform to (43), we get(44)ℒwtDa,t,wα,β,γNht=λℒwt−μℒwtNht.
Using Theorem2 [9], we have(45)ℒwtNht=Nαw0T0sβ−1Nα+μ1−αsβ+αμ+λ1−αsβ+λαNα+μ1−αsβ+αμℒwt.
Letaα=Nα+μ1−α. Then,(46)ℒwtNht=Nαw0Nh0ααsβ−1sβ+αμ/aα+λ1−αaαsβ−1sβ+αμ/aαsℒwt+λαaα1sβ+αμ/aαℒwt.
Therefore,(47)ℒwtNht=Nαw0Nh0ααℒEβ−αμaαtβ+λ1−αaαℒEβ−αμaαtβℒw′t+w0−λμℒddtEβ−αμaαtβℒwt.
Applying inverse Laplace, the preceding equation reduces to(48)wtNht=Nαw0Nh0ααEβ−αμaαtβ+λ1−αaαEβ−αμaαtβw′t−λμddtEβ−αμaαtβwt+λ1−αw0aαEβ−αμaαtβ.
Lettingwt=1 and applying the integration by parts, the preceding equation reduces to(49)Nht=λμ−Nαaαλμ−Nh0Eβ−αμaαtβ.
In the next section, using numerical simulations, we will analyze the impact of the order of the new fractional derivative on the behavior of solutions of the cholera infection model.


## 5. Numerical Simulations

The parameters presented in Table1 are used in the simulation of solutions of the SIR-B model and are either assumed logically or taken from the literature.Table 1
Description of the parameter value and source.
ParameterParameter valueSourceλ15/day[15]βc1.710−8[5]β1.4810−8[5]γ0.025/day[15]η0.004/day[15]μ0.0000548/day[15]φ0.15/day[15]ν0.033/day[15]κ106cells/ml[15]r0.0004EstimatedK108cells/mlEstimated


## 6. Results and Discussion

In this study, the three-dimensional dynamical system is formulated and analyzed, considering the interaction of human to human and human to environment in the transmission of theVibrio cholerae bacterium. The analysis of the SIR-B model of the cholera epidemic shows that we can control cholera more effectively if we apply strong effort on the preventive, sanitation, hygiene, vaccination, and follow-up program. Figure 2 shows that, without control measures, the recovered individuals and others become susceptible. Figure 3 shows that weak control measures do not bring a significant change on the epidemic even though the number of pathogens is reduced. Figure 4 shows that the high contribution of control measures helps in the reduction of the Vibrio cholerae bacterium in freshwater but increases recovered individuals. Also, from Figures 2–4, one can observe that, as control measures increase applied as no control, weak control, and strong control, the recovered individuals increase as 55.02, 67.47, and 674.7, respectively. The SIR-B model of cholera exhibits globally asymptotic stability at disease-free equilibrium if R0>1 and unstable if r>ν. Figure 5 describes that, as the order of fractional derivatives increases, the total population size approaches to disease-free equilibrium of the cholera model. This shows that, as the human population gets free of cholera, the behaviors of the population change along with the order of the new fractional derivative. Furthermore, the memory effect results we obtained in this study are supported with the study done in [9].Figure 2
Simulation of the SIR-B model withu1=0,u2=0,andu3=0.Figure 3
Simulation of the SIR-B model with weak control measures (u1=u2=u3=0.1).Figure 4
Simulation of the SIR-B model with strong control measures (u1=0.95,u2=0.95,andu3=0.95).Figure 5
Simulation of total human population size described in (49) for different values of the fractional order (α=0,α=0.2,α=0.4,α=0.6,andα=0.8) and Nh0=6000.


## 7. Conclusion

The current study pointed out that, as we apply strong control measures at all levels, the concentration ofVibrio cholerae in freshwater will be reduced. Furthermore, the study pointed out that, to significantly reduce the concentration of Vibrio cholerae in reserved freshwater, the intrinsic growth rate of the Vibrio cholerae bacterium must be less than the natural death rate of the Vibrio cholerae bacterium in freshwater. The stability of the disease-free equilibrium SIR-B model is accompanied with the intrinsic growth rate of Vibrio cholerae that needs to be less than its natural death rate and reproduction number less than one. The new fractional derivative results in that the memory effect in the human population has a direct relationship. As the order of the new fractional derivative increases, human memories of cholera infection controlling also increase.



---

*Source: 2104051-2021-11-02.xml*