Abstract
Fortran 90 programs for the computation of real parabolic cylinder functions are presented. The code computes the functions U(a, x), V(a, x) and their derivatives for real a and x (x ≥ 0). The code also computes scaled functions. The range of computation for scaled PCFs is practically unrestricted. The aimed relative accuracy for scaled functions is better than 5 10−14. Exceptions to this accuracy are the evaluation of the functions near their zeros and the error caused by the evaluation of trigonometric functions of large arguments when |a| > x. The routines always give values for which the Wronskian relation for scaled functions is verified with a relative accuracy better than 5 10−14. The accuracy of the unscaled functions is also better than 5 10−14 for moderate values of x and a (except close to the zeros), while for large x and a the error is dominated by exponential and trigonometric function evaluations. For IEEE standard double precision arithmetic, the accuracy is better than 5 10−13 in the computable range of unscaled PCFs (except close to the zeros).
Supplemental Material
Available for Download
Software for "Real parabolic cylinder functions U(a, x), V(a, x)"
- Gil, A., Segura, J., and Temme, N. M. 2004. Integral representations for computing real parabolic cylinder functions. Numer. Math. 98, 105--134. Google Scholar
- Gil, A. Segura, J., and Temme, N. M. Computing the Parabolic Cylinder Functions U(a, x), V(a, x). ACM Trans. Math. Soft., (in this issue). Google Scholar
- Segura, J. and Gil, A. 1998. Parabolic cylinder functions of integer and half-integer orders for non-negative arguments. Comput. Phys. Commun. 115, 69--86.Google Scholar
Index Terms
- Algorithm 850: Real parabolic cylinder functions U(a, x), V(a, x)
Recommendations
Algorithm 914: Parabolic cylinder function W(a, x) and its derivative
A Fortran 90 program for the computation of the real parabolic cylinder functions W(a, ± x), x ≥ 0 and their derivatives is presented. The code also computes scaled functions for a > 50. The functions W(a, ± x) are a numerically satisfactory pair of ...
Computing the real parabolic cylinder functions U(a, x), V(a, x)
Methods for the computation of real parabolic cylinder functions U(a, x), and V(a, x) and their derivatives are described. We give details on power series, asymptotic series, recursion and quadrature. A combination of these methods can be used for ...
Compact Numerical Quadrature Formulas for Hypersingular Integrals and Integral Equations
In the first part of this work, we derive compact numerical quadrature formulas for finite-range integrals $I[f]=\int^{b}_{a}f(x)\,dx$ , where f ( x )= g ( x )| x t | β , β being real. Depending on the value of β , these integrals are defined either in the regular ...
Comments