Table of Contents Author Guidelines Submit a Manuscript
Abstract and Applied Analysis
Volume 2014, Article ID 145398, 13 pages
http://dx.doi.org/10.1155/2014/145398
Research Article

Analysis of a Mathematical Model of Emerging Infectious Disease Leading to Amphibian Decline

Department of Mathematics, Lahore University of Management Sciences, DHA, Lahore Cantt 54792, Pakistan

Received 9 November 2013; Revised 12 February 2014; Accepted 2 March 2014; Published 28 April 2014

Academic Editor: Igor Leite Freire

Copyright © 2014 Muhammad Dur-e-Ahmad 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

We formulate a three-dimensional deterministic model of amphibian larvae population to investigate the cause of extinction due to the infectious disease. The larvae population of the model is subdivided into two classes, exposed and unexposed, depending on their vulnerability to disease. Reproduction ratio has been calculated and we have shown that if , the whole population will be extinct. For the case of , we discussed different scenarios under which an infected population can survive or be eliminated using stability and persistence analysis. Finally, we also used Hopf bifurcation analysis to study the stability of periodic solutions.

1. Introduction

Worldwide catastrophic declines in the amphibian population are perhaps one of the most pressing and discussed problems among the ecologists during the last two decades. Although many of these are attributable due to the habitat loss (see, e.g., [1, 2]), the majority have remained enigmatic till today. Many hypotheses responsible for this decline have been documented in the literature, such as adverse weather patterns [3, 4], acid precipitation [5], environmental pollution [3], increased ultraviolet (UV-B) radiation [6], introduction of predators or competitors [7], infectious disease [3, 8, 9], or a combination of these.

Recently, infectious diseases have become one of the emergent factors behind rapid amphibian decline which often results in extinction of species, for example, the recent extinction of the golden toad in Costa Rica and some species of gastric-brooding frogs in Australia. Investigations reveal that the main cause behind the mass death of these species is two infectious diseases, chytridiomycosis in the rain forest of Australia and Central America and some parts of North America and iridoviral infections in United Kingdom, United States, and Canada [3, 10, 11].

The larval stage of the amphibian is considered to be the most vulnerable towards the spread of infectious disease. Most of the larva population of tropical amphibian species remain alive for 12 to 18 months but some temperate species may also survive as long as 3 years before metamorphosing. Recently it has been observed, both in Australia and Central America, that larval amphibians infected with chytridiomycosis may exhibit disfigurement of their keratinized mouth parts and demonstrate a significant reduction in their growth and development which may eventually cause the complete extinction sometimes [3]. Further, in the case of reduced amphibian population, this infection causes prolonging of the existence of Batrachochytrium [12] and implicates the life-cycle stage as a reservoir host for the pathogen. This kind of larval infection enhancing pathogen-mediated host population extinction has also been reported for invertebrates [13]. Therefore in this work, we keep our focus on the larval stage and investigate the dynamics of spread of disease among them.

This work has been motivated by the recent work of [9], in which only susceptible and infected classes are considered. Since, as discussed earlier, larvae are the key players in the case of density-dependent disease incidence which eventually leads to the host extinction, therefore, it is imperative to keep this component intact while studying the disease based amphibian decline. Since the main focus of the model is the disease based larvae extinction, therefore, in order to look at the complete scenario, we assume that, at the early stage of their lives, the larvae are surrounded by a safe environment which is free of disease, like in a closed shell or some protected holes. For our model, we denote this class by . At the later stage once these unexposed larvae enter into free environment they become vulnerable to disease and we call them susceptibles (denoted by ). These susceptibles can be infected with the disease and enter into the infective class, denoted by . Further, these infected larvae can also recover from the disease and return to the susceptible class. The model is formulated in a way that the number of new infected larvae depends on the present number of susceptible and infective larvae in the sense that the more larvae (susceptible or infective) we have, the more chance of spreading the infections among the healthy or susceptible larvae will there be; that is, the disease can be transferred from one to another.

We organize this paper as follows. First we describe the mathematical model along with its underlying assumptions. Then we introduce many threshold quantities, for example, the reproduction ratios and the critical host density for the disease establishments. We also derive the number of endemic equilibriums and discuss their stabilities. Next we also elaborate our findings using numerical simulations. Finally, we give the conclusion and discuss the potential limitations and impact of our results.

