% Variables zeta=0.5; omegan=1; % Transfer function sys=tf([omegan^2],conv([1],[1 2*zeta*omegan omegan^2])) % Get poles pole(sys)