Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)We present here a Boundary Element Method (BEM) implementation for SH harmonic waves in a class of inhomogeneous anisotropic media The inhomogeneity is assumed to be the same not only for the stiffnesses, but also for the density The implementation is based on a closed form fundamental solution for SH waves derived by Daros [C.H. Daros, A fundamental solution for SH-waves in class of inhomogeneous anisotropic media, Int J Eng Science 46 (2008) 809-817] We show numerical results obtained by the traditional boundary integral equation Moreover, the non-hypersingular traction based BEM is also implemented, allowing the modelling of cracks in inhomogeneous anisotropic media We obtain numerical results for the stress intensity factors (SIF) which are compared to previous published results (C) 2010 WILEY-VCH Verlag GmbH & Co KGaA, Weinheim902113121Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)