function g=GG(u)
% control variate for callopt2.
% this function integrates to 2*(.53)^3+.53^2/2
u=max(0,u-.47);
g=6*u.^2+u;