Bessel-y1-func
evaluates the Bessel function of the second kind, , x > 0.
The approximation is based on Chebyshev expansions.
For more information please review the s17adc function in the NAG document.