Calculation of elastic wave fields has important applications in a variety of engineering fields including NDE (Non-destructive evaluation). Scattering problems have been investigated by numerous authors with different solution schemes. For simple geometries of the scatterers (e.g., cylinders or spheres), the analysis of steady-state elastic wave scattering has been carried out using analytical techniques. For arbitrary geometries and multiple inclusions, numerical methods have been developed. Special finite element methods, e.g., the infinite element method and a hybrid method called the Global-Local finite element method have also been developed for this purpose. Recently, the boundary integral equation method has been used successfully to solve scattering problems. In this paper, a volume integral equation method (VIEM) is proposed as a new numerical solution scheme for the solution of general elasto-dynamic problems in unbounded solids containing multiple inclusions and voids or cracks. A boundary integral equation method (BIEM) is also presented for elastic wave scattering problems. The relative advantage of the volume and boundary integral equation methods for solving scattering problems is discussed.