Integración por partes, se obtiene:
$$\int_{0}^{1}nx^{n-1}\frac{2x}{1+x^2}\,dx = \left.x^n\frac{2x}{1+x^2}\right|_{0}^{1}-2\int_{0}^{1}x^{n}\frac{(1-x^2)}{(1+x^2)^2}\,dx$$
para la realización de dos veces:
$$\int_{0}^{1}nx^{n-1}\frac{2x}{1+x^2}\,dx = 1-2\left.\frac{x^{n+1}}{n+1}\frac{(1-x)^2}{(1+x^2)^2}\right|_{0}^{1}+2\int_{0}^{1}\frac{x^{n+1}}{n+1}\cdot\frac{2x(3-x^2)}{(1+x^2)^3}\,dx$$
pero desde $\frac{2x(3-x^2)}{(1+x^2)^3}$ es no negativo y limitado en $(0,1)$, se sigue que:
$$\int_{0}^{1}nx^{n-1}\frac{2x}{1+x^2}\,dx = 1+O\left(\frac{1}{n^2}\right),$$
por lo que la quería límite es sólo $\color{red}{1}$.
Para ser más exactos, podemos evaluar $I_n = 2n\int_{0}^{1}\frac{x^n}{1+x^2}\,dx$ en términos de la función digamma:
$$ I_n = \frac{n}{2}\left(\psi\left(\frac{n+3}{4}\right)-\psi\left(\frac{n+1}{4}\right)\right)=\sum_{m\geq 1}\frac{4n}{(4m+n-1)(4m+n-3)}$$
o a través de la serie simple representación:
$$ I_n = 2\sum_{m\geq 1}\left(\frac{4m-1}{n+4m-1}-\frac{4m-3}{n+4m-3}\right)$$
para obtener:
$$ I_n \approx 1-\frac{1}{n^2}$$
para un gran $n$s.