ts = 0.5;
al = [1+.3*s, .5*s; .1*s 1+0.7*s];
bl = [.2 .4;.6 .8];
bl = pol(bl);
[a,b,c,d] = lmf2ss(bl,al);
G = ss(a,b,c,d);
