This paper describes the development of a two-dimensional direct boundary integral method for elastotstatics boundary value problems in multilayered materials with cracks. Boundary integral equations for displacements are supplemented with traction boundary integral equations derived by differentiating the dispalcements. The discontinuities over the cracks. The kernel functions for the boundary integral equations are developed form the analytical solution to the wich the interface continuity conditions are fulfilled exactly. The displacement discontinuity kernels are derived by differentiating the point force kernels).