Conical: an extended module for computing a numerically satisfactory pair of solutions of the differential equation for conical functions

Conical functions appear in a large number of applications in physics and engineering. In this paper we describe an extension of our module CONICAL for the computation of conical functions. Specifically, the module includes now a routine for computing the function ${{\rm R}}^{m}_{-\frac{1}{2}+i\tau}(x)$, a real-valued numerically satisfactory companion of the function ${\rm P}^m_{-\tfrac12+i\tau}(x)$ for $x>1$. In this way, a natural basis for solving Dirichlet problems bounded by conical domains is provided.Comment: To appear in Computer Physics Communication

Computing solutions of the modified Bessel differential equation for imaginary orders and positive arguments

We describe a variety of methods to compute the functions $K_{ia}(x)$, $L_{ia}(x)$ and their derivatives for real $a$ and positive $x$. These functions are numerically satisfactory independent solutions of the differential equation $x^2 w'' +x w' +(a^2 -x^2)w=0$. In an accompanying paper (Algorithm xxx: Modified Bessel functions of imaginary order and positive argument) we describe the implementation of these methods in Fortran 77 codes.Comment: 14 pages, 1 figure. To appear in ACM T. Math. Sof

On the computation of moments of the partial non-central chi-squared distribution function

Properties satisfied by the moments of the partial non-central chi-square distribution function, also known as Nuttall Q-functions, and methods for computing these moments are discussed in this paper. The Nuttall Q-function is involved in the study of a variety of problems in different fields, as for example digital communications.Comment: 6 page

Computation of the Marcum Q-function

Methods and an algorithm for computing the generalized Marcum $Q-$function ($Q_{\mu}(x,y)$) and the complementary function ($P_{\mu}(x,y)$) are described. These functions appear in problems of different technical and scientific areas such as, for example, radar detection and communications, statistics and probability theory, where they are called the non-central chi-square or the non central gamma cumulative distribution functions. The algorithm for computing the Marcum functions combines different methods of evaluation in different regions: series expansions, integral representations, asymptotic expansions, and use of three-term homogeneous recurrence relations. A relative accuracy close to $10^{-12}$ can be obtained in the parameter region $(x,y,\mu) \in [0,\,A]\times [0,\,A]\times [1,\,A]$, $A=200$, while for larger parameters the accuracy decreases (close to $10^{-11}$ for $A=1000$ and close to $5\times 10^{-11}$ for $A=10000$).Comment: Accepted for publication in ACM Trans. Math. Soft

Asymptotic approximations to the nodes and weights of Gauss-Hermite and Gauss-Laguerre quadratures

Asymptotic approximations to the zeros of Hermite and Laguerre polynomials are given, together with methods for obtaining the coefficients in the expansions. These approximations can be used as a standalone method of computation of Gaussian quadratures for high enough degrees, with Gaussian weights computed from asymptotic approximations for the orthogonal polynomials. We provide numerical evidence showing that for degrees greater than $100$ the asymptotic methods are enough for a double precision accuracy computation ($15$-$16$ digits) of the nodes and weights of the Gauss--Hermite and Gauss--Laguerre quadratures.Comment: Submitted to Studies in Applied Mathematic

On Non-Oscillating Integrals for Computing Inhomogeneous Airy Functions

Integral representations are considered of solutions of the inhomogeneous Airy differential equation $w''-z w=\pm1/\pi$. The solutions of these equations are also known as Scorer functions. Certain functional relations for these functions are used to confine the discussion to one function and to a certain sector in the complex plane. By using steepest descent methods from asymptotics, the standard integral representations of the Scorer functions are modified in order to obtain non-oscillating integrals for complex values of $z$. In this way stable representations for numerical evaluations of the functions are obtained. The methods are illustrated with numerical results.Comment: 12 pages, 5 figure