Numerical computation of eigenvalues in spectral gaps of Sturm–Liouville operators