Computing solutions of the modified Bessel differential equation for imaginary orders and positive arguments
Amparo Gil, Javier Segura, Nico M. Temme
TL;DR
A variety of methods to compute the functions K, L, IA and X and their derivatives for real and positive real and negative positive values and the implementation in Fortran 77 codes is described.
Abstract
We describe a variety of methods to compute the functions $K_{ia}(x)$, $L_{ia}(x)$ and their derivatives for real $a$ and positive $x$. These functions are numerically satisfactory independent solutions of the differential equation $x^2 w'' +x w' +(a^2 -x^2)w=0$. In an accompanying paper (Algorithm xxx: Modified Bessel functions of imaginary order and positive argument) we describe the implementation of these methods in Fortran 77 codes.
