function h1=GG(m,sig)
%this is the function G of the paper 
% m expectation
% sig standard error 
% computes the expectation of the absolute value  of a N(m,sig^2)
 h1=  2 *sig* dnorm(m/sig) + m.* (2* pnorm(m/sig) -1); 