Computation of autocorrelations of interdeparture times by numerical transform inversion

The generating function of the autocorrelations of the interdeparture times in stationary M/G/1 and GI/M/1 systems involves the probability generating function of the number of customers served in a busy period. The latter function is implicitly determined as a solution to a functional equation. Standard methods for the numerical inversion of generating functions require the values of these functions at many complex arguments. A recently discovered substitution method for contour integrals allows the numerical inversion of implicitly determined generating functions without the numerical solution of the functional equations.