Como una extensión de una pregunta que planteé anteriormente, pensé que sería mejor intentar y final el resultado más general de la forma:
\begin{equation} I_n = \int_{0}^{\infty} \frac{e^{-x^n}}{x^n + 1}\:dx \end{equation} con $n \in \mathbb{R}, n > 1$
Como con la anterior pregunta, estoy interesado en la búsqueda de formas alternativas de solucionar esto que no se basa en el análisis complejo.
Mi Método: se utiliza el mismo método como con mi pregunta anterior. Aquí primero vamos a
\begin{equation} J_n(t) = \int_{0}^{\infty} \frac{e^{-tx^n}}{x^n + 1}\:dx \end{equation}
Vemos que $I_n = J_n(1)$ y que $J_n(0) = \frac{1}{n}\Gamma\left(1 - \frac{1}{n}\right)\Gamma\left(\frac{1}{n}\right)$ (Esto se muestra aquí)
Ahora, tomar la derivada con respecto a '$t$' para lograr \begin{align} J_n'(t) &= \int_{0}^{\infty} \frac{-x^ne^{-tx^n}}{x^n + 1}\:dx = -\int_{0}^{\infty} \frac{\left(x^n + 1 - 1\right)e^{-tx^n}}{x^n + 1}\:dx \\ &= -\left[\int_{0}^{\infty}e^{-tx^n}\:dx - \int_{0}^{\infty}\frac{e^{-tx^n}}{x^n + 1}\:dx \right] \\ &= -\left[ \frac{t^{-\frac{1}{n}}}{n}\Gamma\left(\frac{1}{n}\right) -J_n(t)\right] \end{align}
Que los rendimientos de la ecuación diferencial:
\begin{equation} J_n'(t) - J_n(t) = -\frac{t^{-\frac{1}{n}}}{n}\Gamma\left(\frac{1}{n}\right) \end{equation}
Que los rendimientos de la solución:
\begin{equation} J_n(t) = \frac{1}{n}\Gamma\left(1 - \frac{1}{n}, t\right)\Gamma\left(\frac{1}{n}\right)e^t \end{equation}
Y por último:
\begin{equation} I_n = J_n(1) = \int_{0}^{\infty} \frac{e^{-x^n}}{x^n + 1}\:dx = \frac{e}{n}\Gamma\left(1 - \frac{1}{n}, 1\right)\Gamma\left(\frac{1}{n}\right) \end{equation}
Que para mí, es un buen resultado. Fascinado a ver otros métodos!
Editar - Gracias a spaceisdarkgreen para la camioneta de mi escribir mal de la Función Gamma Incompleta.