besselK_app_ser(z, lambda, details = TRUE)
z
is sufficiently large. Otherwise it
approximates this value. The function returns: an accurate
approximated value of $K_{\lambda}(z)$ (up to 15 dp) and a
finite number of terms in the summation. By using this function
the mathematical difficulties due to the presence of
$K_{\lambda}(z)$ in the density function of the
generalized hyperbolic and related distributions
can be dealt with in a much simpler way. If
z >=705
then approximated value is on an
exponential scale to avoid underflow.besselK_app_ser(73, 9.12, details = TRUE)
besselK_app_ser(1, 5, details = TRUE)
Run the code above in your browser using DataLab