Numerical methods for solving moment equations in kinetic theory of neuronal network dynamics