In this article, He's variational iteration method is applied to linear Sturm-Liouville eigenvalue and boundary value problems, including the harmonic oscillator. In this method, solutions of the problems are approximated by a set of functions that may include possible constants to be determined from the boundary conditions. By computing variations, the Lagrange multipliers are derived and the generalised expressions of variational iterations are constructed. Numerical results show that the method is simple, however powerful and effective. (c) 2009 Elsevier Ltd. All rights reserved.