jdzo Subroutine

This procedure computes the zeros of Bessel functions Jn(x) and
Jn'(x), and arrange them in the order of their magnitudes.

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.

01 August 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.

Output, integer ( kind = 4 ) N(*), the  order of Jn(x) or Jn'(x) associated
with the L-th zero.

Output, integer ( kind = 4 ) M(*), the serial number of the zeros of Jn(x)
or Jn'(x) associated with the L-th zero ( L is the serial number of all the
zeros of Jn(x) and Jn'(x) ).

Output, character ( len = 4 ) P(L), 'TM' or 'TE', a code for designating 
the zeros of Jn(x)  or Jn'(x).  In the waveguide applications, the zeros
of Jn(x) correspond to TM modes and those of Jn'(x) correspond to TE modes.

Output, real ( kind = 8 ) ZO(*), the zeros of Jn(x) and Jn'(x).