clear all;
close all;

s = tf('s');

for c=0.8:0.05:1.2
  Ptilde = 1/((s - 1)*(c*s + 1));

  figure(1);
  bode(Ptilde, 'b');
  hold on;

  figure(2);
  pzmap(Ptilde);
  hold on;
end

P0 = 1/((s - 1)*(s + 1));
%P0 = 1/((s - 1)*(1.2*s + 1));
figure(1);
bode(P0, 'r');
