About this Journal Submit a Manuscript Table of Contents
Abstract and Applied Analysis
Volume 2014 (2014), Article ID 576365, 13 pages
http://dx.doi.org/10.1155/2014/576365
Research Article

Analysis of a Patch Model for the Dynamical Transmission of Echinococcosis

1Department of Mathematics, North University of China, Taiyuan, Shanxi 030051, China
2Department of Medical Engineering and Technology, Xinjiang Medical University, Urumqi 830011, China
3College of Mathematics and System Sciences, Xinjiang University, Urumqi 830046, China

Received 18 January 2014; Accepted 4 February 2014; Published 18 March 2014

Academic Editor: Weiming Wang

Copyright © 2014 Kai Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

A patch model for echinococcosis due to dogs migration is proposed to explore the effect of dogs migration among patches on the spread of echinococcosis. We firstly define the basic reproduction number . The mathematical results show that the dynamics of the model can be completely determined by . If , the disease-free equilibrium is globally asymptotically stable. When , the model is permanence and endemic equilibrium is globally asymptotically stable. According to the simulations, it is shown that the larger diffusion of dogs from the lower epidemic areas to the higher prevalence areas can intensify the spread of echinococcosis. However, the larger diffusion of dogs from the higher prevalence areas to the lower epidemic areas can reduce the spread and is beneficial for disease control.

1. Introduction

Echinococcosis, which is often referred to as hydatid disease, is a parasitic disease that affects both humans and other mammals, such as sheep, dogs, rodents, and horses [1]. The two most clinically relevant species are Echinococcus granulosus and Echinococcus multilocularis, which cause cystic and alveolar echinococcosis respectively. Humans are incidental hosts and, in most cases, do not contribute to continuance of the parasite life cycle, except under unique circumstances [2].

The prevalent scope of echinococcosis in China is approximately 420 square kilometers, accounting for about 41.7% of the territory. The rate of incidence of echinococcosis has increased in the past decade. The operability of echinococcosis exceeds 10/100000 in each year. High-risk group subject to echinococcosis reaches up to 50 million, and the number of domestic animal amount being faced with the infection of echinococcosis is more than one hundred million, in which the amount of dogs is at least 5 million [3].

Mathematical modeling has become an important tool in analyzing the epidemiological characteristics of infectious disease and can provide useful control measures. Various models have been used to study different aspects of echinococcosis [416]. The models included varied primarily on the basis of six key features that were differentially incorporated in their design [17]. These are the inclusion of a “latent” class (with time delay from host exposure to infectiousness); an age structure for definitive and/or intermediate hosts; the presence of density dependent constraints; accounting for seasonality; stochastic parameters; inclusion of a spatial and risk structures.

In [18], in order to explore effective control and prevention measures authors proposed a deterministic model to study the transmission dynamics of echinococcosis in Xinjiang. The results showed that the dynamics of the model was completely determined by the basic reproductive number . The model provided an approximate estimate of the basic reproduction number .

Many epidemic models with population dispersal among patches have been proposed and studied (see [1928]). Wang and Zhao [19] proposed an epidemic model to describe the dynamics of disease spread among patches due to population dispersal. The effect of population dispersal among patches on the spread of a disease was investigated by Jin and Wang in [20]. To understand the effect of transport-related infection on disease spread, an epidemic model for several regions which are connected by transportation of individuals has been proposed by Cui et al. in [21]. In [23], an SIS patch model with nonconstant transmission coefficients was formulated to investigate the effect of media coverage and human movement on the spread of infectious diseases among patches. Qiu [26] developed a mathematical model to explore the effect of host migration between two patches on the spread of a vector-host disease.

To date, few scholars have researched the echinococcosis transmission models with dogs migration among patches. Considering an increasing number of stray dogs, the dispersal is an essential trait for dogs population. Therefore, we expect to explore the effect of dogs migration among patches on the spread of echinococcosis.

