% Part (a)

b1 = [ 4 ];
a1 = [ 1 2 17 ];
poles1 = roots(a1);

pzplot(b1,a1,-0.5);
title('H_1(s)');
print H1_poles.eps -deps
