Fix expansion of polygamma(n, p/q) for n in Z and 1 <= p < q