num = [5 3.5 21.58 -4.6 18];
den = [5 15.5 31.7 22.52 4.8];

factorize(num);
factorize(den);
k = num(1)/den(1);

zplane(num ,den)