Mathematical Problems in Engineering

Volume 2013, Article ID 469072, 8 pages

http://dx.doi.org/10.1155/2013/469072

## Global Dynamics of the Hastings-Powell System

Departamento de Ingeniería Eléctrica y Electrónica, Instituto Tecnológico de Tijuana, Boulevard Alberto Limón Padilla s/n, Mesa de Otay, 22454 Tijuana, BCN, Mexico

Received 4 September 2013; Revised 16 November 2013; Accepted 17 November 2013

Academic Editor: Sebastian Anita

Copyright © 2013 Luis N. Coria. 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

This paper studies the problem of bounding a domain that contains all compact invariant sets of the Hastings-Powell system. The results were obtained using the first-order extremum conditions and the iterative theorem to a biologically meaningful model. As a result, we calculate the bounds given by a tetrahedron with excisions, described by several inequalities of the state variables and system parameters. Therefore, a region is identified where all the system dynamics are located, that is, its compact invariant sets: equilibrium points, periodic-homoclinic-heteroclinic orbits, and chaotic attractors. It was also possible to formulate a nonexistence condition of the compact invariant sets. Additionally, numerical simulations provide examples of the calculated boundaries for the chaotic attractors or periodic orbits. The results provide insights regarding the global dynamics of the system.

#### 1. Introduction

Chaos is an interesting phenomenon that has been studied extensively over recent years. Moreover, it continues to be an active research topic, particularly in the area of analysis of nonlinear dynamic systems. Chaotic nature can be found in a huge variety of disciplines; interesting examples can be found in financial systems [1], cancer tumor growth models [2], power systems [3], weather models [4], electromechanical systems [5], and many more. For instance, complex dynamics are detected in several biological models, with many works published in recent years [6–10]. In particular, this paper analyzes a chaotic biological model which describes a tritrophic food chain interaction within an ecosystem, commonly referred to as the Hastings-Powell (HP) system [11]. Dynamic properties of this model have been studied in some contributions; for instance, see papers [12–14].

The main interest in bounding a dynamical system is due to the fact that such bounds can provide relevant information regarding its global dynamics. The bounds are domains given by inequalities, expressed in terms of the parameters and state variables of a system. Within these bounds, it is possible to find every compact invariant set of the system, which characterize its most interesting dynamics which include periodic, homoclinic, or heteroclinic orbits; equilibrium points, and chaotic behavior. Besides the purely theoretical interest of obtaining such bounds, the results may provide valuable insights regarding the evolution and future of the considered system.

The bounds in this work are obtained by the application of first-order extremum conditions, proposed initially by Krishchenko [15] and improved in recent years by Krishchenko and Starkov: they propose the iterative method of localization of compact invariant sets in [16], ellipsoidal estimates in [17], and localization for time-varyng systems in [18], for Hamiltonian systems [19], and they continue to work on the development of the general theory and applications of the method. One of the most interesting improvements of the method is the iterative theorem, which can be used to sharpen the bounds based on previous results of localization [16]. This work includes some interesting results obtained by applying the iterative theorem (Theorem 3 described below). The localization method has been successfully applied recently on dynamic system in diverse areas such as permanent-magnet motor system [5], Lanford system [20], Lorenz system [21], the amplitude of a plasma instability [22], and Nonetheless, recently an interest has arisen to understand the global dynamics of a biologically meaningful model. Analyzing the HP system contributes with an alternative way to understanding its behavior.

Similar analytic methods (see [8, 23–25] and references therein) calculate the geometric approximation of attractors or special cases of periodic dynamics. We calculate a region that includes all compact invariant sets of the system. Our results can be used to analyze the global dynamics and combined with nonlinear analysis theories may provide important information for engineering applications [26, 27]. Regarding previous results obtained for the global dynamics analysis of the HP system [28–30], we can comment that our approach gives a novel perspective of the boundaries of the system. Besides, various essential properties of dynamics like nonchaoticity and nonexistence of repeating behavior may be analyzed with the help of this approach; see [31, 32].

The paper is organized as follows. Some useful assertions and notations are provided in Section 2. Section 3 presents results concerning the bounds of the Hastings-Powell model. This section is divided into several subsections, which presents diverse applications of the localization method. Section 4 includes numerical simulation examples of the bounds with respect to some compact invariant sets. Finally Section 5 contains conclusions.

