# Approximation of Fresnel Integrals with Applications to Diffraction Problems

Authors: M. Sandoval-Hernandez; H. Vazquez-Leal; L. Hernandez-Martinez; U. A. Filobello-Nino; V. M. Jimenez-Fernandez; A. L. Herrera-May; R. Castaneda-Sheissa; R. C. Ambrosio-Lazaro; G. Diaz-Arango
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

This article introduces two approximations that allow the evaluation of Fresnel integrals without the need for using numerical algorithms. These equations accomplish the characteristic of being continuous in the same interval as Fresnel. Both expressions have been determined applying the least squares method to suitable expressions. Accuracy of equations improves asx increases; as for small values of x, it is possible to achieve an absolute error less than 8×10-5. To probe the efficiency of the equations, two case studies are presented, both applied in the optics field. The first case is related to the semi-infinite opaque screen for Fresnel diffraction. In this case study Fresnel integrals are evaluated with the proposed equations to calculate the irradiance distribution and the Cornu spiral for diffraction computations of the Fresnel diffraction; obtained results show a good accuracy. The second case is related to the double aperture problem for Fresnel diffraction.


---

## Body

## 1. Introduction

Fresnel integrals are two transcendental functions named after Augustin-Jean Fresnel; these are widely used in physics, specially in optics and electromagnetic theory [1, 2]. For instance, in optics, solution of the Fresnel diffraction for a rectangular aperture is calculated, usually, in terms of integrals that cannot be obtained by analytical methods [3–6]. These integrals are known as Fresnel integrals and contain arguments in terms of the rectangular coordinate system (Cartesian coordinate system).In [7] the iterative Fresnel integral method is introduced for numerical calculation of the Fresnel diffraction patterns from rectangular apertures tilted at an arbitrary angle to the optical axis. The work presented in [8] demonstrated that the propagation rate of truncated, nondiffracting, and accelerating beams could be accurately calculated using a method that represents the beams by a finite superposition of Gaussian [1] to be in turn propagated by means of the Fresnel diffraction integral. In [9], the iterative Fresnel integrals method (IFIM) has been applied, for the first time, to the simulation and generation of the complete near-field Fresnel diffraction images crated by N-apertures. The article [2] provides a development for the double Fresnel integral based on adequate variable changes of the problem under study to predict electromagnetic reduction of the double razor edge with reflection on the earth surface. In [10], Fresnel integrals were applied to obtain the derivation of linear FM Pulse Compression Spectra.For some years proposals to evaluate these integrals have been provided. The work from [11] introduces polynomial approximations to evaluate Fresnel integrals. These approximations are applied in diffraction analysis [12]. An algorithm to evaluate Fresnel integrals based on the continued fractions method is provided in [13]. On the other hand, to evaluate Fresnel integrals for high values of |x|, the solution is resort to the expansion in terms of Bessel functions [14, 15]. Besides, in [16], one of the performed approximations was made by the sum of the first term of the Fresnel integrals asymptotic expansion and an exponential function which approaches infinity at the zero of Fresnel function argument.Likewise, several different numerical alternatives have been proposed to evaluate Fresnel integrals. [17] presented the Chebyshev approximations for the Fresnel integrals which employ different values of the argument at different evaluation intervals, depending on its magnitude. In [18], a method for spreadsheet computations of Fresnel integrals to six significant figures is introduced; it is based on successive improvements of known relational approximations which are accurate to just three figures.This work presents two expressions to calculate, analytically, the Fresnel sine integral and Fresnel cosine integral; with our approximations it is possible to evaluate Fresnel integrals in a simple way without the need for using numerical algorithms like Simpson’s rule, among others. These expressions can be evaluated using any calculus software or standard scientific calculator. In this article we assume that exact Fresnel functions are the ones provided by Maple17.The paper is organized as follows. Section2 provides an introduction to the Fresnel integrals. In Section 3, the least squares method is presented. Deduction of equations that approximate Fresnel integrals is performed in Section 4. Study cases are introduced in Section 5. Section 6 presents the discussion about the obtained results and efficiency of the equations. Finally, conclusions of this work are provided in Section 7.


