In certain epidemiological modeling, the population is generally divided into three classes which are susceptible represented by S, infected individual represented by I and recovered individual represented by R. This kind of mathematical model is noted SIR. Recently, many authors have studied the dynamical behavior of epidemic models (see    and references therein). There are two kinds of mathematical models: The continuous-time models described by differential equations, and the discrete-time models described by difference equations. The simplest forms of these models are Ordinary Differential Equations (ODEs)   . In  , a discrete delay model is given to account for transmission by vectors (e.g. mosquitoes), where the delay is used to account for a latent period in the vector. Allowing the vectors latency periods to vary according to some distribution gives a model with a distributed delay  .
The delay appears in the incidence term which is typically the only non-linearity, and is therefore the “cause” of all “interesting behavior”. Various forms have been used for the incidence term, both for ODEs and for delay equations. Common forms include mass action    , saturating
incidence   , and standard (or proportional) incidence
 . Changing the form of the incidence function can potentially change the behavior of the system.
In this paper we study the discrete mathematical model which result from the continuous-time model presented and study in  by C. Connell McCluskey.
From this we use the general incidence term , where is a time delay. We choose the constant so that . The discrete
model is obtained by using the backward Euler method.
The studied of discrete epidemic models is motivate by the fact that there occur situations such that constructing discrete epidemic models is more appropriate approach to understand disease transmission dynamics and to evaluate eradication policies because they permit arbitrary time-step units, preserving the basic features of corresponding continuous-time models  . Furthermore, this allows better use of statistical data for numerical simulations due to the reason that the infection data are compiled at discrete given time intervals. For a discrete epidemic model with immigration of infectives, Jang and Elaydi  showed the global asymptotic stability of the disease-free equilibrium, the local asymptotic stability of the endemic equilibrium and the strong persistence of susceptible class by means of the nonstandard discretization method. In he’s recent work, using a discretization called “ mixed type” formula in Izzo and Vecchio  and Izzo et al.  , Sekiguchi  obtained the permanence of a class of SIR discrete epidemic models with one delay and SEIRS (Susceptible-Latent-Infected-Recovered-Susceptible) discrete epidemic model with two delay if an endemic equilibrium of each model exists.
This paper is organized as follows. In Section 2, we give the discrete model, the equilibrium point and the reproduction rate . In Section 3, the positivity and boundedness of the solution of system (3) are obtained. In addition we proved the existence and uniqueness of disease-free equilibrium and endemic equilibrium . In Section 4, we study the stability of disease-free equilibrium point for . In Section 5, we study the global stability of the endemic equilibrium point for . In section 6, we give the numerical result and their comment. In the last part we give the conclusion.
2. Discrete Mathematical Model
In this section we describe the discrete mathematical model derived by the continuous time model study in  , by C. Connell McCluskey. This continuous time model is given by:
1) A population is divided into susceptible, infectious and recovered classes with sizes , and respectively.
2) B is the recruitment of new individuals, it is into the susceptible class.
3) , and denote respectively the death rates of susceptible, infectious and recovered class.
4) The total exit rate for infectious is , which, for biological reasons we assume is at least as large as ; that is, .
5) The incidence at time t is where the maximum delay
, k is a Lebesgue integral function which gives the relative infectivity of vectors of different infection ages. We choose so that . It is assumed that the support of k has positive measure in any open interval having supremum h so the interval of integration is not artificially extended by concluding with an interval for which the integral is automatically zero.
The form of the function f is of fundamental importance. In this paper we use a general incidence function used in one of he’s discrete version. So we use assumption:
H1 f is a non-negative differentiable function on the non-negative quadrant. Furthermore, f is positive if and only if both arguments are positive.
H2 for all .
The partial derivative of f are denoted by and from the first and second variable.
H3 for all n.
H4 for all n.
Now, we use the backward Euler difference scheme to discretize the model (1). The time step size of our discretization is one. Thus, we obtain the following discrete SIR epidemic model with nonlinear general incidence given by:
Since R does not appear in the first and second equations of system above, it is sufficient to analyses the behavior of solutions of the following system:
The constants and the relation between this constants are given above. In the discrete model the incidence function at time t is
where the maximum delay . Let be a equilibrium point model of (3) so we have,
By adding the equations of system above we get
Let and be respectively disease-free equilibrium and endemic equilibrium point of model (3). The disease-free equilibrium correspond to the case where the infectious class is nil . Thus, we have ; with
. The endemic equilibrium is given by: ; with
Proposition 2.1. The basic reproduction number is given by .
Proof: The Jacobian matrix of system (3) at equilibrium is define by:
Thus, we have:
3. Basic Properties
We suppose that initial condition of system (3) satisfy:
where the space of continuous functions from to equipped with the sup norm: . Standard theory of functional differential equations  can be used to show that the solution of (1) exist and are differentiable for all . We assume any initial condition for which the disease is initially present satisfies for all .
Lemma 3.1. Let be a solution of system (3) with initial condition (11), then we have and for all n.
Proof: Assume that and . From system (3) we have the following system:
By using second equation of system above and the fact that , we have, . So . From the non-negativity of we used the assumption H4. Thus, for all n.
Lemma 3.2. Any solution of system (3), with initial condition (11) satisfy
Proof: Let and ; for biological reasons we assume is at least as large as ; that is, .
By adding the different equation of system (3) we get:
Hence, we have
1) When , then model (3) has only a unique disease-free equilibrium .
2) When , then model (3) has only a unique endemic equilibrium .
Proof: Any equilibrium point of system (3) verified the following system:
By using the second equation of system above we have:
So we can consider the function G defined by,
Hence we have
And also we have,
when , we have . Consequently, there is not any
such that . Therefore, model (3) has a unique disease-free equilibrium .
When, , we have . Therefore, there exists a unique such that .
This implies that model (3) has unique endemic equilibrium .
Remark 3.1. The space is positively invariant and attracting domain for system (3).
Now, let us analyze the behavior of system (3) when the basic reproduction rate is less than one.
4. Stability of the Disease-Free Equilibrium
In this section, we study the stability of diseases-free equilibrium , with .
Theorem 4.1. If , then the diseases-free equilibrium of system (3) is locally asymptotically stable.
Proof: The linearization of system (3) at diseases-free equilibrium point is given by:
Thus, we have:
The matrix M associate of the linearization (23) is given by:
and the linearization system can be rewrite by:
The model (3) is locally asymptotically stable at diseases-free equilibrium point if all eigenvalue of matrix M is greater than one.
Let be characteristic polynomial associate of matrix M. We have,
Let be a eigenvalue of matrix M, thus . This implies that:
From the Equation (27) we have
By using the second member of (28), the fact that and we assume that the matrix M have a eigenvalue which is less than one. So we have:
as a result of, the Equation (28) cannot have roots which is less than one. Hence, is locally asymptotically stable according to the theorem 2 in  .
Theorem 4.2. When , the disease-free equilibrium of system (3) is globally asymptotically stable in K.
Proof: In this proof we used the comparison theorem  . By using the second equation of system (3) and the assumption H3; we get:
Hence we have:
By using the fact we have . So the constant M is greater than one. we conclude that the linearized Equation (32) is stable whenever . By a standard comparison theorem  , as for Equation (32) and substituting in system (3) we get , as . Thus, as for system (3), when . Therefore is globally asymptotically stable in the positively invariant set K if .
5. Global Stability of the Endemic Equilibrium
In this section, we study the stability the stability of endemic equilibrium , with
Theorem 5.1. If , then the endemic equilibrium of system (3) is globally asymptotically stable.
Proof: From the equation of system (3), at endemic equilibrium , we have:
which will be used as substitutions in the calculations below. Let and
We will study the behavior of the Lyapunov functional
which satisfies with equality if and only if
for all . For clarity, the difference , and will be calculated separately and then combined to obtain .
Calculation of the variation : in this calculation, we used the value theorem and we assume that . Note that we have the same result when .
Let us calculate of the variation : in this calculation we used the mane value theorem and we assume that . Note that we have the same result when .
Let now evaluate the variation :
By adding Equations (42)-(44) we obtain
By adding and subtracting to (45) we obtain:
By using the assumption H4 and the fact that the function g is increasing on , we have
so we have . This implies that . Hence, by the Lyapunovs theorems on the global asymptotical stability for difference equations  , we obtain that the endemic equilibrium is globally asymptotically stable.
6. Simulation and Comments
In this section, we presented a numerical result of continuous-time model (1) and the discrete one (2) study above. From this we used a particular incidence
function define by , which is the saturating incidence. In
addition we discuss from the different value of the basic reproduction number . We have the case and . The parameters values used in the simulation are:
from this value we have .
When we change the value of by , we get . It’ is important to notice that the software used is Scilab and the time is in term of weeks or months. In our graphic the red curve give the evolution of the class in the discrete model and the dashed ones give the evolution of the class in continuous-time model.
Figure 1 present the evolution of the susceptibles population through the time, the dashed cuve represent the discrete model and the red one the continuous model when .
Figure 2 give dynamic of the infected population along the time, the dashed cuve represent the discrete model and the red one the continuous model when .
Figure 3 show the evolution of the recovered population through the time, the dashed cuve represent the discrete model and the red one the continuous model when .
Figure 4 (Susceptibles population), Figure 5 (Infected population) and Figure 6 (Recovered population) represent the evolution through the time of the population when . The dashed cuves represent the discrete model and the red one the continuous model.
For all these cuves, we can see the convergence of the red cuves (the discrete model) and the dashed ones (the continuous model)
In this paper, we have studied a discrete SIR epidemic model with general
Figure 1. Graphic of susceptibles class, when .
Figure 2. Graphic of infectious class, when .
Figure 3. Graphic of recovered class, when .
Figure 4. Graphic of susceptibles class, when .
Figure 5. Graphic of infectious class, when .
Figure 6. Graphic of recovered class, when .
incidence. We have proved the global stability of discrete SIR epidemic model by using the comparison theorem from the global stability of disease free equilibrium, when , on the positive invariant set K and we have also proved the local stability of disease free equilibrium. The technique of Lyapunov function is used to proved the global stability of endemic equilibrium, when . We have made the numerical simulation to corroborate theoretical results. From the results obtained in this paper, we can conclude that the Euler backward difference scheme, that is, the discrete dynamical model (2), is obtained with excellent dynamical properties for the step size in the local and global stability of equilibra. These properties are nearly the same as the corresponding continuous-time model (1). In our future work, it shall be important for us to study the same model, but with general positive step size and see how bifurcation can happen.
The authors want to thank the anonymous referee for his valuable comments on the paper.
Aboudramane Guiro provide the subject, wrote the introduction and the conclusion and verified some calculation. Dramane Ouédraogo conceived the study and computed the equilibria and their local stabilities. Harouna Ouédraogo rote mathematical formula, bring up the Lyapunov functional and did all the calculus with the other authors. All the authors read and approved the final manuscript.