2. Formulation of Mathematical Model

In this mathematical model we assume that the disease only affects the larval population. So we divide the larvae into two categories, the susceptibles larvae and the infected larvae . Although this disease is transferable from one to another, the recovery is also possible and on recovering from the disease, an infected larva will reenter into the class of susceptible larvae. Further, for our model we also assume that only susceptibles can contribute to the reproduction process. Thus once the disease is spread, the whole population will go towards extinctions not only due to the illness but also due to the lack of reproduction process.

The model is as follows: Here represents the unexposed class of the larval population which is not yet entered into the susceptible class . Here we assume that the transition from stage to stage requires a certain minimum size. This assumption is valid for many species, for example, Daphnia magna, a water flea, where the individuals typically have a length of 0.8 mm at birth and a length of 2.5 mm once they enter into the exposed class [14]. Further, in the case of some amphibians, the body size at metamorphosis is quite flexible (see, e.g., [15]); therefore, a certain minimum size still seems to be required in these cases. Thus based on these evidences, it is reasonable to assume that the scarcity of resources will prolong the length of the larval stage . So, if there are more larvae in stage and less resources, then it will take longer for a single larva in stage to complete the transition into stage . Therefore, in our model, we assume that is nondecreasing in in order to represent a strong negative feedback from the number of larval population . It is clear from the model that this class will enter into the exposed susceptible class with the rate which further can become infected with a rate . Since, in the case of closed population, the disease can transfer from one individual to the other, therefore by using the law of mass action, we express this by , where is disease transmission rate. Finally the last term in the third equation represents the recovery from disease with the rate .

The detailed description of parameters of our model is given in Table 1.

tab1
Table 1: Description of parameters in model (1).

Assumption 1. All the parameters are positive except , which can also be zero. is a strictly decreasing nonnegative function of such that as .

3. Existence, Positivity, and Boundedness of the Solution

Theorem 2. Assume that Assumption 1 is satisfied; then for nonnegative initial data, there are unique nonnegative solutions of the system which are defined for all nonnegative times. Further all the solutions to the system are uniformly eventually bounded.

Proof. Let and , where Notice that all the partial derivatives , are continuous. Further So , where for all for which it is defined and whenever . Thus all the solutions of the given system are nonnegative by proposition A.1 of [14]. We conclude that there exists a unique solution to the above system with the values in and it is defined in the interval , (Theorem A.4 [14]) and if then To prove the first part of the theorem, we need to show only that . Assume that ; then from model (1) we get Let ; then we get Since , therefore This is a contradiction, so .
Now we show that the solutions are uniformly eventually bounded. Let us define , where is to be determined later. Then by using the system (1), we have By choosing , we get Since we know that as , therefore for each there exist some such that, for all , . So we can divide the proof into two cases.
Case 1  . From (9) we get So if we take , then we get .
Case 2  . Using , again from (9), we have Simplification and neglecting the negative terms (involving ), yields We find some such that whenever .
Therefore in either case . This implies that . Otherwise since for all sufficiently large time and would become negative in finite time, contradiction. Further since and we have shown that our system is weakly dissipative, so it is dissipative by Proposition 3.18 of [14].

4. Existence and Stability of Steady States

In this section we will discuss the existence and stability of the steady states of our model (1). We will have three steady states, the trivial steady states , the disease-free steady state of the form , and the interior steady state of the form . The trivial steady state always exists. To find the disease-free steady state, we proceed as follows. From the equation of our main model (1) we have From (13) we can say that either or . Let us suppose the disease-free case that is, ; then by first and second equations of model (1) we have Now by adding these two equations we get To find the value of , we substitute the value of from (15) and in equation of model (1). We have Notice that will give us trivial state . In order to find the nontrivial state, we will consider the case when . In this case we will have Equivalently Thus we have two disease-free steady states, that is, and , where is given by (17). Observe that there is a threshold condition We can rewrite the condition as . We refer this quantity as reproduction number and it is defined as Now we find the interior equilibrium with . From the equation of model (1), we get .

Adding the first two equations of our original model (1) gives Now by using the values of , we get Therefore we have Now we need to find the value of . By substituting the values of and in the first equation of model (1), we will have

