Numerical simulations of systems of PDEs by variational iteration method

In this Letter, a general framework of the variational iteration method (VIM) is presented for solving systems of linear and nonlinear partial differential equations (PDEs). In VIM, a correction functional is constructed by a general Lagrange's multiplier which can be identified via a variational theory. VIM yields an approximate solution in the form of a rapid convergent series. Comparison with the exact solutions shows that VIM is a powerful method for the solution of linear and nonlinear systems of PDEs.