Abstract
In this paper, we construct a nonstandard finite difference (NSFD) scheme for an SIR epidemic model of childhood disease with constant strategy. The dynamics of the obtained discrete model is investigated. First we show that the discrete model has equilibria which are exactly the same as those of the continuous model. Furthermore, we prove that the conditions for those equilibria to be globally asymptotically stable are consistent with the continuous model for any size of numerical timestep. The analytical results are confirmed by some numerical simulations.
Keywords:
mathematical model; transmission dynamics; basic reproduction number; sensitivity analysis; control strategies1 Introduction
Childhood diseases are the most common form of infectious diseases. These are the diseases such as measles, mumps, chicken pox, rubella, poliomyelitis, etc. to which children are born susceptible and usually contract within five years. Because young children are in frequent contact with each other at school or other place, such a disease can be spread very quickly. Meanwhile, the development of vaccines against infectious children diseases has been booming and protecting children from the diseases. Hence vaccination is considered to be the most effective strategy against childhood diseases, it is essential for us to predict the optimal vaccine coverage level to prevent the spread of theses diseases. A universal effort to extend vaccination coverage to all children began in 1974, when the World Health Organization (WHO) founded the Expanded Program on Immunization (EPI). Mathematical models (see [17]) of deterministic type have often been used to provide deeper insights into the transmission dynamics of a childhood disease and to evaluate control strategies.
In this paper, the total population that is involved in the spread of infection is split into three epidemiological classes: a susceptible class (S), an infected class (I) and a removed class (R) denoting vaccinated as well as recovered people with permanent immunity. We assume that the efficacy of vaccine is 100%, and the natural death rates μ in the classes remain unequal to births, so that the total population N is realistically not constant. Citizens are born into the population at a constant birth rate A with extremely low childhood disease mortality rate. We denote the fraction of citizens vaccinated at birth each year as p (with ) and assume the rest are susceptible. A susceptible individual will move into the infected group through contact with an infected individual, approximated by an average contact rate β. An infected individual recovers at rate γ, and enters the removed class. The removed class also contains people who are vaccinated. The differential equations for the SIR (see [1,4,7]) epidemic model of childhood diseases with constant vaccination strategy are as follows:
The biological background requires that all parameters be nonnegative. Makinde [4] employed the Adomian decomposition method to compute an approximate nonperturbative solutions of model (1.1). Yildirim and Cherruault [7] by qualitative analysis revealed the vaccination reproductive number for disease control and eradication.
However, for practical purposes, it is often necessary to discretize the continuous model. The discrete dynamical system obtained from the discretization should contain as many qualitative properties of the continuous problem as possible. It is shown that many standard methods such as Euler method, RungeKutta method and some other standard finite schemes implemented in a dynamical system can lead to negative solutions for spurious dynamical behaviors such as converging to wrong equilibrium point or wrong periodic cycle or numerical instabilities [810]. In this paper, we propose a numerical scheme to solve model (1.1) by implementing a nonstandard finite difference (NSFD) scheme. This method was originally developed by Mickens [1116]. The nonstandard scheme relied on the following important rules: the standard denominator h in standard discrete derivative is replaced by a denominator function , where ; the nonlinear terms are approximated in a nonlocal way using more than one mesh point. Here, h is the timestep size of numerical integration. Moreover, the fundamental principle for constructing NSFD scheme for differential equations is dynamic consistency, that is, the discretized model maintain essential dynamical properties such as positivity of solutions, boundedness of solutions, monotonicity of solutions, correct number and stability of fixedpoints and other special solutions of the continuous model. This method has been applied to various problems in which the resulting discrete systems preserve dynamical properties of the related continuous models [5,1719]. In [19], the NSFD scheme has been implemented in a special class of SIR epidemic models. Mickens [5] considered a SIR epidemic model with squareroot dynamics.
This paper is organized as follows. In the next section, we present several important properties of solutions to the continuous model. A particular discretization is constructed in Section 3. We illustrate the global asymptotic stability of diseasefree equilibrium and endemic equilibrium in Sections 4 and 5, and provided numerical examples to verify our results in Section 6. Finally, we provide a summary of the obtained results and present a possible extension of this work.
2 Dynamical properties of the continuous model
In this section, we review some dynamical properties of the SIR epidemic model with constant vaccination strategy (1.1). It should be noted that the system of continuous equations, if added together, satisfies the conservation law [16]
The exact solution of equation (2.1) is
where . Any solution , , of model (1.1) satisfies
All valid epidemic models must have this feature since negative population numbers cannot exist as a physical reality. Moreover, it is easy to verify that the domain
is a compact, positively invariant set for model (1.1).
Define the basic reproduction number as follows:
Then the following results can be summarized in Li et al.[3].
Theorem 1Model (1.1) always has a diseasefree equilibriumand has a unique endemic equilibriumwhen, where
Theorem 2For model (1.1), the following results hold.
(i) If, the diseasefree equilibriumis globally asymptotically stable inD. On the other hand, if, is unstable;
(ii) If, then the unique endemic equilibriumis globally asymptotically stable in D.
3 The NSFD scheme
To derive the NSFD scheme, we define the following notation [16]:
where is the constant timestep size. By applying the NSFD scheme to model (1.1), we can obtain the following discretetime SIR epidemic model with constant vaccination strategy
where the denominator function (see [11,12]) is
It is noted that discretized model (3.1) can be considered as using the standard Euler method for the first derivative and a nonlocal expression for nonlinear terms. It can be easily rearranged to get its explicit version
Since all the parameters in model (3.2) are positive, it is clear that if the initial values , and are positive, then the numerical solutions will also be positive for all , namely
Defining and adding the three equations of model (3.1), we get
which is the exact finite scheme for the conservation law as expressed by equation (2.1) (see [16]). A straightforward calculation gives
Similar to the continuoustime case, discrete model (3.1) or equivalent (3.2) also has a compact, positivity invariant set
It is easy to verify that discrete model (3.1) or equivalent (3.2) has the same equilibrium as model (1.1) which is independent of h. It can be described as the following theorem.
Theorem 3For model (3.1) or equivalent (3.2), there always exists a diseasefree equilibriumand has a unique endemic equilibriumwhen, where
4 Global asymptotic stability of diseasefree equilibrium
In this section, we mainly discuss the global asymptotic stability of diseasefree equilibrium . First, we consider the local stability of for model (3.1). In order to obtain the local stability of equilibria for discrete model (3.1) or equivalent (3.2), we notice that the first two equations in model (3.1) or (3.2) do not depend on the third equation, and therefore the third equation can be omitted without loss of generality properties. For the sake of simplicity, we define the following functions [19]:
Obviously, the Jacobian matrix [20,21] at the equilibrium point is given by
Now, we give the following theorem about the local stability of the diseasefree equilibrium for model (3.1).
Theorem 4If, the diseasefree equilibrium pointof discrete model (3.1) or equivalent (3.2) is locally asymptotically stable in. On the other hand, if, is unstable.
Proof Substituting the diseasefree equilibrium to Jacobian matrix (4.1) yields
The eigenvalues of are as follows:
Obviously for all h. From the definition of basic reproductive number (2.4), we can easily conclude that is equivalent to . Thus, if , then the magnitude of eigenvalue is also strictly less than unity irrespective of h. This completes the proof. □
Next, let us consider the global asymptotic stability of diseasefree equilibrium for model (3.1). Since the variable R does not appear in the first and the second equations, it is sufficient to consider the following 2dimensional system:
Set
Similar to the Izzo et al. [[22], proof of Lemma 3.3], we obtain the following basic lemma.
Lemma 1For any solutionof model (4.2), with the initial conditions, , we have that
Furthermore, we define the function as follows:
then we easily get the following lemma.
Lemma 2is strictly monotone increasing function onwith
and
Moreover, if, then there exists a unique solutionofsuch that
By applying techniques in Izzo et al.[22], we now prove the global stability of the diseasefree equilibrium for .
Theorem 5If, then the diseasefree equilibriumof model (3.1) is globally asymptotically stable.
Proof From (4.3) in Lemma 1, for any , there exists an integer such that
Construct the following sequence defined by
Since ϵ is arbitrary, we conclude that if , then
which yields that the sequence is monotone decreasing. Therefore, there exists a nonnegative constant such that . We will prove that for . In fact, if , then we have . Consequently, from (4.5), we have that . Meanwhile, by Lemma 1, we can conclude that . Suppose that , that is, . We first transform model (4.2) into the following form:
We claim that there exists a sequence such that if , the claim is evident. Now we consider the case . By applying the first inequality in (4.4), we obtain . Therefore, there exists a sequence such that . From the first equation of (4.6), we have
that is,
from which it is not difficult to obtain
and
We easily find that (4.7) implies . Meanwhile, combining (4.8) with (4.5), we have . Therefore, we have , which yields . Since , we have .
Finally, we will prove that if , then is uniformly stable. First, we consider the case that there exists a nonnegative integer such that for any . From the first equation of (3.1), we have that for any ,
Meanwhile, from the second equation of model (3.1), we have that for any ,
From the third equation of model (3.1), it is not difficult to obtain
which implies that for any such that
Next, we consider the case that there exists a nonnegative integer such that . Then, by the first equation of (3.1), we obtain
Thus, we have, for , . Then, by the second equation of (3.1) and , we have that for any ,
Moreover, by the first equation of model (3.1), we have, for ,
By applying (4.11), we conclude that (4.10) holds for any . Thus, from (4.9)(4.12), we conclude that is uniformly stable. Hence, if , is globally asymptotically stable. □
5 Global asymptotic stability of the endemic equilibrium
In this section, we mainly discuss the global dynamics of the endemic equilibrium of model (3.1). Before we prove the stability of the endemic equilibrium, we first give the following lemma.
The quadratic equationhas two roots that satisfy, , if and only if the following conditions are satisfied:
Theorem 6If, then the endemic equilibrium pointof discrete model (3.1) or equivalent (3.2) is locally asymptotically stable in.
Proof Assuming that and substituting the endemic equilibrium point to Jacobian matrix (4.1) lead to
The characteristic equation of is given by , straightforward calculation gives
Obviously,
where , that is, . Finally, using the fact that , it is easy to obtain
Hence all the conditions in Lemma 3 are satisfied when . This proves that when , then the endemic equilibrium point is locally asymptotically stable for any h. □
Next, we will prove the permanence of model (3.1) for . Similar to the result of McCluskey in [25,26], we first give the following lemma.
Lemma 4Ifholds, then; inversely, ifholds, then, where
Proof By the second equation of model (3.1), we obtain
This completes the proof of Lemma 4. □
Theorem 7If, then for any solution of model (3.1) with the initial conditions that, ,
where the constantis sufficiently large such that.
Proof For any positive constant ϵ, there exists a sufficiently large positive integer such that
By the first equation of models (3.1) and (5.3), we have that for any ,
Since ϵ is arbitrary, we conclude that (5.1) holds. Now we prove that (5.2) holds. In fact, for any positive constant , it is seen that
We first claim that any solution of model (3.1) does not have the following property: there exists a nonnegative integer such that for all . Suppose on the contrary that there exists a solution of model (3.1) and a nonnegative such that for , it can be seen that for ,
Since , therefore, there exists a positive such that for all ,
We hence set
Thus, by applying Lemma 1, we have that for all .
Furthermore, for the sequence defined by (4.5), we have that for ,
Since , this leads to , which yields a contradiction. Hence the claim is proved.
By the claim above, we are left to consider two possibilities. First, for all n sufficiently large. If this case holds, we get the conclusion of the proof. Second, we investigate the case that oscillates about for all n sufficiently large. Let be sufficiently large such that
By the second equation of model (3.1), we obtain, for ,
If , then by applying a similar discussion above, we obtain for all . We hence prove that for all . Since the interval is arbitrarily chosen, we conclude that for all n sufficiently large. Meanwhile, since q is also arbitrary, we conclude that . This completes the proof. □
By Theorem 6, we easily obtain the permanence of model (3.1) for . Next, by constructing the Lyapunov function, we will prove that is globally asymptotically stable for . Consider the Lyapunov function as follows:
where the function . First, by calculating, we have
In the same way, we have
From (5.4) and (5.5), we have
which implies that is a monotone decreasing sequence. We easily get
Then , which combined with (5.6) implies that
From the first equation of models (3.1) and (5.7), we have
Notice that , which implies that is uniformly stable. Finally, we hence obtain the theorem as follows.
Theorem 8If, then the endemic equilibriumfor model (3.1) or equivalent (3.2) is globally asymptotically stable.
6 Numerical simulations
In this section, numerical simulations will be given to verify theoretical results obtained in the previous section. The simulation is performed using MATLAB software.
(i) We choose , , , , and . By calculation, we have that and the endemic equilibrium . According to Theorem 5, the diseasefree equilibrium of discrete model (3.1) or equivalent (3.2) is globally stable, which is shown in Figure 1.
Figure 1. The solutionsof model (3.1) or equivalent (3.2) are globally asymptotically stable and converge to the diseasefree equilibrium, when.
(ii) Assuming the following parameter values: , , , , and , by calculation, we have and the endemic equilibrium . According to Theorem 8, the endemic equilibrium of discrete model (3.1) or equivalent (3.2) is globally stable, which is depicted in Figure 2.
Figure 2. The solutionsof model (3.1) or equivalent (3.2) are globally asymptotically stable and converge to the diseasefree equilibrium, when.
7 Conclusion
In this paper, we have proposed a discretetime analogue of the continuous SIR epidemic model of childhood diseases with constant vaccination strategy which is derived by the NSFD scheme of Michens. In order to obtain the permanence of model (3.1) for , we offer Lemma 4. Applying the discrete Lyapunov functional technique (see [25,26]) for both cases and , it shown that the global dynamics of this discretetime analogue of the continuous SIR epidemic model is fully determined only by the basic reproduction number . This shows dynamical consistency between the discrete SIR epidemic model and its corresponding continuous model. The NSFD scheme constructed in this paper is for the SIR epidemic model with constant vaccination strategy. For our future work, we will consider an epidemic model with varying vaccination strategy.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
The authors declare that the study was realized in collaboration with the same responsibility. All authors read and approved the final manuscript.
Acknowledgements
The work was supported by the National Natural Science Foundation of P.R. China (11201399, 11301451) and the Natural Science Foundation of Shihezi University (2013ZRKXYQYD05).
References

