next up previous index
Next: C341 Modified Bessel Up: CERNLIB Previous: C339 Dawson's Integral

C340 Modified Bessel Functions I and K of Order 1/3 and 2/3

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

Function subprograms BSIR3, BSKR3 and DBSIR3, DBSKR3 calculate the modified Bessel functions

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

for real arguments x>0 and ν= -2,-1,1,2 . 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 DBSIR3 etc. are not available.

Structure:

FUNCTION subprograms
User Entry Names: BSIR3, BSKR3, EBSIR3, EBSKR3, DBSIR3, DBSKR3, DEBIR3, DEBKR3
Files Referenced: Unit 6
External References: MTLMTR (N002), ABEND (Z035)

Usage:

In any arithmetic expression, BSIR3(X,NU) or DBSIR3(X,NU) has the value INU/3(X) ,
BSKR3(X,NU) or DBSKR3(X,NU) has the value KNU/3(X) ,
EBSIR3(X,NU) or DEBIR3(X,NU) has the value exp(-X) * INU/3(X) ,
EBSKR3(X,NU) or DEBKR3(X,NU) has the value exp(X) * KNU/3(X) , where BSIR3 etc. are of the type REAL, DBSIR3 etc. are of the type DOUBLE PRECISION, and X has the same type as the function name. NU (INTEGER) has one of the values -2,-1,1,2.

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.

Accuracy:

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

Error handling:

Error C340.1: X ≤0 or X<0 , repectively, or NU ≠-2,-1,1,2 .
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) 352, 355, 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.

C341

K.S. Kölbig



next up previous index
Next: C341 Modified Bessel Up: CERNLIB Previous: C339 Dawson's Integral


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