E. Cojocaru
Abstract:Routines for computation of Weber's parabolic cylinder functions and their derivatives are implemented in Matlab for both moderate and great values of the argument. Standard, real solutions are considered. Tables of values are included.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to implement and optimize the calculation of Parabolic Cylinder Functions (PCFs) and their derivatives in Matlab for the cases of medium and large parameter values. Specifically, the author aims to provide an efficient and accurate calculation method to facilitate the use of these special functions in scientific research and engineering applications.
### Background and Problem Description
Parabolic Cylinder Functions are a class of special functions in mathematical physics, originally introduced by Weber in 1869. They are widely used in fields such as quantum mechanics, electromagnetic theory, and statistical physics. However, for larger parameter values or independent variable values, direct calculation of these functions may encounter numerical stability problems, so special algorithms are required to ensure the accuracy and efficiency of the calculation.
### Main Contributions of the Paper
1. **Implementation of Standard Solutions**: The paper details how to implement the standard solutions \( U(a, x) \) and \( V(a, x) \) of the Parabolic Cylinder Function and their derivatives in Matlab.
2. **Handling of Large Parameter Values**: For the case of larger parameter values, the paper provides asymptotic expansion formulas to ensure the stability of numerical calculations.
3. **Handling of Complex Parameters**: The paper also involves the calculation of the gamma function in the case of complex parameters and provides the corresponding Matlab code implementation.
4. **Tabular Data**: To verify the correctness of the implementation, the paper provides multiple tables listing the function values for different parameter and independent variable values, which are convenient for users to compare and verify.
### Summary of Mathematical Formulas
The Parabolic Cylinder Function satisfies the following differential equations:
\[
\frac{d^2 y}{dx^2}-\left(\frac{1}{4}x^2 + a\right)y = 0\quad(1)
\]
\[
\frac{d^2 y}{dx^2}+\left(\frac{1}{4}x^2 - a\right)y = 0\quad(2)
\]
The standard solutions are:
\[
U(a, x)=D_{-a-\frac{1}{2}}
\]
\[
V(a, x)=\frac{1}{\pi\Gamma(\frac{1}{2}+a)}\left[\sin(\pi a)U(a, x)+U(a, -x)\right]
\]
For larger independent variable \( x \), there are asymptotic expansions:
\[
U(a, x)\sim x^{-a-\frac{1}{2}}e^{-\frac{x^2}{4}}\left[1-\frac{(a + \frac{1}{2})(a+\frac{3}{2})}{2x^2}+\cdots\right]
\]
\[
V(a, x)\sim\sqrt{\frac{2}{\pi}}x^{a-\frac{1}{2}}e^{\frac{x^2}{4}}\left[1+\frac{(a-\frac{1}{2})(a-\frac{3}{2})}{2x^2}+\cdots\right]
\]
Through these formulas and implementation methods, the paper successfully solves the problem of efficient calculation of Parabolic Cylinder Functions in Matlab, providing a powerful tool for research in related fields.