Arafa, AAM, Rida, SZ, Khalil, M: Solutions of fractional order model of childhood diseases with constant vaccination strategy. Math. Sci. Lett.. 1, 17–23 (2013)

Cui, Q, Yang, X, Zhang, Q: An NSFD scheme for a class of SIR epidemic model with vaccination and treatment. J. Differ. Equ. Appl.. 20, 416–422 (2014). Publisher Full Text

Li, J, Zhang, J, Ma, Z: Global analysis of some epidemic models with general contact rate and constant immigration. Appl. Math. Mech.. 4, 396–404 (2004)

Makinde, OD: Adomian decomposition approach to a SIR epidemic model with constant vaccination strategy. Appl. Math. Comput.. 184, 842–848 (2007). Publisher Full Text

Mickens, RE: A SIRmodel with squareroot dynamics: an NSFD scheme. J. Differ. Equ. Appl.. 16, 209–216 (2010). Publisher Full Text

Wang, L, Cui, Q, Teng, Z: Global dynamics in a class of discretetime epidemic models with disease courses. Adv. Differ. Equ.. 2013, (2013) Article ID 57. http://www.advancesindifferenceequations.com/content/2013/1/57

Yildirim, A, Cherruault, Y: Analytical approximate solution of a SIR epidemic model with constant vaccination strategy by homotopy perturbation method. Kybernetes. 38, 1566–1575 (2009). Publisher Full Text

