1/2+ (sum 1.0*(-1)^(j+1)*j!*StirlingS2[4+1,j+1]*(zeta(j+1)-1-1/2^(j+1)) from j= 1 to 4)