## 2. Fresnel Integrals

The Fresnel sine and cosine integrals are represented byS(x) and C(x), respectively; these are two integral functions that originate by applying the analysis of Fresnel diffraction phenomena, which is defined by the following integrals:(1)CP=∫0Pcos⁡πx22dx,(2)SP=∫0Psin⁡πx22dx,where P≥0 is a real number and x is a real variable. When P tends to infinite, C=S=1/2 [19].Fresnel integrals (C(x) and S(x)) are odd functions of x. They can be extended to the complex numbers domain; thus, it is possible to obtain analytic functions for one variable [20].(3)Cx=π4ierf⁡ix+-ierf⁡-ix,(4)Sx=π4-ierf⁡ix+ierf⁡-ix.Figure1 shows the Fresnel integrals graph.Figure 1
Fresnel integrals.The Cornu spiral, also known as clothoid or Euler spiral, is the parametric curve generated by the Fresnel integralsS(x) and C(x) from (1) and (2). The spiral is a tangent curve to the abscissa axis at the origin and its radius of curvature decreases inversely in proportion to the distance travelled over it (see Figure 2).Figure 2
Cornu spiral.In the optics field, the Cornu spiral can be employed to obtain a quantitative image from a diffraction pattern. In this way, when the diffraction is observed from a position close to the diffractor obstacle, it is referred to as near-field diffraction or Fresnel zone. Therefore, in order to perform the diffraction pattern analysis, the semi-periodic Fresnel zones are calculated from the observation point. In consequence, the contribution from each zone is a phasor that adds up. At the limit of the infinitesimal width zones the Cornu spiral is obtained.The Cornu spiral [1] has the property that its curvature at any point is proportional to the distance along the curvature from the origin. This property allows using it as transition curve for tracing highways and railroads, given the fact that a mobile that has this displacement at constant speed will have a constant angular acceleration [21, 22]. The transition curve has infinite radius at the tangent point to the straight part of the path, radius R at the tangency point to the uniform circular curve; therefore, the most common type of curves on highways is straight stretch-clothoid-circular-clothoid-straight stretch [21, 22].


## 3. Least Squares Approximations

Let (x0,y0), (x1,y1), (x2,x3), …, (xi,yi) be the coordinates of a data set, such that i=0,1,2,…,n, and the adjustment curve be y=f(x,a0,a1,a2,…,aj), where aj (j=0,1,2,…) are adjustment constants. The least squares approximations attempt to minimize the sum of the squares from the vertical distances of yi values to the ideal model f(x) and obtain the model function S(a0,a1,a2,…,aj), which minimizes the square error defined by(5)Sa0,a1,a2,…,aj=∑i=1nyi-fxi,a0,a1,a2,…,aj2.Therefore, (5) is designed to minimize the error with respect to the sample set. To do this, partial derivatives are defined with respect to every adjustment constant, creating a system of nonlinear equations given as(6)∂S∂a0=∂∂a0∑i=1nyi-fxi,a0,a1,a2,…,aj2=0∂S∂a0=∂∂a0∑i=1nyi-fxi,a0,a1,a2,…,aj2=0∂S∂a1=∂∂a1∑i=1nyi-fxi,a0,a1,a2,…,aj2=0∂S∂a2=∂∂a2∑i=1nyi-fxi,a0,a1,a2,…,aj2=0⋮∂S∂an=∂S∂an∑i=1nyi-fxi,a0,a1,a2,…,aj2=0.Solving system (6) for adjustment constants aj, it is possible to obtain a model that provides the best fit to the data set. This method allows creating a continuous function in the space with a minimized error with respect to the samples [23]. In consequence, it is possible to test various models and select the one that provides the best fit. Figure 3 depicts the methodology. Notice that in Figure 3 the asterisk symbol represents the data set. For every data element a value (xi,yi) is given, where i=1,2,3,…. It is possible to obtain the best function that fits the data set, represented by the blue line.Figure 3
Least squares method.


## 4. Deduction of Approximations

