We have studied scattered pulse shapes by modeling inhomogeneities as a sequence of infinitesimally thin homogeneous layers. With oblique incidence of plane P or SV waves, the reflected-converted-transmitted waves are obtained by taking the calculus limit for the sum of primary interactions of the incident wave with all layer boundaries. The resulting scattered waves thus present themselves naturally in the time domain. For an incident impulse, the scattered pulse shape is merely an analytic function of the depth from which scatter has taken place within the inhomogeneity. The direct application of this simple method appears to be new, and we have found it remarkably accurate when compared with methods in which higher-order boundary interactions are also retained (i.e., Haskell methods and an adaptation in the time domain which also keeps all multiples).In specific studies of P-waves incident (up to 30 degrees away from the vertical) upon a 5 km thick crust-mantle transition, between materials having impedance ratio 1:2.8, we find the scattered pulse shapes are given adequately by our theory, for the passband of short-period seismometers. Indeed, the theory remains remarkably accurate even for long periods, being in error by only 8 per cent at zero frequency.