>>>>> "David" == CALCRTS writes:
David> In src/bessel.lisp:
David> - there is a definition of $airy, which I assume is the
David> Airy function Ai
David> - code to numerically evaluate $airy using code in slatec
Yes, I think I did that because bessel.lisp already defined an $airy
function.
David> In src/numerical/slatec;
David> - there is numerical code for the Airy Bi function, but I can't find
David> this accessibel at the maxima level
Most likely because bessel.lisp didn't have a Bi function.
David> - which numerical code should we use?
David> The slatec code is pretty solid
I would suggest the slatec code. At least it's double-precision. The
code in share seems to be single-precision at best.
Ray