Numerical solution of vector Sturm–Liouville problems with Dirichlet conditions and nonlinear dependence on the spectral parameter

A numerical-analytical iterative method is proposed for solving generalized self-adjoint regular vector Sturm-Liouville problems with Dirichlet boundary conditions. The method is based on eigenvalue (spectral) correction. The matrix coefficients of the equations are assumed to be nonlinear functions of the spectral parameter. For a relatively close initial approximation, the method is shown to have second-order convergence with respect to a small parameter. Test examples are considered, and the model problem of transverse vibrations of a hinged rod with a variable cross section is solved taking into account its rotational inertia.