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