wn = 2;

dr = .5; % damping ratio



wh = wn*sqrt(1-2*dr^2+sqrt(2-4*dr^2+4*dr^4));