How does one calculate these four following integrals?
$$ \int_0^1\frac{\ln(1\pm\varepsilon x)\ln(x)^3}{1\pm \varepsilon x}\,dx,\;\forall\varepsilon\in\{-1,1\}. $$ CONTEXT: Our teacher asks us to to calculate these four integrals using only changes of variables, integrations by parts and the following known result: $$\int_0^1 \frac{\ln^n(x)}{1-x} \; dx=(-1)^n n!\zeta(n+1),\quad \int_0^1 \frac{\ln^n(x)}{1+x} \; dx=(-1/2)^n (-1 + 2^n) \Gamma(1 + n) \zeta (1 + n)$$ without using complex analysis, series, differentiation under the integral sign, double integrals or special functions.
For calculate $ U =\int_0^1 \frac{\ln(1+x)\ln^3 x}{1-x}\,dx\\$ By IBP
$ U =\left[\left(\int_0^x \frac{\ln^3t}{1-t}\,dt\right)\ln(1+x)\right]_0^1-\int_0^1 \frac{1}{1+x}\left(\int_0^x\frac{\ln^3t}{1-t}\,dt\right)\,dx\\ =-6\zeta(4)\ln 2+\int_0^1\int_0^1\left(\frac{\ln^3(tx)}{(1+t)(1+x)}-\frac{\ln^3(tx)}{(1+t)(1-tx)}\right)\,dt\,dx\\ =-6\zeta(4)\ln 2+6\left(\int_0^1\frac{\ln^2 t}{1+t}\,dt\right)\left(\int_0^1\frac{\ln x}{1+x}\,dx\right)+2\left(\int_0^1\frac{\ln^3 t}{1+t}\,dt\right)\left(\int_0^1\frac{1}{1+x}\,dx\right)-\int_0^1 \frac{1}{t(1+t)}\left(\int_0^t \frac{\ln^3 u}{1-u}\,du\right)\,dt\\$ $=-\frac{33}{2}\zeta(4)\ln 2-\frac{9}{2}\zeta(2)\zeta(3)-\int_0^1 \frac{1}{t(1+t)}\left(\int_0^t \frac{\ln^3 u}{1-u}\,du\right)\,dt\\ \overset{\text{IBP}}=-\frac{33}{2}\zeta(4)\ln 2-\frac{9}{2}\zeta(2)\zeta(3)-\left[\ln\left(\frac {t}{1+t}\right)\left(\int_0^t \frac{\ln^3 u}{1-u}\,du\right)\right]_0^1+\int_0^1 \frac{\ln\left(\frac{t}{1+t}\right)\ln^3 t}{1-t}\,dt\\ =-\frac{45}{2}\zeta(4)\ln 2-\frac{9}{2}\zeta(2)\zeta(3)+\int_0^1 \frac{\ln\left(\frac{t}{1+t}\right)\ln^3 t}{1-t}\,dt\\ =-\frac{45}{2}\zeta(4)\ln 2-\frac{9}{2}\zeta(2)\zeta(3)+24\zeta(5)-U\\ U =\boxed{-\frac{45}{4}\zeta(4)\ln 2-\frac{9}{4}\zeta(2)\zeta(3)+12\zeta(5)}$
Precisely, i can't see how to calculate $V=\int_0^1 \frac{\ln(1+x)\ln^3 x}{1+x}\,dx\\$ edit I am also interested in how to calculate the two other integrals
$$\mathcal{I}=\int_0^1\frac{\ln^3x\ln(1+x)}{x(1+x)}dx=\int_0^\infty\frac{\ln^3x\ln(1+x)}{x(1+x)}dx-\underbrace{\int_1^\infty\frac{\ln^3x\ln(1+x)}{x(1+x)}dx}_{x\mapsto 1/x}$$
$$\mathcal{I}=\int_0^\infty\frac{\ln^3x\ln(1+x)}{x(1+x)}dx+\color{blue}{\int_0^1\frac{\ln^3x\ln(1+x)}{1+x}dx}-\int_0^1\frac{\ln^4x}{1+x}dx$$
By adding $\ \mathcal{I}=\int_0^1\frac{\ln^3x\ln(1+x)}{x(1+x)}dx=\int_0^1\frac{\ln^3x\ln(1+x)}{x}dx-\color{blue}{\int_0^1\frac{\ln^3x\ln(1+x)}{1+x}dx}\ $ to both sides, the blue integral nicely cancels out and we get
$$2\mathcal{I}=\int_0^\infty\frac{\ln^3x\ln(1+x)}{x(1+x)}dx-\int_0^1\frac{\ln^4x}{1+x}dx+\underbrace{\int_0^1\frac{\ln^3x\ln(1+x)}{x}dx}_{IBP}$$
$$2\mathcal{I}=\underbrace{\int_0^\infty\frac{\ln^3x\ln(1+x)}{x(1+x)}dx}_{\text{Beta function:}\ 6\zeta(2)\zeta(3)+6\zeta(5)}-\frac54\underbrace{\int_0^1\frac{\ln^4x}{1+x}dx}_{\frac{45}2\zeta(5)}$$
or
$$\mathcal{I}=3\zeta(2)\zeta(3)-\frac{177}{16}\zeta(5)\tag1$$
But
$$\mathcal{I}=\int_0^1\frac{\ln^3x\ln(1+x)}{x}dx-\int_0^1\frac{\ln^3x\ln(1+x)}{1+x}dx$$
$$=-\frac{45}{8}\zeta(5)-\int_0^1\frac{\ln^3x\ln(1+x)}{1+x}dx\tag2$$
Subtracting (1) and (2) yields
$$\int_0^1\frac{\ln^3x\ln(1+x)}{1+x}dx=\frac{87}{16} \zeta(5)-3\zeta(2)\zeta(3)$$
The integral $\int_0^\infty\frac{\ln^3x\ln(1+x)}{x(1+x)}dx$ can be calculated without using beta function:
With $\frac1{1+x}=y$ we have
$$\int_0^\infty\frac{\ln^3x\ln(1+x)}{x(1+x)}\ dx=\int_0^1\frac{\ln^3\left(\frac{x}{1-x}\right)\ln x}{1-x}\ dx$$
$$=\int_0^1\frac{\ln^4x}{1-x}-3\int_0^1\frac{\ln^3x\ln(1-x)}{1-x}+3\underbrace{\int_0^1\frac{\ln^2x\ln^2(1-x)}{1-x}}_{IBP}-\underbrace{\int_0^1\frac{\ln x\ln^3(1-x)}{1-x}\ dx}_{IBP}$$
$$=\int_0^1\frac{\ln^4x}{1-x}-3\int_0^1\frac{\ln^3x\ln(1-x)}{1-x}+2\underbrace{\int_0^1\frac{\ln^3(1-x)\ln x}{x}}_{\large 1-x\to x}-\frac14\underbrace{\int_0^1\frac{\ln^4(1-x)}{x}\ dx}_{\large 1-x\to x}$$
$$=\frac34\int_0^1\frac{\ln^4x}{1-x}\ dx-\int_0^1\frac{\ln^3x\ln(1-x)}{1-x}\ dx$$ $$=\frac34(4!\zeta(5))+\sum_{n=1}^\infty H_n\int_0^1 x^n \ln^3x\ dx$$
$$=18\zeta(5)-6\sum_{n=1}^\infty\frac{H_n}{(n+1)^4}$$
$$=18\zeta(5)-6\sum_{n=1}^\infty\frac{H_n}{n^4}+6\zeta(5)$$
$$=18\zeta(5)-6[3\zeta(5)-\zeta(2)\zeta(3)]+6\zeta(5)$$
$$=6\zeta(2)\zeta(3)+6\zeta(5)$$
Your integral can be related to harmonic series:
$$\int_0^1\frac{\ln^3x\ln(1+x)}{1+x}dx=-\sum_{n=1}^\infty (-1)^nH_n\int_0^1 x^n \ln^3xdx$$
$$=6\sum_{n=1}^\infty\frac{(-1)^nH_n}{(n+1)^4}=-6\sum_{n=1}^\infty\frac{(-1)^nH_{n-1}}{n^4}$$ $$=-6\sum_{n=1}^\infty\frac{(-1)^nH_{n}}{n^4}+6\sum_{n=1}^\infty\frac{(-1)^n}{n^5}$$
$$=-6\sum_{n=1}^\infty\frac{(-1)^nH_{n}}{n^4}-\frac{45}{8}\zeta(5)$$