Numerical search of surface waves in a periodic grating

A two-dimensional reflection grating described by the Neumann problem for the Helmhotz equation is considered. The grating is obtained by adding a periodic sequence of rectangles to the upper half-plane. A method for approximate computation of the scattering matrix in the grating is implemented and its convergence is studied. A modification of the method is used for the numerical search of surface waves.