On a method for solving the inverse Sturm-Liouville problem (Q2316692)

Denote by $\varphi (\lambda, x)$ a solution of the Cauchy problem \[ \begin{aligned} -\varphi^{\prime \prime}+q(x) \varphi &=\lambda \varphi \\ \varphi(\lambda, 0) &=1, \quad \varphi^{\prime}(\lambda, 0)=h, \end{aligned} \] which can also be rewritten in the form \[ \varphi(\lambda, x)=\cos \sqrt{\lambda}x+\int_{0}^{x}G(x, t) \cos \sqrt{\lambda} t d t. \] It is known that $G$ admits the representation: \[ G(x, t)=\sum_{n=0}^{\infty} \frac{a_{n}(x)}{x} P_{2n}\left(\frac{t}{x}\right), \quad 0 \leq t \leq x, \] where $P_{m}$ stands for the Legendre polynomial of order $m$. In this paper, the authors prove that the coefficients $a_{n}(x)$ satisfy an infinite system of linear algebraic equations. They also give a method for recovering the potential based on the above linear algebraic equations.