J. Segura

T. M. Dunster, A. Gil, J. Segura

Linear second order differential equations having a large real parameter and turning point in the complex plane are considered. Classical asymptotic expansions for solutions involve the Airy function and its derivative, along with two infinite series, the coefficients of which are usually difficult to compute. By considering the series as asymptotic expansions for two explicitly defined analytic functions, Cauchy's integral formula is employed to compute the coefficient functions to high order of accuracy. The method employs a certain exponential form of Liouville-Green expansions for solutions of the differential equation, as well as for the Airy function. We illustrate the use of the method with the high accuracy computation of Airy-type expansions of Bessel functions of complex argument.

A. Gil, J. Segura, N. M. Temme

A Fortran 90 module (GammaCHI) for computing and inverting the gamma and chi-square cumulative distribution functions (central and noncentral) is presented. The main novelty of this package are the reliable and accurate inversion routines for the noncentral cumulative distribution functions. Additionally, the package also provides routines for computing the gamma function, the error function and other functions related to the gamma function. The module includes the routines cdfgamC, invcdfgamC, cdfgamNC, invcdfgamNC, errorfunction, inverfc, gamma, loggam, gamstar and quotgamm for the computation of the central gamma distribution function (and its complementary function), the inversion of the central gamma distribution function, the computation of the noncentral gamma distribution function (and its complementary function), the inversion of the noncentral gamma distribution function, the computation of the error function and its complementary function, the inversion of the complementary error function, the computation of: the gamma function, the logarithm of the gamma function, the regulated gamma function and the ratio of two gamma functions, respectively.

A. Gil, J. Segura, N. M. Temme

The computation and inversion of the noncentral beta distribution $B_{p,q}(x,y)$ (or the noncentral $F$-distribution, a particular case of $B_{p,q}(x,y)$) play an important role in different applications. In this paper we study the stability of recursions satisfied by $B_{p,q}(x,y)$ and its complementary function and describe asymptotic expansions useful for computing the function when the parameters are large. We also consider the inversion problem of finding $x$ or $y$ when a value of $B_{p,q}(x,y)$ is given. We provide approximations to $x$ and $y$ which can be used as starting values of methods for solving nonlinear equations (such as Newton) if higher accuracy is needed.

A. Gil, J. Segura, N. M. Temme

An efficient algorithm and a Fortran 90 module (LaguerrePol) for computing Laguerre polynomials $L^{(α)}_n(z)$ are presented. The standard three-term recurrence relation satisfied by the polynomials and different types of asymptotic expansions valid for $n$ large and $α$ small, are used depending on the parameter region. Based on tests of contiguous relations in the parameter $α$ and the degree $n$ satisfied by the polynomials, we claim that a relative accuracy close or better than $10^{-12}$ can be obtained using the module LaguerrePol for computing the functions $L^{(α)}_n(z)$ in the parameter range $z \ge 0$, $-1 < α\le 5$, $n \ge 0$.

A. Gil, J. Segura, N. M. Temme

We describe methods for computing the Kummer function $U(a,b,z)$ for small values of $z$, with special attention to small values of $b$. For these values of $b$ the connection formula that represents $U(a,b,z)$ as a linear combination of two ${}_1F_1$-functions needs a limiting procedure. We use the power series of the ${}_1F_1$-functions and consider the terms for which this limiting procedure is needed. We give recursion relations for higher terms in the expansion, and we consider the derivative $U^\prime(a,b,z)$ as well. We also discuss the performance for small $\vert z\vert$ of an asymptotic approximation of the Kummer function in terms of modified Bessel functions.