dist_poisson_inverse_gaussian function

The Poisson-Inverse Gaussian distribution