% Return two-sided normal Pvalue Prob(|N(0,1)| >= |xx|)
% The constant is sqrt(2)
function pval = normpval(xx)
  pval = erfc(abs(xx)/(1.41421356237309504880));

