airyzo Subroutine

Compute the first NT zeros of Airy functions Ai(x) and Ai'(x), 
a and a', and the associated values of Ai(a') and Ai'(a); and 
the first NT zeros of Airy functions Bi(x) and Bi'(x), b and
b', and the associated values of Bi(b') and Bi'(b).

This routine is copyrighted by Shanjie Zhang and Jianming Jin.  However, 
they give permission to incorporate this routine into a user program 
provided that the copyright is acknowledged.

14 March 2012

Shanjie Zhang, Jianming Jin

Shanjie Zhang, Jianming Jin,
Computation of Special Functions,
Wiley, 1996,
ISBN: 0-471-11963-6,
LC: QA351.C45.

Input, integer ( kind = 4 ) NT, the number of zeros.

Input, integer ( kind = 4 ) KF, the function code.
1 for Ai(x) and Ai'(x);
2 for Bi(x) and Bi'(x).

Output, real ( kind = 8 ) XA(m), a, the m-th zero of Ai(x) or
b, the m-th zero of Bi(x).

Output, real ( kind = 8 ) XB(m), a', the m-th zero of Ai'(x) or
b', the m-th zero of Bi'(x).

Output, real ( kind = 8 ) XC(m), Ai(a') or Bi(b').

Output, real ( kind = 8 ) XD(m), Ai'(a) or Bi'(b)