Return to Revisions

We have $\phi(x)=\frac 1{\sqrt{2\pi}}\exp\left(-\frac{†^2}2\right)$ and $\Phi(x)=\int_{-\infty}^x\phi(t)dt$. We try to compute $$ I(a,b):=\int\phi(x)\Phi\left(\frac{x-b}a\right)dx.$$ Using the dominated convergence theorem, we are allowed to take the derivative with respect to $b$ inside the integral. We have $$\partial_bI(a,b)=\int\phi(x)\left(-\frac 1a\right)\phi\left(\frac{x-b}a\right)dx$$ and \begin{align} 2\pi\phi(x)\phi\left(\frac{x-b}a\right)&=\exp\left(-\frac 12\left(x^2+\frac{x^2}{a^2}-2\frac{bx}{a^2}+\frac{b^2}{a^2}\right)\right)\\\ &=\exp\left(-\frac 12\frac{a^2+1}{a^2}\left(x^2-2\frac b{a^2+1}x+\frac{b^2}{a^2+1}\right)\right)\\\ &=\exp\left(-\frac 12\frac{a^2+1}{a^2}\left(x-\frac b{a^2+1}\right)^2-\frac 12\frac{a^2+1}{a^2}\left(\frac{b^2}{a^2+1}-\frac{b^2}{(a^2+1)^2}\right)\right)\\\ &=\exp\left(-\frac 12\frac{a^2+1}{a^2}\left(x-\frac b{a^2+1}\right)^2\right)\exp\left(-\frac{b^2}{2a^2}\frac{a^2+1-1}{a^2+1}\right)\\\ &=\exp\left(-\frac 12\frac{a^2+1}{a^2}\left(x-\frac b{a^2+1}\right)^2\right)\exp\left(-\frac{b^2}{2(a^2+1)}\right). \end{align} Integrating with respect to $x$, we get that $$\partial_b I(a,b)=-\frac 1{\sqrt{a^2+1}}\phi\left(\frac b{\sqrt{a^2+1}}\right).$$ Since $\lim_{b\to +\infty}I(a,b)=0$, we have \begin{align}I(a,b)&=\int_b^{+\infty}\frac 1{\sqrt{a^2+1}}\phi\left(\frac s{\sqrt{a^2+1}}\right)ds\\\ &=\int_{b/\sqrt{a^2+1}}^{+\infty}\phi(t)dt = 1 - \Phi(b/\sqrt{a^2+1}). \end{align} This can be expressed with the traditional erf function.