## Unrestricted algorithms for elementary and special functions

52. R. P. Brent,
Unrestricted algorithms for elementary and special functions,
in * Information Processing 80 * (edited by S. H. Lavington),
North-Holland, Amsterdam, 1980, 613-619.
CR 22#38728,
MR 81i:68009.
Retyped with minor corrections 1999.
arXiv:1004.3621v1
Abstract:
dvi (2K),
pdf (75K),
ps (25K).

Paper:
dvi (23K)
pdf (176K)
ps (144K).

## Abstract

We describe some "unrestricted" algorithms which are useful for the
computation of elementary and special functions when the precision
required is not known in advance. Several general classes of algorithms
are identified and illustrated by examples.
The topics include:
power series methods,
use of halving identities,
asymptotic expansions,
continued fractions,
recurrence relations,
Newton's method,
numerical contour integration,
and the arithmetic-geometric mean.
## Comments

Most of the algorithms discussed are implemented in the
MP package.
For related work see Brent
[28, 34].
## Errata for original version

Page 615, remove the absolute value signs in
equation (13) and the following paragraph (*x* should be large
and positive here).
Page 616, second half of equation (26):
insert a minus sign after the equals sign.

Page 617, equation (39):
delete " / j ! ".

Page 617, equation (42):
the assumption "*j* < *k*" should be added. Also, the
contour *C* needs to be enlarged slightly.

Page 617, left-hand-side of equation (44):
replace "S_{j,k}" by "S_{2j,k}".

Page 617, ten lines after equation (44):
replace "O(jn^{2})" by "O(j^{2}n)".

Go to next publication

Return to Richard Brent's index page