Integración de una función exponencial con un argumento trigonométrico

No podemos esperar una antiderivada explícita pero

$$x^{\sin(x)}=e^{\sin(x)\log(x)}$$ Ahora, componga series de Taylor alrededor$x=0$tener $$x^{\sin(x)}=\sum_{n=0}^\infty \frac{a_n}{n!} x^n$$ donde el$a_n$son polinomios en$t=\log(x)$.

Calculando los primeros

$$\left( \begin{array}{cc} n & a_n \\ 0 & 1 \\ 1 & t \\ 2 & t^2 \\ 3 & t^3-t \\ 4 & t^4-4 t^2 \\ 5 & t^5-10 t^3+t \\ 6 & t^6-20 t^4+16 t^2 \\ 7 & t^7-35 t^5+91 t^3-t \\ 8 & t^8-56 t^6+336 t^4-64 t^2 \\ 9 & t^9-84 t^7+966 t^5-820 t^3+t \\ 10 & t^{10}-120 t^8+2352 t^6-5440 t^4+256 t^2 \end{array} \right)$$ lo que significa que para$n \ge 1$, te enfrentas a integrales $$I_{m,n}=\int_0^1 x^m \,\log ^n(x)\,dx=(-1)^n \frac {n!}{(m+1)^{n+1} }$$

Usando la tabla anterior, en función del orden de la expansión deberíamos tener las aproximaciones sucesivas para

$$\int_0^1 x^{\sin(x)}\,dx$$ $$\left( \begin{array}{ccc} 0 & 1 & 1.000000 \\ 1 & \frac{3}{4} & 0.750000 \\ 2 & \frac{85}{108} & 0.787037 \\ 3 & \frac{5485}{6912} & 0.793547 \\ 4 & \frac{17089937}{21600000} & 0.791201 \\ 5 & \frac{461505799}{583200000} & 0.791334 \\ 6 & \frac{380113636145857}{480290277600000} & 0.791425 \end{array} \right)$$ mientras que la integración numérica da$0.791403$.

Por diversión, usé todos los valores dados en la primera tabla y obtuve

$$\frac{11803390040762292561031266684259753655627}{14914487726878692033020558868480000000000}=\color{red}{0.7914043}22$$ ser comparado con$\color{red}{0.791404300}$

para el caso de

$$\int_0^{\frac \pi 2} x^{\sin(x)}\,dx$$ lo mismo podría aplicarse pero es más complicado ya que $$I_{m,n}=\int_0^{\frac \pi 2} x^m \,\log ^n(x)\,dx=(-1)^n \frac{\Gamma \left(n+1,(m+1) \log \left(\frac{2}{\pi }\right)\right) } {(m+1)^{n+1} }$$ donde aparece la función gamma incompleta.

Usando todos los términos dados, deberíamos tener$1.52012$mientras que la integración numérica da$1.52029$

Editar (por su curiosidad)

Hubiera sido agradable trabajar con

$$x^{\sin(x)}=\sum_{n=0}^\infty \frac{1}{n!} \big[ \sin(x)\log(x) \big]^n$$

$$\int x^{\sin(x)}\,dx=\sum_{n=0}^\infty \frac{1}{n!}\int \big[ \sin(x)\log(x) \big]^n\,dx$$ las integrales $$I_n=\int \big[ \sin(x)\log(x) \big]^n\,dx$$ son conocidas pero involucran un montón de funciones hipergeométricas (incluso con coeficientes imaginarios).

La excepción es la primera que es fácil de calcular usando una integración por partes

$$I=\int \sin(x)\log(x) \,dx=\text{Ci}(x)-\log (x) \cos (x)$$ a partir del cual $$I=\int_0^1 \sin(x)\log(x) \,dx=\text{Ci}(1)-\gamma$$ daría para la integral definida$0.760188$. $$I=\int_0^{\frac \pi 2} \sin(x)\log(x) \,dx=\text{Ci}\left(\frac{\pi }{2}\right)-\gamma$$ daría para la integral definida$1.46558$.