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-5zbMATH Open1311.65024arXiv1403.7927OpenAlexW2139077355WikidataQ115381616 ScholiaQ115381616MaRDI QIDQ2017615FDOQ2017615

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)

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.

Full work available at URL: https://arxiv.org/abs/1403.7927

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

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)

Cites Work

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

Cited In (3)

• 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
• An improved algorithm and a Fortran 90 module for computing the conical function \(P^m_{-1/2+i\tau}(x)\)

Uses Software

• Digital Library of Mathematical Functions
• 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)