included upper bound for modified bessel function