A multi-scale perturbation method for studying nonlinear oscillations and their stability in higher order systems is developed. The technique is applied to the two-temperature reactor model Stability boundaries involving the various parameters are established and easily interpretable analytical expressions for nonlinear oscillations are presented. The method requires neither phase plane systems nor autonomous equations and thus provides an analytical tool for investigating higher order and space-dependent reactor models.