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 . 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  . 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 . 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 . Also in 2009, Nigeria reported 13,691 cases and 431 deaths  (WHO). Tian and Wang  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,  . 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      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. . 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:
a = rate of exposure to contaminated water;
cw = 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.
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.
, and in the absence of disease, so that:
Since , Equation (4) gives and Equation (1) gives .
Hence DFE is
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.
Let then the solutions of of the system Equation (2)-(5) are positive for all .
Proof: Consider the first equation of the system;
We consider the solution of
Employing the integrating factor method, we obtain
(C = constant of integration)
With the initial condition ,
We now consider the second equation;
We consider also the third equation
We observe that as .
Finally from the fourth equation
1.3.4. Invariant Region
The system Equations (2)-(5) have solutions which are contained in the feasible region:
where the model makes biological sense, and can be shown to be positively invariant and globally attracting in 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
Computing the respective partial derivatives, we obtain;
Evaluating the Jacobian at the disease free equilibrium state
The eigenvalues of are given by
Giving the characteristic equation;
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, R0
The basic reproduction ratio (or number) denoted by R0 is the average number of secondary infections caused by an infectious individual during his or her entire period of infectiousness (Diekmann et al.  ). 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 generalR0 is found through the study and computation of the eigenvalues of the Jacobianmatrix at the disease- or infectious-free equilibrium. Diekmann et al.  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.
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, R0. Furthermore, stability of equilibria can be analyzed usingR0. If , it means that every infectious individual will cause less than one secondary infection and hence the disease will die out and when , every infectious individual will cause more than one secondary infection and hence the disease will invade the population.
A large value of R0 may indicate the possibility of a major epidemic. For the case of a model with a single infected class, R0 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
where is the rate of appearance of new infection in compartment i, is the transfer of infections from one compartment i to another and is the disease-free equilibrium. From the system of differential equations we identify the equations with infectious classes as
We thus have
where . The eigenvalues of are given by
Giving the characteristic equation;
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 Cw, Cs, and τ at seven different control strategies as explained below. The seven control strategies are:
1) No control: This is a situation where there is no control ( ).
In this case
2) Low without treatment: Low efficiency level of the rate of compliance with water hygiene and rate of compliance with environmental sanitation at 25% ( , ).
3) Moderate without treatment: Moderate effectiveness level of rate of compliance with water hygiene and rate of compliance with environmental sanitation at 50% ( , ).
4) High without treatment: High effectiveness level of rate of compliance with water hygiene and rate of compliance with environmental sanitation at 75% ( , ).
5) Low universal: Low level of all the three control strategies at 25% ( ).
6) Moderate universal: Moderate level of all the three control strategies at 50% ( ).
7) High universal: High level of all the three control strategies at 75% ( ).
1.7. MathCAD Simulation
, , , , , , , , , , , .
Define a function that determines a vector of derivative values at any solution point (t, Y):
Define additional arguments for the ODE solver:
: Initial value of independent variable.
: Vector of initial function values.
: Final value of independent variable.
: Number of solution values on [t0, t1].
: Independent variable values.
: First solution function values.
: Second solution function values.
: Third solution function values.
: 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.
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 R0 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.
 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.
 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.
 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.
 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.
 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.
 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.