A Mathematical Model for the Control of Cholera Epidemic without Natural Recovery

Show more

Introduction

Cholera is a contagious infectious disease that is characterized by extreme vomiting, profuse watery diarrhoea and leg pain. It has been found that transmission transpires mostly via absorption of contaminated drinking water or food. Worldwide, almost every year, there are estimated 3 - 5 million cholera cases and 100,000 - 130,000 deaths due to cholera a year as of 2010 [1]. It has a very short incubation period which starts from a few hours to five days. The health of an infected person disintegrates rapidly and death may occur if treatment is not promptly given. Cholera was first discovered in the Indian subcontinent in 1817. The disease reaches all the way through Asian continent in the 1960s, getting into Africa in 1970 and Latin America in 1991 [2] [3]. In many parts of Africa and Asia, the disease is still endemic. Cholera is a disastrous water-borne infectious disease that is caused by the bacterium vibrio cholera. It is a very serious problem in many developing countries due to inadequate access to safe drinking water supply, improper treatment of reservoirs and improper sanitation. In 2012, WHO reported 245, 393 cholera cases and 3034 death cases across 48 countries in which 67% cases occurred in African countries [4]. In 2005, Nigeria had 4477 cases and 174 deaths. There were reported cases of cholera in 2008 in Nigeria in which there were 429 deaths out of 6330 cases. Furthermore, 2304 cases were reported in Niger State in which 114 were death cases [5]. Also in 2009, Nigeria reported 13,691 cases and 431 deaths [6] (WHO). Tian and Wang [7] evidenced that recent years have seen a strong trend of cholera outbreak in developing countries, such as in India (2007), Iraq (2008), Congo (2008), Zimbabwe (2008-2009), Haiti (2010), Kenya (2010) and Nigeria (2010).

In Nigeria, outbreaks of the disease have been taking place with ever-increasing occurrence ever since the earliest outbreak in recent times in 1970, [8] [9]. In summary, the United Nation (UN) unit reports: “despite Nigeria’s oil wealth, more than 70% of the country’s 126 million people live below the poverty line and cholera outbreaks are common in poor urban areas which lack proper sanitation and clean drinking water” (UN Office for the Coordination of Humanitarian Affairs Integrated Regional Information Networks (IFIN) 2005). In the last few decades, several researchers [10] [11] [12] [13] [14] have designed mathematical models to explore the transmission dynamics and control of the disease. However, in the literature, not many studies on cholera in Nigeria using bio-mathematical models are prevalent, hence, the need for this research study.

1.1. Model Formulation

A mathematical model for cholera transmission is developed by making a slight modification to the previous model by Sulayman Fatima *et al*. [15]. Our model is devoid of natural recovery, this is because from the literature we find that cholera is a fatal disease, and a large percentage of those who recover, do so because of treatment. Our model on the other hand incorporates treatment, water hygiene and environmental sanitation. The model contains five variables which are susceptible, infected, recovered, total human population, and the concentration of vibrio cholera in water. The susceptible population is generated either through birth or through immigration. They acquire infection and move to the infected class at the rate:

$\frac{\alpha \left(1-{c}_{w}\right)B}{K+B}$ (1)

*a *= rate of exposure to contaminated water;

*c _{w}* = rate of compliance with water hygiene;

*K *= concentration of vibrio cholera in water that yields 50% chance of catching cholera;

*B *= concentration of vibrio cholera in contaminated water.

The number of infected individuals decreases through recovery from the disease at the rate *τ*, being recovery due to treatment. Bacterial populations in the aquatic environment grow in the water at the rate determined by the environmental factors such as temperature. The variables and parameters of the model are described in Table 1 and Table 2: A compartmental representation of our model is given in Figure 1: The model consists of the following system of nonlinear ordinary differential equations given in (2) to (5):

Table 1. Variables of the cholera model.

Table 2. Parameters of the cholera model.

Figure 1. Compartmental model of controlled cholera transmission model without natural recovery.

$\frac{\text{d}S}{\text{d}t}=\beta N-\frac{a\left(1-{c}_{w}\right)B}{K+B}S+\omega R-\mu S$ (2)

