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