Antiderivative of inverse error function

Theorem

The following formula holds: $$\displaystyle\int \mathrm{erf}^{-1}(x) \mathrm{d}x = -\dfrac{\exp \left( {-[\mathrm{erf}^{-1}(x)]^2} \right) }{\sqrt{\pi}}+C,$$ where $\mathrm{erf}^{-1}$ denotes the inverse error function, $\exp$ denotes the exponential, and $\pi$ denotes pi.