- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- 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

Discrete Dynamics in Nature and Society

Volume 2013 (2013), Article ID 358930, 15 pages

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

## Complex Dynamical Behaviors in a Predator-Prey System with Generalized Group Defense and Impulsive Control Strategy

Department of Mathematics, Qiannan Normal College for Nationalities, Duyun, Guizhou 558000, China

Received 1 November 2012; Accepted 20 May 2013

Academic Editor: Patricia J. Y. Wong

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

#### Abstract

A predator-prey system with generalized group defense and impulsive control strategy is investigated. By using Floquet theorem and small amplitude perturbation skills, a local asymptotically stable prey-eradication periodic solution is obtained when the impulsive period is less than some critical value. Otherwise, the system is permanent if the impulsive period is larger than the critical value. By using bifurcation theory, we show the existence and stability of positive periodic solution when the pest eradication lost its stability. Numerical examples show that the system considered has more complicated dynamics, including (1) high-order quasiperiodic and periodic oscillation, (2) period-doubling and halving bifurcation, (3) nonunique dynamics (meaning that several attractors coexist), and (4) chaos and attractor crisis. Further, the importance of the impulsive period, the released amount of mature predators and the degree of group defense effect are discussed. Finally, the biological implications of the results and the impulsive control strategy are discussed.

#### 1. Introduction

In population dynamics, a functional response of the predator to the prey density refers to the change in the density of prey attached per unit time per predator as the prey density changes and it is assumed to be monotonically increasing in most predator-prey systems. For example, Holling type I, II, and III functional response [1] and the sigmoidal type response function [2] and Ivlev type response function [3] The previous functional responses are prey dependent. But, both predator and prey densities have an effect on the response, such as Beddington-DeAngelis functional response [4, 5] and modified Holling type II and type III response functions [6]

However, some experimental and observational evidence shown that the functional response is not always monotonically increasing, such as Holling type IV [7] and [8]. Group defense is a term used to describe the phenomenon whereby predation is decreased, or even prevented altogether, due to the increased ability of the prey to better defend or disguise itself when it exists in enough large numbers [9–11]. The buffalo group defense was modeled using a generalized group defense in [12], where is a positive integer whose value determines the degree of antipredator behavior and group defense.

Recently, it is of great interest to investigate complex dynamics for impulsive perturbations in populations dynamics. In particular, the impulsive prey-predator population models have been investigated by many researchers. The results of studies of the dynamics of a predator-prey model with nonmonotonic functional response, such as Holling type IV functional response with respect to an impulsive control strategy, were presented in [13–24]. To the best of our knowledge, there are few papers studying the group defense predator-prey with impulsive effect, where the antipredator behavior and group defense effect described by nonmonotonic functional response. Zhang et al. [25] considered a predator-prey system with defensive ability of prey by Holling type IV functional response and impulsive perturbations on the predator: The conditions for the local asymptotically stable prey-eradication periodic solution and permanence of the system are obtained; a series of complex phenomena are displayed by numerical simulation. Furthermore, based on this work, Pei et al. [26] investigated a one-prey multi-predator model with defensive ability of the prey by introducing impulsive biological control strategy: And it shown that the multi-predator impulsive control strategy is more effective than the classical one and makes the dynamical behaviors of the system more complex. Recently, a predator-prey system with impulsive effect and group defense with the nonmonotone function was studied by Li et al. [27], They proved that there exists a locally stable pest-eradication periodic solution when the impulsive period is less than certain critical values; otherwise, the system is permanent. Some complicated dynamics, such as quasiperiodic oscillation, bifurcation, and attractor crisis, were shown by numerical simulations.

