# Dawson(x_):=exp(-x^2)*erfi(x)*pi^(1/2)/2;