# Dawsonminus(x_):=exp(x^2)*erf(x)*pi^(1/2)/2, Dawson_(x)=D_(x);