erf(x) dx = exp(-y^2) dy