Hu, Z, Teng, Z, Jiang, H: Stability analysis in a class of discrete SIRS epidemic models. Nonlinear Anal., Real World Appl.. 13, 2017–2033 (2012). Publisher Full Text

Suryanto, A: A dynamically consistent nonstandard numerical scheme for epidemic model with saturated incidence rate. Int. J. Math. Comput.. 13, 112–123 (2011)

Suryanto, A: Stability and bifurcation of a discrete SIS epidemic model with delay. Proceedings of the 2nd International Conference on Basic Sciences. 1–6 (2012) Indonesia

Mickens, RE: Nonstandard Finite Difference Model of Differential Equations, World Scientific, Singapore (1994)

Mickens, RE: Application of Nonstandard Finite Difference Schemes, World Scientific, Singapore (2000)

Mickens, RE: Dynamic consistency: a fundamental principle for constructing nonstandard finite difference schemes for differential equations. J. Differ. Equ. Appl.. 11, 645–653 (2005). Publisher Full Text

Mickens, RE: Calculation of denominator functions for nonstandard finite difference schemes for differential equations satisfying a positivity condition. Numer. Methods Partial Differ. Equ.. 3, 528–534 (2012)

Mickens, RE: Nonstandard finite difference schemes for differential equations. J. Differ. Equ. Appl.. 8, 823–847 (2002). Publisher Full Text

