scipy.special.erf

scipy.special.erf(x[, out])

y=erf(z) returns the error function of complex argument defined as as 2/sqrt(pi)*integral(exp(-t**2),t=0..z)