Now we need to solve (24) numerically (depending on the structure of the function ) to find the value of . Therefore, in this case, the interior equilibrium is given by The and are nonzero. From the equation of model (1) we get If we solve this for , we get This implies that if This gives the condition for the existence of infected steady state. Observe that we can write (24) as Since is the average time a larva spends in the infected stage and is the rate at which one average infected larva infects a susceptible larva, if they are at density , therefore we can express the left side of (29) as the replacement ratio of the infectious disease from the infected larvae to the susceptible larvae when they have density . That is, The replacement ratio at the susceptible larval density is given by

From the above results we have the following theorem for the existence of the steady states.

Theorem 3. (a) The trivial steady state exists all the time.
(b) The disease-free state exists if and only if .
(c) The infected state exists if and only if and .
(d) Let . Then there is at least one steady state with the infected larvae.
(e) Let ; then there is no or two steady states with the infective larvae.

Proof. Proof of part (b) and part (c) follows directly from the discussion above this theorem.
(d) The condition can be written as The intermediate value theorem implies that there is some which satisfies (24).
(e) The condition can be written as If is strictly increasing function for , then there is no steady state in . Now assume that this function is not increasing and that there is at least one solution of (24). We can choose . If then we have two solutions counting multiplicities. They are steady state with the infectives if ; otherwise that steady state does not exist with infectives. If , then we can choose some . This implies Therefore by the intermediate value theorem, there exists a solution of (20). By the previous theorem, and are components of equilibria with positive infective population.

Theorem 4. Let be a solution of (24); then satisfy (29) if and only if with .

Proof. Let . Since is strictly decreasing function, this is equivalent to Since , so we have the equivalent expression Multiplying both sides with and rearranging the terms imply At infected steady state the equation of model (1) can be written as So we get Similarly if then (31) will not be satisfied.

It is also clear that (24) has one or three solutions depending on the choice of the parameters , , , , and .

Now we will discuss the stability of the trivial steady state and the persistence of the population. The persistence of a population is defined as follows.

Definition 5. The population is called uniformly strongly persistent if there exists , independent of initial data, such that for all solutions of model (1) satisfying , , and we have for all sufficiently large . It is robust uniform persistence, if model (1) depends continuously on a parameter . If is a parameter and uniform persistence holds with replaced by , where is semiflow, for all in the neighborhood of  , then we say that the system (1) is robust uniformly persistent.

Theorem 6. (a) The trivial state is locally asymptotically stable if and only if .
(b) Assume that and (or equivalently ). Then the steady state is locally asymptotically stable if and only if   or .
(c) Assume that and (or equivalently ). Then the endemic equilibrium is locally asymptotically stable.
(d) If , then the disease is robustly uniformly persistent: that is, for every parameter vector for the system (1), there exists a neighborhood of and such that for all solutions of model (1) corresponding to with .

Proof. The Jacobian matrix associated with our model (1) is given by
(a) This Jacobian matrix at steady state is
The three eigenvalues are , and the eigenvalues of the matrix The trace and determinant are given by So the steady state is locally asymptotically stable if and only if .
(b) The Jacobian matrix at is One of the eigenvalues is , which is negative if . The other two eigenvalues come from Note that where .
Notice that is positive if .
(c) The endemic equilibrium is given by , where is the solution of .
The Jacobian matrix at this steady state is given by Since and , so we get where .
Now we have Further we have where is the determinant of the matrix of size 2 resulting from by removing the row and the column: Finally Simplifying and rearranging the terms we get Since , so this implies by the assumption of the structure of (i.e., the case when we have three interior equilibria and this is the most right one). Thus the endemic equilibrium is locally asymptotically stable.
(d) Let be a fixed parameter and . Define , where is the same region defined in Theorem 2 in which all the solutions are uniformly eventually bounded. Notice that both and are positively invariant and that attracts all the solutions in . Considering as a periodic orbit (say, e.g., of period ), then our results will follow by applying the theorems discussed in [16] (see, in particular, Corollary (4.7), Proposition (4.1), and Theorem (3.2) of [16]). Now if we denote the last two equations of model (1) in the form of a matrix , then it is clear that has a positive spectral radius provided . So condition (9) of Corollary 4.7 in [17] is satisfied. Also, is irreducible which, using Theorem 1.1 in [18], implies that has all positive entries for all . This implies that condition (1) of the abovementioned corollary is also satisfied.

