An efficient numerical method for the generalised Kolmogorov equation