Computation of a numerically satisfactory pair of solutions of the differential equation for conical functions of non-negative integer orders

DOI10.1007/S11075-014-9857-5MaRDI QIDQ2017615zbMATH OpenOpenAlexWikidataFDO

Authors T. M. Dunster, Amparo Gil, Javier Segura, Nico M. Temme

Publication date 23 March 2015

Published in Numerical Algorithms (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1403.7927, http://hdl.handle.net/10902/10202

zbMATH Keywords

algorithmLegendre functionsthree-term recurrence relationsconical functions

Mathematics Subject Classification ID

Computation of special functions and constants, construction of tables (65D20) Orthogonal polynomials and functions of hypergeometric type (Jacobi, Laguerre, Hermite, Askey scheme, etc.) (33C45) Numerical approximation and evaluation of special functions (33F05)

Abstract: We consider the problem of computing satisfactory pairs of solutions of the differential equation for Legendre functions of non-negative integer order $mu$ and degree $-frac12+iau$, where $au$ is a non-negative real parameter. Solutions of this equation are the conical functions ${mP}_{-frac12+iau}^{mu}(x)$ and $Q_{-frac12+iau}^{mu}(x)$, $x>-1$. An algorithm for computing a numerically satisfactory pair of solutions is already available when $-1<x<1$ (see cite{gil:2009:con}, cite{gil:2012:cpc}).In this paper, we present a stable computational scheme for a real valued numerically satisfactory companion of the function ${mP}_{-frac12+iau}^{mu}(x)$ for $x>1$, the function $Releft{e}^{-ipimu}{Q}_{-frac12+iau}^{mu}(x)ight$. The proposed algorithm allows the computation of the function on a large parameter domain without requiring the use of extended precision arithmetic.

Recommendations

• Conical: an extended module for computing a numerically satisfactory pair of solutions of the differential equation for conical functions
• Computing the conical function \(P^{\mu}_{-1/2+i\tau}(x)\)
• An improved algorithm and a Fortran 90 module for computing the conical function \(P^m_{-1/2+i\tau}(x)\)
• On numerical algorithms for real functions of confluent hypergeometric type
• Computing solutions of the modified Bessel differential equation for imaginary orders and positive arguments

Cites work

• scientific article; zbMATH DE number 884971 (Why is no real title available?)
• An improved algorithm and a Fortran 90 module for computing the conical function \(P^m_{-1/2+i\tau}(x)\)
• Computing the conical function \(P^{\mu}_{-1/2+i\tau}(x)\)
• Conical functions with one or both parameters large
• Numerical Methods for Special Functions
• The numerical computation of the confluent hypergeometric function U(a,b,z)

Cited in

(3)

• An improved algorithm and a Fortran 90 module for computing the conical function \(P^m_{-1/2+i\tau}(x)\)
• Computing the conical function \(P^{\mu}_{-1/2+i\tau}(x)\)
• Conical: an extended module for computing a numerically satisfactory pair of solutions of the differential equation for conical functions

Describes a project that uses

Uses Software

• DLMF
• Conical

This page was built for publication: Computation of a numerically satisfactory pair of solutions of the differential equation for conical functions of non-negative integer orders

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2017615)