This article proposes two equations which allow the numerical evaluation of the sine and cosine integrals (C(x) and S(x)) given by (1) and (2), respectively.(7)Cx=-sin⁡πx2/2πx+20πexp⁡-200πx+8251-exp⁡-Aπx3+2251-exp⁡-92πx2+1101-exp⁡-Bπxsgnx(8)Sx=-cos⁡πx2/2πx+Aπexp⁡Bπx+8251-exp⁡-Cπx3+2251-exp⁡-Dπx2+1101-exp⁡-92πxsgnxNotice that in (7) and (8) the first term decreases at the rate of 1/πx [13]. In consequence, the first term in these equations dictates the behaviour for each one since it acts as the characteristic term of a solution just like the solution of a differential equation [5]. The following terms represent the stationary state, when approaching infinity tends to 1/2. Using as adjustment functions (7) and (8) and the NonlinearFit command in Maple 17, model parameters are adjusted to the data set reducing the error by least squares. In the procedure of obtaining the approximate models, the data set (xi, yi) consists of 300 elements, obtained from the exact solutions C(x) and S(x) for the 0≤x≤3 interval; therefore, the resultant equations are(9)Cx=-sin⁡πx2/2πx+20πexp⁡-200πx+8251-exp⁡-69100πx3+2251-exp⁡-92πx2+1101-exp⁡-1.55294068198794πxsgnx,(10)Sx=-cos⁡πx2/2πx+16.7312774552827πexp⁡-1.57638860756614πx+8251-exp⁡-0.608707749430681πx3+2251-exp⁡-1.71402838165388πx2+1101-exp⁡-910πxsgnx.In Figures4 and 5 the absolute error is presented for (9) and (10), respectively. The maximum error for (9) is about 0.0297629 when x=1.4. It can be noticed in Figure 4 that accuracy improves as x value increases.Figure 4
Absolute error forCx in the interval 0≤x≤7.Figure 5
Absolute error forS(x) in the interval 0≤x≤7.For (10), the maximum absolute error occurs when x=1.62 and is 0.0207910487. Likewise, in Figure 5 it can be noticed that the error for (10) tends to decrease as x value increases.In [18] two approximations are proposed to evaluate the sine and cosine integrals; there are(11)Cx=12+1+0.926x2+1.792x+3.104x2+εxsin⁡πx22-12+4.142x+3.492x2+6.67x3+εxcos⁡πx22(12)Sx=12-1+0.926x2+1.792x+3.104x2+εxcos⁡πx22-12+4.142x+3.492x2+6.67x3+εxsin⁡πx22where εx≤2×10-3 and x≥0.As a way to compare the accuracy of our approximations, Figure6 depicts the behaviour of the absolute error between (9) and (11), which approximate the cosine integral.Figure 6
Absolute error comparison between (9) and (11).
(a)
Absolute error comparison in the interval0≤x≤8 (b)
Absolute error comparison in the interval8≤x≤12 (c)
Absolute error comparison in the interval30≤x≤30.5Initially the absolute error in (9) is higher than that in (11), but tends to decrease considerably. At x<6.5, the error is lower than (11) (see Figure 6(a)). By x=8, the absolute error for (11) is about 10-3, while that for (9) is lower. For instance, at x=9, the absolute error for (11) has kept constant, but for (9) the absolute error has decreased to 10-4 (see Figure 6(b)). For higher x values the accuracy improves for (11) and (9); for example, at x=30.5, the absolute error is 10-3 and 10-5, respectively (see Figure 6(c)). From this analysis, it can be concluded that (9) provides better accuracy than (11). Figure 7 depicts the absolute error of approximations for the sine integral (10) and (12). When x<3.5, absolute error for (10) is higher than that for (12), with a maximum error of 10-3 (see Figure 7(a)).Figure 7
Absolute error comparison between (10) and (12).
(a)
Absolute error comparison in the interval0≤x≤8 (b)
Absolute error comparison in the interval8≤x≤12 (c)
Absolute error comparison in the interval30≤x≤30.5Whenx>3.5, the accuracy improves notoriously for (10). For instance, in Figure 7(b), notice that at x=11 the absolute error is about 1×10-4, while (12) keeps its absolute error at about 10-3. The accuracy of (10) continues to improve for larger x values. At x=30.4, the absolute error approaches 10-6, while (12) keeps its error magnitude close to 1×10-3 (see Figure 7(c)). Nevertheless, although (9) and (10) exhibit an improved accuracy as x increases, they just handle one term that describes the behaviour of the Fresnel integrals while (11) and (12) handle two terms.Figure8 shows a comparison between approximations, (9) and (10), and exact Fresnel integrals (1) and (2). On the other hand, Figure 9 shows the Cornu graphic generated by the exact equations (9) and (10). It can be noticed that inside the spiral, as x value increases, accuracy improves notoriously.Figure 8
Comparison between exact equations and approximate equations (9) and (10).
(a)
Comparison between (1) and (9) (b)
Comparison between (2) and (10)Figure 9
Comparison between Cornu graph and exact equations.


