This paper shows the application of GFDM to the problem of seismic wave propagation in 2-D. To use this method in unbounded domains one must truncate the computational grid-point avoiding reflection from the edges. PML absorbing boundary condition has then been included in the numerical model proposed in this work.