Nonparaxial propagation of a cylindrical beam with Lanczos reduction

We discuss a numerical method based on Lanczos reduction of modeling nonparaxial propagation of a cylindrical symmetric beam. To illustrate the performance and demonstrate the significant difference between nonparaxial and paraxial beams, we consider Gaussian beam propagation in two different settings.