Bifurcations were studied numerically in the system of partial differential equations, which is a one variant of predator-prey models. The mathematical model takes into account spatial distribution in habitat, active directed predator movements, birth and death process in prey population. The analysis of possible population dynamics development was performed by two qualitatively different discrete sampling techniques (Bubnov–Galerkin’s method and grid method).