Seismic waves in anisotropic media are more complex than in isotropic media. Here we derive the propagating matrices for seismic waves in 2-D transversely isotropic medium (TIM). With eigen-decomposition, eigenvalues and eigenvectors are given in analytical forms, therefore, calculation of propagators are simple and accurate. For a 2-D model of layered media, we compute the seismic responses to an impulse in the f-k domain, and then do a 2-D inverse Fourier transformation. Clear qP and qSV waves can be recognized from the resultant sections.