global A; global B;
A1=[1 -.5;1 0]; B1=[1;1];
A2=[-1 -.5;1 0]; B2=[-2;1];
setlmis([]);
M1=lmivar(2,[1 2]) ;
M2=lmivar(2,[1 2]) ;
X=lmivar(1,[2 1]) ;
lmiterm([1 1 1 X],1,A1','s');
lmiterm([1 1 1 M1],B1,-1,'s');
lmiterm([2 1 1 X],1,A2','s');
lmiterm([2 1 1 M2],B2,-1,'s');
lmiterm([3 1 1 X],1,A1','s');
lmiterm([3 1 1 X],1,A2','s');
lmiterm([3 1 1 M1],B2,-1,'s');
lmiterm([3 1 1 M2],B1,-1,'s');
lmiterm([4 1 1 X],1,-1);
lmis=getlmis;
[tmin,xfeas]=feasp(lmis)
m1=dec2mat(lmis,xfeas,M1)
m2=dec2mat(lmis,xfeas,M2)
x=dec2mat(lmis,xfeas,X)
P=inv(x)
F1=m1*P
F2=m2*P