A basis of Hermite splines is used in conjunction with the collocation method to solve the orbital equations for diatomic molecules. Accurate solutions of the Hartree-Fock equations are obtained using iterative methods over most regions of space, while solving the equations by Gaussian elimination near the nuclear centres. In order to improve the speed and accuracy of our iterative scheme, a new self-adjoint form of the Hartree-Fock equation is derived. Using this new equation, our iterative subroutines solve the Hartree-Fock equations to one part in 106. The Gaussian elimination routines are accurate to better than one part in 108.