$\frac{\text{d}I}{\text{d}t}=\frac{a\left(1-{c}_{w}\right)B}{K+B}S-\left(\tau +\mu +\delta \right)I$ (3)

$\frac{\text{d}B}{\text{d}t}=e\left(1-{c}_{s}\right)I-\rho B$ (4)

$\frac{\text{d}R}{\text{d}t}=\tau I-\left(\omega +\mu \right)R$ (5)

With

$N=S+I+R$ (6)

1.2. Remark

The variables and parameters of the model are presented in Table 1 and Table 2, while their respective baseline values are given in Table 3 and Table 4.

1.3. Model Analysis

Our Cholera model represented by the system of nonlinear differential Equations (2)-(5) is now analyzed qualitatively to obtain insights into its features, which will give a better understanding of the control features on the transmission dynamics of the cholera virus.

1.3.1. Disease Free Equilibrium (DFE)

The disease free equilibrium of model system (1)-(4) is obtained by setting

Table 3. Baseline values for variables of the cholera model in Nigeria as at 2016.

Table 4. Baseline values for parameters of the cholera model.

$\frac{\text{d}S}{\text{d}t}=\frac{\text{d}I}{\text{d}t}=\frac{\text{d}B}{\text{d}t}=\frac{\text{d}R}{\text{d}t}=0$, and in the absence of disease, $I=B=0$ so that:

$\beta N-\frac{a\left(1-{c}_{w}\right)B}{K+B}S+\omega R-\mu S=0$ (7)

$\frac{a\left(1-{c}_{w}\right)B}{K+B}S-\left(\tau +\mu +\delta \right)I=0$ (8)

$e\left(1-{c}_{s}\right)I-\rho B=0$ (9)

$\tau I-\left(\omega +\mu \right)R=0$ (10)

Since $I=B=0$, Equation (4) gives $R=0$ and Equation (1) gives $S=\frac{\beta N}{\mu}$.

Hence DFE is

$\left({S}^{0},{I}^{0},{B}^{0},{R}^{0}\right)=\left(\frac{\beta N}{\mu},0,0,0\right)$ (11)

1.3.2. Positivity of Solutions

Since the model monitors human population, we need to show that all the state variables remain non-negative for all times.

1.3.3. Theorem

Let $\Omega =\left\{\left(S,I,R,B\right)\in {\mathbb{R}}_{+}^{4}:S\left(0\right)>0,I\left(0\right)>0,R\left(0\right)>0,B\left(0\right)>0\right\}$ then the solutions of $\left\{S\left(t\right),I\left(t\right),R\left(t\right),B\left(t\right)\right\}$ of the system Equation (2)-(5) are positive for all $t\ge 0$.

*Proof*: Consider the first equation of the system;

$\frac{\text{d}S}{\text{d}t}=\beta N-\frac{a\left(1-{c}_{w}\right)B}{K+B}S+\omega R-\mu S$

$\Rightarrow \text{}\frac{\text{d}S}{\text{d}t}\ge \beta N-\mu S$

We consider the solution of

$\frac{\text{d}S}{\text{d}t}+\mu S=\beta N$

Employing the integrating factor method, we obtain

$S\left(t\right)=\frac{\beta N}{\mu}+C{\text{e}}^{-\mu t}$ (*C* = constant of integration)

With the initial condition $S\left(0\right)={S}_{0}\Rightarrow C={S}_{0}-\frac{\beta N}{\mu}$,

$\therefore \text{}S\left(t\right)=\frac{\beta N}{\mu}+\left({S}_{0}-\frac{\beta N}{\mu}\right){\text{e}}^{-\mu t}\ge 0$

We now consider the second equation;

$\frac{\text{d}I}{\text{d}t}=\frac{a\left(1-{c}_{w}\right)B}{K+B}S-\left(\tau +\mu +\delta \right)I$

$\therefore \text{}\frac{\text{d}I}{\text{d}t}\ge -\left(\tau +\mu +\delta \right)I$

Integration yields

