- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Journal of Applied Mathematics

Volume 2014 (2014), Article ID 803764, 11 pages

http://dx.doi.org/10.1155/2014/803764

## Periodic Solutions and Homoclinic Bifurcations of Two Predator-Prey Systems with Nonmonotonic Functional Response and Impulsive Harvesting

^{1}Research Institute of Forest Resource Information Techniques, Chinese Academy of Forestry, Beijing 100091, China^{2}College of Mathematics and Information Science, Xinyang Normal University, Xinyang 464000, China

Received 23 October 2013; Accepted 24 February 2014; Published 27 March 2014

Academic Editor: Shiping Lu

Copyright © 2014 Mingzhan Huang and Xinyu Song. 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

Two predator-prey models with nonmonotonic functional response and state-dependent impulsive harvesting are formulated and analyzed. By using the geometry theory of semicontinuous dynamic system, we obtain the existence, uniqueness, and stability of the periodic solution and analyse the dynamic phenomenon of homoclinic bifurcation of the first system by choosing the harvesting rate as control parameter. Besides, we also study the homoclinic bifurcation of the second system about parameter on the basis of the theory of rotated vector field. Finally, numerical simulations are presented to illustrate the results.

#### 1. Introduction

Predator-prey interaction is one of the most important relationships in the ecosystem, so it has long been a focus of study in mathematical ecology. The functional response of the predator to the prey describes how the predator density changes with the prey density. The well-known response functions, which have appeared in a lot of literatures, include linear function, Holling type-II function, Holling type-III function, and so on. The above mentioned response functions are all monotonic, and they are really accurate to describe the predator-prey interactions in many cases. But in the recent decades, scholars in different fields found through experiments that nonmonotonic functional response occurs in some predator-prey interactions. For example, if the prey exhibits group defense [1, 2], the predator growth rate will be inhibited when the prey density reaches a high level. This phenomenon is known to exist widely in nature and considerable work on it has been studied [3–8]. To study the predator-prey interaction when the prey exhibits group defense, Ruan and Xiao in [4] proposed the following model: where and represent the population densities of prey and predator, respectively; and represent the carrying capacity and the intrinsic birth rate of the prey, respectively, and is the death rate of the predator. The function denotes the predator functional response which can model the phenomenon of group defense because there exists such that for and for .

Harvesting strategy of biological resources is also a focus topic in mathematical bioeconomics because it relates to the optimal management of renewable resources. In many literatures, impulsive differential equations are used to model the human action of harvesting. Predator-prey systems with periodic impulsive harvesting have been studied extensively and important results have been achieved [9–12]. In recent years, more and more scholars begin to consider state dependent impulsive harvesting instead of the periodic impulsive form [13–17]. In point of ecology, we would better know some information about the amount of the species when we harvest them, thus we can avoid excessive exploitation and resource exhaustion. To this end, we introduce a reliable real-time monitoring system to estimate the number of the species. Such monitoring systems exist in many cases. According to the feedback information from the monitoring system, we can manage our resources better. For predator-prey system (1), we assume the predator has high commercial value for human beings, and people increase its production mainly through replenishing its prey. In such an ecosystem, we suppose the amount of the prey can be estimated by a monitoring system, and the monitoring data can help us to decide if we harvest the predator or not. To model this harvesting behavior, we can propose the following state-dependent impulsive differential equations: where is a threshold. When the amount of the prey is large than , which means the food of the predator is abundant, and the development of the system coincide with our economic interest. When the amount of the prey drops to the threshold , which means the nutrition of the predator will be deficient, we harvest the predator at rate and replenish some prey at the same time. We denote the recruitment of the prey as .

A lot of research has shown that the population size of the predator is influenced not only by the prey populations but also by the relative rate of prey population growth. If we take the relative growth rate effect into consideration, the above model needs some changes. To this end, we propose the following system: where the term represents the effect of the prey relative growth rate on the predator.

