function y = G0(x)
%
% Compute int_{-\infty}^x t p(t) d t
%
% where p(x) is the N(0,1) density.
%
y = - exp(-x.*x/2)/sqrt(2*pi);
