This work is a generalization of the algorithm proposed by Kuzmak for
approximating solutions of slowly non-autonomous oscillatory systems.
It specifies a method for routine calculation of the integrals which d
etermine the slow evolution of amplitude and phase of the underlying p
eriodic oscillation, and interprets the algorithm from the viewpoint o
f classical, qualitative systems theory.