queue. The method makes use of the inverse discrete Fourier transform, computed by means of the fast Fourier transform. It avoids the direct evaluation of modified Bessel functions and sidesteps difficulties due to the computation of very large and very small intermediate quantities.