function q = qfunction(x);
q = 0.5*erfc(x/sqrt(2));
return
