CODE EXAMPLE FOR MATLAB

matlab tf get poles

% Variables
zeta=0.5;
omegan=1;

% Transfer function
sys=tf([omegan^2],conv([1],[1 2*zeta*omegan omegan^2]))

% Get poles
pole(sys)