next up previous index
Next: C328 Whittaker Function Up: CERNLIB Previous: C326 Clausen Function

C327 Modified Bessel Functions I and K of Order 1/4, 1/2 and 3/4

Routine ID: C327
Author(s): K.S. KölbigLibrary: MATHLIB
Submitter: Submitted: 15.05.1987
Language: FortranRevised: 15.03.1993

Function subprograms BSIR4, BSKR4 and DBSIR4, DBSKR4 calculate the modified Bessel functions

Iν/4(x)&quad;and&quad;Kν/4(x)

for real arguments x>0 and ν=-3,-2,-1,1,2,3 . The value x=0 is permitted for the functions I if ν> 0 . Note that the functions K are even with respect to ν .

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

Structure:

FUNCTION subprograms
User Entry Names: BSIR4, BSKR4,EBSIR4, EBSKR4, DBSIR4, DBSKR4, DEBIR4, DEBKR4
Files Referenced: Unit 6
External References: MTLMTR (N002), ABEND (Z035)

Usage:

In any arithmetic expression, BSIR4(X,NU) or DBSIR4(X,NU) has the value INU/4(X) ,
BSKR4(X,NU) or DBSKR4(X,NU) has the value KNU/4(X) ,
EBSIR4(X,NU) or DEBIR4(X,NU) has the value exp(-X) * INU/4(X) ,
EBSKR4(X,NU) or DEBKR4(X,NU) has the value exp(X) * KNU/4(X) , where BSIR4 etc. are of the type REAL, DBSIR4 etc. are of the type DOUBLE PRECISION, and X has the same type as the function name. NU is of type INTEGER and must have one of the values -3,-2,-1,1,2,3.

Method:

Approximation by rational functions (I for |x|<8, K for 1 ≤x ≤5 ), by an algorithm based on power series (K for 0 < x < 1), or else by truncated Chebyshev series. The cases |ν|=2 are elementary.

Accuracy:

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

Error handling:

Error C327.1: X ≤0 , or X<0 , respectively, or NU -3,-2,-1,1,2,3. The function value is set equal to zero, and a message is written on Unit 6, unless subroutine MTLSET (N002) has been called.

References:

  1. Y.L. Luke, Mathematical functions and their approximations (Academic Press, New York 1975) 350, 357, 363, 366.
  2. N.M. Temme, On the numerical evaluation of the modified Bessel function of the third kind, J. Comp. Phys. 19 (1975) 324--337.

C328



next up previous index
Next: C328 Whittaker Function Up: CERNLIB Previous: C326 Clausen Function


Janne Saarela
Mon Apr 3 15:06:23 METDST 1995