%close all,clear;
figure(1)
t=0:0.01:6;
num=[1];%input('Type in the coefficient vector of the numerator polynomial of H(s):num=');
den=[1 5 8 4];%input('Type in the coefficient vector of the denominator polynomial of H(s):den=');
z=roots(num);
p=roots(den);
z=z';p=p';

h=impulse(num,den,t);
disp('The zeros of the system are:')
z
disp('The poles of the systme are:')
p
[r,p,k]=residue(num,den);
disp('The expression of H(s) in partial fraction expansion form')
r=r'
p=p'
k
subplot(211)
pzmap(num,den)
subplot(212)
plot(t,h)
title('The impulse response h(t) of the system')
axis([0,5,min(h),max(h)])
xlabel('Time in second')
grid on