clear all; close all; s = tf('s'); for c=0.8:0.05:1.2 P = 1/((s - 1)*(c*s + 1)); figure(1); bode(P, 'b'); hold on; figure(2); pzmap(P); hold on; end P0 = 1/((s - 1)*(s + 1)); figure(1); bode(P0, 'r');