  A new numerical method is developed for solution of the Gel'fand - 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.