The purpose of this paper is to model the transmission dynamics of echinococcosis spread between two patches due to dogs migration and describe the dynamics of the model. The remaining part of this paper is organized as follows. The model is presented in Section 2. The basic properties on the positivity and boundedness of solutions computing the basic reproduction number are in Section 3. In Section 4, we establish the global stability of the disease-free equilibrium for the model. In Section 5, we will apply the theory of permanence to obtain the permanence of the model. The global stability theorem of endemic equilibrium is stated and proved in Section 6. In Section 7, we give some examples to illustrate how the dogs migration affects the dynamics of echinococcosis. A brief discussion is given in Section 8.

2. Model Formulation

In this section, we mainly formulate an epidemic model to describe the transmission dynamics of echinococcosis spread between two discrete patches due to dogs diffusion.

We firstly formulate a model for the spread of echinococcosis in the th patch. It follows from [18] that the parameters of humans do not affect dynamical behaviors of echinococcosis model. Hence in the paper we only consider dogs, livestock, and Echinococcus eggs in our model. We divide the dogs population in the th patch into two classes: the susceptible population and the infected population denoted by and , respectively. For livestock population, we divide the total livestock population in the th patch into two classes: susceptible and infectious denoted by and , respectively. The density of Echinococcus eggs in the th patch is denoted by . Our assumptions on the dynamical transmission of echinococcosis in the th patch are demonstrated in the flowchart (Figure 1).

576365.fig.001
Figure 1: Transmission diagram for echinococcosis among dogs, livestock.

If there is no dogs migration among patches, that is, the patches are isolated, we suppose that the echinococcosis dynamics in th patch is governed by All parameters are assumed positive. For the dog population in the th patch, describes the annual recruitment rate; is the natural death rate; denotes the recovery rate of transition from infected to noninfected dogs, including natural recovery rate and recovery due to anthelmintic treatment; describes the transmission of echinococcosis between susceptible dogs and infectious livestock after the ingestion of cyst-containing organs of the infected livestock. For the livestock population in the th patch, is the annual recruitment rate; is the death rate; describes the transmission of echinococcosis to livestock by the ingestion of Echinococcus eggs in the environment. For Echinococcus eggs in the th patch, denotes released rate from infected dogs; is the mortality rate of eggs.

When two patches are connected, we assume that susceptible and infected dogs of every patch leave for patch at a per capita rate . Then the dynamics of echinococcosis is governed by the following model:

Motivated by biological background of model (2), we always assume that all solutions of model (2) satisfy the following positive initial conditions:

We can easily prove that the solution of model (2) with initial conditions (3) satisfies , , , and for all . Here, we omit the proof.

3. Basic Properties and Basic Reproduction Number of the Model

In this section, we mainly present the preliminary results and derive reproduction number for model (2). In order to investigate the dynamics of model (2), we begin with stating some results on model (1). Model (1) has been analyzed in [18]. Model (1) admits a disease-free equilibrium and a unique positive equilibrium , where The reproduction number of model (1) is established in [18], which can be expressed as From Theorems 3 and 5 in [18], we can obtain the following lemma.

Lemma 1. Considering model (1), one has that(a)if , then disease-free equilibrium is globally asymptotically stable;(b)if , then positive equilibrium is globally asymptotically stable.

In order to obtain our main results, we need the following lemma. Consider the following linear equation: We have the following result on system (6).

Lemma 2. System (6) has a unique equilibrium which is globally stable, where

Proof. The Jacobian matrix of (6) at is By simple calculations, the corresponding characteristic equation is where Therefore, all roots of have negative real parts, and hence is globally stable.

For any , we define region as follows: On the ultimate boundedness of solutions for model (2), we have the following result.

Lemma 3. All solutions of model (2) with initial condition (3) ultimately turn into region as .

Proof. Let be any solution of model (2) with initial conditions (3) and let , . From model (2) we have and then from Lemma 2 we have . Hence, for any , there is a such that From the third and fourth equations of model (2), we have and therefore, there exists a such that Finally, from the fifth equation of model (2), we have and then there is a such that Let , and then for all we have This completes the proof of Lemma 3.

According to Lemma 3, all feasible solutions of model (2) enter or remain in the region as becomes large enough. In what follows, the dynamics of model (2) can be considered only in .

