We obtain numerical results to help the detection and characterization of subsurface cracks in solids by the application of P and Rayleigh elastic waves. The response is obtained from boundary integral equations, which belongs to the field of elastodynamics. Once the implementation of the boundary conditions has been done, a system of Fredholm integral equations of the second kind and order zero is found. This system is solved using the method of Gaussian elimination. Resonance peaks in the frequency domain allow us to infer the presence of cracks using spectral ratios. Several models of cracked media were analyzed, where effects due to different crack orientations and locations were observed. The results obtained are in good agreement with those published in the references.