A formulation extending the Haskell-Thompson matrix method to include the effects of anelastic attenuation is presented. The formulation is exact in that no low-loss approximations are made. Consideration is given to nonparallel propagation and attenuation directions with corresponding velocity anisotropy. Examples are presented for models representing soils, the crust, and the core-mantle boundary.