In the next link: Computation of Hankel Transform using FFT (Fourier)
Rainer implemented a great solution given in the next reference: Manuel Guizar-Sicairos and Julio C. Gutiérrez-Vega, "Computation of quasi-discrete Hankel transforms of integer order for propagating optical wave fields," J. Opt. Soc. Am. A 21, 53-58 (2004).
However he does not comment on how to calculate the propagation of an optical field (for example, how to reproduce Fig.3 from the reference he cited). So it would be nice if he or someone else could provide that example.
My main interest is to calculate the Fresnel Diffraction pattern of a circular aperture, however, I do not want to use a rectangular grid, but to solve the problem in polar coordinates (that is why I quote the other question and the author that posted the question)