In this paper, we study a predator-prey system with impulsive effect and generalized group defense with the nonmonotone function : where and represent the prey and the predator populations at time , respectively; are positive. is the intrinsic rate of increase of the prey and is the death rate of the predator, is the carrying capacity of the prey, is the degree of anti-predator behavior and group defense, and is the rate of conversing prey into predator. , is the periodic of the impulse for predator in order to eradicate target pests, protect nontarget pest (or harmless insect) from extinction and drive target pest to extinction, or control target pest at acceptably low level to prevent an increasing pest population from causing an economic loss. is the proportionality constant which represents the rate of mortality due to the applied pesticide; for example, impulsive reduction of the population is possible by harvesting or by poisoning with chemicals used in agriculture. is the number of predators released each time, for example, by artificial breeding of the species or release of some species.

The paper is arranged as follows. In Section 2, some notations and Lemmas are given. In Section 3, using the Floquet theory of impulsive equation and small amplitude perturbation skills, we will prove the local stability of prey-eradication periodic solution when the impulsive period is less than some critical value and give the condition of permanence. In Section 4, by using bifurcation theory, the existence and stability of positive periodic solution are studied when is close to the critical value . In Section 5, the results of numerical examples are shown, and some rich dynamic behaviors are obtained; the effects of the impulsive period, the released amount of mature predators and the coefficient of group defense effect are discussed. Finally, the conclusions are discussed briefly in Section 6.

#### 2. Preliminaries

In this section, we will give some definitions, notations, and lemmas which will be useful for our main results.

Let . Denote by the map defined by the right hand of the first two equations of system (11), and denote by the set of all nonnegative integers. Let , then is said to belong to class if(1)is continuous in and for each , exists,(2) is locally Lipschitzian in .

*Definition 1. *Let ; then for , the upper right derivative of with respect to the impulsive differential system (11) is defined as

*Definition 2. *System (11) is said to be permanent if there exist two positive constants and such that each positive solution of the system (11) satisfies , for all .

The solution of system (11) is a piecewise continuous function , is continuous on , and exists; the smoothness properties of guarantee the global existence and uniqueness of solutions of system (11); for details see [28, 29]. The following lemma is obvious.

Lemma 3. * Let be a solution of system (11) with ; then for all and further for all if .*

And we will use the following important comparison theorem on impulsive differential equation [29].

Lemma 4. *Suppose . Assume that
**
where is continuous in , and for , exists; is nondecreasing. Let be the maximal solution of the scalar impulsive differential equation
**
existing on . Then implies that where is any solution of system (11).*

Finally, we give some basic properties about the following subsystem of system (11): Clearly, when , is a positive periodic solution of system (15). Since is the solution of system (15) with initial value , where ; then one can get the following.

Lemma 5. * Let be a positive periodic solution of system; (15) and every solution of system (15) with , one has , when .**Therefore, one obtains the pest-eradication periodic solution
**
for .*

#### 3. Extinction and Permanence

Firstly, we study the stability of prey-eradication periodic solution.

Theorem 6. * Let be any solution of system (11); then is locally asymptotically stable provided that
*

*Proof. *The local stability of periodic solution may be determined by considering the behavior of small amplitude perturbations of the solution. Consider
There may be written
where satisfies
and , the identity matrix. The linearization of the third and fourth equations of system (11) becomes
Hence, if both eigenvalues of
have absolute values less than one, then the periodic solution is locally stable. Since all eigenvalues of are
if and only if
According to Floquet theory [28] of impulsive differential equation, the prey-eradication solution is locally stable. This completes the proof.

Theorem 7. *There exists a constant , such that for each solution of system (11) with all being large enough.*

*Proof. *Let . It is clear that . We calculate the upper right derivative of along a solution of system (11) and get the following impulsive differential equation:
Let ; then is bounded. Select and such that
where and are two positive constants. According to Lemma 4, we have
where . Hence
Therefore, is ultimately bounded. We obtain that each positive solution of system (11) is uniformly ultimately bounded. This completes the proof.

In the following, we investigate the permanence of system (11).

Theorem 8. *System (11) is permanent if
*

