samer@4: function y=erfratio(x)
samer@4: % erfratio - returns gausspdf(x)./erfp(x) but using a better computation
samer@4: %
samer@4: % erfratio :: real -> nonneg.
samer@4: y = 2./(sqrt(2*pi)*erfcx(x/sqrt(2)));