## 5. Study Cases

This section provides two study cases applied to the optics field, specifically, the diffraction phenomenon. The first addresses the semi-infinite opaque screen subject; the second is about the double aperture problem for Fresnel diffraction.


### 5.1. The Semi-Infinite Opaque Screen for Fresnel Diffraction

For this study case we will consider the semi-infinite opaque screen defined by theΣ plane (see Figure 10), where IP is the intensity at point P in units Jm2/s or W/m2 and I0 is the intensity of these incident parallel wavefronts. From [1], intensity at P is proportional to(13)IP=Io2=12-Sx2+12-Cx2.Figure 10
Semi-infinite plane.For simplicity, we assume that the semi-infinite opaque screen for Fresnel diffraction is illuminated by a plane wave given by(14)vs=s2rλ,where s is the source, r the observation point, and λ the wavelength of the projected light beam. By substituting (14) in (9) and (10), equations are now in terms of s; thus, it is possible to calculate S(v(s)) and C(v(s)). Now, substituting in (13), we obtain(15)IP=Io2=12-Svs2+12-Cvs2.Proposed numerical values arer=4.0m and λ=632.8nm, 540nm, 390nm, 100nm. To obtain the quantitative values of the intensities from the Cornu spiral, the length from point (0.5,0.5) to points P1, P2, P3, P4, and P5 is measured, calculated using (9) and (10).Likewise, Figures11–14 depict the diagrams of the irradiance distribution and Cornu spiral for a semi-infinite opaque screen, obtained using (9) and (10), against the exact equations for different values of λ. The origin of the Cornu spiral diagram corresponds to the edge of the semi-infinite opaque screen, point (P3), where the intensity has decreased a quarter for negative values of v; in numerical terms, we have that C(0)=S(0)=0, and Ip=I0/4 because half of the waveform is obstructed. The amplitude of the perturbation has been reduced by half and the irradiance drops by a quarter. As already mentioned, this occurs at point P3 of the plane (see Figure 10), which can be located at the corresponding calculated irradiance distribution and the Cornu spiral in Figures 11–14.Figure 11
The semi-infinite opaque screen withr=4.0m and λ=632.8nm.
(a)
The corresponding calculated irradiance distribution withr=4.0m and λ=632.8nm (b)
The Cornu spiral for a semi-infinite opaque screen withr=4.0m and λ=632.8nmFigure 12
The semi-infinite opaque screen withr=4.0m and λ=540nm.
(a)
The corresponding calculated irradiance distribution withr=4.0m and λ=540nm (b)
The Cornu spiral for a semi-infinite screen withr=4.0m and λ=540nmFigure 13
The semi-infinite opaque screen withr=4.0m and λ=390nm.
(a)
The corresponding calculated irradiance distribution withr=4.0m and λ=390nm (b)
The Cornu spiral for a semi-infinite screen withr=4.0m and λ=390nmFigure 14
The semi-infinite opaque screen withr=4.0m and λ=100nm.
(a)
The corresponding calculated irradiance distribution withr=4.0m and λ=100nm (b)
The Cornu spiral for a semi-infinite screen withr=4.0m and λ=100nmObservation pointsP1 and P2 are on the shade of the Σ plane. The arc cords decrease monotonically within the region; in consequence, the irradiance abruptly declines. In fact, as λ decreases, the bands move closer to the edge and become thinner [1]. Observation points were calculated using (9) and (10) and compared with the ones obtained by exact equations. These points are depicted in the magnitude diagrams for λ=632.8nm, λ=540nm, λ=390nm, and λ=100nm (see Figures 11(a), 12(a), 13(a), and 14(a), respectively). The points obtained by the approximate method in these diagrams keep a relationship with the approximated Cornu spiral in Figures 11(b), 12(b), 13(b), and 14(b), respectively.