Mickens, RE: Numerical integration of population models satisfying conservation laws: NSFD methods. J. Biol. Dyn.. 1, 427–436 (2007). PubMed Abstract  Publisher Full Text

Ding, X: A nonstandard finite difference scheme for an epidemic model with vaccination. J. Differ. Equ. Appl.. 19, 179–190 (2013). Publisher Full Text

Mickens, RE, Washington, T: A note on an NSFD scheme for a mathematical model of respiratory virus transmission. J. Differ. Equ. Appl.. 8, 525–529 (2012)

Suryanto, A, Kusumawinahyu, WM, Darti, I, Yanti, I: Dynamically consistent discrete epidemic model with modified saturated incidence rate. Comput. Appl. Math.. 32, 373–383 (2013). Publisher Full Text

Strogatz, SH: Nonlinear Dynamics and Chaos, AddisonWesley, Reading (1994)

Izzo, G, Muroya, Y, Vecchio, A: A general discrete time model of population dynamics in the presence of an infection. Discrete Dyn. Nat. Soc.. 2009, (2009) Article ID 143019
Article ID 143019
Publisher Full Text 
Brauer, F, CastilloChavez, C: Mathematical Models in Population Biology and Epidemiology, Springer, New York (2001)

Elaydi, S: An Introduction to Difference Equations, Springer, New York (1992)

McCluskey, CC: Global stability for an SEIR epidemiological model with varying infectivity and infinite delay. Math. Biosci. Eng.. 6, 603–610 (2009). PubMed Abstract

McCluskey, CC: Complete global stability for an SIR epidemic model with delay  distributed or discrete. Nonlinear Anal., Real World Appl.. 10, 55–59 (2010)