Abstract
In this paper, we study the qualitative behavior of a discretetime hostpathogen model for spread of an infectious disease with permanent immunity. The timestep is equal to the duration of the infectious phase. Moreover, the local asymptotic stability, the global behavior of unique positive equilibrium point, and the rate of convergence of positive solutions is discussed. Some numerical examples are given to verify our theoretical results.
MSC: 39A10, 40A05.
Keywords:
difference equations; local stability; global character1 Introduction
It is a wellknown fact that in the population growth, the disease is an important agent controlling the population dynamics. Many experiments show that parasites can reasonably reduce the host population and even take the host population to complete annihilation. This natural phenomenon is successfully modeled by many simple SI type hostparasite models. The most interesting properties of such models are their ability of generating host annihilation dynamics with the ideal parametric values and initial conditions. This is possible, because such models naturally contain the proportion transmission term, which is often referred to as ratiodependent functional response in the case of predatorprey models. In the SI model, the population is subdivided into two classes, susceptibles S and infectives I. The notation SI means that there is a transfer from the susceptible to infective class, susceptibles become infective and do not recover from the infection. Thus, the transfer continues until all individuals become infected. This type of model is very simple, but may represent some complicated dynamical properties. Most of the SI type models consist of the mass action principle, i.e., the assumption that the new cases arise in a simple proportion to the product of the number of individuals which are susceptible and the number of which are infectious. However, this principle has a limited validity and in the discrete models, this principle leads to biologically irrelevant results, unless some restrictions are suggested for the parameters. It is more appropriate for discrete epidemic models to include an exponential factor in the rate of transmission. Exponential difference equations can be used to study the models in population dynamics [13]. We consider here a simple exponential discretetime hostpathogen model for spread of an infectious disease with permanent immunity. The timestep is equal to the duration of the infectious phase. The state variables are , the number of susceptible individuals at time n, and representing the number of individuals, getting the disease (new cases) between times and n,
where β is the number of births between n and , all added to the susceptible class and assumed to be constant over time. So, the difference equation is just a ‘conservation of mass’ for the susceptible class. The first part of the model is just like NicholsonBailey; it comes from assuming that each susceptible escapes infection with probability ; the more infectives there are, the lower the chance of escape. The model ignores mortality in the susceptible class, on the assumption that everyone gets the disease while young, and mortality occurs later in life.
Many ecological competition models are governed by differential and difference equations. We refer to [4,5] and the references therein for some interesting results, related to the global character and local asymptotic stability. As it is pointed out in [6,7], the discrete time models governed by difference equations are more appropriate than the continuous ones when the populations are of nonoverlapping generations. The study of the discretetime models described by difference equations has now been given a great attention, since these models are more reasonable than the continuous time models when populations have nonoverlapping generations. Discretetime models give rise to more efficient computational models for numerical simulations and also show rich dynamics compared to the continuous ones. In recent years, many papers have been published on the mathematical models of biology that discussed the system of difference equations generated from the associated system of differential equations as well as the associated numerical methods. Mathematical models of epidemics have created a major area of research interest during the last few decades. Recently, theory on the effects of parasites on host population dynamics has received much attention and epidemiological models are often used to explain empirical results for hostparasites interaction system. For more details of such biological models, one can see [810].
More precisely, our aim is to investigate local asymptotic stability of unique positive equilibrium point, the global asymptotic character of equilibrium point, and the rate of convergence of positive solutions of system (1). For more results for the systems of difference equations, we refer the reader to [1114].
2 Boundedness and persistence
The following theorem shows that every positive solution of system (1) is bounded and persists.
Theorem 1Every positive solutionof system (1) is bounded and persists.
Proof Let be any positive solution of system (1). It is easy to see that and for all . Then, it follows that and . Consider the following difference equation
with initial condition , then its solution is given by
Then, by comparison we have . Hence,
and
Theorem 2Letbe a positive solution of the system (1). Then, is an invariant set for system (1).
Proof Let be a positive solution of system (1) with initial conditions and . Then, from system (1)
and
Similarly, we have
and
Hence, and . Suppose that the result is true for , i.e., and . Then from system (1), one can easily obtain
and
Hence, the proof is completed. □
3 Linearized stability
Let us consider twodimensional discrete dynamical system of the form
where and are continuously differentiable functions and I, J are some intervals of real numbers. Furthermore, a solution of system (2) is uniquely determined by initial conditions . An equilibrium point of (2) is a point that satisfies
Let be an equilibrium point of a map , where f and g are continuously differentiable functions at . The linearized system of (2) about the equilibrium point is
where and is Jacobian matrix of system (2) about the equilibrium point .
Let be the equilibrium point of system (1), then one has
Then, it follows that is the unique positive equilibrium point of system (1). Moreover, the Jacobian matrix of system (1) about the equilibrium point is given by
The characteristic polynomial of is given by
Lemma 1[15]
Consider the seconddegree polynomial equation
wherepandqare real numbers. Then, the necessary and sufficient condition for both roots of Equation (4) to lie inside the open diskis
Lemma 2[16]
Assume that, , is a system of difference equations andis the fixed point ofF. If all eigenvalues of the Jacobian matrixaboutlie inside the open unit disk, thenis locally asymptotically stable. If one of them has a modulus greater than one, thenis unstable.
Theorem 3Assume that. Then, the unique positive equilibrium pointis locally asymptotically stable.
Proof The characteristic polynomial of about positive equilibrium point is given by
Let
Assume that , and . Then, one has
Then, by Rouche’s theorem and have the same number of zeroes in an open unit disk . Hence, both roots
and
of (5) lie in an open disk , and it follows from Lemma 2 that the equilibrium point is locally asymptotically stable. □
The following theorem shows the necessary and sufficient condition for the local asymptotic stability of a unique positive equilibrium point of system (1).
Theorem 4The unique positive equilibrium pointof system (1) is locally asymptotically stable if and only if.
Proof Let and , then (5) can be written as
Then, and if and only if . Hence, from Lemma 1, the unique positive equilibrium point of system (1) is locally asymptotically stable if and only if . □
4 Global character
The following lemma is similar to Theorem 1.16 of [15].
Lemma 3Letandbe real intervals, and letandbe continuous functions. Consider system (2) with initial conditions. Suppose that the following statements are true:
(i) is nondecreasing in both arguments.
(ii) is nonincreasing inx, and nondecreasing iny.
(iii) Ifis a solution of the system
Then, there exists exactly one equilibrium pointof the system (2) such that.
Theorem 5The unique positive equilibrium pointof system (1) is a global attractor.
Proof Let , and . Then, it is easy to see that is nondecreasing in both x and y. Moreover, is nonincreasing in x, and nondecreasing in y. Let be a solution of the system
Then, one has
and
From system (6), one has
From (7), one has
Furthermore, assuming as in the proof of Theorem 1.16 of [15], it suffices to suppose that
Using the fact that , one has from (8) and (9)
It follows from (10) that . Then, (7) implies that
Using (6) in (11), we obtain
Following the same technique as in the proof of Proposition 4.1 of [17], one has , i.e., . Thus, and, similarly, one can show that . Hence, from Lemma 3, the equilibrium point of system (1) is a global attractor. □
Lemma 4The unique positive equilibrium pointof system (1) is globally asymptotically stable if and only if.
Proof The proof follows from Theorem 4 and Theorem 5. □
5 Rate of convergence
In this section, we determine the rate of convergence of a solution that converges to the unique positive equilibrium point of system (1). Similar methods can be found in [18] and [19].
The following result gives the rate of convergence of solutions of a system of difference equations
where is an mdimensional vector, is a constant matrix, and is a matrix function satisfying
as , where denotes any matrix norm, which is associated with the vector norm
Proposition 1 (Perron’s theorem [20])
Suppose that condition (13) holds. Ifis a solution of (12), then eitherfor all largenor
exists and is equal to the modulus of one the eigenvalues of matrixA.
Proposition 2[20]
Suppose that condition (13) holds. Ifis a solution of (12), then eitherfor all largenor
exists and is equal to the modulus of one the eigenvalues of matrixA.
Let be any solution of system (1) such that , and . To find the error terms, one has from system (1)
and
and
where
Moreover,
Now, the limiting system of error terms can be written as
which is similar to linearized system of (1) about the equilibrium point .
Using Proposition 1, one has the following result.
Theorem 6Assume thatis a positive solution of system (1) such thatand, whereis a unique positive equilibrium point of (1). Then, the error vectorof every solution of (1) satisfies both of the following asymptotic relations
6 Examples
In order to verify our theoretical results and to support our theoretical discussions, we consider several interesting numerical examples in this section. These examples represent different types of qualitative behavior of solutions to the system of nonlinear difference equations (1). All plots in this section are drawn with Mathematica.
Example 1 Let , and let . Then, system (1) can be written as
In this case, the unique equilibrium point . Moreover, in Figure 1, the plot of is shown in Figure 1(a), the plot of is shown in Figure 1(b), and an attractor of system (16) is shown in Figure 1(c). The basic reproductive number of system (16) is .
Figure 1. Plots for system (16).
Example 2 Let , and let . Then, system (1) can be written as
In this case, the unique equilibrium point . Moreover, in Figure 2, the plot of is shown in Figure 2(a), the plot of is shown in Figure 2(b), and an attractor of system (16) is shown in Figure 2(c). The basic reproductive number of system (17) is .
Figure 2. Plots for system (17).
Example 3 Let , and let . Then, system (1) can be written as
In this case, the unique equilibrium point . Moreover, in Figure 3, the plot of is shown in Figure 3(a), the plot of is shown in Figure 3(b), and an attractor of system (18) is shown in Figure 3(c). The basic reproductive number of system (18) is .
Figure 3. Plots for system (18).
Example 4 Let , and let . Then, system (1) can be written as
with the initial conditions , .
In this case, the unique equilibrium point . Moreover, in Figure 4, the plot of is shown in Figure 4(a), the plot of is shown in Figure 4(b), and an attractor of the system (19) is shown in Figure 4(c). The basic reproductive number of system (19) is .
Figure 4. Plots for system (19).
Example 5 Let , and let . Then, system (1) can be written as
with the initial conditions , .
In this case, the unique equilibrium point . Moreover, in Figure 5, the plot of is shown in Figure 5(a), the plot of is shown in Figure 5(b), and an attractor of system (20) is shown in Figure 5(c). The basic reproductive number of system (20) is .
Figure 5. Plots for system (20).
Conclusion and future work
This work is related to the qualitative behavior of an exponential discretetime hostpathogen model for spread of an infectious disease with permanent immunity. We proved that system (1) has a unique positive equilibrium point, which is locally asymptotically stable. The main objective of dynamical systems theory is to predict the global behavior of a system based on the knowledge of its present state. An approach to this problem consists of determining the possible global behaviors of the system and determining which initial conditions lead to these longterm behaviors. In the paper, a general result for global character for such type of systems is proved. Due to the simplicity of our SItype model, we have carried out a systematic local and global stability analysis of it. The most important finding here is that the unique positive equilibrium point can be a global asymptotic attractor for model (1). Moreover, the rate convergence of positive solutions has also been investigated. In such models, there is a threshold parameter that might tell whether a population will increase or die out, or whether an infectious disease will persist or die out within a population. This parameter is commonly known as the basic reproductive number and is denoted by . In epidemiology, this number is defined as the number of newly infected individual, produced by a single infected individual in its period of infectivity. In case of system (1), the basic reproductive number is given by . From our investigations, it is obvious that the unique positive equilibrium point of system (1) is globally asymptotically stable if , and unstable if . Some numerical examples are provided to support our theoretical results. These examples are experimental verifications of theoretical discussions. The qualitative behavior of the general model, where there is host mortality at some constant rate, will be our next aim to study.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All authors contributed equally and approved the final version of the manuscript.
Acknowledgements
The authors would like to thank the anonymous referees for their valuable comments and suggestions leading to improvement of this paper. This work was supported by the Higher Education Commission of Pakistan.
References

