We propose a method of computing certain double integrals that are encountered in dynamic problems of the theory of elasticity for semi-infinite bodies with cracks. Application of this method makes it possible to reduce a system of boundary integral equations for the interaction of a crack with the boundary of the half-space to a form that contains only integrals over the regions occupied by the cracks.In reducing dynamic problems of the theory of elasticity for a half-space with cracks to a system of boundary integral equations it is necessary to overcome significant mathematical difficulties connected with the fact that these equations contain integrals over the unbounded region So that coincides with the boundary of the half-space. The occurrence of such an integral results from satisfying the boundary conditions on the boundary of the half-space. This greatly complicates the application of numerical methods for solving such equations. For that reason lowering the dimension of the system of integral equations by eliminating the integral equations that contain an infinite limit of integration is important.Consider an unbounded body weakened by N + l planar cracks, arbitrarily situated, whose opposite surfaces S~, n = 0, N, are subject to the action of self-balancing forces, varying harmonically in time with frequency w.Choose local coordinate systems Onxln, x2m x3n, n = O, N, in such a way that the domain Sn occupied by the nth crack lies in the xl,~x2n-coordinate plane, and the values x3n = :t:0 correspond to the opposite surfaces S~ of the nth crack. Then the location of the cracks in the body is completely determined by giving the direction cosines of the coordinate axes and the distances between the cracks. We denote by ljkn, mjk,~, and njkm k,n = 0, N, the direction cosines of the Xjn-aXiS in the Okxlkx2kxak-coordinate system. We denote by dkn the distance between the origins of the kth and nth coordinate systems, and by ejkn the direction cosines of the vector dkn in the kth coordinate system. Then an arbitrary point x* with coordinates (xl~, x2n, x3~) in the nth coordinate system has the following coordinates (xlkn, x2kn, x3kn) in the kth coordinate system: 3 3 3 j=l j=l j=1(1)To reduce the system of boundary integral equations of a dynamic problem for a half-space with cracks to a form not containing integrals over an unbounded domain So, it is necessary to compute the following double integrals:where x~n and 4 are the points with coordinates (XlOn, X2On, ::g3On) and (