“…The HAM can guarantee the convergence of the series solutions by auxiliary parameters especially the so-called convergence-controller parameter ℏ [5][6]. In recent years, this method has been successfully employed to solve many types of nonlinear problems in science and engineering such as Troesch's problem [7], the Fitzhugh-Nagumo equation [8], heat radiation equations [9], MHD viscoelastic fluid flow [10], the Coupled nonlinear schrödinger equations [11], the Continuous population models for single and interacting species [12], the Boussinesq problem [13,14], the Sturm-Liouville problems [15], the wave propagation problems [16], Laplace equation with Dirichlet and Neumann boundary conditions [17], differential-difference equation [18], fractional equations [19]. Many authors trying to solve equation (1) by homotopy analysis method [7,20] but for λ=3 only.…”