Recovery of a potential on a quantum star graph from Weyl's matrix (Q6087987)

The paper is concerned with the Sturm-Liouville problem \[ -u''(x) + q(x) u(x) = \lambda u(x) \] on the star-shaped graph \(\Omega\) with the standard matching conditions at the internal vertex: \[ u \in C(\Omega), \quad \sum_{j = 1}^M \partial u_j(v) = 0. \] The authors develop a numerical method for the reconstruction of the potential \(q(x)\) on the graph from the values of the Weyl matrix, which is also called the Dirichlet-to-Neumann map, at a finite number of points: \(\mathbf{M}(\lambda_1), \mathbf{M}(\lambda_2), \dots, \mathbf{M}(\lambda_m)\). The approach is based on the Neumann series of Bessel functions for the solutions of the Sturm-Liouville equations on the graph edges. First, by using the given values of the Weyl matrix, the authors construct a system of linear algebraic equations and, by solving it, find the Neumann series coefficients for separate edges. Then, using the obtained coefficients, eigenvalues of two spectra are found for each edge. Finally, the potential \(q_i(x)\) on each edge is recovered from two spectra. The authors provide numerical examples, which illustrate the efficiency of their method.