A method for studying the dynamical properties of liquids by molecular dynamics simulation is described. Its basis is the measurement of the response to a weak applied field of appropriate character. The explicit form of the mechanical perturbation is worked out in several cases, and details are given of the numerical techniques used in implementing the method.