Abstract
A numerical evaluator for the confluent hypergeometric function for complex arguments with large magnitudes using a direct summation of Kummer's series is presented. Extended precision subroutines using large arrays to accumulate a single numerator and denominator are ultimately used in a single division to arrive at the final result. The accuracy has been verified through a variety of tests and they show the evaluator to be consistently accurate to 13 significant figures, and on rare occasion accurate to only 9 for magnitudes of the arguments ranging into the thousands in any quadrant in the complex plane. Because the evaluator automatically determines the number of significant figures of machine precision, and because it is written in FORTRAN 77, tests on various computers have shown the evaluator to provide consistently accurate results, making the evaluator very portable. The principal drawback is that, for certain arguments, the evaluator is slow; however, the evaluator remains valuable as a benchmark even in such cases.
Supplemental Material
Available for Download
confluent hypergeometric function Gams: c1
- 1 NADrN, M., PERGER, W. F., AND BHALLA, A. J. of Comput. Appl. Math. 39, (1992), 193 200. Google Scholar
- 2 ABRAMOW~TZ, M., AND STEOUN, I.A. Handbook of Mathematical Functions With Formulas, Graphs. and Mathematical Tables. U.S. Government Printing Office, Washington, D.C., 1972 390-413, 504-535, 546 553. Google Scholar
- 3 SLATER, L.J. Confluent Hypergeometric Functions. Cambridge University Press, London, 1960, 58-60.Google Scholar
Index Terms
- Algorithm 707: CONHYP: a numerical evaluator of the confluent hypergeometric function for complex arguments of large magnitudes
Recommendations
An algorithm for the rapid numerical evaluation of Bessel functions of real orders and arguments
We describe a method for the rapid numerical evaluation of the Bessel functions of the first and second kinds of nonnegative real orders and positive arguments. Our algorithm makes use of the well-known observation that although the Bessel functions ...
Orthogonal confluent hypergeometric lattice polynomials
Here is a method of solving the difference-differential equations of the confluent hypergeometric differential equation using a generalized Pochhammer matrix product. This method provides a convenient analytical way to relate various solutions of the ...
Comments