$I\left(t\right)\ge {I}_{0}{\text{e}}^{-\left(\tau +\mu +\delta \right)t}\ge 0$

We consider also the third equation

$\frac{\text{d}B}{\text{d}t}=e\left(1-{c}_{s}\right)I-\rho B$

$\therefore \text{}\frac{\text{d}B}{\text{d}t}\ge -\rho B$

$\Rightarrow \text{}B\left(t\right)\ge {B}_{0}{\text{e}}^{-\rho Bt}$

We observe that as $t\to \infty ,B\left(t\right)\to {B}_{0}$.

Finally from the fourth equation

$\frac{\text{d}R}{\text{d}t}=\tau I-\left(\omega +\mu \right)R$

$\frac{\text{d}R}{\text{d}t}\ge -\left(\omega +\mu \right)R$

$\Rightarrow \text{}R\left(t\right)\ge {R}_{0}{\text{e}}^{-\left(\omega +\mu \right)t}\ge 0$

1.3.4. Invariant Region

The system Equations (2)-(5) have solutions which are contained in the feasible region:

$\Omega =\left\{\left(S,I,R,B\right)\in {\mathbb{R}}_{+}^{4}:0\le S\left(0\right),0\le I\left(0\right),0\le R\left(0\right),0\le B\left(0\right),S+I+R\le \frac{\beta N}{\mu}\right\}$

where the model makes biological sense, and can be shown to be positively invariant and globally attracting in ${\mathbb{R}}_{+}^{4}$ with respect to the system of ODE describing our model.

1.4. Model Analysis

The model system (2)-(5) is analyzed qualitatively to get insights into its dynamical features which give better understanding of the impact control strategies on the transmission dynamics of the cholera disease. First, we study the stability properties of the cholera model.

1.4.1. Stability Analysis of the Cholera Model

Here we consider the Jacobian matrix associated with the system of differential equations for our Cholera model

$J\left(S,I,B,R\right)=\left(\begin{array}{cccc}\frac{\partial {f}_{1}}{\partial S}& \frac{\partial {f}_{1}}{\partial I}& \frac{\partial {f}_{1}}{\partial B}& \frac{\partial {f}_{1}}{\partial R}\\ \frac{\partial {f}_{2}}{\partial S}& \frac{\partial {f}_{2}}{\partial I}& \frac{\partial {f}_{2}}{\partial B}& \frac{\partial {f}_{2}}{\partial R}\\ \frac{\partial {f}_{3}}{\partial S}& \frac{\partial {f}_{3}}{\partial I}& \frac{\partial {f}_{3}}{\partial B}& \frac{\partial {f}_{3}}{\partial R}\\ \frac{\partial {f}_{4}}{\partial S}& \frac{\partial {f}_{4}}{\partial I}& \frac{\partial {f}_{4}}{\partial B}& \frac{\partial {f}_{4}}{\partial R}\end{array}\right)$ (12)

where

${f}_{1}\left(S,I,B,R\right)=\beta N-\frac{a\left(1-{c}_{w}\right)B}{K+B}S+\omega R-\mu S$

${f}_{2}\left(S,I,B,R\right)=\frac{a\left(1-{c}_{w}\right)B}{K+B}S-\left(\tau +\mu +\delta \right)I$

${f}_{3}\left(S,I,B,R\right)=e\left(1-{c}_{s}\right)I-\rho B$

${f}_{4}\left(S,I,B,R\right)=\tau I-\left(\omega +\mu \right)R$

Computing the respective partial derivatives, we obtain;

$J\left(S,I,B,R\right)=\left(\begin{array}{cccc}-\frac{a\left(1-{c}_{w}\right)B}{K+B}-\mu & 0& -\frac{aK\left(1-{c}_{w}\right)}{{\left(K+B\right)}^{2}}& \omega \\ \frac{a\left(1-{c}_{w}\right)B}{K+B}& -\left(\tau +\mu +\delta \right)& \frac{aK\left(1-{c}_{w}\right)}{{\left(K+B\right)}^{2}}& 0\\ 0& e\left(1-{c}_{s}\right)& -\rho & 0\\ 0& \tau & 0& -\left(\omega +\mu \right)\end{array}\right)$ (13)

