A mean extinction-time estimate for a stochastic Lotka-Volterra predator-prey model

The investigation of interacting population models has long been and will continue to be one of the dominant subjects in mathematical ecology; moreover, the persistence and extinction of these models is one of the most interesting and important topics, because it provides insight into their behavior. The mean extinction-time depends on the initial population size and satisfies the backward Kolmogorov differential equation, a linear second-order partial differential equation with variable coefficients; hence, finding analytical solutions poses severe problems, except in a few simple cases, so we can only compute numerical approximations (an idea already mentioned in Sharp and Allen (1998) 1]).In this paper, we study a stochastic Lotka-Volterra model (Allen, 2007) 2, p. 149]; we prove the nonnegative character of its solutions for the corresponding backward Kolmogorov differential equation; we propose a finite element method, whose Matlab code is offered in the Appendix; and, finally, we make a direct comparison between predictions and numerical simulations of stochastic differential equations (SDEs).

[1]  E. J. Allen,et al.  Numerical solution of first passage problems using an approximate Chapman-Kolmogorov relation , 1998 .

[2]  Atsushi Yagi,et al.  Dynamic of a stochastic predator-prey population , 2011, Appl. Math. Comput..

[3]  Linda J. S. Allen,et al.  Disease emergence in deterministic and stochastic models for host and pathogen , 2005, Appl. Math. Comput..

[4]  Yongsheng Ding,et al.  Asymptotic behavior and stability of a stochastic model for AIDS transmission , 2008, Appl. Math. Comput..

[5]  E. Allen.  Modeling with Itô Stochastic Differential Equations , 2007 .

[6]  Mark S. Gockenbach,et al.  Understanding and implementing the finite element method , 1987 .

[7]  Ke Wang,et al.  Persistence and extinction in stochastic non-autonomous logistic systems , 2011 .

[8]  James D. Murray.  Mathematical Biology: I. An Introduction , 2007 .

[9]  L. Allen,et al.  A comparison of three different stochastic population models with regard to persistence time. , 2003, Theoretical population biology.

[10]  Andrew M. Stuart,et al.  Strong Convergence of Euler-Type Methods for Nonlinear Stochastic Differential Equations , 2002, SIAM J. Numer. Anal..

[11]  D. Gilbarg,et al.  Elliptic Partial Differential Equa-tions of Second Order , 1977 .

[12]  Khachik Sargsyan,et al.  Extinction Times for Birth-Death Processes: Exact Results, Continuum Asymptotics, and the Failure of the Fokker-Planck Approximation , 2004, Multiscale Model. Simul..

[13]  H. Schurz,et al.  A comparison of persistence-time estimation for discrete and continuous stochastic population models that include demographic and environmental variability. , 2005, Mathematical biosciences.

[14]  Eric Renshaw,et al.  Asymptotic behaviour of the stochastic Lotka-Volterra model , 2003 .

[15]  L. Allen.  An introduction to stochastic processes with applications to biology , 2003 .

[16]  Claes Johnson.  Numerical solution of partial differential equations by the finite element method , 1988 .

[17]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .

[18]  Xinzhu Meng,et al.  Stability of a novel stochastic epidemic model with double epidemic hypothesis , 2010, Appl. Math. Comput..