Solving integral: $\int_0^1\ln^2{\left(\frac{1+x}{1-x}\right)} dx$ I want to show that the solution of a BVP: $\ u(x) = \ln{\frac{1+x}{1-x}}$ is in $L^2(0,1)$, so I need to show that the integral
$$\int\limits_0^1\ln^2|\frac{x+1}{x-1}|dx < \infty$$
Just looking at the function, however, it's not even defined in the interval $[0,1]$, right? So can this not actually be a solution to the BVP? Even more generally about the integral itself, does that make it just $0$? And if not, can someone explain how an integral of a function can be defined in a region when the function itself is not?
Also, I know that there exist integral calculators online and have looked this up, specifically using:
https://www.integral-calculator.com/
If you put in the given integral it says the integral could not be found and then gives a complex number approximation. What should be made of this? I have taken a complex variables course but don't really see how you could solve this using Residue Calculus since we don't have any symmetries.
EDIT: Note I had originally forgot the absolute value signs - my mistake.

Actually, this integral is very straightforward to evaluate using Cauchy's Theorem in the complex plane.
Consider the complex integral
$$\oint_C dz \, \log^3{\left (\frac{z+1}{z-1} \right )} $$
where $C$ is the following contour in the complex plane:

The outer arc of $C$ has a radius $R$ and the small circular pieces around the branch points at $z=\pm 1$ have radius $\epsilon$.  The contour integral is then equal to
$$e^{i \pi} \int_R^{1+\epsilon} dx \, \log^3{\left (\frac{x-1}{x+1} \right )} + i \epsilon \int_{\pi}^0 dx \, e^{i \phi} \, \log^3{\left (\frac{-1+\epsilon e^{i \phi}+1}{-1+\epsilon e^{i \phi}-1} \right )}  \\ + \int_{-1+\epsilon}^{1-\epsilon} dx \, \left [\log{\left (\frac{1+x}{1-x} \right )} -i \pi\right ]^3 + i \epsilon \int_{\pi}^{-\pi} d\phi \, e^{i \phi} \, \log^3{\left (\frac{1+\epsilon e^{i \phi}+1}{1+\epsilon e^{i \phi}-1} \right )}  \\ - \int_{-1+\epsilon}^{1-\epsilon} dx \, \left [\log{\left (\frac{1+x}{1-x} \right )} +i \pi \right ]^3 + i \epsilon \int_{2 \pi}^{\pi} dx \, e^{i \phi} \, \log^3{\left (\frac{-1+\epsilon e^{i \phi}+1}{-1+\epsilon e^{i \phi}-1} \right )}  \\ - e^{-i \pi} \int_R^{1+\epsilon} dx \, \log^3{\left (\frac{x-1}{x+1} \right )} + i R \int_{-\pi}^{\pi} d\theta \, e^{i \theta} \, \log^3{\left (\frac{R e^{i \theta}+1}{R e^{i \theta}-1} \right )} $$
As $\epsilon \to 0$ and $R \to \infty$, all integrals vanish or cancel except the third and fifth; the contour integral is then equal to
$$\int_{-1}^1 dx \, \left \{\left [\log{\left (\frac{1+x}{1-x} \right )} -i \pi\right ]^3 -\left [\log{\left (\frac{1+x}{1-x} \right )} +i \pi\right ]^3  \right \} $$
By Cauchy's Theorem, the contour integral is equal to zero.  Expanding the integrand, we immediately determine the value of the integral without further computation, because the cube and linear powers of the log cancel.  That is, we are left with the log squared term sought, and a constant term that is trivially dealt with.  The result is, using symmetry of an even function:
$$\int_0^1 dx \, \log^2{\left (\frac{1+x}{1-x} \right )} = \frac{\pi^2}{3}$$