Theorem 7. , and converge to zero as if and only if .

Proof. We will use “La Salle’s Invariance Principle” to prove the theorem.
Case 1  . Consider the Lyapunov function . Since , so is a continuously differentiable function with . Now consider Now define Further, as we have Since and are nonnegative and also , so we have . But from the condition that should be nonnegative we get . So we have As and are positive, therefore we must have .
So is the only subset of which is also one of the steady states and so it is the largest invariant subset of . Since we have shown above that all the solutions of the above system are bounded, thus, by the above result, all the solutions to the given system converge to .
Case 2  . Again consider . In this case we have Define Since we have Thus from the nonnegativity of , we conclude that and so we get and since, at , we get from the first equation of our main model that .
Thus is the only invariant subset of ; therefore all the bounded solutions of the given system will converge to .
Case 3  . Consider the Lyapunov function . Then Define Since
so we have . Therefore, from the nonnegativity of , we conclude that . So we will have either or .
Notice that if , then , so from the strictly decreasing property of , , and therefore Thus and the only possibility is and hence is the only subset of . Again from the first equation of our model (1), , which implies if and only if . Hence is the only invariant subset of ; therefore all the bounded solutions of the given system will converge to .

5. The Disease Transmission Model Special Case

In this section we will assume that with , . Also assume that for . is strictly decreasing as long as it is strictly positive. Let us define , , and . This implies Now by using these substitutions, our main model will become On simplifying these equations, we get First, we will make this system dimensionless. Assume that the , , , , , and are measured in . , and are measured in the units of population size, say . So is measured in . Also is measured in . Now dividing by both sides our system becomes where and vice versa. Rewriting the above system after renaming the parameters, we will have

Here we measured time in and renamed the parameters , and as , and , as , and , and as , and . The above system (70) is dimensionless. Note that for the system .

Again we will have three types of steady states, that is, , and .

Theorem 8. (a) The trivial state always exists.
(b) The disease-free steady state of the form exists and is unique if and only if .
(c) The infected state exists if and only if and .
(d) If , then for the system (70) there exist one or more than one interior steady state.
If then for the system (70) there exists only one interior steady state with the component given by .

Proof. (b) Since we know that the “” component of disease-free steady state comes from this implies which gives and the “” component is given by This steady state will exist if and only if
(c) It directly follows from part (c) of Theorem 3.
(d) From the third equation of the system (70), we have As calculated earlier, the “” component is given by It is clear that the feasibility condition for the existence of that steady state is where as proved in Theorem 3. This is equivalent to Observe that this holds only if , since for , we get , which is not true. So we must have . In this case we get We can write the left-hand side of (80) as where Notice that and both of them will exist only if . Further from (81) it is clear that and have the opposite signs. So we have the following two cases.
Case 1  ( and). This means and or equivalently .
Case 2   ( and ). This means and . Since , so this case is not possible. So we can have only one choice for ; that is, (with and given in (82)) for the existence of . The “” component of the endemic steady state is the solution of the equation By using the value of , we get This implies that for which is a quadratic equation in and the roots of this equation are given by It is clear that . Also will exist only if . That is, we have Simplifying and rearranging the terms yield We have the following two cases.
Case i  . For this case inequality (89) is automatically satisfied.
Case ii   . This is equivalent to . Then on squaring both sides of inequality (89), we get This implies that Now consider ; this will also exist only if ; from (89), we get which can be rewritten as From the above equation notice that will exist only if . Taking square of both sides of this equation we get which implies that .

Theorem 9. (i) The steady states of the form are stable if and only if and .
(ii) The endemic state is locally stable and is unstable whenever they exist.

Proof. (a) The proof directly follows from Theorems 3 (b) and 6 (b).
(b) The Jacobian matrix is given by At the steady state we get Assume that , so we have the trace and determinant are given by Observe that This implies Similarly Now we evaluate the determinant for and . At , and also .
At , we have , so at least one of the eigenvalues is positive and thus this equilibrium is unstable.
As calculated earlier This implies Since the only point of interest here is (because is unstable) and , therefore we have Now for the local asymptotic stability, we want simplification yields It can be easily shown that we can find such a which satisfies (106). Further this is clear from (101) that

