pub fn Bi_deriv(x: f64, mode: Mode) -> f64
This routine computes the Airy function derivative Bi’(x) with an accuracy specified by mode.