“…The HAM [19,20] was first proposed by Liao in 1992 to solve lots of nonlinear problems. This method has been successfully applied to many nonlinear problems, such as physical models with an infinite number of singularities [21], nonlinear eigenvalue problems [22], fractional Sturm-Liouville problems [23], optimal control problems [24,25], Cahn-Hilliard initial value problem [26], semi-linear elliptic boundary value problems [27] and so on [28]. The HAM contains a certain auxiliary parameter which provides us with a simple way to adjust and control the convergence region and rate of convergence of the series solution.…”