Numerical Realization of Taxis-Reaction-Diffusion Model Describing Population Dynamics in Predator-Prey System
Difference schemes for numerical solution of partial differential equations, constructed to describe predator movements caused by trophotaxis, are presented for one-dimensional (segment) and two-dimensional (rectangle) closed domains. Predator and prey population densities and predator satiety are the model variables. It is shown that schemes keep general integral of predator density; schemes satisfy to the nonzero equilibrium stability condition, obtained analytically for initial continuous system. The emergence of spatially heterogeneous patterns is experimentally demonstrated