Gelfand-Levitan-Krein method in one-dimensional elasticity inverse problem

In this article we propose the numerical solution of the one dimensional inverse coefficient problem for seismic equation. We use a dynamical version of Gelfand-Levitan-Krein approach for reducing a nonlinear inverse problem for recovering the shear wave’s velocity and the density of the medium to two sequences of the linear integral equations. We propose numerical algorithm for solving these equations based on a fast inversion of a Toeplitz matrix. The proposed numerical methods base on the structure of the problem and therefore improve the efficiency of the algorithms, compared with standard approaches. We present numerical results for solving considered integral equations.