Research

# An NSFD scheme for SIR epidemic models of childhood diseases with constant vaccination strategy

Qianqian Cui12, Jiabo Xu3, Qiang Zhang2 and Kai Wang4*

Author Affiliations

1 Department of Applied Mathematics, Nanjing University of Science and Technology, Nanjing, 210094, People’s Republic of China

2 College of Sciences, Shihezi University, Shihezi, 832000, People’s Republic of China

3 Xinjiang Institute of Engineering, Urumqi, 830091, People’s Republic of China

4 Department of Medical Engineering and Technology, Xinjiang Medical University, Urumqi, 830011, People’s Republic of China

For all author emails, please log on.

Advances in Difference Equations 2014, 2014:172  doi:10.1186/1687-1847-2014-172

 Received: 10 December 2013 Accepted: 4 June 2014 Published: 24 June 2014

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

### 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 time-step. The analytical results are confirmed by some numerical simulations.

##### Keywords:
mathematical model; transmission dynamics; basic reproduction number; sensitivity analysis; control strategies

### 1 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 [1-7]) 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:

(1.1)

The biological background requires that all parameters be nonnegative. Makinde [4] employed the Adomian decomposition method to compute an approximate non-perturbative 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, Runge-Kutta 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 [8-10]. 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 [11-16]. 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 time-step 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 fixed-points 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,17-19]. In [19], the NSFD scheme has been implemented in a special class of SIR epidemic models. Mickens [5] considered a SIR epidemic model with square-root 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 disease-free 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]

(2.1)

The exact solution of equation (2.1) is

(2.2)

where . Any solution , , of model (1.1) satisfies

(2.3)

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:

(2.4)

Then the following results can be summarized in Li et al.[3].

Theorem 1Model (1.1) always has a disease-free equilibriumand has a unique endemic equilibriumwhen, where

Theorem 2For model (1.1), the following results hold.

(i) If, the disease-free 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 time-step size. By applying the NSFD scheme to model (1.1), we can obtain the following discrete-time SIR epidemic model with constant vaccination strategy

(3.1)

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

(3.2)

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

(3.3)

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 continuous-time 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 disease-free equilibriumand has a unique endemic equilibriumwhen, where

### 4 Global asymptotic stability of disease-free equilibrium

In this section, we mainly discuss the global asymptotic stability of disease-free 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

(4.1)

Now, we give the following theorem about the local stability of the disease-free equilibrium for model (3.1).

Theorem 4If, the disease-free 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 disease-free 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 disease-free 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 2-dimensional system:

(4.2)

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

(4.3)

Furthermore, we define the function as follows:

then we easily get the following lemma.

Lemma 2is strictly monotone increasing function onwith

and

(4.4)

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 disease-free equilibrium for .

Theorem 5If, then the disease-free 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

(4.5)

Then, for , we have

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:

(4.6)

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,

As , we obtain

from which it is not difficult to obtain

(4.7)

and

(4.8)

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 ,

then, for any , we obtain

which implies that for any ,

(4.9)

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

(4.10)

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 ,

(4.11)

Moreover, by the first equation of model (3.1), we have, for ,

which implies that for ,

(4.12)

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.

Lemma 3[23,24]

The quadratic equationhas two roots that satisfy, , if and only if the following conditions are satisfied:

(i) ,

(ii) ,

(iii) .

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

where , , .

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

If holds, then we have that

This completes the proof of Lemma 4. □

Theorem 7If, then for any solution of model (3.1) with the initial conditions that, ,

(5.1)

(5.2)

where the constantis sufficiently large such that.

Proof For any positive constant ϵ, there exists a sufficiently large positive integer such that

(5.3)

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 ,

Consequently, for , we have

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 ,

Thus, we have that for all ,

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

(5.4)

In the same way, we have

(5.5)

From (5.4) and (5.5), we have

(5.6)

which implies that is a monotone decreasing sequence. We easily get

Then , which combined with (5.6) implies that

(5.7)

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 disease-free 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 disease-free 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 disease-free equilibrium, when.

### 7 Conclusion

In this paper, we have proposed a discrete-time 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 discrete-time 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 (2013ZRKXYQ-YD05).

### References

1. 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)

2. 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

3. 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)

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

5. Mickens, RE: A SIR-model with square-root dynamics: an NSFD scheme. J. Differ. Equ. Appl.. 16, 209–216 (2010). Publisher Full Text

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

7. 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

8. 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

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

10. 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

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

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

13. 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

14. 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)

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

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

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

18. 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)

19. 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

20. Ross, SL: Differential Equations, Blaisdell, Waltham (1964)

22. 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

23. Brauer, F, Castillo-Chavez, C: Mathematical Models in Population Biology and Epidemiology, Springer, New York (2001)

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

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

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