Boundary integral method