The C-L method was generalized from Liapunov-Schmidt reduction method, comb
ined with theory of singularities, for study of non-autonomous dynamical sy
stems to obtain the typical bifurcating response curves in the system param
eter spaces. This method has been used, ar an example, to analyze the engin
eering nonlinear dynamical problems by obtaining the bifurcation programs a
nd response curves which are useful in developing techniques of control to
subharmonic instability of large rotating machinery.