“…In [8], the authors proved convergence of their method to a bounded variation weak solution for Equation (2); these results can readily be extended to the model and scheme from [9]. Here, we develop for the model in [8] a novel second-order high-resolution scheme in the spirit of those developed in [5,7,23,35,40]. To establish convergence of this method, a substantial modification of the techniques used in [5,7] is necessary due to the explicit modelling of the genetic or phenotypic heterogeneity between individuals, nonlocal nonlinearities, and quadratic nonlinearities.…”