“…The stable regions correspond to finite solutions of the Mathieu–Hill equations, while the unstable regions correspond to solutions that increase indefinitely in time [ 3 ]. The stability charts of the Mathieu–Hill equations published in the literature are usually calculated using various methods, i.e., harmonic balance method [ 4 , 6 , 7 , 8 , 9 , 10 , 11 ], Galerkin method [ 11 , 12 ], multiscale methods [ 13 ], averaging method [ 14 , 15 ], discrete singular convolution [ 16 ], homotopy perturbation method [ 17 ], and small parameter method (Poincaré method) [ 3 , 11 ]. There is a wealth of literature available on parametric vibrations that comprehensively addresses all fundamental aspects of the subject [ 3 , 4 , 6 , 7 , 8 , 18 , 19 ].…”