Abstract
In the present paper numerical solution of the coupled system of sine-Gordon equations are studied in the weak sense. A first order and two second order of accuracy unconditionally stable difference schemes corresponding to the system of sine-Gordon equations are considered. Solutions of these difference schemes are obtained in the space of distributions by variational methods. The fixed point theory and finite difference method is used together in the implementations of numerical experiments which verify theoretical statements.