use asymptotic behavior for small arguments if possible and the usual
calculation fails
* add tests for modified bessel functions I and K