Abstract
In this paper, the sufficient conditions that guarantee the convergence of the variational iteration method when applied to solve a coupled system of nonlinear partial differential equations are presented. Especial attention is given to the error bound of the nth term of the resultant sequence. Numerical examples to show the efficiency of the method are presented.