#### 2. Preliminaries of Bounding

This section is devoted to describing the general theorems used to obtain the bounds of the system. This method is commonly referred to as the localization of compact invariant sets; however, in this paper we refer to it as the “localization method.” The general localization method for a nonlinear system was described by Starkov and Krishchenko in [15, 20–22]. This section reviews useful results from these works.

Consider a nonlinear system where is a -differentiable vector field and is the state vector. Let be a -differentiable function such that is not the first integral of (1). Take the restriction of on a set . denotes the set , where is a Lie derivative with respect to vector field corresponding to the system (1). Now define ; .

Theorem 1. *Each compact invariant set of (1) is contained in the localization set [15] .*

If the location of all compact invariant sets is considered inside the domain , then the localization set is valid, with defined in Theorem 1. Let be a subset in .

Proposition 2. *If , then the system (1) has no compact invariant sets located in [15].*

A refinement of the localization set is realized with the help of the iterative theorem, which is stated as follows.

Theorem 3. *Let , be a sequence of functions from [16]. Sets
**
with
**
contain any compact invariant set of the system (1) and
*

#### 3. The Hastings-Powell System

The model of interest describes a tritrophic food chain. This system was obtained originally by Hogeweg and Hesper [33]. Later, Hastings and Powell [11] found chaotic dynamics, and is now known as the Hastings-Powell (HP) system in many papers from the field of ecology [12, 34, 35]. The model is given by the following system of equations: with describing a type II functional response of both of the consumer species. Here, is time, is the top-predator, is the prey and is a predator that consumes the prey . The constant is the intrinsic growth rate. is the carrying capacity of species . Constants and are conversion rates of prey to predator of species and , respectively. and are constant death rates for species and respectively. Constants and , for , parametrize the saturation functional response. is the prey population level where the predation rate per unit prey is half its maximum value.

Model (5) includes a set of 10 parameters that imply a difficult analysis. As is presented in [11], the parameters can be reduced to 6 with the linear transformation given by Then, it is possible to have a simpler nondimensional form written as where the new parameters are given by , , , , , ; finally with denotes the derivatives of with respect to .

This model was used to describe the dynamics of classical fishpond management for the tilapia fish culture [12, 36]. These works consider , like the corresponding density of young tilapia (prey), developed tilapia (predator), and variable describes the density of tucunare fish which is considered the top-predator population. Gomes et al. [12] and Varriale and Gomes [36] analyzed the model and found a rich variety of dynamics. Parameters , , , , and , for system (8), are used traditionally; see [11]. Figure 1 shows the named “tea-cup” chaotic attractor obtained for specific parameters of the system (8). From an engineering perspective, the main interest of the HP system is to study secondary effects of the exploitation of biological resources and the harvest in an ecosystem.

In the following section, the bounds for the system (8) are calculated with the application of first-order extremum conditions and the iterative theorem. Numerical simulations are performed to show the bounds of the tea-cup attractor and other interesting compact invariant sets.

#### 4. Bounding the Hastings-Powell System

This section is devoted to calculating the bounds of system (8). The bounds are calculated with the help of coordinate functions, linear combinations, and rational functions. Therefore, a tetrahedral localization with excisions (cuts over its geometric body) is derived as a result. Additionally, some nonexistence conditions are derived.

Since variables , in (8), have a biological interpretation we examine compact invariant sets that lie within the positive orthant . In addition, all parameters of the model are supposed to be positive and is an invariant set. In what follows, denotes the vector field of the system under consideration. Also, for the sake of simplicity of notations, and .

##### 4.1. Tetrahedral Localization

Firstly, a compact localization set in will be derived. This means that the region of localization delimits all variables of the system. This localization set is given by the following theorem.

Theorem 4. *Let ; if then all compact invariant sets are located in the tetrahedron
*

*Proof. *We have . Therefore, the set is defined by
taking . Now by using the last expression for
This means that
and the localization set is given by .

