Diatomic molecule symmetry eigenfunctions by direct diagonalization