En mi curso, tengo que demostrar la fórmula siguiente $$I=\int_0^\infty \frac{e^{\cos(ax)}\cos\left(\sin (ax)+bx\right)}{c^2+x^2}dx =\frac{\pi}{2c}\exp\left(e^{-ac}-bc\right)$$ para $a,b,c>0.$
Sé que esta integral se puede resolver fácilmente con el análisis complejo utilizando $$f(z)=\frac{1}{2} \ \mathbb{R} \left(\int_{-\infty}^\infty \frac{\exp\left(e^{iaz}+ibz\right)}{c^2+z^2}dz\right)$$ pero ahora mismo estoy en un curso que trata del análisis real. He intentado utilizar el método integral de parametrización $$I'(a)=-\int_0^\infty \frac{xe^{\cos(ax)}\sin(\sin(ax)+(a+b)x)}{c^2+x^2}dx $$ pero no parece más fácil de manejar. Intenté diferenciarlo de nuevo, pero me salió una forma horrible. Se me ocurrió una idea para diferenciar con respecto al parámetro $b$ y establecer una ecuación diferencial $$I''(b)+x^2I(b)=0$$ Si introduzco esta EDO en W|A, obtengo $$I(b)=c_1\cos(bx^2)+c_2\sin(bx^2)$$ ¡Definitivamente está mal! Después de ver la respuesta de Samrat, intenté conectarme de nuevo a W|A y obtuve $$I(b)=c_1 D_{-1/2}((i+1)b)+c_2 D_{-1/2}((i-1)b)$$
donde $D_n(z)$ es la función cilíndrica parabólica, pero no tengo ni idea de lo que significa.
¿Alguna idea? Gracias de antemano.