On a stochastic differential equation modeling of prey-predator evolution

Abstract

We study a stochastic differential equation model of prey-predator evolution. To keep in line with a known deterministic model we include the social and interaction terms with the drift in our model, and the randomness arises as fluctuations in the ecosystem. The notion of equilibrium population level and special types of fluctuations force us to work with degenerate elliptic operators. We consider the propagation of the population system in a sufficiently large but bounded domain. This enables us to look at not only the population extinction, but also the explosion beyond a certain level. Both extinction and explosion are possible; and, when they are not, we show that the population asymptotically reaches the equilibrium level. We show that the extinction, explosion and saturation probabilities satisfy, as functions of the initial population size, an integral equation arising out of a Dirichlet problem for a non-degenerate elliptic equation; and these probabilities are also smooth solutions of the Dirichlet problems. They are also used to express the solution of another Dirichlet problem for a degenerate elliptic equation.