clear,close all,
num=input('Type in the numerator coefficients =');
den=input('Type in the denominator coefficients =');
[z,p,k]=tf2zp(num,den);
m=abs(p);
disp('Zeros are at');disp(z);
disp('Poles are at');disp(p);
disp('Gain constant');disp(k);
disp('Raduis of poles');disp(m);
sos=zp2sos(z,p,k);
disp('Second-order sections');disp(real(sos));
zplane(num,den)