In this article, a new technique is proposed to simulated two-dimensional Zakharov–Kuznetsov equation with the initial condition. The idea of this technique is based on Taylors' series in its derivation. Two test problems are presented to illustrate the performance of the new scheme. Analytical approximate solutions that we obtain are compared with variational iteration method (VIM) and homotopy analysis method (HAM). The results show that the new scheme is efficient and better than the other methods in accuracy and convergence