In this paper, we mainly discuss the dynamics properties of the systems (2) and (3). The paper is organized as follows. In Section 2, some notations and definitions of the geometric theory of semicontinuous dynamical systems are provided. In Section 3, by using the geometry theory of semicontinuous dynamic system, we firstly study the existence, uniqueness, and orbital stability of periodic solutions and analyse the dynamic phenomenon of homoclinic bifurcation of system (2) by choosing the harvesting rate as control parameter. Then, we study the homoclinic bifurcation of system (3) about parameter on the basis of the theory of rotated vector field. The paper ends with a brief discussion and some numerical simulations.

#### 2. Preliminaries

In this section, we give some notations and definitions of the geometric theory of semicontinuous dynamical systems which will be useful for the following discussion.

*Definition 1 (see [9]). *Consider the state-dependent impulsive differential equations

We define the dynamic system consisting of the solution mapping of the system (4) a semicontinuous dynamical system, denoted as . We require that the initial point of the system (4) should not be in the set , that is, , and the function is a continuous mapping that satisfies . Here is called the impulse mapping, where and represent the straight lines or curves in the plane , is called the impulse set, and is called the phase set.

*Remark 2. *For the systems (2) and (3), , , and for any , we have .

*Definition 3 (see [9]). *For the semicontinuous dynamical system defined by the state-dependent impulsive differential equations (4), the solution mapping consists of two parts.(1)Let denote the poincaré mapping with the initial point of the following system:
If , then .(2)If there exists a time point such that , , and , then (see Figure 1(a)).

*Remark 4. *For (2) in Definition 3, if , and there exists a time point such that , , and , then .

If and , then we can repeat the above steps and have the following form:

*Definition 5 (see [9]). *If there exists a point and a time point such that and , then is called an order one periodic solution of the system (4) whose period is (see Figure 1(b)). The orbit of the order one periodic solution is called an order one cycle. If there exists a singularity in the order one cycle, we call it an order one singular cycle. If the singularity is a saddle, we call it an order one homoclinic cycle.

*Definition 6 (see [9]). *If there exists a point and a time point such that , and , and there also exists a time point such that and , then is called an order two periodic solution of the system (4) whose period is (see Figure 1(c)). Analogously, we can define the order periodic solution of the system (4).

*Definition 7. *Suppose is an order one periodic solution of the system (4). If for any , there must exist and , such that for any point , we have for ; then we call the order one periodic solution orbitally asymptotically stable.

*Definition 8 (see [9]). *Suppose the impulse set and phase set of the system (4) are straight lines and a coordinate system can be defined in the phase set . Let point and its coordinate is . Assume that the trajectory from the point intersects the impulse set at a point , and, after impulsive effect, the point is mapped to the point with the coordinate ; then we call point the order one successor point of point , and the order one successor function of point is .

*Remark 9. *For Definition 8, if the trajectory from the point intersects the impulse set again at a point , and, after impulsive effect, the point is mapped to the point with the coordinate , then the point is obviously the order one successor point of point ; we also call point the order two successor point of point , and the order two successor function of point is . If the process can be repeated over and over again, then we can define the order successor point of point (which we denote as and its coordinate is ) and the order successor function of which we denote as .

*Remark 10. *For the systems (2) and (3), for any point , we define the directed distance between point and the -axis as the coordinate of point . In this paper, we denote the coordinate of point as .

Lemma 11 (see [9]). *Successor function is continuous.*

Lemma 12. *For the systems (2) and (3), if there exists two points , such that , then there must exist a point which is between the points and such that ; thus, the system must have an order one periodic solution which passes through the point .*

*Proof. *By Lemma 11, we can easily get that there must exist a point which is between the points and such that . According to Definition 6, we know is an order one periodic solution. That completes the proof.

Lemma 13. *For the systems (2) and (3), suppose point and , then the system does not have an order periodic solution which passes through the point .*

#### 3. Periodic Solution and Homoclinic Bifurcation

In this section, we firstly discuss the existence, uniqueness, and stability of the periodic solution of the system (2) by using differential equation geometry theory then study the dynamic phenomena of homoclinic bifurcation of both systems (2) and (3). We now list some results about the system (1) that are given by Ruan and Xiao in [4].

