Implementation of the Polya-Aeppli distribution