Simple algebraic calculation shows that model (2) always has a unique disease-free equilibrium . According to the concepts of next generation matrix and reproduction number presented in [29, 30], we define Noting that the disease-free equilibrium of model (2) is , then where Denote . After extensive algebraic calculations, we can obtain whereFrom the proof of Theorem 2 in [30], it follows that whereand is the maximum real part of the eigenvalues of matrix .

Using Theorem 2 in [30], we can easily obtain the following stability result.

Theorem 4. For model (2), one has that(a)if , then disease-free equilibrium is locally asymptotically stable;(b)if , then disease-free equilibrium is unstable.

4. Global Stability of the Disease-Free Equilibrium

We start by considering the global stability of disease-free equilibrium when .

Theorem 5. The disease-free equilibrium of model (2) is globally asymptotically stable in if .

Proof. From Theorem 4 we find that disease-free equilibrium is locally asymptotically stable if . In the following we only need to prove the global attractiveness of . From (24) we can see that if , then . Hence, there is a small enough number such that , where and
Let be any solution of model (2) in , then From model (2), it follows that Define an auxiliary linear system: Since system (29) is a linear system, the globally stability of origin is determined by the stability of matrix . Since , then all the eigenvalues of matrix have negative real parts. It then follows that each solution of (29) satisfies By the comparison principle we have Then the limiting system of model (2) is By Lemma 2 we find that there is a unique equilibrium of system (32), which is globally asymptotically stable. Thus, according to the theory of asymptotic autonomous systems [31], we finally obtain that disease-free equilibrium is globally asymptotically stable for model (2) when . This completes the proof of Theorem 5.

5. Permanence

We now turn to the case where . We first establish the permanence for model (2).

Theorem 6. Let , . If , then model (2) is permanent. Furthermore, model (2) also has at least one positive equilibrium .

Proof. Define satisfies model (2),
In order to prove Theorem 6, it suffices to show that repels uniformly the solutions of .
Firstly, by the form of model (2), it is easy to see that both and are positively invariant. Clearly, is relatively closed in . Furthermore, model (2) is point dissipative (see Lemma 3).
We now show that if , , then Assume It suffices to show that Suppose not, then there exists a such that at least one of , , , , , or is greater than zero. Here we only consider the case , , , , , and , . The other case can be deduced in the same way. Since it follows that there is an small enough such that , , , and , , for all . If , then we have This means that for all . If , it then follows from model (2) that It then follows that there exists an such that By the same way we can obtain that there exists an such that If , then we have This means that for all ; if , it then follows from model (2) that It then follows that there exists an such that By the same way we can obtain that there exists an such that Thus for all we have , , , , and , . This contradicts the assumption that . This proves (35).
From (24) we can see that if , then . Hence, there is a small enough number such that , where and is given by (26). Let and we can see the fact that . Hence we can choose small enough such that Suppose is a solution of model (2) with . We now claim that Suppose, for the sake of contradiction, that there exists a such that , , and , , for all . Then by model (2) we have for . Consider the following auxiliary system: As in our analysis in Lemma 2, system (51) has a unique positive equilibrium which is globally stable. By (48) and comparison principle, there is a such that , , , and for all . Consequently, for , we have Consider an auxiliary system The coefficient matrix of the right hand of (53) is . Since matrix has a positive eigenvalues with a positive eigenvector, it follows from a comparison principle that , , and as , , which leads to a contradiction. This proves (49). Hence . Clearly, every forward orbit in converges to . By Theorem 4.6 in [32] we are able to conclude that model (2) is uniformly persistent with respect to . Thus, by a well-known result in persistence theory in [33] we know that model (2) has at least one positive equilibrium . This completes the proof of Theorem 6.

6. Global Stability of

We further have the following result on the stability of the endemic equilibrium.

Theorem 7. If , then model (2) admits a unique equilibrium , which is globally asymptotically stable.

