This paper is devoted to applying a numerical method for solution of the
2D sine-Gordon equation. The bivariate multiple quadratic
quasi-interpolation ( MQQI ) method is adopted to simulate this
equations, which the first order spatial derivative is approximated by
MQQI, the second spatial and time derivative are approximated by forward
difference. One of the merit of this scheme is its simple structure and
easy implementation. In the meanwhile, we present truncation and total
error of this scheme, high accuracy and efficiency of the method are
verified by numerical experiments. In addition, the optimal value of
parameters are investigated in this article based on Luh [10, 11].