Application of the incremental harmonic balance method to cubic non-linearity systems