Proof. In Lemma 3, we have proved that and as , . Therefore, in model (2) we can represent and by and , , respectively, and the model (2) will degenerate into the following system with six equations: By Lemma 3, the dynamics of system (54) can be focused on the following region: We will use the theory of cooperate system to prove the global stability of system (54). Therefore, we only verify the assumption in Corollary 3.2 [34] for system (54). Let and then is a continuously differentiable map. Clearly and for all with , . Since for , we have that is cooperative on . For every and , we have Using the same argument, we can show that is strictly sublinear on . By computing , we have whereClearly, is irreducible for . From (24) we can see that if , then . Since , we have . By Corollary 3.2 in [34], one can conclude that system (54) admits a unique positive equilibrium , which is globally asymptotically stable. According to the theory of asymptotic autonomous systems [31], we further obtain that endemic equilibrium is globally attractive for model (2).

7. Simulations

To complement the mathematical analysis carried out in the previous section, we now investigate some of the numerical properties of the two-patch model (2).

Example 8. Take parameters in model (2) as follows:
, , , , , , , , , , , , , , , , , and . If the two patches are isolated, by simple calculations we have , .

From Lemma 1 we have that the disease will die out in the first patch and will be endemic in the second patch (see Figure 2). From Figure 3 we can easily see that will be larger than 1 under the condition of a larger and a smaller . This means that the larger diffusion of dogs from the lower epidemic areas to the higher prevalence areas can intensify the spread of echinococcosis (see Figure 4). However, when is small and is large, will be smaller than 1. This indicates that the larger diffusion of dogs from the higher prevalence areas to the lower epidemic areas can reduce the spread and is beneficial for disease control (see Figure 5).

fig2
Figure 2: Time series of echinococcosis disease , , when the two patches are isolated for the parameters given in Example 8.
576365.fig.003
Figure 3: Surface plot of as a function of and for the parameters given in Example 8.
fig4
Figure 4: Time series of echinococcosis disease , , when the two patches are connected with , .
fig5
Figure 5: Time series of echinococcosis disease , , when the two patches are connected with , .

Example 9. We use the parameters given in Example 8 except that , , , , and . If the two patches are isolated, by simple calculations we have , .

It follows from Lemma 1 that the disease will die out in both two patches when they are isolated (see Figure 6). However, from Figure 7 we can see that is not always less than 1. This suggests that dogs diffusion can cause the spread of echinococcosis in two patches (see Figure 8).

fig6
Figure 6: Time series of echinococcosis disease , , when the two patches are isolated for the parameters given in Example 9.
576365.fig.007
Figure 7: Surface plot of as a function of and for the parameters given in Example 9.
fig8
Figure 8: Time series of echinococcosis disease , , when the two patches are connected with , .

8. Discussion

In this paper, in order to model the transmission dynamics of echinococcosis spread between two patches due to dogs migration a patch model for echinococcosis is proposed. We define the basic reproduction number . The mathematical results show that the dynamics of the model is completely determined by . If , the disease-free equilibrium is globally asymptotically stable. When , the model is permanence and endemic equilibrium is globally asymptotically stable. According to the simulation we have that the larger diffusion of dogs from a low epidemic area to the high prevalence area can intensify the disease spread. However, the larger diffusion of dogs from the high prevalence area to a low epidemic area can reduce the disease spread and is beneficial to disease control. Additionally, the model presented in this paper can be extended to describe the dynamical transmission of echinococcosis with dogs migration among more than two patches. We leave these in our future work.

Conflict of Interests

The authors declare that they have no financial and personal relationships with other people or organizations that can inappropriately influence their work and there are no professional or other personal interests of any nature or kind in any product, service, and/or company that could be construed as influencing the position presented in, or the review of, this paper.

Acknowledgments

This work is State Key Laboratory Basic Medicine Base of Xinjiang Hydatid Disease Research Fund (XJDX0202-2013-2), Cultivation Project Young Scientific Innovative Talents (2013731017), the National Natural Science Foundation of China (11201399, 11271312), the Scientific Research Programmes of Colleges in Xinjiang (XJEDU2012S20), and Academic Discipline Project of Xinjiang Medical University Health Measurements and Health Economics (XYDXK50780308).

