Computation of Sturm-Liouville eigenvalues via Lagrangian interpolation