A new numerical approach for the efficient computation of Floquet multipliers within the harmonic balance technique