QCD evolution equations: Numerical algorithms from the Laguerre expansion