### 5.2. The Double Aperture Problem for Fresnel Diffraction

Light diffraction coming through a rectangular aperture is analysed in [1, 25]. This study case considers the double aperture for Fresnel diffraction presented and studied in [24] (see Figure 15). We assume that the light source with wavelength λ is S, at a distance P0, and the aperture is centred on the Cartesian system x,y; the origin of the coordinate system is located at the exact centre of the double aperture (see Figure 15). Upper and lower edges of the aperture are located at z=c and z=-c, respectively. For this system, the width for each aperture is (b-a) and separation between centres is (a+b). Besides, it is considered that diffracted light is observed on a screen located at qo away from it. Reference [24] can be consulted for further detailed analysis of this study case.Figure 15
Configuration of the double aperture problem for Fresnel diffraction.Simulations performed in Matlab for this study case (based on [24]) are presented in Figure 16(a). This work considers the same aperture proportions given in [24], these are 2a=2mm, 2b=4mm, 2c=6mm, q0=400mm, image area 10mm×10mm, and step size of 0.01mm. Additionally, the illumination wavelength is λ=632mm.Figure 16
Computer simulated Fresnel diffraction from double rectangular apertures.
(a)
Computer simulated Fresnel diffraction using approximate equations and iterative Fresnel integrals method [24] (b)
Computer simulated Fresnel diffraction using approximate equations (9) and (10)To generate the simulation shown in Figure16(b), built-in functions of Matlab were replaced by the proposed functions (9) and (10) (see the Appendix). The Pearson correlation coefficient, r, is wildly employed in statistical analysis, pattern recognition, and image processing [26, 27]. For the monochromatic digital images, the Person correlation coefficient is defined as(16)r=∑i=1xi-xmyi-ym∑i=1xi-xm2∑i=1yi-ym2,where xi and yi are intensity values of ith pixel in 1st and 2nd image, respectively; xm and ym are mean intensity values of first and second image, respectively. The correlation coefficient is r=1 when both images are completely identical. For this study case, correlation value is very high between the images generated with a value of 0.999964736.


## 5.1. The Semi-Infinite Opaque Screen for Fresnel Diffraction

