Aquí es un poco más largo de método (en comparación con @Robert Z hábil de la serie enfoque de solución) que utiliza Feynman el truco de la diferenciación bajo el signo integral.
Vamos
$$I(a) = \int_0^\infty \frac{\sin (ax)}{x} e^{-x^2} \, dx, \quad a > 0.$$
Tenga en cuenta que $I(0) = 0$ y estamos obligados a encontrar $I(1)$.
La diferenciación $I(a)$ con respecto al parámetro de $a$ da
$$I'(a) = \int_0^\infty \cos (ax) e^{-x^2} \, dx. \tag1$$
En la integración de (1) por las partes conduce a la
$$I'(a) = \frac{2}{a} \int_0^\infty x e^{-x^2} \sin (ax) \, dx.$$
También, la diferenciación de (1) con respecto al parámetro de $a$ rendimientos
$$I''(a) = -\int_0^\infty x e^{-x^2} \sin (ax) \, dx = -\frac{a}{2} I'(a).$$
Si establecemos $u(a) = I'(a)$ por encima de la de segundo orden de la ecuación diferencial se puede reducir a la siguiente ecuación diferencial de primer orden
$$u'(a) = -\frac{a}{2} u(a).$$
La solución de los rendimientos
$$u(a) = I'(a) = K e^{-a^2/4}, \tag1$$
donde $K$ es una constante a determinar. Para encontrar esta constante ajuste de $a = 0$ $I'(a)$ conduce a
$$I'(0) = \int_0^\infty e^{-x^2} \, dx = \frac{\sqrt{\pi}}{2} \cdot \frac{2}{\sqrt{\pi}} \int_0^\infty e^{-x^2} \, dx = \frac{\sqrt{\pi}}{2} \cdot \text{erf} (\infty) = \frac{\sqrt{\pi}}{2}.$$
Así que en la configuración de $a= 0$ en (1) nos encontramos con $K = \sqrt{\pi}/2$. Así
$$I'(a) = \frac{\sqrt{\pi}}{2} e^{-a^2/4}.$$
Ahora como $I(0) = 0$, tenemos
$$I(1) = \int_0^1 I'(a) \, da = \frac{\sqrt{\pi}}{2} \int_0^1 e^{-a^2/4} \, da.$$
El cumplimiento de una sustitución de $a \mapsto 2a$ conduce a
$$I(1) = \sqrt{\pi} \int_0^{1/2} e^{-a^2} \, da = \frac{\pi}{2} \cdot \frac{2}{\sqrt{\pi}} \int_0^{1/2} e^{-a^2} \, da.$$
Y puesto que a partir de la representación integral de la función de error , que está dada por
$$\text{erf}(x) = \frac{2}{\sqrt{\pi}} \int_0^x e^{-t^2} \, dt,$$
uno tiene
$$\int_0^\infty \frac{\sin x}{x} e^{-x^2} \, dx = \frac{\pi}{2} \text{erf} \left (\frac{1}{2} \right ),$$
como era de esperar.