sum(n=1:inf,HarmonicNumber[n]-log(n+1/2)-1/(24n^2)-gamma)