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