Algorithm 588. Fast Hankel transforms using related and lagged convolutions - Anderson, W. L. (ACM Trans. Math. Software, 1982)
Single-precision Fortran, maximum accuracy 8S.