In the computation of a three–dimensional steady creeping flow around a rigid body, the total body force and torque are well predicted using a boundary integral equation (BIE) with a single concentrated pair Stokeslet- Rotlet located at an interior point of the body. However, the distribution of surface tractions are seldom considered. Then, a completed indirect velocity BIE of Fredholm type and second-kind is employed for the computation of the pointwise tractions, and it is numerically solved by using either collocation or Galerkin weighting procedures over flat triangles. In the Galerkin case, a full numerical quadrature is proposed in order to handle the weak singularity of the tensor kernels, which is an extension for fluid engineering of a general framework (Taylor, 2003, “Accurate and Efficient Numerical Integration of Weakly Singulars Integrals in Galerkin EFIE Solutions,” IEEE Trans. on Antennas and Propag., **51**(7), pp. 1630–1637). Several numerical simulations of steady creeping flow around closed bodies are presented, where results compare well with semianalytical and finite-element solutions, showing the ability of the method for obtaining the viscous drag and capturing the singular behavior of the surface tractions close to edges and corners. Also, deliberately intricate geometries are considered.