G=abv(sbs(nd2sys(1,conv([1 0 100],[0.01 1])),nd2sys(1,[1 0 100])),...
      sbs(nd2sys([1 0.1],conv([1 0 100],[1 1])),nd2sys(1,[1 0 100])));
w=logspace(-2,2,100);
[A,B,C,D]=unpck(G);
[a,b,c,d]=minreal(A,B,C,D);
g=pck(a,b,c,d);
p=spoles(g);
z=tzero(a,b,c,d);
gf=frsp(g,w);
rga=vrga(gf);
[u,s,v]=vsvd(gf);
figure(1);
vplot('liv,lm',vdiag(rga),'-',sel(rga,1,2),'--',sel(rga,2,1),'--');
legend('-','L11L22','--','L12L21');
figure(2);
vplot('liv,lm',s);