*Proof. *Suppose is a solution of system (11) with . From Theorem 7 we may assume that , and . Let
According to Lemmas 4 and 5, we have for all large enough. In the following, we want to find such that for all large enough. We will do it in the following two steps for convenience.*Step *1. Since
we can select small enough such that , , and

We will prove there exists such that . Otherwise, according to the above assumption, we get , and by Lemmas 4 and 5, we have and , where
and is the solution of the following equation:
Therefore, there exists a such that
Let and let . We can get
Integrating (38) on , we have
Then as , which is a contradiction to the boundedness of . Hence there exists a such that .*Step *2. If for all , then our aim is obtained.

Hence we only need to consider those solutions which leave the region and reenter again. Let . Then is impulsive point or nonimpulsive point.*Case 1 . * If is impulsive point, there exist a such that . Then for and
Choose such that
where . Let . Then, there exists a such that . Otherwise , . Consider (36) with ; we have
where . Then
for , which implies (39) holds for . As in Step 1, we have
The first and third equations of system (11) given
Integrating the above equation on , we can get
and thus
a contradiction.

Let ; then . For , we have For , the same arguments can be continued since .

*Case 2*If is nonimpulsive point, then for and ; suppose . There are two possible cases for .

*.**Case 2.1.*for all . As in step 1, we can prove that there must be a such that ; Let , then and . For , we have For , the same arguments can be continued since .

*Case 2.2.*There exists a such that . Let ; then for and . For , (45) holds true; integrating (45) on , we have Since for , the same arguments can be continued.

Hence for all . The proof is completed.

*Remark 9. *Let
Since , as , and
so has a unique positive root, denoted by . From Theorems 6 and 8 we know that is a threshold. If , then pest-eradication periodic solution is asymptotically stable; if , then system (11) is permanent.

*Remark 10. * If ; that is, there are without taking any pest-management strategy, large numbers of preys (pest) would coexisting with predators (natural enemy). If , that is, there is periodic spraying pesticide (or harvesting) only. Thus, we can easily obtain that is the threshold. If ; that is, there is periodic releasing of predator (natural enemy) only, without periodic spraying pesticide (or harvesting). We can easily get that is the threshold. Comparing with the classic methods (such as biological control or chemical control), the integrated pest management (IPM) is a better one, since and . Some numerical examples will be given in Section 5.

#### 4. Bifurcation and Existence of Positive Periodic Solution

In this section, we deal with the existence of a nontrivial periodic solution to system (11) near the prey-eradication periodic solution via bifurcation.

Let ; system (11) becomes as follows: All notations used in this section are the same as those in [30]. Let be the flow associated to (53); we have , where , where is the root of , where .

Lemma 11 (see [30]). * If and , then one has the following.*(a)*If , then one has a bifurcation. Moreover, one has a bifurcation of a nontrivial periodic solution of (53) if and a subcritical case if .*(b)*If , then one has an undetermined case.*

In order to apply Lemma 11, we compute the following: If , this corresponds to satisfying Further, we can get Note that then Since it is easy to verify that and In order to determine the sign of , let We have Thus, we can conclude that , since and is strictly increasing. Therefore, we have . In view of and according to Lemma 11, we obtain the following result.

Theorem 12. *System (11) has a positive periodic solution if and is close to , where satisfies
**
and the nontrivial periodic solution is supercritical case via bifurcation, which means that the positive periodic solution is stable.*

#### 5. Numerical Analysis

In this section, we will study the impulsive effect on system (11) and show that the impulsive perturbations cause complicated dynamical behavior for system (11). The influence of , and may be documented by stroboscopically sampling one of the variables over a range of their values. Stroboscopic map is a special case of the Poincaré map for periodically forced system or periodically pulsed system. Fixing points of the stroboscopic map correspond to periodic solutions of system (11) having the same period as the pulsing term; periodic points of period about stroboscopic map correspond to entrained periodic solutions of system (11) having exactly times the period of the pulsing; invariant circles correspond to quasi-periodic solutions of system (11); system (11) possibly appear chaotic (or strange) attractors.

