rgsl::dawson
pub fn dawson(x: f64) -> f64
This routine computes the value of Dawson’s integral for x.