Calculate modified Bessel function of the first kind and order 0
i0kappa calculates individual values of I_0(), the modified Bessel function of the first kind and order 0, so long as the argument is at least -3.75. Polynomial approximations are used for the separate cases in which the argument is more than or not more than 3.75. See Abramowitz and Stegun (1965), in particular formulae (9.8.1) and (9.8.2) on p.378.
Usually this program will be used within other programs. For general background on this and related functions, see texts such as Andrews, Askey and Roy (1999). The particular problem motivating this program was supporting calculations for the von Mises distribution in circular statistics, on which see for example Fisher (1993).
. i0kappa 5
Nicholas J. Cox, University of Durham, U.K. email@example.com
Abramowitz, M. and Stegun, I.A. (eds) 1965. Handbook of mathematical functions. New York: Dover.
Andrews, G.E., Askey, R. and Roy, R. 1999. Special functions. Cambridge: Cambridge University Press.
Fisher, N.I. 1993. Statistical analysis of circular data. Cambridge: Cambridge University Press.