Bounds and algorithms for the K-Bessel function of imaginary order

Using the paths of steepest descent, we prove precise bounds with numerical implied constants for the modified Bessel function of imaginary order and its first two derivatives with respect to the order. We also prove precise asymptotic bounds on more general (mixed) derivatives without working out numerical implied constants. Moreover, we present an absolutely and rapidly convergent series for the computation of and its derivatives, as well as a formula based on Fourier interpolation for computing with many values of . Finally, we have implemented a subset of these features in a software library for fast and rigorous computation of .