Fully analytical method for calculation of the molecular integrals over Slater−type orbitals with non−integer principal quantum numbers are proposed. These integrals are expressed through relativistic molecular auxiliary functions derived in our previous paper [Phys. Rev. E 91, 023303 (2015)]. The procedure for computation of the molecular auxiliary functions is detailed. It applies both in relativistic and non-relativistic electronic structure theory. It is capable of yielding highly accurate molecular integrals for all ranges of orbital parameters and quantum numbers.