*Example 13. *Let , , ,,,,,,,and with initial value .

From Remark 10, large numbers of preys (pest) could coexist with predators (natural enemy) with periodic oscillations, if we are not taking pest-management strategy () (Figure 1(a)). ; that is, there are periodic spraying pesticide (or harvesting) only, without releasing of predator (natural enemy), large numbers of preys (pest) coexist with periodic oscillation, but predators (natural enemy) rapidly decrease to zero when (Figure 1(b)). ; that is, there is periodic releasing of predator (natural enemy) only; without spraying pesticide (or harvesting), a few of preys (pest) coexist with predators (natural enemy) when (Figure 1(c)). We cannot make the prey population eradicate when . From Theorem 6, we know that the prey-eradication periodic solution is asymptotically stable provided that . A typical prey-eradication periodic solution of the system (11) is shown in Figure 1(d), where we observe how the variable oscillates in a stable cycle. In contrast, the prey population rapidly decreases to zero when . Hence, the integrated pest management (IPM) is better than the classic methods (such as biological control or chemical control).

According to Theorem 12, if the impulsive periodic and is close to , the prey eradication solution becomes unstable, there is a supercritical bifurcation, then the prey and predator can coexist on a stable positive periodic solution when (Figure 2). Therefore, in order to control the pest populations, we would choose an appropriate impulsive periodic . and close to would be a better one.

Let and fix other parameter sets of values; we have displayed bifurcation diagrams for the pest population and the predator population for impulsive period over and . We find that by increasing the impulsive period , system (11) undergoes a process of period-doubling cascade chaos crisis and high-order periodic oscillations (Figure 3). When increases from to , there is a cascade of period-doubling bifurcations leading to chaos (Figure 4). When , the chaos suddenly disappears and a -periodic solution appears, then the -periodic solution abruptly disappears and the chaos abruptly appears again when , these constituting several types of crises (Figure 5). However, when and , it appears that attractors are nonunique, coexistence of stranger attractor with -periodic solution (Figure 6). Obviously, which one of the attractors is reached depends on the initial values.

*Example 14. * Let ,,,,,,,,, and with initial value . We investigate the effect of on the system (11). Figure 7 showed bifurcation diagrams obtained by stroboscopically sampling the pest population and the predator population for over . The resulting bifurcation diagrams clearly showed that system (11) has rich dynamics, including period-doubling bifurcation, period-halving bifurcation, and chaos. When increases from to , there is a period-halving bifurcation leading to a -periodic solution (Figure 8).

*Example 15. * Let ,,,,,,,,,,and with initial value . We consider the effect of on the system (11). The resulting bifurcation diagrams (Figure 9), the pest population , and the predator population for over clearly showed that system (11) has complex dynamics, such as period-doubling bifurcation, high-order periodic oscillation, and chaos. In Figure 10, the typical high-order oscillation of system (11) is shown: periodic solutions when , respectively. Further, Figure 11 shown the maximin and mean amount of prey population and predator population of system (11) with over .

From bifurcation diagrams in Figures 3, 7, and 9, we can easily see that the dynamical behavior of these three cases is very complicated, which includes (1) high-order quasi-periodic and periodic oscillations, (2) period-doubling bifurcation, (3) period-halving bifurcations, (4) nonunique dynamics (meaning that several attractors coexist), and (5) cries (the phenomenon of “crisis” in chaotic attractors can suddenly appear or disappear, or change size discontinuously as a parameter smoothly varies).

#### 6. Conclusion

In this paper, we have investigated a predator-prey system with generalized group defense and concerning impulsive control strategy for pest control in detail. We have shown that there exists an asymptotically stable pest-eradication periodic solution if the impulsive period is less than the critical value . If we choose our impulsive control strategy, in order to drive the pest to extinction, we can determine the impulsive period according to the effect of the chemical pesticides on the populations and the cost of releasing natural enemies such that .

