$\int_{a}^{b}e^{x^2}dx=\frac{1}{2}\sqrt{\pi}(erfi(b )-erfi(a))$ fascinating how this can be done using erfi functions defined in the whole complex plane