The multi-pole PML (MPML) is tested on models that simulate seismic waves traveling through the subsurface. Using a recursive integration technique a stretching function consisting of the sum of multiple stretching functions is implemented in the velocity-stress finite difference time domain wave equations. The MPML is implemented in both the rotated staggered grid (RSG) and the Virieux grid. The performance of the MPML is tested on a square model, rectangular model and a rectangular model with a free-surface and compared to other types of PML’s implemented in these models. The main result is that the MPML can be implemented in the velocity stress wave equations giving stable results similar to other PML types.