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.