Quasi-periodic motions and their stability are addressed from the point of view of different harmonic balance-based approaches. Two numerical methods are used: a generalized multidimensional version of harmonic balance and a modification of a classical solution by harmonic balance. The application to the case of a nonlinear response of a Duffing oscillator under a bi-periodic excitation has allowed a comparison of computational costs and stability evaluation results. The solutions issued from both methods are close to one another and time marching tests showing a good agreement with the harmonic balance results confirm these nonlinear responses. Besides the overall adequacy verification, the observation comparisons would underline the fact that while the 2D approach features better performance in resolution cost, the stability computation turns out to be of more interest to be conducted by the modified 1D approach.