We present analytical formulas for the calculation of the two-center two-electron integrals in the basis of Slater geminals and products of Slater orbitals. Our derivation starts with establishing a inhomogeneous fourth-order ordinary differential equation that is obeyed by the master integral, the simplest integral with inverse powers of all interparticle distances. To solve this equation it was necessary to introduce a new family of special functions which are defined through their series expansions around regular singular points of the differential equation. To increase the power of the interparticle distances under the sign of the integral we developed a family of open-ended recursion relations. A handful of special cases of the integrals is also analysed with some remarks on simplifications that occur. Additionally, we present some numerical examples of the master integral that validate the usefulness and correctness of the key equations derived in this paper. In particular, we compare our results with the calculations based on the series expansion of the exp(-γ r12) term in the master integral. ; Comment: 28 pages, 0 figures, 7 tables