By introducing irreducible tensorial operators, the recently developed method of solving multi-center molecular integrals by coupling the two-center exchange integrals with a Fourier transform scheme is further generalized to handle any Slater type orbitals of arbitrarily higher angular momentum states. Numerical procedures can be applied to the final expressions to obtain numerical results.