function y = q(x)
y = erfc(x/sqrt(2))/2;