In this paper, the sufficient conditions that guarantee the convergence of the variational iteration method when applied to solve a coupled system of non-linear partial differential equations are presented. Special 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.