But, in a real world, complete eradication of pest populations is generally not possible, nor is it biologically or economically desirable. A good-pest control program should reduce pest population to levels acceptable to the public. When , the stability of the pest-eradication periodic solution is lost, system (11) is permanent, and there exists a nontrivial periodic solution when is close to . The smaller the period, the fewer the pests. Therefore, we can control the pest population below some economic threshold ( is defined as the pest population level that produces damage equal to the costs of preventing damage) by choosing appropriate impulsive period and the number of mature predator released , according to the degree of antipredator behavior and group defense , making an integrated pest-management strategy every period . Then, the periodic releasing of natural enemies and spraying pesticides change the properties of the system without impulses and our results suggest an effective approach in the pest control.

Numerical results show that system (11) can take on various kinds of periodic fluctuations and several types of attractor coexistence and is dominated by high-order periodic oscillations, quasi-periodic oscillations, and chaotic oscillations. These results imply that the presence of pulses destroys equilibria, initiates multiple attractors, quasi-periodic oscillations, and chaos, and makes the dynamical behaviors more complex.

#### Acknowledgment

This work was supported by the Science and Research Project Foundation of Educational Department of Guizhou Province (no. 2010096).

#### References

- C. S. Holling, “The functional response of predator to prey density and its role in mimicry
and population regulation,”
*Memoirs of the Entomological Society of Canada*, vol. 1, no. 45, pp. 1–60, 1965. View at Google Scholar - H. I. Freedman,
*Deterministic Mathematical Models in Population Ecology*, vol. 57 of*Monograph Textbooks Pure and Applied Mathematics*, Marcel Dekker, New York, NY, USA, 1980. View at Zentralblatt MATH · View at MathSciNet - V. S. Ivlev,
*Experimental Ecology of the Feeding of Fishes*, Yale University Press, New Haven, Conn, USA, 1961. - D. L. DeAngelis, R. A. Goldstein, and R. V. O'Neill, “A Model for trophic interaction,”
*Ecology*, vol. 56, pp. 881–892, 1975. View at Google Scholar - J. R. Beddington, “Mutual interference between parasites or predators and its effect on searching efficiency,”
*Journal of Animal Ecology*, vol. 44, pp. 331–340, 1975. View at Google Scholar - D. A. Miller, J. B. Grand, T. F. Fondell, and M. Anthony, “Predator functional response and prey survival: direct and indirect interactions affecting a marked prey population,”
*Journal of Animal Ecology*, vol. 75, no. 1, pp. 101–110, 2006. View at Publisher · View at Google Scholar · View at Scopus - J. F. Andrews, “A mathematical model for the continuous culture of microorganisms
utilizing inhibitory substrates,”
*Biotechnology and Bioengineering*, vol. 10, pp. 707–723, 1968. View at Google Scholar - H. I. Freedman and S. G. Ruan, “Hopf bifurcation in three-species food chain models with group defense,”
*Mathematical Biosciences*, vol. 111, no. 1, pp. 73–87, 1992. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · 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 - H. I. Freedman and H. S. Quan, “Interactions leading to persistence in predator-prey systems with group defence,”
*Bulletin of Mathematical Biology*, vol. 50, no. 5, pp. 517–530, 1988. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. M. Xiao and S. G. Ruan, “Codimension two bifurcations in a predator-prey system with group defense,”
*International Journal of Bifurcation and Chaos in Applied Sciences and Engineering*, vol. 11, no. 8, pp. 2123–2131, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Ddumba, J. Y. T. Mugisha, J. W. Gonsalves, and G. I. H. Kerley, “The role of predator fertility and prey threshold bounds on the global and local dynamics of a predator-prey model with a prey out-flux dilution effect,”
*Applied Mathematics and Computation*, vol. 218, no. 18, pp. 9169–9186, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. W. Zhang, F. Y. Wang, and L. S. Chen, “A food chain model with impulsive perturbations and Holling IV functional response,”
*Chaos, Solitons and Fractals*, vol. 26, no. 3, pp. 855–866, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. L. Xiong, Y. Xue, and S. Y. Li, “A food chain system with Holling IV functional responses and impulsive effect,”
*International Journal of Biomathematics*, vol. 1, no. 3, pp. 361–375, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Baek, “A food chain system with Holling type IV functional response and impulsive perturbations,”
*Computers & Mathematics with Applications*, vol. 60, no. 5, pp. 1152–1163, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. J. Jiao, X. Z. Meng, and L. S. Chen, “A stage-structured Holling mass defence predator-prey model with impulsive perturbations on predators,”
*Applied Mathematics and Computation*, vol. 189, no. 2, pp. 1448–1458, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. W. Zhang, D. J. Tan, and L. S. Chen, “Chaos in periodically forced Holling type IV predator-prey system with impulsive perturbations,”
*Chaos, Solitons and Fractals*, vol. 27, no. 4, pp. 980–990, 2006. View at Publisher · View at Google Scholar · View at MathSciNet - Z. J. Liu and R. H. 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 - Q. Wang, B. X. Dai, and Y. M. Chen, “Multiple periodic solutions of an impulsive predator-prey model with Holling-type IV functional response,”
*Mathematical and Computer Modelling*, vol. 49, no. 9-10, pp. 1829–1836, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Baek and Y. Do, “Stability for a Holling type IV food chain system with impulsive perturbations,”
*Kyungpook Mathematical Journal*, vol. 48, no. 3, pp. 515–527, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C.-Y. Huang, Y.-J. Li, and H.-F. Huo, “The dynamics of a stage-structured predator-prey system with impulsive effect and Holling mass defence,”
*Applied Mathematical Modelling*, vol. 36, no. 1, pp. 87–96, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Z. Wang and M. Zhao, “Dynamic analysis of an impulsively controlled predator-prey model with Holling type IV functional response,”
*Discrete Dynamics in Nature and Society*, vol. 2012, Article ID 141272, 18 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H.-F. Huo and W.-T. Li, “Dynamics of a nonautonomous semiratio-dependent predator-prey system with nonmonotonic functional responses,”
*Discrete Dynamics in Nature and Society*, vol. 2006, Article ID 70656, 19 pages, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X. P. Li and W. S. Yang, “Permanence of a semi-ratio-dependent predator-prey system with nonmonotonic functional response and time delay,”
*Abstract and Applied Analysis*, vol. 2009, Article ID 960823, 6 pages, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. W. Zhang, L. Z. Dong, and L. S. Chen, “The study of predator-prey system with defensive ability of prey and impulsive perturbations on the predator,”
*Chaos, Solitons & Fractals*, vol. 23, no. 2, pp. 631–643, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Z. Pei, C. G. Li, L. S. Chen, and C. Wang, “Complex dynamics of one-prey multi-predator system with defensive ability of prey and impulsive biological control on predators,”
*Advances in Complex Systems*, vol. 8, no. 4, pp. 483–495, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Y. Li, Z. L. Xiong, and X. Wang, “The study of a predator-prey system with group defense and impulsive control strategy,”
*Applied Mathematical Modelling*, vol. 34, no. 9, pp. 2546–2561, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. D. Bainov and D. D. Simeonov,
*Impulsive Differential Equations: Periodic Solutions and Applications*, Longman UK, 1993. - V. Lakshmikantham, D. D. Baĭnov, and P. S. Simeonov,
*Theory of Impulsive Differential Equations*, vol. 6, World Scientific Publishing, Teaneck, NJ, USA, 1989. View at MathSciNet - A. Lakmeche and O. Arino, “Bifurcation of non trivial periodic solutions of impulsive differential equations arising chemotherapeutic treatment,”
*Dynamics of Continuous, Discrete and Impulsive Systems*, vol. 7, no. 2, pp. 265–287, 2000. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet