This paper describes an improved method for computing Fresnel integrals with an accuracy better than 1X10 -9, based on a known approximate formula for a different integral. Outside the range of validity of this formula, Taylor series expansions are used to