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

Dawson(0,x_):=Dawson(x);