A Fortran program is given for calculating wave functions of the molecular hydrogen ion, expanded in terms of single-center Laguerre-type orbitals. Using this program, the radial quantum number has been extended to 203, and accuracy has been attained in energy of the order of 10−6 a.u. The electron-nucleus Coulomb integrals are evaluated numerically by applying Gaussian quadratures.