The error function is defined as
$$ erf (x )= \frac{2}{\sqrt{\pi}}\int_0^x e^{-t^2} dt . $$
But expanding the integrand $e^{-t^2}$ into Taylor series, it is easy to get the Taylor series of erf. We then see that it defines an entire function.
The question is, what is it order of growth?
For this kind of problem with $\text{erf}(x)$, an approximation such as $$\text{erf}(x)\sim \text{sgn}(x)\sqrt{1-e^{-\frac{4 }{\pi } x^2}} $$ could be sufficient.