A new numerical method is developed for solution of the Gelfand - Levitan - Marchenko inverse scattering integral equations. The method is based on the fast inversion procedure of a Toeplitz Hermitian matrix and special bordering technique. The method is highly competitive with the known discrete layer peeling method in speed and exceeds it noticeably in accuracy at high reflectance.
Download