Evaluating the Jacobian at the disease free equilibrium state

$\left({S}^{0},{I}^{0},{B}^{0},{R}^{0}\right)=\left(\frac{\beta N}{\mu},0,0,0\right)$ gives

${J}^{*}\left({S}^{0},{I}^{0},{B}^{0},{R}^{0}\right)=\left(\begin{array}{cccc}-\mu & 0& -\frac{a\left(1-{c}_{w}\right)}{K}& \omega \\ 0& -\left(\tau +\mu +\delta \right)& \frac{a\left(1-{c}_{w}\right)}{K}& 0\\ 0& e\left(1-{c}_{s}\right)& -\rho & 0\\ 0& \tau & 0& -\left(\omega +\mu \right)\end{array}\right)$ (14)

The eigenvalues of ${J}^{*}$ are given by

$\left|{J}^{*}-\lambda I\right|=0$

Giving the characteristic equation;

${\lambda}^{4}+A{\lambda}^{3}+B{\lambda}^{2}+C\lambda +D=0$ (15)

where

$A=3\mu +\rho +\tau +\omega +\delta $

$B=\frac{3K{\mu}^{2}-ae+ae{c}_{s}+ae{c}_{w}+3K\mu \rho +2K\mu \tau +K\rho \tau +2K\mu \delta +2K\mu \omega +K\rho \delta +K\rho \omega +K\tau \omega +K\omega \delta -ae{c}_{s}{c}_{w}}{K}$

$C=\frac{\begin{array}{l}K{\mu}^{3}+3K{\mu}^{2}\rho +K{\mu}^{2}\tau +K{\mu}^{2}\omega +K{\mu}^{2}\delta -2ae\mu -ae\omega +2ae\mu {c}_{s}+2ae\mu {c}_{w}+ae\omega {c}_{s}+ae\omega {c}_{w}\\ \text{}+2K\mu \rho \tau +2K\mu \rho \omega +2K\mu \rho \delta +K\mu \tau \omega +K\rho \tau \omega +K\mu \omega \delta +K\rho \omega \delta -2ae\mu {c}_{s}{c}_{w}-ae\omega {c}_{s}{c}_{w}\end{array}}{K}$

$D=\frac{\begin{array}{l}K{\mu}^{3}\rho -a{\mu}^{2}e-a\mu \omega e+a{\mu}^{2}{c}_{s}e+a{\mu}^{2}{c}_{w}e+K{\mu}^{2}\rho \tau +K{\mu}^{2}\rho \omega +K{\mu}^{2}\rho \delta \\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+a\mu \omega {c}_{s}e+K\mu \rho \tau \omega +K\mu \rho \omega \delta -a{\mu}^{2}{c}_{s}{c}_{w}e-a\mu \omega {c}_{s}{c}_{w}e\end{array}}{K}$

1.4.2. Remark

We have asymptotic stability of the disease free equilibrium provide all the roots of the quartic Equation (15) have negative real parts.

1.5. The Basic Reproduction Ratio, *R*_{0}

The basic reproduction ratio (or number) denoted by *R*_{0} is the average number of secondary infections caused by an infectious individual during his or her entire period of infectiousness (Diekmann *et al**. * [16] ). The basic reproduction ratio is an important non-dimensional quantity in epidemiology because it sets the threshold in the study of any given disease, both for predicting its outbreak and for evaluating its control strategies.

1.5.1. The Next Generation Matrix Approach

In general*R*_{0} is found through the study and computation of the eigenvalues of the Jacobianmatrix at the disease- or infectious-free equilibrium. Diekmann *et al. * [16] follow a different approach which is *the next generation matrix *method*. *This procedure converts a system of ordinary differential equations of a model of infectious disease dynamics to an operator (or matrix) that translate from one generation of infectious individuals to the next. The basic reproductive number is then defined as the spectral radius (dominant eigenvalue) of this operator.

1.5.2. Remark

