It is shown that the classical perturbation procedure for treating nonlinear systems leads to solutions expressed as Fourier-like series with slowly varying coefficients. These slowly varying coefficients contain the information about the long term behavior of the system. Inconsistently, the classical perturbation procedure expresses these coefficients as power series, a mode of expression which has notoriously poor long term validity. An operational procedure is presented for treating oscillations having slowly variable amplitudes and frequencies. An extension of the usual impedance concepts is presented for expressing the frequency characteristics of both linear and nonlinear elements when oscillations with many frequencies are present simultaneously and when these oscillations vary in both frequency and amplitude. From these methods, a perturbation procedure is devised which permits the behavior of systems to be computed with any order of accuracy, using only the algebraic processes which are characteristic of operational procedures. This procedure avoids expressing its results in terms of the local time. Instead, it expresses them in terms of the fundamental characteristics of the oscillations which axe present. As a consequence, the final solutions have the much desired long term validity and they may be used to obtain asymptotic estimates of the behavior of the system. The method is able to treat systems containing nonlinear perturbing elements and elements which we have described as moderately nonlinear. By means of examples it is shown that it is a straightforward process to treat systems to second order accuracy. This level of accuracy covers a large number of the intercoupling effects that characterize the more sophisticated nonlinear phenomena