This localization means that every compact invariant set of the Hastings-Powell system (8) is bounded in the region given by . This tetrahedral localization includes the restriction of , which coincides with the biological feasible web-food existence condition () given by [11]. The upper limit of the localization set given by varies with resoect to (the constant death rate of the top-predator), where the bound of can be obtained from the localization set , and with , this coincides with that given by [11], (three lines before Table 1 in [11], page 897), and the obtained latter in this paper.

##### 4.2. Excisions by Means of Additional Linear Functions

With we have obtained a compact localization; nevertheless, some additional excisions (cuts over the tetrahedron body) can be obtained with the help of coordinate functions and the application of the iterative theorem. The next theorem gives a localization set from a coordinate function. The first result is given by the next theorem

Theorem 5. *All compact invariant sets located in are contained in the localization set
*

*Proof. *If the localizing function is taken, then , and
Thus,
and the localization set is given by .

Now, with the application of the iterative theorem, the next theorem can be stated.

Theorem 6. *All compact invariant sets located in are contained in the localization set , with
*

*Proof. *If the localizing function is used, then and
Hence solving respecting and , then
So because , we have
As a result, one can get

The latter formula of is an alternative solution of the upper bound for derived from . In this case, the bound of was obtained in a similar way by [11], but in this work we obtained a tighter bound by applying the iterative theorem. Here, no restrictions are imposed on the localization set. The variation of the bound w.r.t. (the constant death rate of the middle-predator) implies that the minimal value of the bound is and, as expected, coincides with the bound . Now, an additional bound is given by the following theorem.

Theorem 7. *If then all compact invariant sets are located in
*

*Proof. *Proposing the localizing function , then
Therefore by taking by completing squares in a similar way to Theorem 4, taking , then
and the localization set is given by .

This localization gives an additional excision of the localization set.

##### 4.3. Rational Functions

Here, rational functions are applied to sharpen the bounds; this is possible because of the iterative theorem. This improves previous bounds; however, there are restrictions on them. These restrictions are not fulfilled with the parameters given in the literature, but they are included for knowledge in future research. Then with the rational function it is possible to state the following theorem.

Theorem 8. *Suppose that
**
Then, all compact invariant sets located in are contained in the set
*

*Proof. *The localizing function is applied respecting , then
Therefore, on the formula
holds. The lower bound for the coefficient of in the last formula with respect to the set will be calculated with

In a similar way, it is possible to state the next.

Theorem 9. *Suppose that and
**
Then, all compact invariant sets located in are contained in the set
*

*Proof. *If the localizing function is used respecting , then
It follows from the latter formula that
is valid on the set . Indeed, since , then includes
then the localization set is obtained because of Theorem 3 (iterative theorem). Proof completed.

##### 4.4. Nonexistence Conditions

Non-existence conditions in dynamic systems involve the possibility to find some range of parameters where the system lacks compact invariant sets. Therefore, in the case of system (8), if it is possible to find two nonexistence conditions, which implies that if a compact invariant set exists for the system, then, must be located inside the coordinate planes. This conclusion can be obtained due to the main characteristic of the biological system regarding the positiveness and invariance of the first orthant . For system (8), two non-existence conditions are calculated and are given by the next theorems.

Theorem 10. *If , then there are no compact invariant sets located in .*

*Proof. *If the function is taken, then
and can the following set be obtained
Therefore, if , . The proof is complete because of the Proposition 2.

Now, it is possible to find an additional nonexistence condition.

Theorem 11. *If , then there are no compact invariant sets located in .*

*Proof. *If the function is taken, then , and the following set can be obtained:
Therefore, if , . The proof is complete based on Proposition 2.

##### 4.5. Condition of Survival of Predator Populations

At this point, it is relevant to mention one noteworthy remark concerning dynamics of system (8). The necessary condition for the population to survive is , . Indeed, if then we have that which entails in case . The case can be considered similarly and in case .

The nonexistence conditions coincide with the condition of survival of predators. Here, it is important to mention that they share a similar interpretation but have specific differences between them. A nonexistence condition means that there are no compact invariant sets in . The condition of survival for the predator means that if the condition is fulfilled, then the population of predator will be extinct. The combination of both results allows us to assert the attractivity of coordinate planes and . Therefore, if one of the conditions is fulfilled, the dynamics of the system (8) are trapped inside a specific coordinate plane and a reduced order system is obtained.

