# Dawsonm(x_):=exp(x^2)*erf(x)*pi^(1/2)/2; Dawson2(x_):=Dawsonm(x);