The math is sound. The method is not. We recast the azimuthal mode-coupled surface equations as a block-Toeplitz convolution in the mode index and diagonalize this coupling via a discrete Fourier transform in \(n\), then briefly review the basic definition and core properties of the DFT used in this step.