(Ai,Aip,Bi,Bip)=airy(z) calculates the Airy functions and their derivatives
evaluated at real or complex number z. The Airy functions Ai and Bi
are two independent solutions of y’‘(x)=xy. Aip and Bip are the first derivatives
evaluated at x of Ai and Bi respectively.