Many types of molecular integrals involving Slater functions can be expressed, with the -function method in terms of sets of one-dimensional auxiliary integrals whose integrands contain two-range functions. After reviewing the properties of these functions (including recurrence relations, derivatives, integral representations, and series expansions), we carry out a detailed study of the auxiliary integrals aimed to facilitate both the formal and computational applications of the -function method. The usefulness of this study in formal applications is illustrated with an example. The high performance in numerical applications is proved by the development of a very efficient program for the calculation of two-center integrals with Slater functions corresponding to electrostatic potential, electric field, and electric field gradient.