den=[1,10,8,17,16,5];p=roots(den),
ii=find(real(p)>=0);n=length(ii);                   %ii：不穩(wěn)定根在p向量中的下標(biāo)；
if(n>0),disp('系統(tǒng)不穩(wěn)定，不穩(wěn)定的根的個(gè)數(shù)');n   %n：不穩(wěn)定根在p向量中的下標(biāo)的個(gè)數(shù)
else,disp('系統(tǒng)穩(wěn)定');end