n!
gamma(1 + n)*polygamma(0, 1 + n)
/ 2 \ \polygamma (0, 1 + n) + polygamma(1, 1 + n)/*gamma(1 + n)
/ 3 \ \polygamma (0, 1 + n) + 3*polygamma(0, 1 + n)*polygamma(1, 1 + n) + polygamma(2, 1 + n)/*gamma(1 + n)