# 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;