new function

find the damping ratio...but you have to load data and specify freq
manuly (r=... load())