We deduce that whether a disease becomes persistent or dies out in a community is totally dependent on the value of the basic reproduction ratio, *R*_{0}. Furthermore, stability of equilibria can be analyzed using*R*_{0}. If
${R}_{0}<1$, it means that every infectious individual will cause less than one secondary infection and hence the disease will die out and when
${R}_{0}>1$, every infectious individual will cause more than one secondary infection and hence the disease will invade the population.

A large value of *R*_{0} may indicate the possibility of a major epidemic. For the case of a model with a single infected class, *R*_{0} is simply the product of the infection rate and the mean duration of the infection.

1.5.3. The Basic Reproduction Ratio for the Cholera Model

The basic reproduction ratio is obtained by taking the dominant eigenvalue of

$F{V}^{-1}=\left[\frac{\partial {F}_{i}\left({E}_{0}\right)}{\partial {x}_{j}}\right]{\left[\frac{\partial {V}_{i}\left({E}_{0}\right)}{\partial {x}_{j}}\right]}^{-1}$ (16)

where
${F}_{i}$ is the rate of appearance of new infection in compartment *i*,
${V}_{i}$ is the transfer of infections from one compartment *i* to another and
${E}_{0}$ is the disease-free equilibrium. From the system of differential equations we identify the equations with infectious classes as

$\frac{\text{d}I}{\text{d}t}=\frac{a\left(1-{c}_{w}\right)B}{K+B}S-\left(\tau +\mu +\delta \right)I$ (17)

and

$\frac{\text{d}B}{\text{d}t}=e\left(1-{c}_{s}\right)I-\rho B$ (18)

We thus have

${F}_{i}=\left[\begin{array}{c}\frac{a\left(1-{c}_{w}\right)BS}{K+B}\\ 0\end{array}\right]$ (19)

Hence

$F=\left[\begin{array}{cc}\frac{\partial {F}_{1}}{\partial I}& \frac{\partial {F}_{1}}{\partial B}\\ \frac{\partial {F}_{2}}{\partial I}& \frac{\partial {F}_{2}}{\partial B}\end{array}\right]=\left[\begin{array}{cc}0& \frac{a\left(1-{c}_{w}\right){S}^{0}}{K}\\ 0& 0\end{array}\right]$ (20)

Similarly

${V}_{i}=\left[\begin{array}{c}\left(\tau +\mu +\delta \right)I\\ -e\left(1-{c}_{s}\right)I+\rho B\end{array}\right]$ (21)

Hence

$V=\left[\begin{array}{cc}\frac{\partial {V}_{1}}{\partial I}& \frac{\partial {V}_{1}}{\partial B}\\ \frac{\partial {V}_{2}}{\partial I}& \frac{\partial {V}_{2}}{\partial B}\end{array}\right]=\left[\begin{array}{cc}\left(\tau +\mu +\delta \right)& 0\\ -e\left(1-{c}_{s}\right)& \rho \end{array}\right]$ (22)

$\Rightarrow \text{}{V}^{-1}={\left[\begin{array}{cc}\left(\tau +\mu +\delta \right)& 0\\ -e\left(1-{c}_{s}\right)& \rho \end{array}\right]}^{-1}=\left[\begin{array}{cc}\frac{1}{\left(\tau +\mu +\delta \right)}& 0\\ \frac{e\left(1-{c}_{s}\right)}{\rho \left(\tau +\mu +\delta \right)}& \frac{1}{\rho}\end{array}\right]$ (23)

$\therefore \text{}F{V}^{-1}=\left[\begin{array}{cc}0& \frac{a\left(1-{c}_{w}\right){S}^{0}}{K}\\ 0& 0\end{array}\right]\left[\begin{array}{cc}\frac{1}{\left(\tau +\mu +\delta \right)}& 0\\ \frac{e\left(1-{c}_{s}\right)}{\rho \left(\tau +\mu +\delta \right)}& \frac{1}{\rho}\end{array}\right]=\left[\begin{array}{cc}\frac{ae\left(1-{c}_{w}\right)\left(1-{c}_{s}\right){S}^{0}}{K\left(\tau +\mu +\delta \right)}& \frac{a\left(1-{c}_{w}\right){S}^{0}}{K\rho}\\ 0& 0\end{array}\right]$

