Code covered by the BSD License
function y=q(x) y=0.5.*erfc(x./sqrt(2));
Contact us at files@mathworks.com