Abstract
This article studies the numerical solution of the two-dimensional sine-Gordon equation (SGE) using a split-step Chebyshev Spectral Method. In our method we split the 2D SGE by considering one dimension at a time, first along x and then along y. In each fractional step we solve a 1D SGE. Time integration is handled by a finite difference scheme. The numerical solution is then compared with many of the known numerical solutions found throughout the literature. Our method is simple to implement and second order accurate in time and has spectral convergence. Our method is both fast and accurate.