Rotation method for accelerating multiple-spherical Bessel function integrals against a numerical source function
bessel plot