Algorithm 644

This algorithm is a package of subroutines for Computing Bessel functions H v (1) ( z ), H v (2) ( z ), I v (z) , J v (z) , K v (z) , Y v (z) and Airy functions Ai( z ), Ai′( z ), Bi( z ), Bi′( z ) for orders v ≥0 and complex z in −π<arg z ≤π. Eight callable subroutines and their double-precision counterparts are provided. Exponential scaling and sequence generation are auxiliary options.