format short
num = [5 3.5 21.58 -4.6 18];
den = [5 15.5 31.7 22.52 4.8];
[z,p,k]= tf2zp (num,den)
sos = zp2sos(z,p,k)
zplane(z,p)