Hi I am stuck on showing that
∫∞0ln(1+x)ln(1+x−2)xdx=πG−3ζ(3)8
where G is the Catalan constant and ζ(3) is the Riemann zeta function. Explictly they are given by
G=β(2)=∞∑n=0(−1)n(2n+1)2, ζ(3)=∞∑n=11n3.
I have tried using
ln(1+x)=∞∑n=1(−1)n+1xnn,
but didn't get very far.
Answer
The infinite sum in Chen Wang's answer, that is, ∞∑n=1H4nn2, can be evaluated using contour integration by considering the function f(z)=πcot(πz)[γ+ψ(−4z)]z2,
where ψ(z) is the digamma function and γ is the Euler-Mascheroni constant.
The function f(z) has poles of order 2 at the positive integers, simple poles at the negative integers, simple poles at the positive quarter-integers, and a pole of order 4 at the origin.
The function ψ(−4z) does have simple poles at the positive half-integers, but they are cancelled by the zeros of cot(πz).
Now consider a square on the complex plane (call it CN) with vertices at ±(N+12)±i(N+12).
On the sides of the square, cot(πz) is uniformly bounded.
And when z is large in magnitude and not on the positive real axis, ψ(−4z)∼ln(−4z).
So ∫CNf(z) dz vanishes as N→∞ through the positive integers.
Therefore,
∞∑n=1Res[f(z),n]+∞∑n=1Res[f(z),−n]+Res[f(z),0]+∞∑n=0Res[f(z),2n+14]=0.
To determine the residues, we need the following Laurent expansions.
At the positive integers,
γ+ψ(−4z)=141z−n+H4n+O(z−n)
and
πcot(πz)=1z−n+O(z−n).
At the origin,
γ+ψ(−4z)=14z−4ζ(2)z−16ζ(3)z2+O(z3)
and
πcot(πz)=1z−2ζ(2)z+O(z3).
And at the positive quarter-integers,
γ+ψ(−4z)=141z−2n+14+O(1)
and
πcot(πz)=(−1)nπ+O(z−2n+14).
Then at the positive integers,
f(z)=1z2(141(z−n)2+H4nz−n+O(1)),
which implies
Res[f(z),n]=Res[14z21(z−n)2,n]+Res[1z2H4nz−n,n]=−12n3+H4nn2.
At the negative integers,
Res[f(z),−n]=γ+ψ(4n)n2=H4n−1n2=H4nn2−14n3.
At the origin,
f(z)=1z2(14z2−ζ(2)2−4ζ(2)−16ζ(3)z+O(z2)),
which implies
Res[f(z),0]=−16ζ(3).
And at the positive quarter-integers,
f(z)=π4z2(−1)nz−2n+14+O(1),
which implies
Res[f(z),2n+14]=Res[π4z2(−1)nz−2n+14,2n+14]=4π (−1)n(2n+1)2.
Putting everything together, we have
−12∞∑n=11n3+2∞∑n=1H4nn2−14∞∑n=11n3−16ζ(3)+4π∞∑n=0(−1)n(2n+1)2
=−12ζ(3)+2∞∑n=1H4nn2−14ζ(3)−16ζ(3)+4πG=0.
Therefore,
∞∑n=1H4nn2=678ζ(3)−2πG.
EDIT:
I found the Laurent expansion of ψ(−4z) at the positive integers by using the functional equation of the digamma function to express ψ(4z) as
ψ(4z)=ψ(4z+4n+1)−14z+4n−14z+4n−1−…−14z.
Then I evaluated the limit lim and the limit \lim_{z \to -n} \Big(\psi(4z) + \frac{1}{4} \frac{1}{z+n} \Big) = - \gamma +H_{4n} .
This leads to the expansion \gamma + \psi (-4z) = \frac{1}{4} \frac{1}{z-n} + H_{4n} + \mathcal{O}(z-n) .
I did something similar to find the expansion at the positive quarter-integers.