References

  1. S. A. Berger and J. S. Marr, Human Parasitic Diseases Sourcebook, Jones & Bartlett Learning, Sudbury, Mass, USA, 2006.
  2. C. N. Macpherson, “An active intermediate host role for man in the life cycle of Echinococcus granulosus in Turkana, Kenya,” American Journal of Tropical Medicine and Hygiene, vol. 32, no. 2, pp. 397–404, 1983. View at Scopus
  3. W. Zhenghuan, W. Xiaoming, and L. Xiaoqing, “Echinococcosis in China, a review of the epidemiology of Echinococcus spp.,” EcoHealth, vol. 5, no. 2, pp. 115–126, 2008. View at Publisher · View at Google Scholar · View at Scopus
  4. M. Roberts, J. Lawson, and M. Gemmell, “Population dynamics in Echinococcosis and cysticercosis: mathematical model of the life-cycle of Echinococcus granulosus,” Parasitology, vol. 92, no. 3, pp. 621–641, 1986. View at Scopus
  5. P. Torgerson, B. Shaikenov, A. Rysmukhambetova, A. Ussenbayev, A. M. Abdybekova, and K. K. Burtisurnov, “Modelling the transmission dynamics of Echinococcus granulosus in dogs in rural Kazakhstan,” Parasitology, vol. 126, no. 5, pp. 417–424, 2003. View at Publisher · View at Google Scholar · View at Scopus
  6. P. Torgerson, “Canid immunity to Echinococcus spp.: impact on transmission,” Parasite Immunology, vol. 28, no. 7, pp. 295–303, 2006. View at Publisher · View at Google Scholar · View at Scopus
  7. P. Torgerson, “The use of mathematical models to simulate control options for echinococcosis,” Acta Tropica, vol. 85, no. 2, pp. 211–221, 2003. View at Publisher · View at Google Scholar · View at Scopus
  8. D. Heinzmann, A. D. Barbour, and P. R. Torgerson, “A mechanistic individual-based two-host interaction model for the transmission of a parasitic disease,” International Journal of Biomathematics, vol. 4, no. 4, pp. 443–460, 2011. View at Publisher · View at Google Scholar · View at MathSciNet
  9. M. Roberts and M. Aubert, “A model for the control of echinococcus multilocularis in France,” Veterinary Parasitology, vol. 56, no. 1–3, pp. 67–74, 1995. View at Publisher · View at Google Scholar · View at Scopus
  10. M. Vervaeke, S. Davis, H. Leirs, and R. Verhagen, “Implications of increased susceptibility to predation for managing the sylvatic cycle of Echinococcus multilocularis,” Parasitology, vol. 132, no. 6, pp. 893–901, 2006. View at Publisher · View at Google Scholar · View at Scopus
  11. H. Ishikawa, Y. Ohga, and R. Doi, “A model for the transmission of Echinococcus multilocularis in Hokkaido, Japan,” Parasitology Research, vol. 91, no. 6, pp. 444–451, 2003. View at Publisher · View at Google Scholar · View at Scopus
  12. T. Nishina and H. Ishikawa, “A stochastic model of Echinococcus multilocularis transmission in Hokkaido, Japan, focusing on the infection process,” Parasitology Research, vol. 102, no. 3, pp. 465–479, 2008. View at Publisher · View at Google Scholar · View at Scopus
  13. F. Hansen, K. Tackmann, F. Jeltsch, C. Wissel, and H. Thulke, “Controlling echinococcus multilocularis-ecological implications of field trials,” Preventive Veterinary Medicine, vol. 60, no. 1, pp. 91–105, 2003. View at Publisher · View at Google Scholar · View at Scopus
  14. K. Takumi and J. Van der Giessen, “Transmission dynamics of Echinococcus multilocularis; its reproduction number, persistence in an area of low rodent prevalence, and effectiveness of control,” Parasitology, vol. 131, no. 1, pp. 133–140, 2005. View at Publisher · View at Google Scholar · View at Scopus
  15. K. Takumi, A. de Vries, M. L. Chu, J. Mulder, P. Teunis, and J. van der Giessen, “Evidence for an increasing presence of Echinococcus multilocularis in foxes in The Netherlands,” International Journal for Parasitology, vol. 38, no. 5, pp. 571–578, 2008. View at Publisher · View at Google Scholar · View at Scopus
  16. N. Kato, K. Kotani, S. Ueno, and H. Matsuda, “Optimal risk management of human alveolar echinococcosis with vermifuge,” Journal of Theoretical Biology, vol. 267, no. 3, pp. 265–271, 2010. View at Publisher · View at Google Scholar · View at MathSciNet
  17. J. A. M. Atkinson, G. M. Williams, L. Yakob et al., “Synthesising 30 years of mathematical modelling of Echinococcus transmission,” PLoS Neglected Tropical Diseases, vol. 7, no. 8, p. e2386, 2013.
  18. K. Wang, X. Zhang, Z. Jin, H. Ma, Z. Teng, and L. Wang, “Modeling and analysis of the transmission of Echinococcosis with application to Xinjiang Uygur Autonomous Region of China,” Journal of Theoretical Biology, vol. 333, pp. 78–90, 2013. View at Publisher · View at Google Scholar · View at MathSciNet
  19. W. Wang and X.-Q. Zhao, “An epidemic model in a patchy environment,” Mathematical Biosciences, vol. 190, no. 1, pp. 97–112, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  20. Y. Jin and W. Wang, “The effect of population dispersal on the spread of a disease,” Journal of Mathematical Analysis and Applications, vol. 308, no. 1, pp. 343–364, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  21. J. Cui, Y. Takeuchi, and Y. Saito, “Spreading disease with transport-related infection,” Journal of Theoretical Biology, vol. 239, no. 3, pp. 376–390, 2006. View at Publisher · View at Google Scholar · View at MathSciNet
  22. Y. Yang, J. Wu, J. Li, and Z. Ma, “Global dynamics-convergence to equilibria-of epidemic patch models with immigration,” Mathematical and Computer Modelling, vol. 51, no. 5-6, pp. 329–337, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  23. D. Gao and S. Ruan, “An SIS patch model with variable transmission coefficients,” Mathematical Biosciences, vol. 232, no. 2, pp. 110–115, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  24. J. J. Tewa, S. Bowong, and B. Mewoli, “Mathematical analysis of two-patch model for the dynamical transmission of tuberculosis,” Applied Mathematical Modelling, vol. 36, no. 6, pp. 2466–2485, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  25. R. Lintott, R. Norman, and A. Hoyle, “The impact of increased dispersal in response to disease control in patchy environments,” Journal of Theoretical Biology, vol. 323, pp. 57–68, 2013. View at Publisher · View at Google Scholar · View at MathSciNet
  26. Z. Qiu, “Dynamics of an epidemic model with host migration,” Applied Mathematics and Computation, vol. 218, no. 8, pp. 4614–4625, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  27. J. Liu, “Spread of disease in a patchy environment,” Chaos, Solitons and Fractals, vol. 45, no. 7, pp. 942–949, 2012. View at Publisher · View at Google Scholar · View at Scopus
  28. D. Gao and S. Ruan, “A multipatch Malaria model with logistic growth populations,” SIAM Journal on Applied Mathematics, vol. 72, no. 3, pp. 819–841, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  29. O. Diekmann, J. A. P. Heesterbeek, and J. A. J. Metz, “On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations,” Journal of Mathematical Biology, vol. 28, no. 4, pp. 365–382, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  30. P. Van den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  31. H. R. Thieme, “Convergence results and a Poincaré-Bendixson trichotomy for asymptotically autonomous differential equations,” Journal of Mathematical Biology, vol. 30, no. 7, pp. 755–763, 1992. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  32. H. R. Thieme, “Persistence under relaxed point-dissipativity (with application to an endemic model),” SIAM Journal on Mathematical Analysis, vol. 24, no. 2, pp. 407–435, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  33. V. Hutson and K. Schmitt, “Permanence and the dynamics of biological systems,” Mathematical Biosciences, vol. 111, no. 1, pp. 1–71, 1992. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  34. X. Q. Zhao and Z. J. Jing, “Global asymptotic behavior in some cooperative systems of functional-differential equations,” The Canadian Applied Mathematics Quarterly, vol. 4, no. 4, pp. 421–444, 1996. View at Zentralblatt MATH · View at MathSciNet