s=tf('s')

G1=2/(s*10+1);
G2=0.5/(s*5+1);

G=series(G1,G2)

damp(G);
[r, p, k] = residue(G.num{1}, G.den{1})

opt = stepDataOptions('StepAmplitude', 0.1)
step(G, opt)

K=3;
W0=4;
psy=0.25;
I=(K*W0^2)/(s^2+2*psy*W0*s+W0^2);
step(I)
damp(I)
%Ex 2


%H=4/(s*(s+2*0.5*2));
%zero(H)
%pole(H) %0 ; -2
%P1=s+2;
%P2=s;
%[r, p, k] = residue(H.num{1}, H.den{1}); % -2/s+2 ; 2/s
 %damp(H);

%HU=H*1/s; %s2*(s+2)
%     4
%  -----------
%  s^3 + 2 s^2
%[r, p, k] = residue(HU.num{1}, HU.den{1}); % 2/s^2 + -1/s+ 1/(s+2)
%t=0:0.1:5;                                                                         
%h=2*t-1+exp(-2*t);

