Considere la posibilidad de $\,t:=x^n\,$, entonces la integral se convierte en
\begin{align}
I(n)&:=\int_0^1 \frac{x^n}{1+x} \, dx\\
&=\int_0^1 \frac{t}{1+t^{1/n}} \frac{t^{1/n-1}}n\, dt\\
&=\frac 1n \int_0^1 \frac{t^{1/n}}{1+t^{1/n}}\, dt\\
\end{align}
Usted puede expandir $\,\dfrac{t^{1/n}}{n\;(1+t^{1/n})}\,$ en la serie como $n\to \infty$ (es decir, expandir $e^{\log(t)/n}$) (más rápido con un CAS) :
$$\frac 1{2n}+\frac{\log(x)}{4n^2}-\frac{\log(x)^3}{48 n^4}+\cdots$$
y la integración entre el $0$ $1$ debe devolver :
$$\frac 1{2n}-\frac{1}{4n^2}+\frac{1}{8 n^4}-\cdots$$
(la integral de la $\log(x)^n$ no fue dada hace mucho tiempo como $(-1)^n\,n!$)
Solución explícita:
Los coeficientes pueden ser obtenidos a partir de la generación de la función de los números de Bernoulli , de modo que el resultado final es simplemente :
$$I(n)\sim \sum_{k>0} \frac {1-2^{k}}{k}\frac{B_{k}}{n^k}$$