ElMetwally, E, Grove, EA, Ladas, G, Levins, R, Radin, M: On the difference equation . Nonlinear Anal.. 47, 4623–4634 (2001). Publisher Full Text

Papaschinopoulos, G, Radin, MA, Schinas, CJ: On the system of two difference equations of exponential form: , . Math. Comput. Model.. 54, 2969–2977 (2011). Publisher Full Text

Papaschinopoulos, G, Schinas, CJ: On the dynamics of two exponential type systems of difference equations. Comput. Math. Appl.. 64(7), 2326–2334 (2012)

Ahmad, S: On the nonautonomous LotkaVolterra competition equation. Proc. Am. Math. Soc.. 117, 199–204 (1993). Publisher Full Text

Tang, X, Zou, X: On positive periodic solutions of LotkaVolterra competition systems with deviating arguments. Proc. Am. Math. Soc.. 134, 2967–2974 (2006). Publisher Full Text

Zhou, Z, Zou, X: Stable periodic solutions in a discrete periodic logistic equation. Appl. Math. Lett.. 16(2), 165–171 (2003). Publisher Full Text

Liu, X: A note on the existence of periodic solution in discrete predatorprey models. Appl. Math. Model.. 34, 2477–2483 (2010). Publisher Full Text

Allen, LJS: An Introduction to Mathematical Biology, Prentice Hall, New York (2007)

