Abstract
We describe a numerical method that can be used to calculate the propagation of light in a medium of constant (possibly complex) index of refraction n. The method integrates the Rayleigh–Sommerfeld diffraction integral numerically. After an appropriate change of integration variables, the integrand of the diffraction integral is split into a slowly varying and an (often fast) oscillating quadratic factor. The slowly varying factor is approximated by a spline fit, and the resulting Fresnel integrals are subsequently integrated exactly. Although the method is not as fast as methods involving a fast Fourier transform, such as plane-wave propagation or Fresnel approximation, it is accurate over a greater range than these methods.
© 2005 Optical Society of America
Full Article | PDF ArticleMore Like This
John C. Heurtley
J. Opt. Soc. Am. 63(8) 1003-1008 (1973)
Soheil Mehrabkhani and Thomas Schneider
Opt. Express 25(24) 30229-30240 (2017)
Arvind S. Marathay and John F. McCalmont
J. Opt. Soc. Am. A 21(4) 510-516 (2004)