Nonstandard finite-difference methods for predator-prey models with general functional response

Predator-prey systems with linear and logistic intrinsic growth rate of the prey are analyzed. The models incorporate the mutual interference between predators into the functional response which stabilizes predator-prey interactions in the system. Positive and elementary stable nonstandard (PESN) finite-difference methods, having the same qualitative features as the corresponding continuous predator-prey models, are formulated and analyzed. The proposed numerical techniques are based on a nonlocal modeling of the growth-rate function and a nonstandard discretization of the time derivative. This discretization approach leads to significant qualitative improvements in the behavior of the numerical solution. In addition, it allows for the use of an essentially implicit method for the cost of an explicit method. Applications of the PESN methods to specific predator-prey systems are also presented.

[1]  Murray E. Alexander,et al.  A non-standard numerical scheme for a generalized Gause-type predator–prey model , 2004 .

[2]  J. Beddington,et al.  Mutual Interference Between Parasites or Predators and its Effect on Searching Efficiency , 1975 .

[3]  Robert M. May,et al.  Stability and Complexity in Model Ecosystems , 2019, IEEE Transactions on Systems, Man, and Cybernetics.

[4]  Hristo V. Kojouharov,et al.  NONSTANDARD NUMERICAL METHODS FOR A CLASS OF PREDATOR-PREY MODELS WITH PREDATOR INTERFERENCE , 2007 .

[5]  R. Mickens.  Nonstandard Finite Difference Models of Differential Equations , 1993 .

[6]  Eric Deleersnijder,et al.  A high-order conservative Patankar-type discretisation for stiff systems of production-destruction equations , 2003 .

[7]  F. Brauer,et al.  Mathematical Models in Population Biology and Epidemiology , 2001 .

[8]  Donald L. DeAngelis,et al.  A MODEL FOR TROPHIC INTERACTION , 1975 .

[9]  Hristo V. Kojouharov,et al.  Analysis and numerical simulation of phytoplankton-nutrient systems with nutrient loss , 2005, Math. Comput. Simul..

[10]  C. S. Holling,et al.  The functional response of predators to prey density and its role in mimicry and population regulation. , 1965 .

[11]  S. Hsu,et al.  Global analysis of the Michaelis–Menten-type ratio-dependent predator-prey system , 2001, Journal of mathematical biology.

[12]  J. F. Gilliam,et al.  FUNCTIONAL RESPONSES WITH PREDATOR INTERFERENCE: VIABLE ALTERNATIVES TO THE HOLLING TYPE II MODEL , 2001 .

[13]  Edward H. Twizell,et al.  An unconditionally convergent discretization of the SEIR model , 2002, Math. Comput. Simul..

[14]  Yang Kuang,et al.  Global qualitative analysis of a ratio-dependent predator–prey system , 1998 .

[15]  M. Doebeli.  Genetic Variation and Persistence of Predator-prey Interactions in the Nicholson–Bailey Model , 1997 .

[16]  Roumen Anguelov,et al.  Contributions to the mathematics of the nonstandard finite difference method and applications , 2001 .

[17]  R. May,et al.  Stability and Complexity in Model Ecosystems , 1976, IEEE Transactions on Systems, Man, and Cybernetics.

[18]  Edward H. Twizell,et al.  An unconditionally convergent finite-difference scheme for the SIR model , 2003, Appl. Math. Comput..

[19]  D. Dimitrov,et al.  Stability-preserving finite-difference methods for general multi-dimensional autonomous dynamical systems , 2007 .

[20]  Hristo V. Kojouharov,et al.  Complete mathematical analysis of predator-prey models with linear prey growth and Beddington-DeAngelis functional response , 2005, Appl. Math. Comput..

[21]  Tzy-Wei Hwang,et al.  Global analysis of the predator–prey system with Beddington–DeAngelis functional response , 2003 .

[22]  Alet Roux,et al.  An Improved Theta-method for Systems of Ordinary Differential Equations , 2003 .

[23]  Hristo V. Kojouharov,et al.  Positive and elementary stable nonstandard numerical methods with applications to predator-prey models , 2006 .

[24]  Hristo V. Kojouharov,et al.  Nonstandard finite-difference schemes for general two-dimensional autonomous dynamical systems , 2005, Appl. Math. Lett..

[25]  Donald L. DeAngelis,et al.  A Model for Tropic Interaction , 1975 .

[26]  Tzy-Wei Hwang,et al.  Uniqueness of limit cycles of the predator–prey system with Beddington–DeAngelis functional response , 2004 .

[27]  R Arditi,et al.  Parametric analysis of the ratio-dependent predator–prey model , 2001, Journal of mathematical biology.