Implementation of the Polya–Aeppli distribution
