- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
International Journal of Ecology
Volume 2014 (2014), Article ID 429086, 12 pages
Stabilizing Effect of Prey Refuge and Predator’s Interference on the Dynamics of Prey with Delayed Growth and Generalist Predator with Delayed Gestation
Ecological Modelling Laboratory, Department of Zoology, Visva-Bharati University, Santiniketan 731 235, India
Received 7 February 2014; Accepted 11 April 2014; Published 30 April 2014
Academic Editor: Daniel I. Rubenstein
Copyright © 2014 Debaldev Jana. 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.
In the present paper, I study a prey-predator model with multiple time delays where the predator population is regarded as generalist. For this regard, I consider a Holling-Tanner prey-predator system where a constant time delay is incorporated in the logistic growth of the prey to represent a delayed density dependent feedback mechanism and the second time delay is considered to account for the length of the gestation period of the predator. Predator’s interference in predator-prey relationship provides better descriptions of predator's feeding over a range of prey-predator abundances, so the predator's functional response here is considered to be Type II ratio-dependent. In accordance with previous studies, it is observed that delay destabilizes the system, in general, and stability loss occurs via Hopf bifurcation. There exist critical values of delay parameters below which the coexistence equilibrium is stable and above which it is unstable. Hopf bifurcation occurs when delay parameters cross their critical values. When delay parameters are large enough than their critical values, the system exhibits chaotic behavior and this abnormal behavior may be controlled by refuge. Numerical computation is also performed to validate different theoretical results. Lyapunov exponent, recurrence plot, and power spectral density confirm the chaotic dynamical behaviors.
In ecology, predation describes a biological interaction where a predator (an organism that is hunting) feeds on its prey (the organism that is attacked). Predators may or may not kill their prey prior to feeding on them, but the act of predation often results in the death of its prey and the eventual absorption of the prey’s tissue through consumption. Predators can have profound impacts on the dynamics of their prey that depend on how predator consumption is affected by prey density (the predator’s functional response). A generalist species is able to thrive in a wide variety of environmental conditions and can make use of a variety of different resources (e.g., a heterotroph with a varied diet). Omnivores are usually generalists. Herbivores are often specialists, but those that eat a variety of plants may be considered generalists. Consumption by a generalist predator is expected to depend on the densities of all its major prey species (its multispecies functional response, or MSFR), but most studies of generalists have focused on their functional response to only one prey species [1–4]. Recently, many authors have explored the dynamics of a class of the semi-ratio-dependent [5–8]/Holling-Tanner [9–12]/Leslie-Grower [13, 14] predator-prey systems with functional responses as follows: where and are, respectively, the prey and predator densities at time . The prey population is assumed to grow logistically to its carrying capacity with intrinsic growth rate in absence of predator. is the maximal per capita prey consumption rate, and is the amount of prey necessary for the relative biomass growth rate of the predator to be half its maximum; the predator consumes the prey according to the functional response and grows logistically with growth rate and carrying capacity proportional to the population size of the prey (or prey abundance). is the measure of the food quality that the prey provides, which is converted to predator birth. For more background and biological adjustments of system (1), one can see [7, 15, 16] and the references cited therein.
Predator’s functional response is considered as Type II ratio-dependent [17–25] because a ratio-dependent predator-prey model does not show the so called paradox of enrichment [17, 26, 27] and biological control paradox . The simplification of the food web structure due to predator interference allows qualitative predictions concerning the response of a food web to an external perturbation [25, 29–32]. Ginzburg and Akcakaya  demonstrated a positive relationship between the abundances of all trophic levels and concluded that the ratio-dependent models provide a better representation of predator-prey interactions than prey dependent models. Predator’s functional response, defined as the amount of prey catch per predator per unit of time, is affected by the structure of prey habitat and predator’s hunting ability [34, 35]. Anderson  and Johnson  showed that the dynamics of a local population largely depend upon attributes of the local habitats. In addition, a spatial refuge protects a constant proportion of prey from predation. Mite predator-prey interactions often exhibit spatial refugia which afford the prey some degree of protection from predation and reduce the chance of extinction due to predation. Maynard Smith  shows that the presence of a constant proportion refuge does not alter the dynamical stability of the neutrally stable Lotka-Volterra model, while a constant number refuge of any size replaces the neutrally stable behaviour with a stable equilibrium. Hassel  shows that adding a large refuge to a model, which in the absence of a refuge exhibits divergent oscillations, replaces the oscillatory behaviour with a stable equilibrium. These mathematical models and a number of experiments indicate that refugia have a stabilizing effect on predator-prey interactions, but, as Taylor  has mentioned, it would be an over simplification to assume this is always the case. Kar  and Jana  proposed and studied a predator-prey system where prey population is subjected to refuge from their predator population. These papers lead a model by incorporating a refuge protecting of the prey, where is the density of prey population and is constant measures of the degree or strength of prey refuge. This leaves of the prey available to the predator, so the Holling-Tanner predator-prey model with ratio-dependent functional response with prey refuge is given by
Delay is frequently used in a predator-prey model to make the model biologically more realistic. Recently, researchers are using more than one delay to study the effect of past history of the system populations [43–47]. In the second phase of this study, I consider two delays in the model system (2). One discrete delay is considered in the specific growth rate of prey to incorporate the effect of density dependence feedback mechanism which takes units of time to respond to changes in the prey population . The second delay is incorporated in the negative feedback of the predator density [49–52]. I thus obtain the following multidelayed predator-prey model: The model system (3) has to be studied with the following initial conditions: The organization of the paper is as follows. Section 2 deals with the analysis of the model system (3). Numerical studies are given in Section 3, and, finally, a brief discussion is presented in Section 4.
2. Mathematical Results of System (3)
The system (3) has only one interior equilibrium point given by , where and . The equilibrium point is biologically feasible if . Linearizing the system (3) at , I get where The corresponding characteristic equation is given by where I now discuss the following cases.
Theorem 1. The interior equilibrium of the system (3) exists and becomes locally asymptotically stable in absence of delays if and hold.
Case 2 (). In this case, the characteristic equation (7) becomes
Let be a root of (11). Then I have
This leads to
It follows that (12) has no positive roots if the following conditions are satisfied:
Hence, all roots of (13) will have negative real parts when if conditions of Theorem 1 and (H2) are satisfied.
(H3) If Theorem 1 and (H3) hold, then (13) has a unique positive root . Substituting into (12), I have where is the unique positive root of (13).
(H4) If (H1) and (H4) hold, then (13) has two positive roots and . Substituting into (12), we obtain If is a root of (11) satisfying , resp.) and (or , resp.), I obtain by hypothesis (H4).
Similarly, I can show that
Theorem 2. Assume that and conditions of Theorem 1 are satisfied, then the equilibrium is locally asymptotically stable for and unstable for . Furthermore, the system (3) undergoes a Hopf bifurcation at when , where
Case 3 (). In this case, I allow a gestation time period for the predator and also a constant time delay affecting the density dependent growth rate of the prey. I fix at some value from its stability range and regard as a free parameter. I also assume that the model parameters are such, that hypothesis (H1) holds. Let with be a root of (7). Then,
Note that is always negative. I define
Then it is easy to verify that and , as . Thus, (22) has finitely many positive roots, denoted by . For every fixed , where , there exists a sequence , where
such that (22) holds. Let . When , (7) has a pair of pure imaginary roots for . In the following, I assume that
(H5) Therefore, by the general Hopf bifurcation theorem of functional differential equations, see, for example, , I obtain the following result considering the change of stability of system (3).
Theorem 3. Suppose that the parameters in system (3) are such that hypotheses and ( hold true, and . Then the coexistence equilibrium is locally asymptotically stable when , and it is unstable when . Moreover, Hopf bifurcation occurs when .
Theorem 4. Assume that and the conditions in hypothesis hold true. Then the equilibrium is locally asymptotically stable for and unstable for . Furthermore, the system (3) undergoes Hopf bifurcation when , where and is the unique positive root of the polynomial as follows:
Case 5 ( and is within its stability range). This choice of parameter regime corresponds to the biologically interesting case, when the predator individuals have a fixed gestation period and the growth of the prey population is affected by delayed density dependent mechanisms. The proof of the main result follows similar lines to that in Case 3 and I only summarize the stability results in the following theorem.
Theorem 5. Assume that the model parameters are such that the hypotheses in hold true, and . Then the coexistence equilibrium is asymptotically stable when , and it is unstable when . Hopf bifurcation occurs when , where and with
3. Numerical Computations
I perform numerical computations to observe various dynamics of the coexistence equilibrium for both the systems (2) and (3). I consider the fixed parameter values as , and . With these parameter values . Here, and , and initial value is considered as for each simulation. First I consider the nondelayed system (2). In this case, the system is asymptotically stable (Figure 1).
I now study the dynamic behavior of the delay-induced system (3). When (Case 2), one can compute from Theorem 2 that , , and . Therefore, the coexisting equilibrium is asymptotically stable for (Figure 2(b)) and unstable for (Figure 2(c)). When , the system (3) undergoes a Hopf bifurcation at . System behaviors in Case 2 can be demonstrated more prominently if I plot the bifurcation diagram in the three-dimensional space . Figure 2(a) shows that, when , the coexistence equilibrium is stable for but the instability sets in when .
For Case 3, I take any value of from its stability range , say that , and consider as a parameter. One can compute the value of as and the corresponding critical value of as . Thus, for fixed stable value of , the system (3) exhibits stable behavior around for and unstable oscillatory behavior for (Figure 3). A Hopf bifurcation occurs when takes the critical value .
In Case 4, when , one can compute from Theorem 4, and . One can also verify that and . Therefore, the coexisting equilibrium is asymptotically stable for , unstable for , and when , the system (3) undergoes a Hopf bifurcation at . These results can be observed from the bifurcation diagram Figure 4.
In the last case (Case 5), I take any value of from its stability range , say that , and consider as a parameter. As before, I calculate and the corresponding critical value . The bifurcation diagram Figure 5 demonstrates that the system (2) is locally asymptotically stable around when and unstable when . These results can be observed from the bifurcation diagram Figure 5.
3.1. Chaotic Dynamics
I have analytically studied the local behavior of the system (3) when the delay parameters are within or slightly above their critical values. One interesting topic in the delay-induced system is to study the dynamical behavior of the system when the delay parameters are far away from their critical values, or they assume large values. To observe the system dynamics, I consider and , both are beyond their stability range, and then the system shows chaotic behavior (Figure 6).
In order to characterize the irregular behavior, I perform the standard numerical diagnostics, for example, Lyapunov exponent, sensitivity of the solutions to initial condition, power spectral density, and recurrence plot. Figure 7(a) shows the Lyapunov exponent of with respect to the delay in a smaller range [0.45, 0.539] for fixed . Positive Lyapunov exponent indicates that the system exhibits chaotic behavior [42, 54, 55]. Sensitivity of the system trajectories to the initial values is shown in Figure 7(b). To show this, I compute the error , where and . This phenomena ensures the chaotic nature of the system [42, 54, 56]. The power spectra of the predator population are presented in Figure 7(c). The irregular broad peaks of this figure are indicative of chaos and randomness [42, 54, 57]. The recurrence plot of the system for the parameter values as in Figure 6 is represented in Figure 7(d). The random points on the time-time plane ensure that more or less the same values of phase trajectories take place without any rhythm, indicating the chaotic nature of the system [42, 54, 58].
3.2. Biological Control of Chaotic Dynamics
Many ecological situations show that their unstable or abnormal behavior may be controlled by any external perturbation [14, 21] or habitat structure [34, 36–40]. Recently Jana  shows that chaotic behavior of a discrete predator-prey system at low level of prey refuge must be controlled to its stable coexistence by increasing the degree of prey refuge. Parameter set of Figure 6 depicts the chaotic behavior of the system (3). Standing at this situation, if I increase the intensity of degree of prey refuge (from to ), then again the system (3) settles down to its stable coexistence (Figure 8).
Effect of prey refuge has an extended impact on the predator-prey models. However, both field and laboratory experiments confirm that intensity of prey refuge reduces predation rates by decreasing encounter rates between predator and prey. On the other hand, a predator-prey model becomes more realistic in presence of different delays which are unavoidable elements in physiological and ecological processes. In this paper, I have studied a multidelayed predator-prey model where the prey species is subject to partial refuge from predator population. Also I have considered here that the predator is a generalist type. A time delay is considered in the growth rate of the prey to represent density dependent feedback mechanism, and a second delay is introduced to account for the gestation time of the predator. The objective was to study the role of the effects of the interaction of multiple constant time delays and the prey refuge on the dynamical behavior of a predator-prey system. To have a better understanding of the complexities in natural systems, I can construct larger systems of differential equations containing more parameters. A second approach that is gaining prominence is the inclusion of time delay terms in the differential equations. Complex dynamical behaviour arises as a consequence of time delay in a biological system (with significant time delay) which may exhibit limited cycle oscillation and chaos. In case of the multidelayed system, I obtained sufficient conditions in terms of the system parameters for the stability of the coexistence equilibrium. I observed that the coexistence equilibrium is locally asymptotically stable if the delay parameters are lower than some critical values and it is unstable when the parameter values exceed these critical values. One important question in the multidelayed predator-prey system is to study the behavior of the system when the delay parameter is significantly large and hence it is far away from its critical value. As the estimated length of delays to preserve stability and the critical length of time delays for Hopf bifurcation are dependent on the system’s parameters, it is possible to impose some control, which will prevent the possible abnormal oscillation in the population density. To explore the behavior of the system in this case, I performed extensive numerical simulations. My simulation results show that a refuged predator-prey system with multiple delays may exhibit different interesting (e.g., chaotic) behavior when the delay is large enough. Finally I showed that the chaotic characteristic of the multidelayed system was controlled by increasing the intensity of the prey refuge.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
This research is supported by UGC (Dr. D. S. Kothari Postdoctoral Fellowship), India; no. F.4-2/2006(BSR)/13-1004/2013(BSR).
- J. Matthiopoulos, K. Graham, S. Smout et al., “Sensitivity to assumptions in models of generalist predation on a cyclic prey,” Ecology, vol. 88, no. 10, pp. 2576–2586, 2007.
- T. M. Temesgen, “Bifurcation analysis on the dynamics of a genralist predator-prey system,” International Journal of Ecosystem, vol. 2, no. 3, pp. 38–43, 2012.
- Y. Kang and L. Wedekin, “Dynamics of a intraguild predation model with generalist or specialist predator,” Journal of Mathematical Biology, vol. 67, no. 5, pp. 1227–1259, 2013.
- R. K. Upadhyay and S. R. K. Iyengar, Introduction to Mathematical Modelling and Chaotic Dynamics, Taylor and Francis, 2013.
- Q. Wang, M. Fan, and K. Wang, “Dynamics of a class of nonautonomous semi-ratio-dependent predator-prey systems with functional responses,” Journal of Mathematical Analysis and Applications, vol. 278, no. 2, pp. 443–471, 2003.
- X.-X. Liu, “A note on periodic solutions for semi-ratio-dependent predator-prey systems,” Applied Mathematics, vol. 25, no. 1, pp. 1–8, 2010.
- L. Hongying and W. Weiguo, “Dynamics of a delayed discrete semiratiodependent predator-prey system with Holling type IV functional response,” Advances in Difference Equation, vol. 7, pp. 2–19, 2011.
- C. Lu and L. Zhang, “Permanence and global attractivity of a discrete semi-ratio dependent predator-prey system with Holling II type functional response,” Journal of Applied Mathematics and Computing, vol. 33, no. 1-2, pp. 125–135, 2010.
- R. M. May, Stability and Complexity in Model Ecosystems, Princeton University Press, Princeton, NJ, USA, 1974.
- J. D. Murray, Mathematical Biology, Springer, Berlin, Germany, 1989.
- P. A. Braza, “The bifurcation structure of the Holling-Tanner model for predator-prey interactions using two-timing,” SIAM Journal on Applied Mathematics, vol. 63, no. 3, pp. 889–904, 2003.
- S.-B. Hsu and T.-W. Hwang, “Hopf bifurcation analysis for a predator-prey system of Holling and Leslie type,” Taiwanese Journal of Mathematics, vol. 3, no. 1, pp. 35–53, 1999.
- H. C. Tsai and C. P. Ho, “Global stability for the leslie-gower predator- prey system with time-delay and holling’s type functional response,” Tunghai Science, vol. 6, pp. 43–72, 2004.
- N. Zhang, F. Chen, Q. Su, and T. Wu, “Dynamic behaviors of a harvesting Leslie-Gower predator-prey model,” Discrete Dynamics in Nature and Society, vol. 2011, Article ID 473949, 14 pages, 2011.
- J. T. Tanner, “The stability and the intrinsic growth rates of prey and predator populations,” Ecology, vol. 56, pp. 855–867, 1975.
- S. B. Hsu and T. W. Hwang, “Uniqueness of limit cycles for a predator-prey system of Holling and Lesile type, Canad,” Applied Mathematics Quarterly, vol. 6, pp. 91–117, 1998.
- Y. Kuang and E. Beretta, “Global qualitative analysis of a ratio-dependent predator-prey system,” Journal of Mathematical Biology, vol. 36, no. 4, pp. 389–406, 1998.
- A. Martin and S. Ruan, “Predator-prey models with delay and prey harvesting,” Journal of Mathematical Biology, vol. 43, no. 3, pp. 247–267, 2001.
- C. Jost, O. Arino, and R. Arditi, “About deterministic extinction in ratio-dependent predator-prey models,” Bulletin of Mathematical Biology, vol. 61, no. 1, pp. 19–32, 1999.
- D. Xiao and S. Ruan, “Global dynamics of a ratio-dependent predator-prey system,” Journal of Mathematical Biology, vol. 43, no. 3, pp. 268–290, 2001.
- D. Jana, S. Chakraborty, and N. Bairagi, “Stability, nonlinear oscillations and bifurcation in a delay-induced predator-prey system with harvesting,” Engineering Letters, vol. 20, no. 3, pp. 238–246, 2012.
- D. Xiaq and L. S. Jennings, “Bifurcations of a ratio-dependent predator-prey system with constant rate harvesting,” SIAM Journal on Applied Mathematics, vol. 65, no. 3, pp. 737–753, 2005.
- D. Xiao, W. Li, and M. Han, “Dynamics in a ratio-dependent predator-prey model with predator harvesting,” Journal of Mathematical Analysis and Applications, vol. 324, no. 1, pp. 14–29, 2006.
- S. Chakraborty, S. Pal, and N. Bairagi, “Predator-prey interaction with harvesting: mathematical study with biological ramifications,” Applied Mathematical Modelling, vol. 36, no. 9, pp. 4044–4059, 2012.
- H. R. Akcakaya, R. Arditi, and L. R. Ginzburg, “Ratio-dependent predation: an abstraction that works,” Ecology, vol. 76, no. 3, pp. 995–1004, 1995.
- N. G. Hairston, F. E. Smith, and L. B. Slobodkin, “Community structure, population control and competition,” American Naturalist, vol. 94, pp. 421–425, 1960.
- M. L. Rosenzweig, “Paradox of enrichment: destabilization of exploitation ecosystems in ecological time,” Science, vol. 171, no. 3969, pp. 385–387, 1971.
- R. F. Luck, “Evaluation of natural enemies for biological control: a behavioral approach,” Trends in Ecology and Evolution, vol. 5, no. 6, pp. 196–199, 1990.
- D. Kesh, D. Mukherjee, A. K. Sarkar, and A. B. Roy, “Ratio dependent predation: a bifurcation analysis,” The Korean Journal of Computational & Applied Mathematics, vol. 5, no. 2, pp. 295–305, 1998.
- F. Bartumeus, D. Alonso, and J. Catalan, “Self-organized spatial structures in a ratio-dependent predator-prey model,” Physica A: Statistical Mechanics and its Applications, vol. 295, no. 1-2, pp. 53–57, 2001.
- L. R. Ginzburg and R. X. J. Jensen, “From controversy to consensus: the indirect interference functional response,” Verhandlungen des Internationalen Verein Limnologie, vol. 30, no. 2, pp. 297–301, 2008.
- C. Cosner, D. L. Deangelis, J. S. Ault, and D. B. Olson, “Effects of spatial grouping on the functional response of predators,” Theoretical Population Biology, vol. 56, no. 1, pp. 65–75, 1999.
- L. R. Ginzburg and H. R. Akcakaya, “Consequences of ratio-dependent predation for steady-state properties of ecosystems,” Ecology, vol. 73, no. 5, pp. 1536–1543, 1992.
- D. Alstad, Basic Populas Models of Ecology, Prentice Hall, Upper Saddle River, NJ, USA, 2001.
- O. Anderson, “Optimal foraging by largemouth bass in structured environments,” Ecology, vol. 65, no. 3, pp. 851–861, 1984.
- T. W. Anderson, “Predator responses, prey refuges, and density-dependent mortality of a marine fish,” Ecology, vol. 82, no. 1, pp. 245–257, 2001.
- D. W. Johnson, “Predation, habitat complexity, and variation in density-dependent mortality of temperate reef fishes,” Ecology, vol. 87, no. 5, pp. 1179–1188, 2006.
- J. Maynard Smith, Models in Ecology, Cambridge University Press, Cambridge, UK, 1974.
- M. P. Hassel, The Dynamics of Arthropod Predator-Prey Systems, Princeton University Press, Princeton, NJ, USA, 1978.
- R. J. Taylor, Predation, Chapman & Hall, New York, NY, USA, 1984.
- T. K. Kar, “Stability analysis of a prey-predator model incorporating a prey refuge,” Communications in Nonlinear Science and Numerical Simulation, vol. 10, no. 6, pp. 681–691, 2005.
- D. Jana, “Chaotic dynamics of a discrete predator-prey system with prey refuge,” Applied Mathematics and Computation, vol. 224, pp. 848–865, 2013.
- M. Liao, X. Tang, and C. Xu, “Bifurcation analysis for a three-species predator-prey system with two delays,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 1, pp. 183–194, 2012.
- X.-Y. Meng, H.-F. Huo, and X.-B. Zhang, “Stability and global Hopf bifurcation in a delayed food web consisting of a prey and two predators,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 11, pp. 4335–4348, 2011.
- G.-P. Hu, W.-T. Li, and X.-P. Yan, “Hopf bifurcations in a predator-prey system with multiple delays,” Chaos, Solitons and Fractals, vol. 42, no. 2, pp. 1273–1285, 2009.
- Y. Song, Y. Peng, and J. Wei, “Bifurcations for a predator-prey system with two delays,” Journal of Mathematical Analysis and Applications, vol. 337, no. 1, pp. 466–479, 2008.
- X.-P. Yan and Y.-D. Chu, “Stability and bifurcation analysis for a delayed Lotka-Volterra predator-prey system,” Journal of Computational and Applied Mathematics, vol. 196, no. 1, pp. 198–210, 2006.
- H. I. Freedman, Deterministic Mathematical Models in Population Ecology, HIFR Consulting, Edmonton, Canada, 1987.
- Y. Kuang, Delay Differential Equations with Applications in Population Dynamics, Academic Press, New York, NY, USA, 1993.
- N. Bairagi and D. Jana, “On the stability and Hopf bifurcation of a delay-induced predator-prey system with habitat complexity,” Applied Mathematical Modelling, vol. 35, no. 7, pp. 3255–3267, 2011.
- S. Ruan, “On nonlinear dynamics of predator-prey models with discrete delay,” Mathematical Modelling of Natural Phenomena, vol. 4, no. 2, pp. 140–188, 2009.
- S. Ruan and J. Wei, “On the zeros of transcendental functions with applications to stability of delay differential equations with two delays,” Dynamics of Continuous, Discrete and Impulsive Systems Series A: Mathematical Analysis, vol. 10, no. 6, pp. 863–874, 2003.
- J. Hale, Theory of Functional Differential Equations, Springer, Berlin, Germany, 1977.
- S. Mandal, D. Jana, A. B. Roy, and N. C. Majee, “Chaotic behavior of a class of neural network with discrete delays,” International Journal of Modern Nonlinear Theory and Application, vol. 2, no. 1A, pp. 97–101, 2013.
- S. H. Strogatz, Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering, Westview Press, 2009.
- K. Wang, W. Wang, H. Pang, and X. Liu, “Complex dynamic behavior in a viral model with delayed immune response,” Physica D: Nonlinear Phenomena, vol. 226, no. 2, pp. 197–208, 2007.
- A. H. Nayfeh and B. Balachandran, Applied Nonlinear Dynamics, Wiley, New York, NY, USA, 1995.
- N. Marwan, M. Carmen Romano, M. Thiel, and J. Kurths, “Recurrence plots for the analysis of complex systems,” Physics Reports, vol. 438, no. 5-6, pp. 237–329, 2007.