Algebraic theory is used to fit the set of bound vibrational energy levels of H-2 and its isotopomers, up to the dissociation limit and to an accuracy of within 1 cm(-1) in each case. Coordinate space potential functions are constructed from these algebraic Hamiltonians, and the limitations of the algebraic approach are assessed.