Abstract: As an excellent numerical method, boundary element method (BEM) has been widely applied in various scientific and engineering problems. In this paper, a new boundary integral method is obtained based on Somigliana's equation and the properties of Green's function by referring to the idea of boundary element method. It can be used to find the analytic solution of linear elastic problems. The boundary integral method can also be obtained from Betti's reciprocity theorem. By using this new method, the classical problem of elastic circular inclusion under a uniform tensile field at infinity is solved. Firstly, the perfect bonding between inclusion and matrix is set up, and the displacement and stress at interface are expanded according to Fourier series. According to the symmetry of the problem and the orthogonality of trigonometric function, the hypothesis is simplified and the number of undetermined coefficients is reduced. Secondly, the appropriate trial functions are selected (these trial functions satisfy the condition of displacement single value and the control equation of linear elasticity without body force). And the boundary integral method is used to calculate the displacement and stress at the interface. Then the displacement and stress in the domain are solved using similar tricks. The exact analytical solution of the problem is obtained, which is exactly the same with the results in literatures. When the elastic modulus of the inclusion is zero or tends to infinity, it degenerates to the analytical solution of the problem of circular hole or rigid inclusion. The solution process shows that if the problem has boundary conditions at infinity, trial functions should meet the boundary condition at infinity. If the domain of the problem contains the coordinate origin, the displacement and stress of trial functions at the origin should be limited. The results show that the method is effective.