Most seismic imaging algorithms assume single-scattering, because this linearizes the inverse problem. Multiply-scattered energy can be significant, however, especially if there are strong reflectors in the model. We setup a series expansion that allows for the inclusion of multipy-scattered waves through a sequence of linear problems rather than a single nonlinear problem. This results in an iterative algorithm for imaging with multiply-scattered waves that allows the imaging of structures not illuminated by singly scattered waves by constructing a sequence of images of the subsurface. A downside of this approach is the dependence of subsequent images on previous ones; extensions to mitigate this problem will also be discussed.