# What is $\int_{0}^{z} e^{-a^{2} x^{2}} {\rm erf}(bx)\, dx$?

The integral $$\int_{0}^{z} e^{-a^{2} x^{2}} {\rm erf}(bx)\, dx$$ is related to the convolution of two half-normal distributions. This can be inferred from this question on MSE. The following expression for the definite version of this integral is known. See, for instance, "A Table of Integrals of the Error Functions" by Edward W. Ng and Murray Geller, equation 4.3.2 (link here). It is as follows: $$\int_{0}^{\infty} e^{-a^{2} x^{2} } {\rm erf} (bx) \,dx = \frac{\sqrt{\pi}}{2a} - \frac{1}{a \sqrt{\pi}} \tan^{-1} \big{(} \frac{a}{b} \big{)} .$$

However, I don't know about a similar expression for the indefinite integral. Do you know about a paper or other source that adresses the indefinite integral? Or do you know how to compute it yourself?

The reason I ask about this integral is because I'd like to compute the convolution of two half-normal distributions. Please let me know if you have a reference in which the convolution of two half-normal distributions with unequal variances is calculated.

This indefinite integral is a special function, called Owen's T:

$$\int_0^z e^{-a^2 x^2}{\rm erf}\,(bx)\,dx=\frac{\arctan(b/a)}{a\sqrt\pi}-\frac{2\sqrt\pi}{a} T\left(\sqrt{2} az,b/a\right)$$

Here is the requested derivation:

$$\int_0^z e^{-a^2x^2}{\rm erf}\,(bx)\,dx=\frac{2}{\sqrt\pi}\int_0^z e^{-a^2x^2}\left(\int_0^{bx}e^{-t^2}dt\right)dx=$$ $$\frac{2}{\sqrt\pi}\int_0^z e^{-a^2x^2}\left(\int_0^{b}e^{-(yx)^2}x\,dy\right)dx=$$ $$\frac{2}{\sqrt\pi}\int_0^b\left(\int_0^z e^{-(a^2+y^2)x^2}xdx\right)dy=\\ \frac{1}{\sqrt\pi}\int_0^b\frac{e^{-(a^2+y^2)z^2}-1}{a^2+y^2}\,dy=$$ $$\frac{1}{a\sqrt\pi}\left[ \arctan(b/a)-2\pi T(\sqrt{2}az,b/a)\right]$$ with the definition $T(z,b)=\frac{1}{2\pi}\int_0^b(1+y^2)^{-1}\exp[-\tfrac{1}{2}z^2(1+y^2)]\,dy$ of Owen's T-function.

• Thank you for your answer! Do you have a source in which this expression for the indefinite integral is derived? Feb 9, 2018 at 15:36
• I don't know of a published reference (it's not in Wikipedia's list of Gaussian integrals), but I have added a derivation. Feb 9, 2018 at 20:28

Maple does not evaluate this integral. So I suspect there is no known evaluation.

I did find that integration by parts can switch the $a$ and $b$:

$$\int \!{{\rm e}^{-{a}^{2}{x}^{2}}}{\rm erf} \left(bx\right)\,{\rm d}x= {\frac {{\rm erf} \left(bx\right)\sqrt {\pi }\;{\rm erf} \left(ax \right)}{2 a}}-\frac{b}{a}\int \!{ {{{\rm e}^{-{b}^{2}{x}^{2}}}{\rm erf} \left(ax\right)}}\,{\rm d}x$$