Sturm-Liouville problems with parameter dependent potential and boundary conditions (Q2468133)

The author considers the Sturm-Liouville eigenvalue problem \[ -y''+ q(x, \mu) y = \mu^2 y,\quad x \in [x_0, x_1], \] with the boundary constraint \[ A(y(x_0), y'(x_0), y(x_1), y'(x_1))^T = 0, \] where \(q(x, \mu) = q_0 (x) + q_1 (x) + q_2 (x) \mu^2\), \(q_0, q_1, q_2\) are complex valued integrable functions over \((x_0, x_1)\) and \(A = (a_{ij} (\mu))\) is a \(2 \times 4\) matrix with rank 2. Denote by \(y_c (x, \mu), y_s (x, \mu)\) the (fundamental) solutions with the boundary conditions \(y_c (x_0) = 1, y_c'(x_0) = 0\); \(y_s (x_0) = 0, y_s'(x_0) = 1\), respectively. Then each eigenvalue \(\lambda = \mu^2\) can be obtained as a root of \(B(\mu)\), where \(B(\mu) = \) det \((A \omega_1 | A \omega_2)\) with \(\omega_1 = (1, 0, y_c (x_1, \mu), y_c' (x_1, \mu))^T\); \(\omega_2 = (0, 1, y_s (x_1, \mu), y_s' (x_1, \mu))^T\). To obtain \(y_c (x_1, \mu), y_s (x_1, \mu)\), a modified sampling theorem is used [\textit{C. E. Shin, S.-Y. Chung} and \textit{D. Kim}, J. Math. Anal. Appl. 291, 50--65 (2004; Zbl 1044.94518)] since the classical Whittaker-Shannon-Kotel'nokov theorem can not be applied here due to the regularity insufficiency. In sum, just as the numerical examples in this article shows, the method applied consists of the following steps: {\parindent=6mm \begin{itemize}\item[(1)] Use integral representations to get a control of the (exponential) growth of \(y_s (x_1, \mu), y_s ' (x_1, \mu)\) and \(y_c (x_1, \mu), y_c ' (x_1, \mu)\) (as functions of \(\mu\)), and accordingly to get the sampling points. \item[(2)] Use a numerical method (in these examples the Fehlberg 4--5 order Runge-Kutta method is used) to get the sampling values of the functions considered in (1). \item[(3)] Recover the functions considered in (1) (a truncated version). \item[(4)] Compute the eigenvalues from \(B(\mu)\). \end{itemize}}