We present a fast algorithm to calculate Coulomb/exchange integrals of prolate spheroidal electronic orbitals, which are the exact solutions of the single-electron, two-center Schrödinger equation for diatomic molecules. Our approach employs Neumann’s expansion of the Coulomb repulsion 1/jx yj, solves the resulting integrals symbolically in closed form and subsequently performs a numeric Taylor...