This paper describes the progress that has been made during the decade 1973–1983 in the study of the method of harmonic balance, known to many engineers as the describing function method. It is now possible to give rigorous existence proofs and optimized error bounds for limit cycles that the describing function method has predicted in non-linear systems: even in “highly” non-linear systems such as those with discontinuities and multiple-valued relations.