The Sine-Gordon (SG) equations are very important in that they can accurately model many essential physical phenomena. In this paper, the Jacobi-Gauss-Lobatto collocation (JGL-C) and Generalized Lagrange Jacobi-Gauss-Lobatto collocation (GLJGL-C) methods are adopted and compared to simulate the (2 + 1)-dimensional nonlinear SG equations. In order to discretize the time variable t, the Crank-Nicolson method is employed. For the space variables, two numerical methods based on the aforementioned collocation methods are applied. Furthermore, error estimation for both methods is provided. The present numerical method is truly effective, free of integration and derivative, and easy to implement. The given examples and the results assert that the GLJGL-C method outperforms the JGL-C method in terms of computation speed. Also, the presented methods are very valid, effective, and reliable.