# convert to erf(x); toerf(erfc(x_)) := 1-erf(x); toerf(erfi(x_)) := -i*erf(-i*x); toerf(Phi(x_)):= erf(x/sqrt(2))/2+1/2;