function y = qfunc(x)
  
  y = 0.5*erfc(x/sqrt(2));
  
  endfunction