Numerical study of coexistence of populations in an environmental niche

The predator-prey interactions on the spatial heterogeneous two-dimensional area are described. The model is written as a system of nonlinear parabolic equations for two closely related predator populations and two prey populations competing for the general resource. It is shown that under certain relationships between the parameters and the variable natural habitat resource functions, the model belongs to the class of the cosymmetric dynamical systems. In this case, there is a continuous family of stationary distributions of the coexistent populations. The simulation experiment is based on the method of straight lines, and on the scheme of staggered grids. The balance method is used for the approximation in spatial variables of the task on a rectangular area. The results showing the model capabilities for describing the formation of the population stationary distributions are presented. The formation of the biological structures is studied under the growth parameter heterogeneity; the conditions for the coexistence of closely related types are analyzed.