In this paper, we apply the Modified Variational Iteration Method (MVIM) for solving systems of partial differential equations. The proposed modification is made by introducing He’s polynomials in the correction functional of Variational Iteration Method (VIM). Several examples are given to verify the reliability and efficiency of the method. The fact that the MVIM solves nonlinear problems wit...