Lagrangian Homotopy Analysis Method using the Least Action Principle

The Homotopy Analysis Method (HAM) is a powerful technique which allows to derive approximate solutions of both ordinary and partial differential equations. We propose to use a variational approach based on the Least Action Principle (LAP) in order to improve the efficiency of the HAM when applied to Lagrangian systems. The extremization of the action is achieved by varying the HAM parameter, therefore controlling the accuracy of the approximation. As case studies we consider the harmonic oscillator, the cubic and the quartic anharmonic oscillators, and the Korteweg-de Vries partial differential equation. We compare our results with those obtained using the standard approach, which is based on the residual error square method. We see that our method accelerates the convergence of the HAM parameter to the exact value in the cases in which the exact solution is known. When the exact solution is not analytically known, we find that our method performs better than the standard HAM for the cases we have analyzed. Moreover, our method shows better performance when the order of the approximation is increased and when the nonlinearity of the equations is stronger.
Comment: 26 pages, 7 figures, 15 tables and 3 appendices