  We provide analytical functions approximating $\int e^{-x^2} dx$, the basis
of which is the kink soliton and which are both accurate (error $< 0.2 %$) and
simple. We demonstrate our results with some applications, particularly to the
generation of Gaussian random fields.
