This work is part of a continuing effort to develop a capability for quantifying the scattering of ultrasonic waves by arbitrarily shaped flaws. The general problem of elastodynamic behavior of a homogeneous, isotropic defect in an otherwise homogeneous, isotropic fullspace is cast as a Boundary Integral Equation (BIE). A general scattering model is needed to provide information for probability of detection (POD) models and inversion schemes for cases when low or high frequency approximations are not appropriate. Previously the Boundary Element Method (BEM), a method for solving the BIE, was adapted to NDE and the void problem was investigated [1]. Here we focus on the inclusion problem, experimental verification, and to overall extensions of the capability.