e5p3.m

% Example 5.3 
Kp = [1 10];  w = logspace(-1,1);
for k = 1:2
  Gnum = [1 0.1];  Gden = [1 1+0.25*Kp(k)];
    Gjomega = freqs(Gnum,Gden,w); Gmag = abs(Gjomega);
      plot(w,Gmag)
      title('Sensitivity Functions')
      xlabel('Omega'), ylabel('Sensitivity')
      grid
      hold on
end
hold off