Algorithm for numerical solutions to the kinetic equation of a spatial population dynamics model with coalescence and repulsive jumps (Q2028046)

The purpose of this paper is development of an efficient algorithm to solve the repulsion-jump coalescence kinetic equation in a high-dimensional space. Different techniques including time-space discretization (composite Simpson and trapezoidal rules, Runge-Kutta methods, and adjustable system-size schemes), periodic, Dirichlet, and asymptotic boundary conditions are combined for studying infinite systems on the basis of their finite samples. The algorithm is applied to population dynamics simulations of one-dimensional systems with various initial spatially inhomogeneous densities and forms of the jump, coalescence, and repulsion kernels. A comprehensive discussion of the obtained numerical results is provided for some specific choices of the model parameters and initial densities.