We consider Sturm–Liouville differential operator with potential having a finite number of simple discontinuities. This paper is devoted to the numerical solution of such inverse spectral problems. The main result of this work is a procedure that is able to recover both the points of discontinuities as well as the heights of the jumps. Following, using these results, we may apply a suitable numerical method (for example, the generalized Rundell–Sacks algorithm with a special form of the reference potential) to reconstruct the potential more precisely.