Spectrum completion and inverse Sturm–Liouville problems

Abstract

Given a finite set of eigenvalues of a regular Sturm–Liouville problem for the equation , the potential of which is unknown. We show the possibility to compute more eigenvalues without any additional information on the potential . Moreover, considering the Sturm–Liouville problem with the boundary conditions and , where are some constants, we complete its spectrum without additional information neither on the potential nor on the constants and . The eigenvalues are computed with a uniform absolute accuracy. Based on this result, we propose a new method for numerical solution of the inverse Sturm–Liouville problem of recovering the potential from two spectra. The method includes the completion of the spectra in the first step and reduction to a system of linear algebraic equations in the second. The potential is recovered from the first component of the solution vector. The approach is based on special Neumann series of Bessel functions representations for solutions of Sturm–Liouville equations possessing remarkable properties and leads to an efficient numerical algorithm for solving inverse Sturm–Liouville problems.