Airy-bi-func
evaluates an approximation to the Airy function Bi(x). It is based on a number of Chebyshev expansions.
For more information please review the s17ahc function in the NAG document.