“…Recently, several methods have been introduced and developed to obtain approximate solutions for (NDEs) due to their complexity and the difficulty of solving them through traditional perturbation techniques. For example, variational iteration method, 5 homotopy perturbation method, 6 max-min approach, [7][8][9] global residue harmonic balance method (GRHBM) for obtaining higher-order approximate solutions, [10][11][12] modified homotopy perturbation method, [13][14][15] energy balance method, 16,17 Hamiltonian approach, [18][19][20] iteration perturbation technique, 21 coupled homotopy-variational approach, [22][23][24] frequency-amplitude formulation, 25,26 multiple scales technique, 27 parameter expansion method, 28 averaging method, 29 iteration method, 30 and Laplace variational iteration method. 31 The harmonic balance method (HBM) is one of the main techniques for obtaining approximate analytical solutions to NDEs describing oscillatory systems.…”