English

A numerical algorithm for computing the zeros of parabolic cylinder functions in the complex plane

Numerical Analysis 2025-03-27 v2 Numerical Analysis Classical Analysis and ODEs

Abstract

A numerical algorithm (implemented in Matlab) for computing the zeros of the parabolic cylinder function U(a,z)U(a,z) in domains of the complex plane is presented. The algorithm uses accurate approximations to the first zero plus a highly efficient method based on a fourth-order fixed point method with the parabolic cylinder functions computed by Taylor series and carefully selected steps, to compute the rest of the zeros. For a|a| small, the asymptotic approximations are complemented with a few fixed point iterations requiring the evaluation of U(a,z)U(a,z) and U(a,z)U'(a,z) in the region where the complex zeros are located. Liouville-Green expansions are derived to enhance the performance of a computational scheme to evaluate U(a,z)U(a,z) and U(a,z)U'(a,z) in that region. Several tests show the accuracy and efficiency of the numerical algorithm.

Keywords

Cite

@article{arxiv.2412.13085,
  title  = {A numerical algorithm for computing the zeros of parabolic cylinder functions in the complex plane},
  author = {T. M. Dunster and A. Gil and D. Ruiz-Antolín and J. Segura},
  journal= {arXiv preprint arXiv:2412.13085},
  year   = {2025}
}

Comments

To appear in BIT Numerical Mathematics

R2 v1 2026-06-28T20:39:07.892Z