The approximation method in a Hilbert space is applied to the analysis of a generalized non-singular Sturm-Liouville eigenvalue problem. The variable coefficients in the system equation and the solution are approximated by a suitable linear independent basis, the expansion coefficients are explicitly expressed as a backward recursive formula and the eigenvalues are obtained by solving an algebraical characteristic equation. The computation is simple and straightforward.