This paper presents the derivation of the non-singular integral representation of stresses in two- and three-dimensional elastostatics. In contrast to the strongly singular and weakly singular integral representations, the numerical computation of the nearly singular integrals is eliminated because all the integrands are made finite in this new formulation even if the internal point approaches the boundary. Thus the method gives accurate numerical results even in that portion of a solid which is very close to a discretized boundary. Three test problems are analysed in which we present a comparison of the accuracies achieved by the numerical computations based on the use of strongly singular, weakly singular and non-singular integral representations of stresses.