clear all;
close all;

s = tf('s');
load freqresp.mat;

figure(1);
nyquist(P0_g);

K = ss(0.15); %
clp = feedback(P0, K);
eig(clp)

figure(2);
bodemag(P0, 'b', clp, 'r');
legend('openloop', 'closedloop');
xlim([0.5 50]);