x=1;
t=x*2^(1/2);
fi=.5*(1+erf(t))