Special Case. Here we will discuss the special case when there is no disease, that is, when . Also assume that for . Under this condition our system given by (70) will be reduced to a two-dimensional system and is given by Observe that is a steady state for this system. We can write the steady state form of this system into matrix form as for any nontrivial steady state, ; the determinant of this matrix should be zero; that is, On rearranging the terms, we get This implies

This steady state will exist if . Now the Jacobian matrix at is given by The trace and determinant are given by Notice that, from (111), . So we get So this internal steady state is locally asymptotically stable if , that is, when , and unstable if , that is, when . Thus we have that is the point of bifurcation.

Theorem 10. Assume that the positive parameters and can be chosen such that ; then becomes the first coordinate of a unique nontrivial steady state that is unstable spiral point. Every solution that does not start at the nontrivial steady state or the origin converges towards a periodic orbit. In particular there exists an orbitally stable periodic orbit.

Proof. Since we have shown earlier that all the solutions are bounded and we have only one interior steady state for this limiting system which is also nondegenerate, that is, all the eigenvalues of the associated Jacobian matrix evaluated at this steady state point are nonzero, so by Theorem A-15 of [19], there exists a locally asymptotically stable periodic orbit.

So far we have shown that, with proper choice of the values of and , we can find which becomes the first coordinate of an unstable spiral point or in particular there exists a periodic solution. Here we will use the approach of the theory of Hopf bifurcation and in particular the concepts of supercritical and subcritical bifurcation to discuss the stability of these periodic orbits. Actually we need to find a Jacobian matrix of the form with .

Again consider the system equation (86). Define . This implies This implies Similarly This implies Let . So we have or Observe that the component of the interior steady state of this system satisfies . The Jacobian matrix is given by Evaluating the Jacobian matrix at the interior steady state we get Since we also know that at (i.e., Hopf bifurcation point) we have , so at this point the Jacobian matrix will become Now we need to choose in such a way that which gives So with this choice of , we have Hopf bifurcation. Now assume that is the vector field associated with (122); that is, The stability of the bifurcating periodic orbit is determined by the sign of the following number: Now we calculate these partial derivatives: So by substituting these values we get Since , so the bifurcating periodic orbits are asymptotically stable; that is, we have the case of supercritical bifurcation.

6. Discussion

In this paper, we have introduced and analyzed a model of disease based amphibian decline. We found basic reproduction number , which guarantees the extinction of disease population in the presence of disease as shown in Figure 1. We proved that only the disease-free steady state is stable and the stability of the other two steady states is conditional and depends on some other environmental factors whenever . We also have shown that not only is the sufficient condition for the survival of infected population but also for the case when disease may or may not persist depending on certain constraints (see Figures 1, 2, and 3). In our model we also have considered a small population size of . The large population case will be relatively easy to handle in terms of the function . In this case, since is decreasing, therefore for large population size in (24) and this will take the form

145398.fig.001
Figure 1: This figure shows the time series of model (1) when . It is clear that all the population will extinct in this case. The parameters used are .
145398.fig.002
Figure 2: The case of . It is clear that, under certain parameter values, only larvae and susceptible population will survive while the infected population will be extinct. All parameter values are the same except .
145398.fig.003
Figure 3: Phase portrait of versus . The existence of a periodic solution for model (1) is evident. Parameters values are the same as of Figure 2.

This is the case when we have large number of susceptibles.

