A procedure is given for the computation of the transient state occupancy probabilities of the

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.