For the system (1), there is always a hyperbolic saddle point at the origin and an equilibrium in the -axis. When , the system (1) may have one or two positive equilibria. If the positive equilibria exist, we denote them as , , where Now, we give some results appeared in [4].

Lemma 14 (see [4]). *If and (or if and ), then the system (1) has three equilibria: two hyperbolic saddles and and a globally asymptotically stable equilibrium in the interior of the first quadrant.*

Lemma 15 (see [4]). *If and (or if and ), then the system (1) has four equilibria: two hyperbolic saddles and , a hyperbolic stable node , and a stable (an unstable) equilibrium , and system (1) has no closed orbits.*

##### 3.1. Periodic Solution and Homoclinic Bifurcation of System (2) about Parameter *β*

According to the model (2), the threshold and the recruitment of the prey should satisfy the condition by ecological significance. For this consideration, we have the following results.

Theorem 16. *If , (or if , ), and , then there must exist a fixed value such that for every , the system (2) has a unique order one periodic solution in region , where region is the region enclosed by the -axis, the impulse set , and the unstable flow of the saddle .*

*Proof. *According to Lemma 14, the system (1) has three equilibria: two hyperbolic saddles and and a globally asymptotically stable equilibrium in the interior of the first quadrant. For convenience, we suppose the -axis intersects the impulse set and the phase set at point and point , respectively. The unstable flow of also passes through the impulse set and the phase set, we denote the intersections by point and point , respectively. Besides, we suppose the vertical isocline intersects the impulse set and the phase set at point and point , respectively. Then the region is the interior of the closed curve , where denotes the equilibrium (see Figure 2).

According to the impulsive conditions of the system (2), there must exist a fixed value , when , point is mapped to the point after impulsive effect, that is to say, .

When , after impulsive effect, point is mapped to a point between points and . Select a point next to in the phase set , the trajectory of the system (2) from point must intersect the impulse set at a point which is next to point , and after impulsive effect, point is mapped to a point which is the order one successor point of ; then, we have . That is to say point is between points and . The trajectory from point must intersect the impulse set again at a point , and after impulsive effect, point is mapped to a point . Since distinct trajectories do not intersect, we can easily have and . Obviously, point is the order one successor point of point , and then we have the following results of the order one successor function:
By Lemma 12, we know that there must exist a point in the phase set which is between the points and such that . Then we know the system (2) has an order one periodic solution which passes through the point .

In the following, we prove the uniqueness of the order one periodic solution. Arbitrarily choose two points and in the phase set , where . The trajectories of the system (2) through points and must intersect the impulse set at some points and , respectively, and after impulsive effect, the points and must be mapped to the phase set at some points and , respectively (see Figure 3). Because distinct trajectories do not intersect, we can easily have and ; then we get the order one successor functions must satisfy
which means the order one successor function is monotone decreasing in the line segment , thus there exists only one point such that .

Besides, for any point , it is easy to know that . So the system (2) has only one order one periodic solution in the region . That completes the proof.

Theorem 17. *Under the conditions of Theorem 16, the order one periodic solution of the system (2) is orbitally asymptotically stable.*

*Proof. *According to Theorem 16, the system (2) has a unique order one periodic solution that passes through the point which is in the phase set , where . The trajectory of the system (2) through point intersects the impulse set at point , and after impulsive effect, the point is mapped to point , where . Besides, the trajectory of the system (2) from point must intersect the impulse set again at a point , and after impulsive effect, the point is mapped to a point in the phase set , where (see Figure 4).

Repeat the above steps, the trajectory from point will come across impulsive effect infinite times. Denote the phase point corresponding to the th impulsive effect by , . Obviously, is also the order successor point of point . Then we have
Thus , is a monotonically increasing sequence, and , is a monotonically decreasing sequence. Because the order one successor function is monotone decreasing in the line segment , we have
and furthermore,
Pick any point different from the point ; without loss of generality, we assume (otherwise, ; the discussions are similar), there must exist a natural number such that . The trajectory from point will also undergo impulsive effect infinite times, and we denote the phase point corresponding to the th impulsive effect as , . For any natural number , we have and . Then , is also a monotonically increasing sequence, , is also a monotonically decreasing sequence, and we get
This indicates that the trajectory of the system (2) from point ultimately tends to the trajectory passing through point , which means the order one periodic solution is orbitally asymptotically stable. That completes the proof.

