sum [ polygamma(0, n + 3)]/(n+1/a) from n = 2 to m