Algorithm xxx: Modified Bessel functions of imaginary order and positive argument
| dc.creator | Gil, Amparo | |
| dc.creator | Segura, Javier | |
| dc.creator | Temme, Nico M. | |
| dc.date | 2004-01-13 | |
| dc.date.accessioned | 2026-07-07T03:20:49Z | |
| dc.date.available | 2026-07-07T03:20:49Z | |
| dc.description | 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'' +x w' +(a^2 -x^2)w=0$. The code also computes exponentially scaled functions. The range of computation is $(x,a)\in (0,1500]\times [-1500,1500]$ when scaled functions are considered and it is larger than $(0,500]\times [-400,400]$ for standard IEEE double precision arithmetic. The relative accuracy is better than $10^{-13}$ in the range $(0,200]\times [-200,200]$ and close to $10^{-12}$ in $(0,1500]\times [-1500,1500]$. | |
| dc.description | 6 pages, 4 figures. To appear in ACM T. Math. Software | |
| dc.identifier | https://arxiv.org/abs/cs/0401008 | |
| dc.identifier | http://arxiv.org/abs/cs/0401008 | |
| dc.identifier.uri | http://salesiana.dossiersoluciones.com/handle/123456789/31961 | |
| dc.subject | Mathematical Software | |
| dc.subject | Numerical Analysis | |
| dc.subject | G.4 | |
| dc.title | Algorithm xxx: Modified Bessel functions of imaginary order and positive argument | |
| dc.type | text |