Calculus · real student question

Evaluate the integral from 0 to 1 of (ln x)^4 * ln(1 - x) / (x(1 - x)) dx.

Question

Evaluate I=01ln4xln(1x)x(1x)dx.I=\int_0^1 \frac{\ln^{4}x\,\ln(1-x)}{x(1-x)}\,dx.

Step-by-step solution

  1. Split the kernel into two standard pieces. Partial fractions give 1x(1x)=1x+11x\frac{1}{x(1-x)}=\frac1x+\frac1{1-x}, so I=I1+I2I=I_1+I_2 with I1=01ln4xln(1x)xdx,I2=01ln4xln(1x)1xdx.I_1=\int_0^1\frac{\ln^{4}x\,\ln(1-x)}{x}\,dx,\qquad I_2=\int_0^1\frac{\ln^{4}x\,\ln(1-x)}{1-x}\,dx . The two halves need different series, which is exactly why the naive symmetry argument that sends x1xx\mapsto 1-x fails here.

  2. Evaluate the first half with the logarithm series. Using ln(1x)=n1xnn\ln(1-x)=-\sum_{n\ge1}\frac{x^{n}}{n} and the standard moment 01xn1ln4xdx=4!n5=24n5\int_0^1 x^{n-1}\ln^{4}x\,dx=\frac{4!}{n^{5}}=\frac{24}{n^{5}}, I1=n11n24n5=24ζ(6).I_1=-\sum_{n\ge1}\frac{1}{n}\cdot\frac{24}{n^{5}}=-24\zeta(6).

  3. Evaluate the second half with the harmonic generating function. Here ln(1x)1x=n1Hnxn\frac{\ln(1-x)}{1-x}=-\sum_{n\ge1}H_n x^{n}, and 01xnln4xdx=24(n+1)5\int_0^1 x^{n}\ln^{4}x\,dx=\frac{24}{(n+1)^{5}}, so I2=24n1Hn(n+1)5.I_2=-24\sum_{n\ge1}\frac{H_n}{(n+1)^{5}} . Shifting the index with m=n+1m=n+1 and Hm1=Hm1mH_{m-1}=H_m-\frac1m turns this into I2=24(m1Hmm5ζ(6)).I_2=-24\left(\sum_{m\ge1}\frac{H_m}{m^{5}}-\zeta(6)\right).

  4. Add the halves. The ζ(6)\zeta(6) terms cancel: I=I1+I2=24ζ(6)24m1Hmm5+24ζ(6)=24m1Hmm5.I=I_1+I_2=-24\zeta(6)-24\sum_{m\ge1}\frac{H_m}{m^{5}}+24\zeta(6)=-24\sum_{m\ge1}\frac{H_m}{m^{5}} . Everything now hinges on one linear Euler sum.

  5. Apply Euler's formula for the linear sum. With m1Hmms=(1+s2)ζ(s+1)12k=1s2ζ(k+1)ζ(sk)\sum_{m\ge1}\frac{H_m}{m^{s}}=\left(1+\frac{s}{2}\right)\zeta(s+1)-\frac12\sum_{k=1}^{s-2}\zeta(k+1)\zeta(s-k) at s=5s=5, m1Hmm5=72ζ(6)ζ(2)ζ(4)12ζ(3)2.\sum_{m\ge1}\frac{H_m}{m^{5}}=\frac72\zeta(6)-\zeta(2)\zeta(4)-\frac12\zeta(3)^{2}. Hence I=84ζ(6)+24ζ(2)ζ(4)+12ζ(3)2I=-84\zeta(6)+24\zeta(2)\zeta(4)+12\zeta(3)^{2}.

  6. Reduce to a closed form and check numerically. Using ζ(6)=π6945\zeta(6)=\frac{\pi^{6}}{945} and ζ(2)ζ(4)=π26π490=π6540\zeta(2)\zeta(4)=\frac{\pi^{2}}{6}\cdot\frac{\pi^{4}}{90}=\frac{\pi^{6}}{540}, the π6\pi^{6} terms combine as 445π6+245π6=245π6-\frac{4}{45}\pi^{6}+\frac{2}{45}\pi^{6}=-\frac{2}{45}\pi^{6}, giving I=12ζ(3)22π645I=12\zeta(3)^{2}-\frac{2\pi^{6}}{45}. Numerically this is 25.389119-25.389119\ldots, and double-exponential quadrature of the original integral returns 25.389119-25.389119\ldots as well.

Answer

I=12ζ(3)22π64525.389119I=12\zeta(3)^{2}-\frac{2\pi^{6}}{45}\approx -25.389119

Need to solve a different problem like this? Open the solver →