den1 = [1 5]; den1 = conv([1 0 0.5],[1 -1.5]); den=conv(den1,[1 1 1]); sys=tf(num,den) rlocus(sys)