% Uncertain transfer function of the first actuator
%
hold off
omega = logspace(1,4,100);
K0 = 1.08; T0 = 0.005; 
Gnom = nd2sys([K0],[T0 1]);
Gnom_g = frsp(Gnom,omega);
for K = 0.9*K0:0.04*K0:1.1*K0
    for T = 0.8*T0:0.05*T0:1.2*T0
        G = nd2sys([K],[T 1]);
        G_g = frsp(G,omega);
        vplot('liv,lm',Gnom_g,'r-',G_g,'c--')
        hold on
    end
end
grid
xlabel('Frequency (rad/s)')
ylabel('Magnitude')
title('Uncertain transfer function of the first actuator')
%legend('|G_{nom}(j\omega)|','|G(j\omega)|',3)
hold off