A symbolic-numerical algorithm for the computation of matrix elements in the parametric eigenvalue problem

A symbolic-numerical algorithm for the computation of the matrix elements in the parametric eigenvalue problem to a prescribed accuracy is presented. A procedure for calculating the oblate angular spheroidal functions that depend on a parameter is discussed. This procedure also yields the corresponding eigenvalues and the matrix elements (integrals of the eigenfunctions multiplied by their derivatives with respect to the parameter). The efficiency of the algorithm is confirmed by the computation of the eigenvalues, eigenfunctions, and the matrix elements and by the comparison with the known data and the asymptotic expansions for small and large values of the parameter. The algorithm is implemented as a package of programs in Maple-Fortran and is used for the reduction of a singular two-dimensional boundary value problem for the elliptic second-order partial differential equation to a regular boundary value problem for a system of second-order ordinary differential equations using the Kantorovich method.