#### 5. Numerical Simulations

This section presents numerical simulations for a few special cases of parameters, given by Kuznetsov and Rinaldi [8]; the localization set is valid, and as expected, the localization set is different for each special case. This boundary depends on the parameters and , the constant death rates of the species, meaning that the localization set is described by the parameters of the system. Compared with numerical methods, localization surface is given by a mathematical function of the parameters and state variables of the system.

Inequalities for the existence of other localization sets are not satisfied for these values of parameters. In particular, and are negative. However, if then the localization by the set exists. Besides, as it was indicated in [11], biologically reasonable food chains are realized when , which entails .

##### 5.1. Periodic Orbits

###### 5.1.1. Periodic Orbit [8]

With the parameters ; ; ; ; ; ; [8] it is a periodic orbit behavior is obtained. The bounds are given by the set ; in this case do not cut the tetrahedron . The periodic orbit and the bounds for these parameters are shown in Figure 2, where , and , and .

###### 5.1.2. Two Periodic Orbits [8]

With the set of parameters ; ; ; ; ; ; [8] is possible to obtain different periodic orbits for the system (8).Figure 3 shows two periodic orbits obtained for this set of parameters and the boundary surfaces. In this case do not cut the tetrahedron . Here .

###### 5.1.3. Periodic Orbit [8]

For parameters ; ; ; ; ; ; [8], an interesting behavior is given by a periodic orbit of the system (8). Figure 4 shows a periodic orbit obtained for this set of parameters and the bounds given by . In this case and do not cut the tetrahedron .Here .

###### 5.1.4. Periodic Orbit [8]

Taking the set of parameters , , , , and , [8], it is possible to obtain another periodic orbit. Similar to the previous localizations, do not cut the tetrahedron . Figure 5 shows the boundary given by , and the periodic orbit for these parameters. Here ; ; and .

###### 5.1.5. Periodic Orbit [8]

Consider , , , , , and [8] and again do not cut the tetrahedron . Figure 6 shows the boundary set and this periodic orbit, localization sets are given by , and .

##### 5.2. Chaotic Attractor

Parameters , , , , , for system (8), are used traditionally, see [7, 8, 11, 12, 36]. In this case, the bifurcation parameter is varied within the segment . A chaotic attractor can be observed at ; see in [12]. Figure 7 provides a localization of this attractor by the set , with ; ; .

The upper side of the plot in Figure 7 is a piece of the face of the tetrahedral localization ; the latter has three excisions corresponding to the set .

#### 6. Conclusions

This paper describes how to obtain a compact boundary containing all compact invariant sets of a biological system which is given by a three-species food model called the Hastings-Powell system. These bounds are given by a tetrahedron with linear excisions. The linear excisions are obtained by using coordinate and a rational localizing function. The bounds depend on the constant death rates for each species.

Our method requires a candidate with some function, and we have to use heuristic methods for its choice. Here we deal with biological systems that lie within the positive orthant ; therefore, our choice of that function is atypical compared with traditional polynomial systems where quadratic-style functions are usually proposed. Here linear, coordinate, and rational functions give better results. Moreover, the iterative method for localization can be applied to sharpen the bounds. Two nonexistence conditions of compact invariant sets are described as well, and two conditions for the survival of predator populations are provided.

Numerical simulations are performed to show the bounds of some compact invariant sets, that is, periodic orbits and chaotic attractors. The bounds provide information regarding the global dynamics of the Hastings-Powell system. The nonexistence conditions allowed us to find an important contribution about the dynamics of the system combined with the conditions of predator survival, about attractivity of coordinate planes and . Moreover, numerical simulations show a graphical interpretation of the mathematical work.

#### Conflict of Interests

The author declare that there is no conflict of interests.

#### Acknowledgments

This work was supported by DGEST Project TIJ-ING-2012-111 and the PROMEP Research Group ITTIJ-CA-6. The author would like to acknowledge Dr. Leonardo Trujillo and Dr. Konstantin E. Starkov for their contributions in various aspects of the preparation in this work. The author is pleased to acknowledge the helpful and insightful comments of the reviewers. A part of this work was presented in the IFAC CHAOS 2009 congress [37].

