SUMMARYA multiharmonic method for analysis of non-linear dynamic systems submitted to periodic loading conditions is presented. The approach is based on a systematic use of the fast Fourier transform. The exact linearization of the resulting equations in the frequency domain allows to obtain full quadratic convergence rate.