Assessment of numerical methods for the evaluation of higher-order harmonics in diffusion theory