A recently proposed flutter prediction method called mu-omega method in frequency domain is generalized to the static divergence analysis. Numerical study is conducted using this method to analyze the stability boundary of complex aeroelastic system. First, the flutter characteristic of a high aspect ratio swept wing is calculated. The critical flutter speed and flutter frequency obtained by the mu-omega method have a relative error of 0.83% and 0.39% compared with those solved by the p-k method, respectively. After then, static divergence analysis is carried out for a delta wing model. The critical divergence speed is obtained by the mu-omega method with a relative error of 2.02% compared with that given by the p-k method. It is indicated that the mu-omega method has good numerical accuracy and is applicable in the stability analysis of complex aeroelastic systems.