Bessel-k0-func
evaluates an approximation to the modified Bessel function of the second kind,. The function is based on Chebyshev expansions.
For more information please review the s18acc function in the NAG document.