Many numerical methods were developed for this type of nonlinear partial
differential equations such as the Adomian Decomposition Method (ADM)
[3], the Homotopy Perturbation Method (HPM) [4], the Homotopy Analysis
Method (HAM) [5] and Variational Iteration Method (VIM) [6].