Theorem 18. *Under the conditions of Theorem 16, the system (2) has no order periodic solution in region , where .*

*Proof. *For any point , , the trajectory from point will undergo impulsive effect infinite times. Denote the impulsive point and phase point corresponding to the impulsive effect by and , respectively. Similar to the discussions in Theorem 17, we can easily know , , then we have . Obviously, point is the order successor point of point , then the order successor function . By Lemma 13, there does not exist an order periodic solution through point .

For any point , according to the proof of Theorem 17, there must exist a natural number such that . Denote the order successor point of point by point , then we have and . Then , is amonotonically increasing sequence where , , is a monotonically decreasing sequence, and
When , we have and the order successor function . When , we have and . That is to say there does not exist an order periodic solution passing through point .

Analogously, for any point , we can prove there does not exist an order periodic solution passing through point .

From the above discussion, we know that the order one periodic solution is the unique periodic solution of the system (2) and there is no order periodic solution in region . That completes the proof.

Theorem 19. *If , (or if , ), and , then there must exist two fixed values and which satisfy such that*(1)*when , the system (2) has an order one homoclinic cycle;*(2)*if , then the homoclinic cycle of system (2) disappears and bifurcates an order one periodic solution in region , where region is the region enclosed by the -axis, the impulse set , and the unstable flow of the saddle . Furthermore, the order one periodic solution is unique and is orbitally asymptotically stable;*(3)*if , then the homoclinic cycle of system (2) also disappears and there is no periodic solution in region .*

*Proof. *For , it is easy to know that . According to Lemma 15, the system (1) has four equilibria: two hyperbolic saddles and , a hyperbolic stable node , and a node or focus . Obviously, the -axis, the vertical isocline , and the unstable and stable flow of are all intersected with the impulse set and phase set . For convenience, we denote the intersections by points and , and , and , and and , respectively (see Figure 5). Then the region is the interior of the closed curve , where and denote the equilibria and , respectively. According to the impulsive conditions of the system (2), there must exist a fixed value , when , point is mapped to the point after impulsive effect. Also, there must exist a fixed value , when , point is mapped to the point after impulsive effect. Thus we have and .

When , we have ; we know the curve is an order one circle which has the saddle in it; that is to say, the system (2) has an order one homoclinic cycle.

If , after impulsive effect, point is mapped to a point between points and . Similar to the discussion of Theorems 16 and 17, we can prove that the system (2) has an order one periodic solution which passes through a point and is orbitally asymptotically stable. Besides, for any point , the trajectory of the system (2) from point must ultimately tend to the equilibrium and does not undergo any impulsive effect, that is to say the system (2) has no periodic solution passing through the point ; so we can know the system (2) has a unique order one periodic solution in the region which is orbitally asymptotically stable.

When , we have . Then the trajectory of the system (2) through point must tend to the equilibrium after coming across once impulsive effect and the trajectory passing through point must tend to the equilibrium without undergoing any impulsive effect. So we get that the system (2) has no periodic solution in region . That completes the proof.

*3.2. Homoclinic Bifurcation of System (3) about Parameter **δ*

*δ*

*By above analysis, we know there exists a such that system (2) has an order one homoclinic cycle and when the parameter is appropriately changed, the homoclinic cycle disappears and bifurcates a unique stable order one periodic solution. In the following, we will choose as a control parameter and study homoclinic bifurcation of system (3) by using the theory of rotated vector fields. For the sake of convenience, we give some properties about rotated vector fields.*

*The perturbed system of system (4) with parameter is as follows:
Let
then we have the following definitions.*

*Definition 20 (see [18]). *For any point on the trajectory of system (15), if , then system (15) constitutes positive rotated vector fields concerning the parameter ; otherwise, if , system (15) constitutes negative rotated vector fields.

