The boundary integral in Boundary Element Method affects the precision and the speed of the method. If the boundary integral with constant element, the nonsingular integrals are popularly calculated by the Gauss numerical integral, and the singular integrals are popularly calculated by the analytical integral. This paper presents an alternative way with Gauss formula to transform the double integral in elastic problem on 3-d into the linear integrals on the boundary of each subdomains, so that all the singular integrals and nonsingular integrals are calculated by analytical method. The example indicates that this method makes the precision and the speed of BEM improve.