FullSimplify[HarmonicNumber[n, 2] - Zeta[2] +PolyGamma[0, n + 1]]