function nf = NatFreq(F,Sys) % Natural Frequencies, Damping Ratios, and Real Roots % =================================================== nf = 1; Eigenvalues = esort(eig(F)) [NaturalFrequencies,DampingRatios] = damp(Sys)