This subroutine calculates the First Kind Bessel Function of
order 1, for any real number X. The polynomial approximation by
series of Chebyshev polynomials is used for 0<X<8 and 0<8/X<1. References:
M.Abramowitz, I.A.Stegun, Handbook of Mathematical Functions, 1965. C.W.Clenshaw, Chebyshev Series for Mathematical Functions, Mathematical Tables, Vol. 5, 1962.

Nodes of different colours represent the following:

Solid arrows point from a procedure to one which it calls. Dashed
arrows point from an interface to procedures which implement that interface.
This could include the module procedures in a generic interface or the
implementation in a submodule of an interface in a parent module.