Added algorithms for Bessel function with real parameter K_r(x)
Added support for small eigenvalues through specifying negative R