Time-dependent wavepacket techniques are used together with the Moller operator formulation of scattering theory to derive an exact expression that relates scattering matrix elements to a correlation function between incoming reactant and outgoing product wavepackets. Using this expression, a computational strategy is developed for the numerical evaluation of scattering matrix elements that is focused on the independent propagation of reactant and product wavepackets. Computational efficiency is achieved through the use of different coordinates and interaction Hamiltonians that are well suited for the reactant and product channels. A simple one-dimensional scattering problem is used to illustrate the essential features of the method.