If a>1
and 0<h<=cut.point
, a series expansion is used,
truncated after jmax
terms.
If a>1
and h>cut.point
, an asymptotic approximation is used.
In the other cases, various reflection properties of the function
are exploited. See the reference below for more information.