For this study case we will consider the semi-infinite opaque screen defined by theΣ plane (see Figure 10), where IP is the intensity at point P in units Jm2/s or W/m2 and I0 is the intensity of these incident parallel wavefronts. From [1], intensity at P is proportional to(13)IP=Io2=12-Sx2+12-Cx2.Figure 10
Semi-infinite plane.For simplicity, we assume that the semi-infinite opaque screen for Fresnel diffraction is illuminated by a plane wave given by(14)vs=s2rλ,where s is the source, r the observation point, and λ the wavelength of the projected light beam. By substituting (14) in (9) and (10), equations are now in terms of s; thus, it is possible to calculate S(v(s)) and C(v(s)). Now, substituting in (13), we obtain(15)IP=Io2=12-Svs2+12-Cvs2.Proposed numerical values arer=4.0m and λ=632.8nm, 540nm, 390nm, 100nm. To obtain the quantitative values of the intensities from the Cornu spiral, the length from point (0.5,0.5) to points P1, P2, P3, P4, and P5 is measured, calculated using (9) and (10).Likewise, Figures11–14 depict the diagrams of the irradiance distribution and Cornu spiral for a semi-infinite opaque screen, obtained using (9) and (10), against the exact equations for different values of λ. The origin of the Cornu spiral diagram corresponds to the edge of the semi-infinite opaque screen, point (P3), where the intensity has decreased a quarter for negative values of v; in numerical terms, we have that C(0)=S(0)=0, and Ip=I0/4 because half of the waveform is obstructed. The amplitude of the perturbation has been reduced by half and the irradiance drops by a quarter. As already mentioned, this occurs at point P3 of the plane (see Figure 10), which can be located at the corresponding calculated irradiance distribution and the Cornu spiral in Figures 11–14.Figure 11
The semi-infinite opaque screen withr=4.0m and λ=632.8nm.
(a)
The corresponding calculated irradiance distribution withr=4.0m and λ=632.8nm (b)
The Cornu spiral for a semi-infinite opaque screen withr=4.0m and λ=632.8nmFigure 12
The semi-infinite opaque screen withr=4.0m and λ=540nm.
(a)
The corresponding calculated irradiance distribution withr=4.0m and λ=540nm (b)
The Cornu spiral for a semi-infinite screen withr=4.0m and λ=540nmFigure 13
The semi-infinite opaque screen withr=4.0m and λ=390nm.
(a)
The corresponding calculated irradiance distribution withr=4.0m and λ=390nm (b)
The Cornu spiral for a semi-infinite screen withr=4.0m and λ=390nmFigure 14
The semi-infinite opaque screen withr=4.0m and λ=100nm.
(a)
The corresponding calculated irradiance distribution withr=4.0m and λ=100nm (b)
The Cornu spiral for a semi-infinite screen withr=4.0m and λ=100nmObservation pointsP1 and P2 are on the shade of the Σ plane. The arc cords decrease monotonically within the region; in consequence, the irradiance abruptly declines. In fact, as λ decreases, the bands move closer to the edge and become thinner [1]. Observation points were calculated using (9) and (10) and compared with the ones obtained by exact equations. These points are depicted in the magnitude diagrams for λ=632.8nm, λ=540nm, λ=390nm, and λ=100nm (see Figures 11(a), 12(a), 13(a), and 14(a), respectively). The points obtained by the approximate method in these diagrams keep a relationship with the approximated Cornu spiral in Figures 11(b), 12(b), 13(b), and 14(b), respectively.


## 5.2. The Double Aperture Problem for Fresnel Diffraction

Light diffraction coming through a rectangular aperture is analysed in [1, 25]. This study case considers the double aperture for Fresnel diffraction presented and studied in [24] (see Figure 15). We assume that the light source with wavelength λ is S, at a distance P0, and the aperture is centred on the Cartesian system x,y; the origin of the coordinate system is located at the exact centre of the double aperture (see Figure 15). Upper and lower edges of the aperture are located at z=c and z=-c, respectively. For this system, the width for each aperture is (b-a) and separation between centres is (a+b). Besides, it is considered that diffracted light is observed on a screen located at qo away from it. Reference [24] can be consulted for further detailed analysis of this study case.Figure 15
Configuration of the double aperture problem for Fresnel diffraction.Simulations performed in Matlab for this study case (based on [24]) are presented in Figure 16(a). This work considers the same aperture proportions given in [24], these are 2a=2mm, 2b=4mm, 2c=6mm, q0=400mm, image area 10mm×10mm, and step size of 0.01mm. Additionally, the illumination wavelength is λ=632mm.Figure 16
Computer simulated Fresnel diffraction from double rectangular apertures.
(a)
Computer simulated Fresnel diffraction using approximate equations and iterative Fresnel integrals method [24] (b)
Computer simulated Fresnel diffraction using approximate equations (9) and (10)To generate the simulation shown in Figure16(b), built-in functions of Matlab were replaced by the proposed functions (9) and (10) (see the Appendix). The Pearson correlation coefficient, r, is wildly employed in statistical analysis, pattern recognition, and image processing [26, 27]. For the monochromatic digital images, the Person correlation coefficient is defined as(16)r=∑i=1xi-xmyi-ym∑i=1xi-xm2∑i=1yi-ym2,where xi and yi are intensity values of ith pixel in 1st and 2nd image, respectively; xm and ym are mean intensity values of first and second image, respectively. The correlation coefficient is r=1 when both images are completely identical. For this study case, correlation value is very high between the images generated with a value of 0.999964736.


