“…According to this method, the amplitude and phase of the first harmonic of v are considered slowly varying quantities, that is, v ( t ) t = w ( t ) + h ( t ), where w ( t ) = W ( t )cos[ ωt + θ ( t )] denotes the first harmonic of v and h ( t ) the sum of the remaining harmonics. Thus, by using the averaging method, equation is reduced to a system of first‐order equations for the amplitude and the phase of w ( t ), which take the form (Appendix) where the dot denotes the time derivative, and I c and I s denote the cosine and sine coefficients of the fundamental harmonic of i ( w + h , v in ), respectively. The quantity h is given by the solution of the linear equation where i h ( w , v in ) denotes the harmonics of i ( v , v in ) that we assume to be slightly dependent on the harmonics of v , so that i h ( w + h , v in ) ≈ i h ( w , v in ).…”