function ps53(n)
if nargin<1, n=5; end
q=poly(-rand(n,1));
p=((-1).^(0:n)).*q;
h=hsvd(tf(p,q));
close(gcf);bar(h);grid