I was attempting to evaluate a series $$\sum_{n=1}^\infty \frac{1}{n} \ln\left(1+\frac{1}{n}\right)$$

Since $$\frac{1}{n}\ln\left(1+\frac{1}{n}\right)=\int_0^1 \frac{1}{n(n+t)}dt,$$

I rewrote it as $$\sum_{n=1}^\infty \int_0^1 \frac{1}{n(n+t)} \; dt$$

and switched the sum and integral:

$$\int_0^1 \sum_{n=1}^\infty \frac{1}{n(n+t)}dt$$

The sum is related to digamma. Specifically,

$$\frac{\gamma-\psi(t+1)}{t}=\sum_{n=1}^\infty \frac{1}{n(n+t)}$$

Now, integrating this is the problem:

$$\int_0^1 \frac{\gamma+\psi(t+1)}{t} \; dt=1.257746887\ldots$$

That $t$ in the denominator causes a fit. I tried integration by parts, to no avail.

I ran it through Maple and it gave me $$1.257746887\ldots$$ which is indeed what the sum converges to.

Does anyone know if we can evaluate the above digamma integral? Perhaps a numerical approximation is the best we can do? After all, Maple would not give me a closed form; just the decimal solution. Can it be related to zeta or some other advanced function somehow?

Thanks very much for any input.


This formula appeared in earlier work of Coffey 'Series of zeta values, the Stieltjes constants, and a sum S_γ(n)' (see (38a) and following expressions),
in his paper of 2011 'Series representations of the Riemann and Hurwitz zeta functions and series and integral representations of the first Stieltjes constant' (1.23 and next)

as well as in Steven Finch's 'Continued Fraction Transformation III' (page 5)
see too his article concerning the related 'Khintchine-Levy constants' (page 62 the expression for $\ln(W)$ at the middle right)

It seems that it may be written in many ways! For example :

$$S=-\sum_{k=2}^\infty \zeta'(k)$$