## 6. Discussion

This section presents numerical results of the previous study cases. For the first case, the comparison is performed using the Cornus spiral and observation points of the irradiance distribution. We calculated the absolute error difference from(0.5,0.5) to every observation points of the diffraction pattern. Table 1 shows data for the case when r=4.0m and λ=632.8nm; the error is 0.0185625807 and 0.006802459 at observation points P1 and P2, respectively.Table 1
Observation points comparison forr=4.0m and λ=632.8nm.
Observation point Exact Approximate Absolute error 1 0.227105905 0.208543324 0.018562581  2 0.375123732 0.381926192 0.006802459  3 0.707106781 0.707106781 0.000000000  4 1.655561529 1.643873584 0.011690945  5 1.247629737 1.245639720 0.001990017Notice how the error decreases on observation pointsP1 and P2 (see Tables 1–4), because they are internally moving within the spiral since its geometry is contracting and are in function of the diffracted light wavelength. For observation points P4 and P5, the error has increased to about 8×10-7 because they have moved to the right, closer to the origin. This happens for two reasons: the wavelength is getting smaller and we approach the region where the error value is higher for (9) and (10). Nevertheless, in Table 4, the absolute error at P5 has diminished due to the accuracy of our proposed expressions (see Figures 4 and 5). For observation point P3, the error is zero becuase the intersection with y axis is always the same. This error is in good agreement with the exact spiral because (9) and (10) evaluated at x=0 are equal to zero.Table 2
Observation points comparison forr=4.0m and λ=540nm.
Observation point Exact Approximate Absolute error 1 0.212072842 0.199156973 0.012915869  2 0.358007985 0.364436149 0.006428163  3 0.707106781 0.707106781 0.000000000  4 1.655525630 1.643873930 0.011651701  5 1.247648595 1.245609700 0.002038895Table 3
Observation point comparison forr=4.0m and λ=390nm.
Observation point Exact Approximate Absolute error 1 0.183070676 0.187027639 0.003956963  2 0.322264740 0.326617422 0.004352681  3 0.707106781 0.707106781 0.000000000  4 1.655510405 1.643827879 0.011682526  5 1.247659193 1.245607630 0.002051563Table 4
Observation point comparison forr=4.0m and λ=100nm.
Observation point Exact Approximate Absolute error 1 0.094716894 0.094879587 0.000162693  2 0.185219606 0.188826768 0.003607161  3 0.707106781 0.707106781 0.000000000  4 1.655561120 1.643873413 0.011687707  5 1.247601422 1.245601596 0.001999826For the second study case, the Pearson correlation coefficient value was0.999964736, which indicates that properties of diffraction patterns from the reported simulation in [24] and the ones obtained from our approximations (see Figure 16) are identical. This result shows the accuracy that can be accomplished by using our proposal on applications for Fresnel diffraction.


## 7. Conclusions

Approximations proposed on this article allow evaluating Fresnel sine and Fresnel cosine integrals with high accuracy. In the first study case, we built the Cornu spiral and the irradiance distribution with good results since the lower wavelength of the diffraction, the higher accuracy on the shading zone, and the zero are kept at the origin of the Cornu diagram. The proposed equations in this work are continuous throughout the range of real numbers and are evaluated at the origin, with an exact value of zero. As for the second study case, simulations showed a correlation close to the unit. Therefore, the simulations of the double aperture problem for Fresnel diffraction produce near identical figures; thus, the simulation in Figure16(a) presented in [24] keeps an almost perfect correlation against the simulation in Figure 16(b). With the obtained results, the proposed equations can be employed in diverse applications, specially for values of x>1.5; as the value of x increases, accuracy improves considerably. Applying the experience of approximate Fresnel functions on this paper, as a future work, we will attempt to improve precision of the results. One possible path to follow would be the use of a different technique, besides least squares method, that allows improving precision closer to the origin region; among other techniques to be applied is the power series extender method [28, 29]. This method may be useful because it provides highly accurate approximations at the origin, while maintaining good prediction for large domains.



---

*Source: 4031793-2018-07-04.xml*