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}