General expressions are obtained for the calculation of the matrix elements of the Coulomb repulsion between electrons in two-electron diatomic molecules using the adiabatic basis. Based on these formulae, the algorithm for calculating these matrix elements for a hydrogen molecule in a wide range of internuclear distances is realised. The authors' results are compared with the asymptotic calculations and the previous calculations (1979) in the vicinity of the equilibrium position.