Numerical solution of inverse spectral problems for Sturm-Liouville operators with discontinuous potentials

We consider Sturm-Liouville differential operators on a finite interval with discontinuous potentials having one jump. As the main result we obtain a procedure of recovering the location of the discontinuity and the height of the jump. Using our result, we apply a generalized Rundell-Sacks algorithm of Rafler and Böckmann for a more effective reconstruction of the potential and present some numerical examples.