Rewrite Airy functions in terms of Bessel functions

Note: Maybe we should not implement these formulae but
rewrite to hyper first and then call hyperexpand? This
will however allow for less explicit control.
Fix mathematical errors in the expansion formulae for Airy(c*(d*z**n)**m)
Documentation of Airy functions
Test cases for Airy functions