BUG: special: switch to AMOS implementation of Airy functions at large arguments

Goes on top of gh-500 (only commits a40d621 and c0d48b1 are for this PR).

Use the Airy function implementations from AMOS instead of Cephes at large real arguments. The Cephes implementation apparently doesn't use a correct asymptotic expansion.