The aim of this work is to develop a compact finite-difference approach for modeling the dynamics of predator and prey based on reaction-diffusion-advection equations with variable coefficients. Methods. To discretize a spatially inhomogeneous problem with nonlinear terms of taxis and local interaction, the balance method is used. Species densities are determined on the main grid whereas fluxes are computed at the nodes of the staggered grid. Integration over time is carried out using the high-order Runge-Kutta method. Results. For the case of one-dimensional annular interval, the finite-difference scheme on the three-point stencil has been constructed that makes it possible to increase the order of accuracy compared to the standard second-order approximation scheme. The results of computational experiment are presented and comparison of schemes for stationary and non-stationary solutions is carried out. We conduct the calculation of accuracy order basing on the Aitken process for sequences of spatial grids. The calculated values of the effective order accuracy for the proposed scheme were greater than the standard two: for the diffusion problem, values of at least four were obtained. Decrease was obtained when directional migration was taken into account. This conclusion was also confirmed for non-stationary oscillatory regimes. Conclusion. The results demonstrate the effectiveness of the derived scheme for dynamics of predator and prey system in a heterogeneous environment.