Affiliation:
1. U. Autónoma de Madrid, Santander, Spain
2. U. de Cantabria, Santander, Spain
3. CWI, Amsterdam, The Netherlands
Abstract
Fortran 77 programs for the computation of modified Bessel functions of purely imaginary order are presented. The codes compute the functions
K
ia
(
x
),
L
ia
(
x
) and their derivatives for real
a
and positive
x
; these functions are independent solutions of the differential equation
x
2
w" +
xw
' + (
a
2
−
x
2
)
w
= 0. The code also computes exponentially scaled functions. The range of computation is (
x
,
a
) ∈ (0,1500] × [−1500, 1500] when scaled functions are considered and it is larger than (0,500] × [−400, 400] for standard IEEE double precision arithmetic. The relative accuracy is better than 10
−13
in the range (0,200] × [−200, 200] and close to 10
−12
in (0, 1500] × [−1500, 1500].
Publisher
Association for Computing Machinery (ACM)
Subject
Applied Mathematics,Software
Cited by
22 articles.
订阅此论文施引文献
订阅此论文施引文献，注册后可以免费订阅5篇论文的施引文献，订阅后可以查看论文全部施引文献