where ${S}^{0}=\frac{\beta N}{\mu}$. The eigenvalues of $F{V}^{-1}$ are given by

$\left|F{V}^{-1}-\lambda I\right|=0$ (24)

or

$\left|\left[\begin{array}{cc}\frac{ae\left(1-{c}_{w}\right)\left(1-{c}_{s}\right){S}^{0}}{K\left(\tau +\mu +\delta \right)}& \frac{a\left(1-{c}_{w}\right){S}^{0}}{K\rho}\\ 0& 0\end{array}\right]-\lambda \left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]\right|=0$ (25)

Giving the characteristic equation;

$\frac{K\rho \left(\tau +\mu +\delta \right){\lambda}^{2}-{S}^{0}ae\left(1-{c}_{w}\right)\left(1-{c}_{s}\right)\lambda}{K\rho \left(\tau +\mu +\delta \right)}=0$

i.e.,

$K\rho \left(\tau +\mu +\delta \right){\lambda}^{2}-{S}^{0}ae\left(1-{c}_{w}\right)\left(1-{c}_{s}\right)\lambda =0$ (26)

$\therefore \text{}\lambda \left[K\rho \left(\tau +\mu +\delta \right)\lambda -{S}^{0}ae\left(1-{c}_{w}\right)\left(1-{c}_{s}\right)\right]=0$

$\Rightarrow \text{}\lambda =0,\text{}\lambda =\frac{ae\left(1-{c}_{w}\right)\left(1-{c}_{s}\right){S}^{0}}{K\rho \left(\tau +\mu +\delta \right)}$ (27)

Hence

${R}_{0}=\frac{ae\left(1-{c}_{w}\right)\left(1-{c}_{s}\right){S}^{0}}{K\rho \left(\tau +\mu +\delta \right)}$ (28)

Substituting ${S}^{0}=\frac{\beta N}{\mu}$.

$\Rightarrow \text{}{R}_{0}=\frac{\beta Nae\left(1-{c}_{w}\right)\left(1-{c}_{s}\right)}{K\rho \mu \left(\tau +\mu +\delta \right)}$ (29)

1.6. Control Measures with Various Strategies

As shown in Equations (2) to (5), our model incorporates three new control measures namely treatment, water hygiene and environmental sanitation. In this study, it is of main interest to see the effect of control measures at different levels of control strategies. Using the baseline values for variables and parameters as in Table 3 and Table 4 we compute the effect of control measures *C** _{w}*,

1) No control: This is a situation where there is no control ( $\tau ={C}_{w}={C}_{s}=0$ ).

In this case

${R}_{0}=\frac{\beta Nae}{K\rho \mu \left(\mu +\delta \right)}$ (30)

2) Low without treatment: Low efficiency level of the rate of compliance with water hygiene and rate of compliance with environmental sanitation at 25% ( $\tau =0$, ${C}_{w}={C}_{s}=0.25$ ).

$\Rightarrow \text{}{R}_{0}=\frac{9\beta Nae}{16K\rho \mu \left(\mu +\delta \right)}$ (31)

3) Moderate without treatment: Moderate effectiveness level of rate of compliance with water hygiene and rate of compliance with environmental sanitation at 50% ( $\tau =0$, ${C}_{w}={C}_{s}=0.50$ ).

$\Rightarrow \text{}{R}_{0}=\frac{\beta Nae}{4K\rho \mu \left(\mu +\delta \right)}$ (32)

4) High without treatment: High effectiveness level of rate of compliance with water hygiene and rate of compliance with environmental sanitation at 75% ( $\tau =0$, ${C}_{w}={C}_{s}=0.75$ ).

$\Rightarrow \text{}{R}_{0}=\frac{\beta Nae}{16K\rho \mu \left(\mu +\delta \right)}$ (33)

5) Low universal: Low level of all the three control strategies at 25% ( $\tau ={C}_{w}={C}_{s}=0.25$ ).

