The homotopy perturbation method [4], proposed first by He in
1998 and was further developed and improved by He [3,5,7,6]. The method yields a very rapid convergence of the solution
series in the most cases. Usually, one iteration leads to high accuracy of the solution. Although goal of Hes homotopy
perturbation method was to find a technique to unify linear and nonlinear