  Several numerical evaluations of the density and distribution of convolution
of independent gamma variables are compared in their accuracy and speed. In
application to renewal processes, an efficient formula is derived for the
probability mass function of the event count.