*Lemma 21 (see [18]). In the positive (negative) rotated vector fields of system (15), the rotated direction of vector fields is counterclockwise (clockwise) when parameter changes from to .*

*Theorem 22. If , (or if , ), , , , , and , then the system (3) has a unique order one periodic solution. Furthermore, it is orbitally asymptotically stable.*

*Proof. *From the discussion in Theorem 19, we know that when and , the system (3) has an order one homoclinic circle .

For the system (3), we denote and . By simple calculation, we have
Divided by the vertical isocline , the rotated direction of vector fields below the vertical isocline is counterclockwise, but above the vertical isocline, the rotated direction of vector fields of system (3) is clockwise. When the parameter changes from 0 (the value in which heteroclinic cycle exists) to , both the unstable and stable flow of move away from their original positions. We denote the new intersections of the unstable and stable flow of and the impulse set and phase set are and , and and , respectively (see Figure 6). According to Theorem 19, we have ; then there must exist a point which is in the impulse set and below the point such that . Obviously, system (3) must have a trajectory passing through the point , and we suppose it intersects with the phase set at point .

Obviously, point is the order one successor point of point , then the order one successor function . Besides, select a point next to in the phase set ; the trajectory of the system (3) from point must intersect the impulse set at a point which is next to point and above point , and after impulsive effect, point is mapped to a point which is the order one successor point of , then we have . That is to say point is above point . Then we have the successor function . By Lemma 12, we know that there must exist a point in the phase set which is between the points and such that . Then we know the system (3) has an order one periodic solution which passes through the point .

Similar to the discussion of Theorems 16 and 17, we can prove the uniqueness and stability of the order one periodic solution. That completes the proof.

*4. Numerical Simulations and Discussions*

*In this paper, we have proposed two predator-prey models with nonmonotonic functional response and state dependent impulsive harvesting. If we take no account of the predator harvesting, Ruan and Xiao (2001) showed that the system (1) has no periodic solution under the conditions listing in Theorem 16. Our numerical simulations also confirm such conclusions (see Figure 7). However, under the same conditions, if we further consider an impulsive predator harvesting, we get that the system (2) has a periodic solution in some cases. We also showed that the existence is mainly dependent on the harvesting rate . We found that there exists a threshold , and the system must have a unique orbitally asymptotically stable order one periodic solution when the harvesting rate is higher than . When the periodic solution exists, our numerical simulations show that, with a bigger harvesting rate, the solution has a longer period and a larger amplitude (see Figure 8).*

*Besides, we also prove that the system (2) does not have order , periodic solutions in Theorem 18; then the order one periodic solution is the unique periodic solution. These results illustrate that if we can satisfy the conditions listing in Theorem 16, we can not only form a sustainable ecological system but also obtain large amounts of valuable biological resources.*

*Under the conditions listing in Theorem 19, the system (1) also does not have any periodic solution (see Figure 9) and we also give the conditions under which the system (2) has an order one periodic solution. Besides, according to the conclusions of Theorem 19, we can choose the parameter as a bifurcation parameter such that the impulsive differential equation system (2) exhibits the phenomenon of homoclinic bifurcation. Similar to many ordinary differential equation systems, we can see that there exists a bifurcation point for the system (2). When , the system (2) has an order one homoclinic cycle which is the unique order one circle in region . When gradually changes from to , the order one homoclinic cycle is broken and a new order one periodic solution which is orbitally asymptotically stable is generated (see Figure 10). When gradually changes from to , the order one homoclinic cycle is also broken but no new order one periodic solution is generated at the same time; that is to say, the system (2) will have no periodic solution (see Figure 11). These results illustrate that it demands reasonable control of the harvest yield to form a good ecological environment.*

*Conflict of Interests*

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

*Acknowledgments*

*This work is supported by the National Natural Science Foundation of China (11171284, 11371306, and 11301453), the Innovation Scientists and Technicians Troop Construction Projects of Henan Province (104200510011), Program for Innovative Research Team (in Science and Technology) in University of Henan Province (2010IRTSTHN006), and Sci-Tech Program Project of Henan Province (122300410034).*

