A harmonic-balance based method for computing Floquet's multipliers in Lur'e systems