In this article, systems of linear and nonlinear partial differential equations solve
with variational iteration method (VIM). The VIM is to construct correction
functional using general Lagrange multipliers identified optimally via the
variational theory and the initial approximations can be freely chosen with
unknown constants. In this article, we used (VIM) to solve partial differential
equation system. This method is very simple and reduces computation and readily
and fast converges to the answer. The procedure is well organized.