In this study a composite numerical method is implemented for solving the nonlinear coupled system of sine Gordon equations. The first and second order of accuracy unconditionally stable difference schemes are used in fixed point iteration for the approximate solution of the nonlinear coupled system. The unique solvability of the system is considered. Some numerical results are presented in tables and error analysis is given. Presented numerical results are carried out using the MATLAB software.