The modified Bessel function KiŅχ, also known as the Macdonald function, finds application in the Kontorovich-Lebedev integral transform when χ and Ņ are real and positive. In this paper, a comparison of three codes for computing this function is made. These codes differ in algorithmic approach, timing, and regions of validity. One of them can be tested independently of the other two through Wronskian checks, and therefore is used as a standard against which the others are compared.
Citation: NIST Interagency/Internal Report (NISTIR) - 6596Report Number:
NIST Pub Series: NIST Interagency/Internal Report (NISTIR)
Pub Type: NIST Pubs
Bessel functions, computation of special functions, Macdonald function