figure(2)
[sv1,w1]=sigma(sys1);

[sv2,w2]=sigma(sys2);

[sv3,w3]=sigma(sys3);


%  plot the maximum singular values


plot(w1,log(sv1(1,:)),'r');
hold on
plot(w2,log(sv2(1,:)),'g');
plot(w3,log(sv3(1,:)),'b');
hold off
norm(sys1,2)
norm(sys1,inf)

norm(sys2,2)
norm(sys2,inf)

norm(sys3,2)
norm(sys3,inf)
return
