We propose an accurate and efficient method to compute vibrational spectra of molecules, based on exact diagonalization of an algebraically calculated matrix based on powers of Morse coordinate. The present work focuses on the 1D potential of diatomic molecules: as typical examples, we apply this method to the standard LennardJones oscillator, and to the ab-initio potential of the H2 molecule. ...