% Return the Gaussian function Prob(N(0,1) <= xx)
% The constant is sqrt(2)
function pval = normint(xx)
  pval = (1+erf(xx/(1.41421356237309504880)))/2;


