La integral $I(a)=\int_0^{\infty}\sin(x^a)=\Im\int_0^{\infty}\exp(ix^a)$ se calcula fácilmente utilizando la analiticidad del integrando: Gira el contorno de integración un ángulo de $\frac{\pi}{2a}$ para obtener $I(a)=\Im\left(e^{\frac{i\pi}{2a}}\int_0^{\infty}e^{-x^a} \right)$ o
$$ I(a)=\frac{1}{a}\Gamma\left(\frac{1}{a}\right)\sin\left(\frac{\pi}{2a}\right) $$
Veamos ahora
$$ P_n=\left(n^{\frac12}\prod_{r=1}^{n-1}I\left(\frac{n}r\right)\right)^{\frac1n}=n^{\frac1{2n}}\exp\left(\frac1n\underbrace{\sum_{r=1}^{n-1}\log\left(I\left(\frac{n}{r}\right)\right)}_{s_{n-1}}\right) $$
Ahora observamos que $I(1)=0$ así que $s_n=s_{n-1}$ así que
$$ \frac{1}{n}s_{n-1}=\frac{1}{n}s_n=\frac{1}{n}\sum_{r=1}^{n}\log\left(\frac{r}{n}\sin\left(\frac{\pi r}{2n}\right)\Gamma\left(\frac{r}{n}\right)\right) $$
que es la definición de una integral de Riemann. Puesto que también el prefactor en $P_n$ es finito ( $n^{\frac{1}{2n}}\sim_{\infty}1$ ) podemos dividir el límite
$$ L=\lim_{n\rightarrow\infty}P_n=\lim_{n\rightarrow\infty}n^{\frac{1}{2n}}\lim_{n\rightarrow\infty}e^{\frac{1}{n}s_n}=1\cdot e^J \,\,\quad \left(\star\right) $$
donde $J$ viene dado por
$$ J=\int_0^1\log(I(x))dx=\\\int_0^1\log(x)dx+\int_0^1\log\left(\sin\left(\frac{\pi }{2}x\right)\right)dx+\int_0^1\log\left(\Gamma(x)\right)dx $$
o
$$ J=-1-\log(2)+\frac{1}{2}\log(2\pi) $$
donde utilizamos un resultado clásico de Raabe así como este favorito de todos los tiempos .
Ahora podemos concluir de $\left(\star\right)$ que
$$ L=\frac{\sqrt{2\pi}}{2e} $$
de acuerdo con @math110
3 votos
Después de jugar con WA+dx) parece que $$ I_{(p/q)} = \cos(\theta) \Gamma((p+q)/p) $$ para algún ángulo normalmente bonito $\theta$ . Aunque todavía no he conseguido el patrón exacto.