Complementary error function ($\mathrm{erfc}z=1-\mathrm{erf}z$)

The complementary error function describes the Gaussian integral area remaining beyond the endpoint. For a real $x$,It is an entire function for complex arguments when defined by $1-\mathrm{erf}(z)$. On the real line it decreases from two to zero, with $\mathrm{erfc}(0)=1$ and $\mathrm{erfc}(-x)=2-\mathrm{erfc}(x)$.^[ref]

For example, $\mathrm{erfc}(1)\approx 0.157299$ is a sizeable tail, whereas $\mathrm{erfc}(8)\approx 1.12243\times 10^{-29}$ is extremely small. In the tested floating-point arithmetic, a direct evaluation of $1-\mathrm{erf}(8)$ gives zero because the computed error function has already rounded to one. Evaluating the complementary error function directly preserves this tail. It is the useful quantity in an upper-tail normal distribution probability and in the decaying side of an error-function diffusion front.

For large positive real $x$, the leading asymptotic expansion isThe first corrections are $1-1/(2x^2)+3/(4x^4)-\cdots$ multiplying the leading term. This is an asymptotic expansion, not an assertion that summing arbitrarily many terms converges at a fixed $x$.^[ref] The complementary error function tail bound provides finite-$x$ inequalities instead.

For every real $x>0$,For a short proof, put $I(x)={\int }_x^{\infty }{e}^{-t^2}dt$. Integration by parts givesThe last integral is positive and strictly smaller than $I(x)/x^2$. Using these two inequalities gives the upper and lower bounds after multiplication by $2/\sqrt{\pi }$. Dividing both bounds by $e^{-x^2}/(\sqrt{\pi }x)$ also proves the leading asymptotic equivalence by the squeeze theorem. For small positive $x$, the bounds remain true but can be loose; they are especially useful when the Gaussian integral tail is already small.

The Faddeeva function is an entire function built from the complementary error function. Its integral representation for $\mathrm{Im}z>0$ is $w(z)=(i\pi {)}^{-1}{\int }_{-\infty }^{\infty }{e}^{-t^2}/(t-z)dt$, as normalized in NIST's integral representation. Its reflection identity is $w(-z)=2e^{-z^2}-w(z)$, directly from the odd function property of the error function.

For a pole off the real axis, substitution $s=2t$ into the Faddeeva function integral gives $G(a)=i\pi w(a/2)$ when $\mathrm{Im}a>0$, and $G(a)=-i\pi w(-a/2)$ when $\mathrm{Im}a<0$. The difference of boundary values is the residue theorem jump $2\pi i{e}^{-a^2/4}$. In particular, for a lower-half-plane pole, $G(a)+2\pi i{e}^{-a^2/4}=i\pi w(a/2)$, which packages a crossed residue and a nearby saddle into one smooth expression.