zbMATH — the first resource for mathematics

Geometry Search for the term Geometry in any field. Queries are case-independent.
Funct* Wildcard queries are specified by * (e.g. functions, functorial, etc.). Otherwise the search is exact.
"Topological group" Phrases (multi-words) should be set in "straight quotation marks".
au: Bourbaki & ti: Algebra Search for author and title. The and-operator & is default and can be omitted.
Chebyshev | Tschebyscheff The or-operator | allows to search for Chebyshev or Tschebyscheff.
"Quasi* map*" py: 1989 The resulting documents have publication year 1989.
so: Eur* J* Mat* Soc* cc: 14 Search for publications in a particular source with a Mathematics Subject Classification code (cc) in 14.
"Partial diff* eq*" ! elliptic The not-operator ! eliminates all results containing the word elliptic.
dt: b & au: Hilbert The document type is set to books; alternatively: j for journal articles, a for book articles.
py: 2000-2015 cc: (94A | 11T) Number ranges are accepted. Terms can be grouped within (parentheses).
la: chinese Find documents in a given language. ISO 639-1 language codes can also be used.

a & b logic and
a | b logic or
!ab logic not
abc* right wildcard
"ab c" phrase
(ab c) parentheses
any anywhere an internal document identifier
au author, editor ai internal author identifier
ti title la language
so source ab review, abstract
py publication year rv reviewer
cc MSC code ut uncontrolled term
dt document type (j: journal article; b: book; a: book article)
Algorithm 609 a portable FORTRAN subroutine for the Bickley functions Ki n (x). (English) Zbl 0545.33001

The author discribes a FORTRAN subroutine for the computation of the so- called Bickley functions

Ki n (x)= x Ki n-1 (t)dt(x0,n=1,2,···)

where Ki 0 (x)=K 0 (x) is the modified Bessel function of the second kind. The algorithm consists of two main parts: For 0x2 it uses a four-term recurrence relation with starting values for n=0,1,2 obtained from a power series. For 2<x< a uniform asymptotic expansion is used for n0. The author also discusses the computation of the exponential integral E n (x), the gamma ratio R(x)=Γ(x)/Γ(x+1 2), and the digamma function ψ (x), all used as auxiliary functions. The description of the procedures is fundamentally mathematical and a careful error analysis as well as extensive tests are presented. The author emphazises that the main feature of these subroutines is that they are portable for an accuracy up to 18 digits. He states that, by using the special initialization routines provided, the user has a choice of 15 machine environments representing more than 19 machines and operating systems.

Reviewer: K.S.Kölbig
33-04Machine computation, programs (special functions)
33B15Gamma, beta and polygamma functions
33C10Bessel and Airy functions, cylinder functions, 0 F 1
33E99Other special functions
41A30Approximation by other special function classes