A non-standard numerical scheme for a generalized Gause-type predator–prey model

Abstract A non-standard finite-difference scheme is constructed to simulate a predator–prey model of Gause-type with a functional response. Using fixed-point analysis, it is shown that the scheme preserves the physical properties of the model and gives results that are qualitatively equivalent to the real dynamics of the model. It is also shown that the scheme undergoes a supercritical Hopf bifurcation for a specific value of the bifurcation parameter (k0). This leads to the existence of a stable limit cycle created by the scheme when the bifurcation parameter passes through k0, as predicted by the continuous model. The scheme is used to simulate the model with the functional responses of Holling-types II and III. The simulation results generated by the non-standard finite-difference scheme are compared with those obtained from the standard methods such as forward-Euler and Runge–Kutta methods. These comparisons show that the standard methods give erroneous results that disagree with the theoretical predictions of the model. However, it is proved that the proposed non-standard finite-difference scheme is consistent with the asymptotic dynamics of the model. Numerical simulations are presented to support these facts.

[1]  P. Ricciardi,et al.  Lyapunov functions for a generalized Gause-type model , 1995 .

[2]  J. Lambert.  Computational Methods in Ordinary Differential Equations , 1973 .

[3]  Mary R. Myerscough,et al.  Stability, persistence and structural stability in a classical predator-prey model , 1996 .

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

[5]  V. Křivan,et al.  Alternative Food, Switching Predators, and the Persistence of Predator‐Prey Systems , 2001, The American Naturalist.

[6]  H. Regan,et al.  The Currency and Tempo of Extinction , 2001, The American Naturalist.

[7]  Jitsuro Sugie.  Two-Parameter Bifurcation in a Predator–Prey System of Ivlev Type , 1998 .

[8]  Jitsuro Sugie,et al.  On a predator-prey system of Holling type , 1997 .

[9]  S. Bacher,et al.  Functional response of a generalist insect predator to one of its prey species in the field , 2002 .

[10]  S. M. Moghadas.  Some Conditions for the Nonexistence of Limit Cycles in a Predator-Prey System , 2002 .

[11]  Paul Glendinning,et al.  Stability, instability and chaos , by Paul Glendinning. Pp. 402. £45. 1994. ISBN 0 521 41553 5 (hardback); £17.95 ISBN 0 521 42566 2 (paperback) (Cambridge). , 1997, The Mathematical Gazette.

[12]  Murray E. Alexander,et al.  A Positivity-preserving Mickens-type Discretization of an Epidemic Model , 2003 .

[13]  Robert M. May,et al.  Theoretical Ecology: Principles and Applications , 1981 .

[14]  Joseph W.-H. So,et al.  Global stability and persistence of simple food chains , 1985 .

[15]  Robert E. Kooij,et al.  A predator-prey model with Ivlev's functional response , 1996 .

[16]  G. Smith,et al.  Numerical Solution of Partial Differential Equations: Finite Difference Methods , 1978 .

[17]  Seyed M. Moghadas,et al.  Existence of limit cycles for predator–prey systems with a class of functional responses , 2001 .

[18]  Ronald E. Mickens,et al.  Discretizations of nonlinear differential equations using explicit nonstandard methods , 1999 .

[19]  Jonathan M. Jeschke,et al.  PREDATOR FUNCTIONAL RESPONSES: DISCRIMINATING BETWEEN HANDLING AND DIGESTING PREY , 2002 .

[20]  Y. Kuznetsov.  Elements of Applied Bifurcation Theory , 2023, Applied Mathematical Sciences.