#### References

- D. A. Hsieh, “Chaos and nonlinear dynamics: application to financial markets,”
*The Journal of Finance*, vol. 46, no. 5, pp. 1839–1877, 1991. View at Google Scholar - A. Dalgleish, “The relevance of non-linear mathematics (chaos theory) to the treatment of cancer, the role of the immune response and the potential for vaccines,”
*Monthly Journal of the Association of Physicians*, vol. 92, no. 6, pp. 347–359, 1999. View at Google Scholar · View at Scopus - A. Ben-Tal, V. Kirk, and G. Wake, “Banded chaos in power systems,”
*IEEE Transactions on Power Delivery*, vol. 16, no. 1, pp. 105–110, 2001. View at Publisher · View at Google Scholar · View at Scopus - E. Tziperman, L. Stone, M. A. Cane, and H. Jarosh, “El Nino chaos: overlapping of resonances between the seasonal cycle and the Pacific ocean-atmosphere oscillator,”
*Science*, vol. 264, no. 5155, pp. 72–74, 1994. View at Google Scholar · View at Scopus - L. N. Coria and K. E. Starkov, “Bounding a domain containing all compact invariant sets of the permanent-magnet motor system,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 14, no. 11, pp. 3879–3888, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. N. K. Anderson, C.-H. Hsieh, S. A. Sandin et al., “Why fishing magnifies fluctuations in fish abundance,”
*Nature*, vol. 452, no. 7189, pp. 835–839, 2008. View at Publisher · View at Google Scholar · View at Scopus - A. Klebanoff and A. Hastings, “Chaos in three-species food chains,”
*Journal of Mathematical Biology*, vol. 32, no. 5, pp. 427–451, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. A. Kuznetsov and S. Rinaldi, “Remarks on food chain dynamics,”
*Mathematical Biosciences*, vol. 134, no. 1, pp. 1–33, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. K. Upadhyay, R. K. Naji, and N. Kumari, “Dynamical complexity in some ecological models: effect of toxin production by phytoplankton,”
*Lithuanian Association of Nonlinear Analysts (LANA). Nonlinear Analysis. Modelling and Control*, vol. 12, no. 1, pp. 123–138, 2007. View at Google Scholar · View at MathSciNet - J. T. Wootton, “Experimental species removal alters ecological dynamics in a natural ecosystem,”
*Ecology*, vol. 91, no. 1, pp. 42–48, 2010. View at Publisher · View at Google Scholar · View at Scopus - A. Hastings and T. Powell, “Chaos in a three-species food chain,”
*Ecology*, vol. 72, no. 3, pp. 896–903, 1991. View at Google Scholar · View at Scopus - A. A. Gomes, E. Manica, and M. C. Varriale, “Applications of chaos control techniques to a three-species food chain,”
*Chaos, Solitons and Fractals*, vol. 35, no. 3, pp. 432–441, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. McCann and P. Yodzis, “Biological conditions for chaos in a three-species food chain,”
*Ecology*, vol. 75, no. 2, pp. 561–564, 1994. View at Google Scholar · View at Scopus - K. McCann and P. Yodis, “Bifurcation structure of a three-species food chains model,”
*Theoretical Population Biology*, vol. 48, no. 1, pp. 93–125, 1995. View at Google Scholar - A. P. Krishchenko, “Estimations of domains with cycles,”
*Computers & Mathematics with Applications*, vol. 34, no. 2–4, pp. 325–332, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. P. Krishchenko and K. E. Starkov, “Iteration method of the localization of periodic orbits,” in
*Proceedings of the International Conference on Physics and Control (PhysCon '05)*, pp. 602–605, St. Petersburg, Russia, August 2005. View at Publisher · View at Google Scholar · View at Scopus - K. E. Starkov and A. P. Krishchenko, “Localization of periodic orbits of polynomial systems by ellipsoidal estimates,”
*Chaos, Solitons and Fractals*, vol. 23, no. 3, pp. 981–988, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. P. Krishchenko and K. E. Starkov, “Localization of compact invariant sets of nonlinear time-varying systems,”
*International Journal of Bifurcation and Chaos in Applied Sciences and Engineering*, vol. 18, no. 5, pp. 1599–1604, 2008. View at Publisher · View at Google Scholar · View at MathSciNet - K. E. Starkov, “Universal localizing bounds for compact invariant sets of natural polynomial Hamiltonian systems,”
*Physics Letters A*, vol. 372, no. 41, pp. 6269–6272, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. P. Krishchenko and K. E. Starkov, “Localization of compact invariant sets of nonlinear systems with applications to the Lanford system,”
*International Journal of Bifurcation and Chaos in Applied Sciences and Engineering*, vol. 16, no. 11, pp. 3249–3256, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. P. Krishchenko and K. E. Starkov, “Localization of compact invariant sets of the Lorenz system,”
*Physics Letters A*, vol. 353, no. 5, pp. 383–388, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Krishchenko and K. Starkov, “Estimation of the domain containing all compact invariant sets of a system modelling the amplitude of a plasma instability,”
*Physics Letters A*, vol. 367, no. 1-2, pp. 65–72, 2007. View at Publisher · View at Google Scholar · View at Scopus - B. Deng, “Food chain chaos with canard explosion,”
*Chaos*, vol. 14, no. 4, pp. 1083–1092, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. A. Kuznetsov, S. Muratori, and S. Rinaldi, “Homoclinic bifurcations in slow-fast second order systems,”
*Nonlinear Analysis. Theory, Methods & Applications*, vol. 25, no. 7, pp. 747–762, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Muratori and S. Rinaldi, “Low- and high-frequency oscillations in three-dimensional food chain systems,”
*SIAM Journal on Applied Mathematics*, vol. 52, no. 6, pp. 1688–1706, 1992. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. E. Starkov and L. N. Coria, “Global dynamics of the Kirschner-Panetta model for the tumor immunotherapy,”
*Nonlinear Analysis. Real World Applications*, vol. 14, no. 3, pp. 1425–1433, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. E. Starkov, L. N. Coria, and L. T. Aguilar, “On synchronization of chaotic systems based on the Thau observer design,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 17, no. 1, pp. 17–25, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. I. Freedman and P. Waltman, “Mathematical analysis of some three-species food-chain models,”
*Mathematical Biosciences*, vol. 33, no. 3-4, pp. 257–276, 1977. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. I. Freedman and J. W.-H. So, “Global stability and persistence of simple food chains,”
*Mathematical Biosciences*, vol. 76, no. 1, pp. 69–86, 1985. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Kuang, “Global stability and persistence in diffusive food chains,”
*The ANZIAM Journal*, vol. 43, no. 2, pp. 247–268, 2001. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. E. Starkov, “Compact invariant sets of the static spherically symmetric Einstein-Yang-Mills equations,”
*Physics Letters A*, vol. 374, no. 15-16, pp. 1728–1731, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. E. Starkov, “Compact invariant sets of the Bianchi VIII and Bianchi IX Hamiltonian systems,”
*Physics Letters A*, vol. 375, no. 36, pp. 3184–3187, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. Hogeweg and B. Hesper, “Interactive instruction on population interactions,”
*Computers in Biology and Medicine*, vol. 8, no. 4, pp. 319–327, 1978. View at Google Scholar · View at Scopus - K. E. Lonngren, B. Er-Wei, and U. Ahmet, “Dynamics and synchronization of the Hastings-Powell model of the food chain,”
*Chaos, Solitons & Fractals*, vol. 20, no. 2, pp. 387–393, 2004. View at Google Scholar - C. T. Perretti, G. Sugihara, and S. B. Munch, “Nonparametric forecasting outperforms parametric methods for a simulated multispecies system,”
*Ecology*, vol. 94, pp. 794–800, 2013. View at Google Scholar - M. C. Varriale and A. A. Gomes, “A study of a three species food chain,”
*Ecological Modelling*, vol. 110, no. 2, pp. 119–133, 1998. View at Publisher · View at Google Scholar · View at Scopus - K. Starkov and L. N. Coria, “Bounding the domain of some three species food systems,” in
*Proceedings of the 2nd IFAC Conference on Analysis and Control of Chaotic Systems*, Queen Mary University of London, London, UK, 2009, http://www.ifac-papersonline.net/.