Moshier (1989, §6.14) provides minimax rational approximations
for calculating , , , .
They are in terms of the variable
, where
when is positive,
when is negative,
and when .
The approximations apply when , that is,
when or .
The precision in the coefficients is 21S.
…
►These expansions are for real arguments and are supplied in sets of four for each function, corresponding to intervals , , , .
…
►
•
Corless et al. (1992) describe a method of approximation based on
subdividing into a triangular mesh, with values of ,
stored at the nodes. and are then
computed from Taylor-series expansions centered at one of the nearest nodes.
The Taylor coefficients are generated by recursion, starting from the stored
values of ,
at the node. Similarly for
, .
…
►For uniform asymptotic approximations for the zeros of in the interval when with
fixed, see Olver (1997b, p. 469).
►
§14.16(iii) Interval
►
has exactly one zero in the interval if either of the following sets of conditions holds:
…
►For all other values of and (with ) has no zeros in the interval .
►
has no zeros in the interval when , and at most one zero in the interval when .
…
►This is a multivalued function of with branch point at
.
…
►
is a single-valued analytic function on and real-valued when ranges over the positive real numbers.
…
►
…
►The principal value is
…This is an analytic function of on , and is two-valued and discontinuous on the cut shown in Figure 4.2.1, unless .
…
…
►Figure 4.3.2 illustrates the conformal mapping of the strip onto the whole -plane cut along the negative real axis, where and (principal value).
…Lines parallel to the real axis in the -plane map onto rays in the -plane, and lines parallel to the imaginary axis in the -plane map onto circles centered at the origin in the -plane.
In the labeling of corresponding points is a real parameter that can lie anywhere in the interval .
…
►In the graphics shown in this subsection height corresponds to the absolute value of the function and color to the phase.
…
►
…
►Table 22.5.1 gives the value of each of the 12 Jacobian elliptic functions, together with its -derivative (or at a pole, the residue), for values of that are integer multiples of , .
…
►Table 22.5.2 gives , , for other special values of .
…
►
…
►If the results agree within significant figures, then it is likely—but not certain—that the truncated asymptotic series will yield at least correct significant figures for larger values of .
…
►In this way we arrive athyperasymptotic expansions.
…
►17408, compared with the correct value
…
►17045, which is much closer to the true value.
…
►Comparison with the true value
…
…
►These constraints guarantee that the orthogonality only involves the integral , as above.
►For other cases there may also be, in addition to a possible integral as in (18.30.10), a finite sum of discrete weights on the negative real -axis each multiplied by the polynomial product evaluated at the corresponding values of , as in (18.2.3).
…
►
…
►Assume that again has the expansion (2.3.7) and this expansion is infinitely differentiable, is infinitely differentiable on , and each of the integrals , , converges at
, uniformly for all sufficiently large .
…
►Assume also that and are continuous in and , and for each the minimum value of in is at
, at which point vanishes, but both and are nonzero.
When Laplace’s method (§2.3(iii)) applies, but the form of the resulting approximation is discontinuous at
.
…
►
being the value of
at
.
We now expand in a Taylor series centered at the peak value
of the exponential factor in the integrand:
…
…
►
…
►
…
►In Figures 19.3.7 and 19.3.8 for complete Legendre’s elliptic integrals with complex arguments, height corresponds to the absolute value of the function and color to the phase.
…
►
…