sum_(n=1,5000)1./(n^2(gamma+digamma(n+1)))