pgdraw.moments function

Compute exact first and second moments for the Polya-Gamma distribution, PG(b, c)