SUGERENCIA:
Como $x^4-1=(x^2)^2-1=(x^2+1)(x^2-1),$
$$\frac{x^4 + 1}{x^2 +1}=\frac{x^4-1+2}{x^2 +1}=x^2-1+\frac2{x^2+1}$$
Generalización :
Si el integrando es $$\frac{a_0+a_1x+a_2x^2+a_3x^3+\cdots}{ax^2+bx+c},$$
solo hay que dividir para obtener el cociente de la forma $b_0+b_1x+b_2x^2+\cdots$ que es fácilmente integrable y la parte restante será de la forma $\frac{px+q}{ax^2+bx+c}$
Si $p=0$ $I_2$ por debajo.
demás que establezca $px+q=r\frac{d(ax^2+bx+c)}{dx}+ s=r(2ax)+rb+s$
La comparación de los coeficientes de $x,2ar=p\implies r=\frac p{2a}$
y la Comparación de las constantes, $rb+s=q\implies s=q-rb=q-\frac{pb}{2a}$
$$\text{So, }\int \frac{px+q}{ax^2+bx+c} dx$$
$$=r \int \frac{d(ax^2+bx+c)}{dx}\frac1{ax^2+bx+c}dx+s\int\frac1{ax^2+bx+c}dx $$
$$=\frac p{2a}\int \frac{d(ax^2+bx+c)} {ax^2+bx+c}+\left(q-\frac{pb}{2a}\right)\int\frac1{ax^2+bx+c}dx $$
$$\text{Now, } \int \frac{d(ax^2+bx+c)} {ax^2+bx+c}=\ln|ax^2+bx+c|+C$$
$$\text{and } I_2=\int\frac1{ax^2+bx+c}dx =\int\frac{4a}{(2ax+b)^2+4ac-b^2}dx $$
Si $4ac-b^2=0,$ puesto $2ax+b=u$ $I_2$
Si $4ac-b^2>0, 4ac-b^2=t^2$(es decir, ) $I_2=4a\int \frac{dx}{(2ax+b)^2+t^2}$ y poner $2ax+u=t\tan\theta$
Si $4ac-b^2<0, 4ac-b^2=-t^2$(es decir, ) $I_2=4a\int \frac{dx}{(2ax+b)^2-t^2}$ y poner $2ax+u=t\sec\theta$