[num, den] = numden(H);
z = sym2poly(expand(num));
p = sym2poly(expand(den));
zplane(z, p);