A preliminary study of the capability of the finite difference and finite element methods (FDM, FEM) to evaluate eigenvalues of one-, two-, and three-dimensional self-adjoint operators is reported with reference to applications dealing with the description of vibrational levels. Results of harmonic oscillator model potentials and ab initio PES for the water molecule are obtained by using the FDM. In spite of the large matrices used, low accuracy, nonvariational results are found. A different method, based on FEM and normal coordinates, is therefore proposed. Two nearly harmonic cases are studied and it is shown that variational results with higher accuracy can be obtained with a moderate cost. The vibrational levels of the water molecule are also calculated in order to compare the results with those of the FDM treatment.