$\Rightarrow \text{}{R}_{0}=\frac{9\beta Nae}{4K\rho \mu \left(1+4\mu +4\delta \right)}$ (34)

6) Moderate universal: Moderate level of all the three control strategies at 50% ( $\tau ={C}_{w}={C}_{s}=0.50$ ).

$\Rightarrow \text{}{R}_{0}=\frac{\beta Nae}{2K\rho \mu \left(1+2\mu +2\delta \right)}$ (35)

7) High universal: High level of all the three control strategies at 75% ( $\tau ={C}_{w}={C}_{s}=0.75$ ).

$\Rightarrow \text{}{R}_{0}=\frac{\beta Nae}{4K\rho \mu \left(3+4\mu +4\delta \right)}$ (36)

1.7. MathCAD Simulation

$\beta :=0.000107$, $a:=0.00000289$, $K:={10}^{6}$, $\tau :=0.2$, $\rho :=0.33$, $\omega :=0.003$, $\mu :=0.0000526$, $\delta :=0.2$, $\Omega :=740$, ${c}_{w}:=0.25$, ${c}_{s}:=0.25$, $e:=10$.

Define a function that determines a vector of derivative values at any solution point (t, Y):

$D\left(t,Y\right):=\left[\begin{array}{c}\beta \cdot \Omega -\frac{a\cdot \left(1-{c}_{w}\right)\cdot {Y}_{2}\cdot {Y}_{0}}{K+{Y}_{2}}+\omega \cdot {Y}_{3}-\mu \cdot {Y}_{0}\\ \frac{a\cdot \left(1-{c}_{w}\right)\cdot {Y}_{2}}{K+{Y}_{2}}\cdot {Y}_{0}-\left(\tau +\mu +\delta \right)\cdot {Y}_{1}\\ e\cdot \left(1-{c}_{s}\right)\cdot {Y}_{1}-\rho \cdot {Y}_{2}\\ \tau \cdot {Y}_{1}-\left(\omega +\mu \right)\cdot {Y}_{3}\end{array}\right]$

Define additional arguments for the ODE solver:

$\text{t0:}=\text{0}$ : Initial value of independent variable.

$\text{Y0:}=\left(\begin{array}{c}640\\ 100\\ {10}^{6}\\ 0\end{array}\right)$ : Vector of initial function values.

$\text{t1:}=1\text{0}$ : Final value of independent variable.

$\text{num:}=1\times {10}^{3}$ : Number of solution values on [t0, t1].

$\text{t:}=\text{S}{1}^{\langle 0\rangle}$ : Independent variable values.

$\text{S:}=\text{S}{1}^{\langle 1\rangle}$ : First solution function values.

$\text{I:}=\text{S}{1}^{\langle 2\rangle}$ : Second solution function values.

$\text{B:}=\text{S}{1}^{\langle 3\rangle}$ : Third solution function values.

$\text{R:}=\text{S}{1}^{\langle 4\rangle}$ : Fourth solution function values.

The table of values for the matrix S is given in Figure 2.

1.8. Results and Discussion

The graphical profiles in Figures 3(a)-(d) depict respectively the variation of the

Figure 2. Cross section of solution matrix.

Figure 3. (a) Variation of susceptible population with time. (b) Variation of infective population with time. (c) Variation of concentration of vibriocholera population with time. (d) Variation of recovered human population with time.

intrinsic variables of the our model, namely the susceptible class *S*, the infected class *I*, the concentration of vibrio cholera population *B*, and the removed class *R*. The variation of the susceptible population shows a sharp rise from the initial value over time. However, the infective population deceases exponentially over the same time period, while the concentration of vibrio cholera decreases exponentially to zero due, to intervention measures. This in turn enables an increase in the recovered population, which reaches a saturation level over time.

1.9. Conclusion