Brauer, F, CastilloChavez, C: Mathematical Models in Population Biology and Epidemiology, Springer, Berlin (2000)

EdelsteinKeshet, L: Mathematical Models in Biology, McGrawHill, New York (1988)

Kalabuŝić, S, Kulenović, MRS, Pilav, E: Dynamics of a twodimensional system of rational difference equations of LeslieGower type. Adv. Differ. Equ. (2011). BioMed Central Full Text

Din, Q: Dynamics of a discrete LotkaVolterra model. Adv. Differ. Equ.. 2013, Article ID 95 (2013)

Din, Q, Donchev, T: Global character of a hostparasite model. Chaos Solitons Fractals. 54, 1–7 (2013)

Din, Q, Qureshi, MN, Khan, AQ: Dynamics of a fourthorder system of rational difference equations. Adv. Differ. Equ.. 2012, Article ID 215 (2012)

Grove, EA, Ladas, G: Periodicities in Nonlinear Difference Equations, Chapman & Hall/CRC Press, Boca Raton (2004)

Sedaghat, H: Nonlinear Difference Equations: Theory with Applications to Social Science Models, Kluwer Academic, Dordrecht (2003)

Papaschinopoulos, G, Radin, MA, Schinas, CJ: Study of the asymptotic behavior of the solutions of three systems of difference equations of exponential form. Appl. Math. Comput.. 218, 5310–5318 (2012). Publisher Full Text

Kulenović, MRS, Nurkanović, M: Asymptotic behavior of a competitive system of linear fractional difference equations. Adv. Differ. Equ.. 2006, Article ID 019756 (2006)

Clark, CA, Kulenović, MRS, Selgrade, JF: On a system of rational difference equations. J. Differ. Equ. Appl.. 11(7), 565–580 (2005). Publisher Full Text

Pituk, M: More on Poincare’s and Perron’s theorems for difference equations. J. Differ. Equ. Appl.. 8, 201–216 (2002). Publisher Full Text