Fix moment calculations for negative c in PG(b, c)

Closes #7