% Plant (2.26), G(s)=3(-2s+1)/(5s+1)(10s+1)
g=nd2sys([-2,1],[50 15 1],3);

% frequency response of G
points=301;
omega=logspace(-3,1,points);
g_f=frsp(g,omega);
