Multistability for a mathematical model of the dynamics of predators and preys in a heterogeneous area

We consider the system of reaction–diffusion–advection equations describing the evolution of the spatial distributions of two populations of predators and two prey populations. This model allows us to consider directed migration, the Holling functional response of the second kind, and the hyperbolic prey growth function. We obtain conditions on the parameters under which cosymmetries exist. As a result, multistability is realized, i.e., the one- and two-parameter families of stationary solutions appear. For a homogeneous environment, we analytically derive explicit formulas for equilibria. With a heterogeneous habitat, we computed distributions of species using the method of lines and the scheme of staggered grids. We present the results of violation of cosymmetry and transformation of the family in the case of invasion of a predator.