Which erf is scipy using? There is one in cephes/ndtr.c and another in cdflib/erf.f . The two don't give the same results for smallish x, and I suspect that erf.f is the more accurate because ndtr computes it by 1.0 - erfc(x) .