“…This behavior has been analyzed in the past by means of harmonic balance methods [6,7], linearized harmonic balance procedures [8], asymptotic perturbation techniques which combine the harmonic balance procedure and the method of multiple time scales [9], the Linstedt-Poincaré procedure [10], the method of averaging of Krylov-Bogoliubov-Mitropolskii [10][11][12], parameter-perturbation Linstedt-Poincaré techniques [13,14] which employ an artificial or book-keeping parameter and expand both the solution and some constants that appear (or are introduced) in the differential equation in terms of this parameter, artificial parameter-Linstedt-Poincaré techniques [15,16] based on the introduction of a linear term proportional to the unknown frequency of oscillation and a new independent variable and the use of either the third-order equation or a system of a first-order and a second-order ordinary differential equations, etc. Parameter-perturbation methods are extensions of the homotopy perturbation technique which introduces an artificial parameter and expands both the solution and the unknown frequency of oscillation in terms of this parameter and has been applied to a variety of problems [17].…”