An efficient algorithm for solving high order Sturm-Liouville problems using variational iteration method (Q2871111)

An efficient algorithm based on the generalized variational iteration method is proposed in order to solve \(2m\)-order Sturm-Liouville problems. After the Lagrange multiplier is identified, a suitable iteration formula is established providing a sequence of successive approximations convergent to the exact solution. Sufficient condition for this convergence is obtained and the a priori error estimate is provided. The proposed method is illustrated and tested on three numerical experiments involving second, fourth, and sixth order Sturm-Liouville problems. The case of stiff equations is discussed in another numerical example. The obtained numerical results show that the variational iteration method is an efficient tool to compute the eigenvalues of high \(2m\)-order Sturm-Liouville problems.