…
►Minimax polynomial approximations (§3.11(i)) for and in terms of with can be found in Abramowitz and Stegun (1964, §17.3) with maximum absolute errors ranging from 4×10⁻⁵ to 2×10⁻⁸.
Approximations of the same type for and for are given in Cody (1965a) with maximum absolute errors ranging from 4×10⁻⁵ to 4×10⁻¹⁸.
Cody (1965b) gives Chebyshev-series expansions (§3.11(ii)) with maximum precision 25D.
…
►The accuracy is controlled by the number of terms retained in the approximation; for real variables the number of significant figures appears to be roughly twice the number of terms retained, perhaps even for near with the improvements made in the 1970 reference.
…
►Figure 22.3.25:
as a function of complex , , .
…
Magnify3DHelp►►►Figure 22.3.26: Density plot of as a function of complex , , .
…
Magnify►►►Figure 22.3.27: Density plot of as a function of complex , , .
…
Magnify
…
A. R. DiDonato and A. H. Morris (1987)Algorithm 654: Fortran subroutines for computing the incomplete gamma function ratios and their inverses.
ACM Trans. Math. Software13 (3), pp. 318–319.
A. R. DiDonato and A. H. Morris (1992)Algorithm 708: Significant digit computation of the incomplete beta function ratios.
ACM Trans. Math. Software18 (3), pp. 360–373.
G. W. Hill (1981)Algorithm 571: Statistics for von Mises’ and Fisher’s distributions of directions: , and their inverses [S14].
ACM Trans. Math. Software7 (2), pp. 233–238.
…
►Inside the cusp, that is, for , the zeros form pairs lying in curved rows.
…
►Just outside the cusp, that is, for , there is a single row of zeros on each side.
…
►Near , and for small and , the modulus
has the symmetry of a lattice with a rhombohedral unit cell that has a mirror plane and an inverse threefold axis whose and repeat distances are given by
…The rings are almost circular (radii close to and varying by less than 1%), and almost flat (deviating from the planes by at most ).
…There are also three sets of zero lines in the plane related by rotation; these are zeros of (36.2.20), whose asymptotic form in polar coordinates is given by
…
…
►Similar definitions apply in the case : we transform into the origin by replacing in (16.8.1) by ; again compare §2.7(i).
…
►Equation (16.8.3) is of order .
…
►
…
►Let with and .
…The integers , , and are characteristics of the machine.
…
►In the case of the normalized binary interchange formats, the representation of data for binary32 (previously single precision) (, , , ), binary64 (previously double precision) (, , , ) and binary128 (previously quad precision) (, , , ) are as in Figure 3.1.1.
…
►
, and
…Then rounding by chopping or rounding down of gives , with maximum relative error .
…