An *SIR-B* model was modified and analyzed. The model incorporates treatment, water hygiene as well as environmental sanitation. The model employs a system of nonlinear ordinary differential equations, which is studied in detail for its stability properties. The basic reproduction ratio *R*_{0} was computed for the various control parameters and it is discovered that with proper combination of control measures, the spread of cholera could be minimized. Numerical simulation of the cholera model was performed using the numerical software MathCAD14, and the graphical profiles of the intrinsic variables are depicted. We conclude that improvement in treatment, water hygiene and the environmental sanitation is indeed effective in eradicating the cholera epidemic.

References

[1] WHO (World Health Organization) (2013) Cholera 2012. The Weekly Epidemiological Record, 88, 321-336.

[2] Codecco, C.T. (2001) Endemic and Epidemic Dynamics of Cholera: The Role of the Aquatic Reservoir. BMC Infectious Diseases, 1, Article No. 1.

https://doi.org/10.1186/1471-2334-1-1

[3] Wearing, H.J., Rohani, P. and Keeling, M.J. (2005) Appropriate Models for Management of Infectious Disease. PLOS Medicine, 2, e174.

https://doi.org/10.1371/journal.pmed.0020174

[4] WHO (World Health Organization) (2013) Media Center: Cholera.

https://www.who.int/news-room/fact-sheets/detail/cholera

[5] NBS (National Bureau of Statistics) (2009) Social Statistics in Nigeria.

https://www.nigerianstat.gov.ng/

[6] WHO (World Health Organization) (2012) Global Task on Cholera Control. Cholera Country Profile: Nigeria.

https://www.who.int/cholera/publications/final%20outbreak%20booklet%20260105-OMS.pdf

[7] Tian, J.P. and Wang, J. (2011) Global Stability for Cholera Epidemic Models. Mathematical Biosciences, 232, 341.

https://doi.org/10.1016/j.mbs.2011.04.001

[8] Epstein, P.R. (1993) Algal Blooms in the Spread and Persistence of Cholera. Biosystems, 31, 209-221.

https://doi.org/10.1016/0303-2647(93)90050-M

[9] Lawoyin, T.O., Ogunbode, W.A., Olumide, E.A.A. and Onadeko, M.O. (2004) Outbreak of Cholera in Ibadan, Nigeria. European Journal of Epidemiology, 15, 365-368.

https://doi.org/10.1023/A:1007547117763

[10] Capasso, V. and Paveri-Fontana, S.L. (1979) A Mathematical Model for the 1973 Cholera Epidemic in the European Mediterranean Region. Revue d’Epidémiologie et de Santé Publique, 27, 121-132.

[11] Pascual, M., Bouma, M.J. and Dobson, A.P. (2002) Cholera and Climate: Revisiting the Quantitative Evidence. Microbes and Infection, 4, 237-245.

https://doi.org/10.1016/S1286-4579(01)01533-7

[12] Hartley, D.M., Morris, J.B. and Smith, D.L. (2006) Hyperinfectivity: A Critical Element in the Ability of V. cholerae to Cause Epidemics. PLOS Medicine, 3, e7.

https://doi.org/10.1371/journal.pmed.0030007

[13] Tien, J.H. and Earn, D.J.D. (2010) Multiple Transmission Pathways and Disease Dynamics in a Waterborne Pathogen Model. Bulletin of Mathematical Biology, 72, 1506-1533.

https://doi.org/10.1007/s11538-010-9507-6

[14] Misra, A.K. and Singh, V. (2012) A Delay Mathematical Model for the Spread and Control of Water Borne Diseases. Journal of Theoretical Biology, 301, 49-56.

https://doi.org/10.1016/j.jtbi.2012.02.006

[15] Fatima, S., Krishnarajah, I., Jaffar, M.Z.A.M. and Adam, M.B. (2014) A Mathematical Model for the Control of Cholera in Nigeria. Research Journal of Environmental and Earth Sciences, 6, 321-325.

https://doi.org/10.19026/rjees.6.5776

[16] Diekman, O., Heesterbeek, J.A.P. and Metz, J.A.P. (1990) On the Definition and Computation of the Basic Reproduction Ratio in the Model of Infectious Disease in Heterogeneous Populations. Journal of Mathematical Biology, 2, 265-382.