This paper presents a boundary element formulation for 3-D linear and viscoelastic bodies subjected to the body force of gravity. The Laplace transformation is first used to suppress the time variable, and solutions of displacements and stresses are found in the transformed domain. The time domain solutions are then found by an accurate and efficient numerical inversion method which requires only real calculations for all quantities. Input and output data, and solutions in the transformed and time domains are connected through an Interactive Data Language code written by the authors. While particular solutions of stresses and displacements related to the body force of gravity (which is applied at time t 0 and is kept constant) are derived, the Green's functions in the Laplace domain are obtained through the correspondence principle. The new formulation has been implemented into an existing 3-D BEM program, and several numerical examples involving 3-D viscoelastic bodies are presented. Although the discussion in this paper focuses on Maxwell viscoelastic and isotropic media, other linear isotropic and even anisotropic viscoelastic models can also be incorporated, without difficulty, into the 3-D viscoelastic BEM program.