Further, in our model, we just consider the simple case when only the larvae in stage are subject to intrastage competition. A future direction for this model is to consider a rather more complicated case when both secure larvae and exposed larvae will compete for the resources. In this case, the term will be replaced by , which may cause some oscillatory effects due to the density dependent development rate for both and .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

  1. R. M. Lehtinen, S. M. Galatowitsch, and J. R. Tester, “Consequences of habitat loss and fragmentation for wetland amphibian assemblages,” Wetlands, vol. 19, no. 1, pp. 1–12, 1999. View at Google Scholar · View at Scopus
  2. S. N. Stuart, J. S. Chanson, N. A. Cox et al., “Status and trends of amphibian declines and extinctions worldwide,” Science, vol. 306, no. 5702, pp. 1783–1786, 2004. View at Publisher · View at Google Scholar · View at Scopus
  3. P. Daszak, L. Berger, A. A. Cunningham, A. D. Hyatt, D. Earl Green, and R. Speare, “Emerging infectious diseases and amphibian population declines,” Emerging Infectious Diseases, vol. 5, no. 6, pp. 735–748, 1999. View at Google Scholar · View at Scopus
  4. S. A. Diamond, G. S. Peterson, J. E. Tietge, and G. T. Ankley, “Assessment of the risk of solar ultraviolet radiation to amphibians. III. Prediction of impacts in selected northern midwestern wetlands,” Environmental Science and Technology, vol. 36, no. 13, pp. 2866–2874, 2002. View at Publisher · View at Google Scholar · View at Scopus
  5. B. A. Pierce, “The effect of acid precipitation on amphibians,” Ecotoxicology, vol. 2, no. 1, pp. 65–77, 1993. View at Google Scholar · View at Scopus
  6. L. E. Licht and K. P. Grant, “The effects of ultraviolet radiation on the biology of amphibians,” American Zoologist, vol. 37, no. 2, pp. 137–145, 1997. View at Google Scholar · View at Scopus
  7. S. J. Hecnar and R. T. M'Closkey, “The effects of predatory fish on amphibian species richness and distribution,” Biological Conservation, vol. 79, no. 2-3, pp. 123–131, 1997. View at Publisher · View at Google Scholar · View at Scopus
  8. K. R. Lips, F. Brem, R. Brenes et al., “Emerging infectious disease and the loss of biodiversity in a Neotropical amphibian community,” Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 9, pp. 3165–3170, 2006. View at Publisher · View at Google Scholar · View at Scopus
  9. H. R. Thieme, T. Dhirasakdanon, Z. Han, and R. Trevino, “Species decline and extinction: synergy of infectious disease and Allee effect?” Journal of Biological Dynamics, vol. 3, no. 2-3, pp. 305–323, 2009. View at Publisher · View at Google Scholar · View at MathSciNet
  10. L. Berger, R. Speare, P. Daszak et al., “Chytridiomycosis causes amphibian mortality associated with population declines in the rain forests of Australia and Central America,” Proceedings of the National Academy of Sciences of the United States of America, vol. 95, no. 15, pp. 9031–9036, 1998. View at Publisher · View at Google Scholar · View at Scopus
  11. S. E. Drury, R. E. Gough, and A. A. Cunningham, “Isolation of an iridovirus-like agent from common frogs (Rana temporaria),” The Veterinary Record, vol. 137, no. 3, pp. 72–73, 1995. View at Google Scholar · View at Scopus
  12. L. J. Rachowicz, R. A. Knapp, J. A. T. Morgan et al., “Emerging infectious disease as a proximate cause of amphibian mass mortality,” Ecology, vol. 87, no. 7, pp. 1671–1683, 2006. View at Google Scholar · View at Scopus
  13. A. Richards, J. Cory, M. Speight, and T. Williams, “Foraging in a pathogen reservoir can lead to local host population extinction: a case study of a Lepidoptera-virus interaction,” Oecologia, vol. 118, no. 1, pp. 29–38, 1999. View at Publisher · View at Google Scholar · View at Scopus
  14. H. R. Thieme, Mathematics in Population Biology, Princeton Series in Theoretical and Computational Biology, Princeton University Press, Princeton, NJ, USA, 2003. View at MathSciNet
  15. J. P. Collins, “Intrapopulation variation in the body size at metamorphosis and timing of metamorphosis in the bullfrog, Rana castesbeiana,” Ecology, vol. 60, pp. 738–749, 1979. View at Google Scholar
  16. P. L. Salceanu and H. L. Smith, “Persistence in a discrete-time stage-structured fungal disease model,” Journal of Biological Dynamics, vol. 3, no. 2-3, pp. 271–285, 2009. View at Publisher · View at Google Scholar · View at MathSciNet
  17. P. L. Salceanu, “Robust uniform persistence in discrete and continuous dynamical systems using Lyapunov exponents,” Mathematical Biosciences and Engineering, vol. 8, no. 3, pp. 807–825, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  18. H. L. Smith, Monotone Dynamical Systems, An Introduction to the Theory of Competitive and Cooperative Systems, vol. 41 of Mathematical Surveys and Monographs, American Mathematical Society, Providence, RI, USA, 1995. View at MathSciNet
  19. F. Brauer and J. A. Nohel, The Qualitative Theory of Ordinary Differential Equations, Dover Books on Mathematics, New York, NY, USA, 1989.