$$f(z)=\frac{\log z}{(1+z^2)^2} = \frac{\log z}{(z+i)^2(z-i)^2}$$
donde utilizamos la rama principal del logaritmo.
$$P.V.\int_{-R}^R f(z)\, dz = \int_{\epsilon}^R f(z)\, dz+\int_{-R}^\epsilon f(z)\, dz = \int_{\epsilon}^R f(z)\, dz + \int_{\epsilon}^R\frac{\log z+i\pi}{(1+z^2)^2}\, dz =$$
$$= 2\int_{\epsilon}^R f(z)\, dz + i\pi\int_{\epsilon}^R \frac{dz}{(1+z^2)^2}$$
A continuación, considerando el contorno $\Gamma$ que atraviesa $-R$ a $R$ (sangrado en $z=0$ en torno a $\epsilon e^{i \theta}$ ) y vuelve a recorrer $Re^{i \theta}$ sólo hay un polo de $f(z)$ en $\Gamma$ . En $R \to \infty$ la integral alrededor del arco desaparece, por lo que
$$2\pi i \operatorname*{Res}_{z = i} f(z) = 2\int_{\epsilon}^\infty f(z)\, dz + \int_{0}^\pi f(\epsilon e^{i \theta})i\epsilon e^{i \theta}\, d\theta + i\pi\int_{\epsilon}^\infty \frac{dz}{(1+z^2)^2}$$
y como $\epsilon \to 0$ vemos que el lado derecho se convierte en
$$2\int_{0}^\infty f(z)\, dz + \frac{i \pi^2}{4}$$
Los polos son de orden 2, por tanto
$$\operatorname*{Res}_{z = i} f(z) = \lim_{z \to i}\,\frac{d}{dz}(z-i)^2f(z) = \frac{\pi+2 i}{8}$$
Así
$$-\frac{\pi}{2}+ \frac{i \pi^2}{4} = 2\int_{0}^\infty f(z)\, dz + \frac{i \pi^2}{4}$$
y finalmente encontramos
$$\int_{0}^\infty \frac{\log z}{(1+z^2)^2}\, dz = -\frac{\pi}{4}$$
$\blacksquare$