1/4 (-n-n^2+2 n HarmonicNumber(n))