Calculation of Fresnel integrals of real and complex arguments up to 28 significant digits

Mofreh R. Zaghloul, Leen Alrawas

Research output: Contribution to journalArticlepeer-review

Abstract

We present numerical algorithms for efficient and highly accurate computations of Fresnel’s sine and cosine integrals S(z) and C(z) for real (z = x) and complex (z = x + iy) arguments. The algorithms are based on series expansion for small values of |z|, expansion in Chebyshev subinterval polynomial approximations for intermediate real values (y = 0), toggled with sum of half-integer order Bessel functions for complex arguments, together with asymptotic series expansion for large values of |z|. The present algorithms, implemented in a Fortran elemental module, can be run using any of the single, double, or quadruple precision arithmetic. Results from the present code have been benchmarked versus comprehensive data tables, generated using Maple, Mathematica, and Matlab symbolic toolbox. Compared to a maximum of 16 significant digits in the literature, the present algorithms can calculate S(x) and C(x) up to 28 significant digits in the range xϵ[0, 106] and S(z) and C(z) to the same accuracy in the domain where |S(z)| and |C(z)| each is less than the largest finite floating point number in the precision under consideration with ∣argz∣<π2 or more specifically 3×10-4<yx<2×103 for S(z) and 4×10-4<yx<2.2×103 for C(z).

Original languageEnglish
Pages (from-to)489-506
Number of pages18
JournalNumerical Algorithms
Volume96
Issue number2
DOIs
Publication statusPublished - Jun 2024

Keywords

  • Fortran
  • Fresnel integrals
  • Special functions

ASJC Scopus subject areas

  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Calculation of Fresnel integrals of real and complex arguments up to 28 significant digits'. Together they form a unique fingerprint.

Cite this