A linearly scaling grid method for generalized Langevin and other integrodifferential equations

A grid-based method for numerical solution of systems of integrodifferential dynamical equations is introduced that has computational costs which scale linearly with the grid size. Grid parameters are automatically generated from a user supplied error tolerance. We illustrate the accuracy of the method by solving a variety of representative sets of equations. The method is much faster than comparable quadratically scaling algorithms.