function res = exact1(x)
  C1  = (3*exp(-0.5) - 1) / erf( 1/sqrt(2)) ;
  res = -1 + ( 1 + C1 * erf( x ./ sqrt(2) ) ) .* exp( (x.^2) ./ 2 ) ;
end