% hw 09_p2

clear all; close all

s = tf( 's' )

K = 1;

sysG = (1)/((s-1));
sysH = (s+2)/((s+1)^2+1);

%sysRY = ( K*sysG ) / ( 1 + K*sysG*sysH );
%sysRY0 = (s^2+2*s+2) / ( s^3 + s^2 + s );

figure(11)

rlocus( sysG*sysH );
axis( [ -5 5, -5, 5 ])
axis('equal')
set(findall(gcf,'type','line'),'linewidth',3);

figure(12)

bode( sysG*sysH );
grid
%axis('equal')
set(findall(gcf,'type','line'),'linewidth',3);
grid


figure(13)

nyquist( sysG*sysH );
axis('equal')
set(findall(gcf,'type','line'),'linewidth',3);
axis( [ -1.5 0.5, -0.5, 0.5 ])