*References*

- J. S. Tener,
*Muskoxen*, Queens Printer, Ottawa, Canada, 1965. - J. C. Holmes and W. M. Bethel, “Modification of intermediate host behavior by parasites,”
*Zoological Journal of the Linnean Society*, vol. 51, pp. 123–149, 1972. View at Google Scholar - G. S. K. Wolkowicz, “Bifurcation analysis of a predator-prey system involving group defence,”
*SIAM Journal on Applied Mathematics*, vol. 48, no. 3, pp. 592–606, 1988. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Ruan and D. Xiao, “Global analysis in a predator-prey system with nonmonotonic functional response,”
*SIAM Journal on Applied Mathematics*, vol. 61, no. 4, pp. 1445–1472, 2001. View at Publisher · View at Google Scholar · View at MathSciNet - H. I. Freedman and G. S. K. Wolkowicz, “Predator-prey systems with group defence: the paradox of enrichment revisited,”
*Bulletin of Mathematical Biology*, vol. 48, no. 5-6, pp. 493–508, 1986. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. Mischaikow and G. Wolkowicz, “A predator-prey system involving group defense: a connection matrix approach,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 14, no. 11, pp. 955–969, 1990. View at Publisher · View at Google Scholar · View at MathSciNet - J. B. Collings, “The effects of the functional response on the bifurcation behavior of a mite predator-prey interaction model,”
*Journal of Mathematical Biology*, vol. 36, no. 2, pp. 149–168, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. H. Edwards, “The influence of high substrate concentrations on microbial kinetics,”
*Biotechnology and Bioengineering*, vol. 12, no. 5, pp. 679–712, 1970. View at Publisher · View at Google Scholar · View at Scopus - L. S. Chen, “Pest control and geometric theory of Semi-continuous dynamical system,”
*Journal of Beihua University*, vol. 12, pp. 1–9, 2011. View at Google Scholar - P. Yongzhen, L. Changguo, and C. Lansun, “Continuous and impulsive harvesting strategies in a stage-structured predator-prey model with time delay,”
*Mathematics and Computers in Simulation*, vol. 79, no. 10, pp. 2994–3008, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. Liu and R. Tan, “Impulsive harvesting and stocking in a Monod-Haldane functional response predator-prey system,”
*Chaos, Solitons and Fractals*, vol. 34, no. 2, pp. 454–464, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. Negi and S. Gakkhar, “Dynamics in a Beddington-DeAngelis prey-predator system with impulsive harvesting,”
*Ecological Modelling*, vol. 206, no. 3-4, pp. 421–430, 2007. View at Publisher · View at Google Scholar · View at Scopus - G. Zeng, L. Chen, and L. Sun, “Existence of periodic solution of order one of planar impulsive autonomous system,”
*Journal of Computational and Applied Mathematics*, vol. 186, no. 2, pp. 466–481, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. Jiang, Q. Lu, and L. Qian, “Complex dynamics of a Holling type II prey-predator system with state feedback control,”
*Chaos, Solitons & Fractals*, vol. 31, no. 2, pp. 448–461, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. Nie, J. Peng, Z. Teng, and L. Hu, “Existence and stability of periodic solution of a Lotka-Volterra predator-prey model with state dependent impulsive effects,”
*Journal of Computational and Applied Mathematics*, vol. 224, no. 2, pp. 544–555, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Huang, S. Liu, X. Song, and L. Chen, “Periodic solutions and homoclinic bifurcation of a predator-prey system with two types of harvesting,”
*Nonlinear Dynamics*, vol. 73, no. 1-2, pp. 815–826, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. Dai, M. Zhao, and L. Chen, “Homoclinic bifurcation in semi-continuous dynamic systems,”
*International Journal of Biomathematics*, vol. 5, no. 6, Article ID 1250059, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - Y. Q. Ye,
*Limit Cycle Theory*, Shanghai Science and Technology Press, Shanghai, China, 1984.

*
*