An Easily Implemented Generalised Inverse Gaussian Generator

Using the ‘ratio’ method an easily implemented algorithm is derived for the generalised inverse Gaussian distribution. Computer timings and efficiency calculations show that the procedure is fast over a wide range of distribution parameter values.