Sé el OP dijo que "sin complejos de análisis," pero yo voy a ofrecer una manera de hacer esta integral usando el Teorema de los Residuos de todos modos porque es agradable ver el poder de métodos complejos, con estas integrales.
El análisis de espejo que se calcula en esta respuesta para un caso más complicado. Empezaremos considerando el siguiente contorno de la integral:
$$\oint_C dz \frac{\log{(1+z)} \log{z}}{z (1+z^2)} $$
donde $C$ es un ojo de la cerradura de contorno de radio exterior $R$ y un radio interior de $\epsilon$ como se define en la siguiente figura:
Es decir, $C$ evita los puntos de ramificación en$z=0$$z=-1$. Como se demuestra en el vinculado a la respuesta, cada logaritmo proporciona un salto de $-i 2 \pi$ a través de sus respectivas rama cortada. (La razón es debido a que el recorrido de las agujas del reloj alrededor del punto de ramificación como el contorno se recorre en sentido positivo.) Debido a que la integral sobre los arcos se desvanecen como $R \to \infty$$\epsilon \to 0$, el contorno de la integral es igual a
$$-i 2 \pi \int_0^{\infty} dx \frac{\log{(1+x)}}{x (1+x^2)} - i 2 \pi \int_{e^{i \pi}}^{\infty \, e^{i \pi}} dt \frac{\log{t}}{t (1+t^2)}$$
Tenga en cuenta que la rama de corte para $\log{z}$ $[0,\infty]$ y, para los argumentos $z$ fed en $\log{z}$, $\operatorname{arg}{z} \in [0,2 \pi)$. También, la rama de corte para $\log{(1+z)}$ $(-\infty,-1]$ y, para los argumentos $z$ fed en $\log{(1+z)}$, $\operatorname{arg}{z} \in (-\pi,\pi]$.
La segunda integral se evalúa el uso de $t=u e^{i \pi}$, y es igual a
$$\begin{align}\int_1^{\infty} du \frac{i \pi +\log{u}}{u (1+u^2)} &= i \pi \int_1^{\infty} du \frac{1}{u (1+u^2)} + \int_1^{\infty} du \frac{\log{u}}{u (1+u^2)} \\ &= i \pi \int_0^1 du \frac{u}{1+u^2} - \int_0^1 du \frac{u \log{u}}{1+u^2}\\ &= i \frac{\pi}{2} \log{2} - \left [ \frac{d}{d\alpha}\sum_{k=0}^{\infty} (-1)^k \int_0^1 du \, u^{2 k +1+\alpha} \right ]_{\alpha=0} \\ &= i \frac{\pi}{2} \log{2} + \frac14 \sum_{k=0}^{\infty} \frac{(-1)^k}{(k+1)^2} \\ &= i \frac{\pi}{2} \log{2} + \frac{\pi^2}{48} \end{align}$$
El contorno de la integral es también igual a $i 2 \pi$ veces la suma de los residuos en los polos $z_1=e^{i \pi/2}$$z_2=e^{i 3 \pi/2}$. Esta suma es igual a
$$\begin{align}\sum_{k=1}^2 \operatorname*{Res}_{z=z_k} \frac{\log{(1+z)} \log{z}}{z (1+z^2)} &= \frac{\log{(1+i)} \log{\left ( e^{i \pi/2}\right )}}{i (2 i)} + \frac{\log{(1-i)} \log{\left ( e^{i 3 \pi/2}\right )}}{(-i) (-2 i)}\\&= -\frac12 \left [\log{\left (\sqrt{2} e^{i \pi/4} \right )} i \frac{\pi}{2} + \log{\left (\sqrt{2} e^{-i \pi/4} \right )} i \frac{3\pi}{2} \right ]\\ &= -i \frac{\pi}{2} \log{2} - \frac{\pi^2}{8} \end{align}$$
Por lo tanto, podemos encontrar la deseada integral de
$$\int_0^{\infty} dx \frac{\log{(1+x)}}{x (1+x^2)} + i \frac{\pi}{2} \log{2} + \frac{\pi^2}{48} = i \frac{\pi}{2} \log{2} + \frac{\pi^2}{8} $$
o
$$\int_0^{\infty} dx \frac{\log{(1+x)}}{x (1+x^2)} = \frac{5 \pi^2}{48} $$
como iba a ser mostrado.