Merge pull request #501 from pv/special-fresnel

BUG: special: add a missing term in the asymptotic expansion of Fresnel C/S
BUG: special: switch to AMOS implementation of Airy functions at large arguments

The Cephes implementation gives wrong results (relative tolerance >
1e-12) after |x| > 10.
GEN: special: regenerate files