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 $$\mid \mathit{\arg}\ z\mid <\frac{\pi }{2}$$
∣
arg
z
∣
<
π
2
or more specifically $$3\times {10}^{-4}<\left|\frac{y}{x}\right|<2\times {10}^3$$
3
×
10
-
4
<
y
x
<
2
×
10
3
for S(z) and $$4\times {10}^{-4}<\left|\frac{y}{x}\right|<2.2\times {10}^3$$
4
×
10
-
4
<
y
x
<
2.2
×
10
3
for C(z).