A high order numerical scheme for a nonlinear nonlocal reaction–diffusion model arising in population theory

This paper provides a numerical method for nonlinear equation arising in mathematical biology. It is an extension of another one recently proposed for the linear, less realistic, situation. The main novel result is the proof that the convergence of the numerical method is of order four, as to our knowledge no similar high accuracy results exist yet in the current literature for usually employed simulation schemes for nonlocal equations.