function y = my_normcdf(x)
% computes the cumulative normal distribution

y = 0.5*erfc(-x/sqrt(2));
end

