C312: Bessel Functions J and Y of Orders Zero and One

Author(s): K.S. Kölbig	Library: MATHLIB
Submitter:	Submitted: 18.10.1967
Language: Fortran	Revised: 15.03.1993

Function subprograms BESJ0, BESJ1, BESY0, BESY1 and DBESJ0, DBESJ1, DBESY0, DBESY1 calculate the Bessel functions

for real arguments x, where x>0 for and .

On CDC and Cray computers, the double-precision versions DBESJ0 etc. are not available.

Structure:

FUNCTION subprograms
User Entry Names: BESJ0, BESJ1, BESY0, BESY1, DBESJ0, DBESJ1, DBESY0, DBESY1
Files Referenced: Unit 6
External References: MTLMTR, ABEND

Usage:

In any arithmetic expression,

where BESJ0 etc. are of type REAL, DBESJ0 etc. are of type DOUBLE PRECISION, and X has the same type as the function name.

Method:

Approximation by truncated Chebyshev series.

Accuracy:

BESJ0 etc. (except on CDC and Cray computers) have full single-precision accuracy. For most values of the argument X, DBESJ0 etc. (and BESJ0 etc. on CDC and Cray computers) have an accuracy of approximately one significant digit less than the machine precision.

Error handling:

Error C312.1: for or . The function value is set equal to zero, and a message is written on Unit 6 unless subroutine MTLSET has been called.

References:

1. Y.L. Luke, Mathematical functions and their approximations (Academic Press, New York 1975) 322-324.