Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Subroutine package for Bessel functions of a complex argument and nonnegative order

Technical Report ·
OSTI ID:5249727
The machine implementation of a Bessel function package, which was outlined in mathematical form in two previous documents, is described for Bessel functions H/sub ..nu..//sup (1)/(z), H/sub ..nu..//sup (2)/(z), I/sub ..nu../(z), J/sub ..nu../(z), K/sub ..nu../(z), Y/sub ..nu../(z) and Airy Functions Ai(z), Ai'(z), Bi(z), Bi'(z) with ..nu.. greater than or equal to 0 and /sup -/..pi.. < arg z less than or equal to ..pi... The less obvious aspects of the implementation are the main topics of interest. These topics include estimation of function magnitudes for underflow and overflow tests, scaling near underflow and overflow limits, the use of recurrence in connection with scaling, the use of machine constant routines and an error package for portability.
Research Organization:
Sandia National Labs., Albuquerque, NM (USA)
DOE Contract Number:
AC04-76DP00789
OSTI ID:
5249727
Report Number(s):
SAND-85-1018; ON: DE86000047
Country of Publication:
United States
Language:
English