setlmis([]);
epsilon=lmivar(1,[1 1]);
P=lmivar(1,[n-m 1]);
lmiterm([1 1 1 P],A',1,'s'); % LMI #1: A'*P+P*A
lmiterm([1 1 1 P],.5*alpha,1,'s'); % LMI #1: alpha*P (NON SYMMETRIC?)
lmiterm([1 1 1 epsilon],.5*1,eye(n-m),'s'); % LMI #1: epsilon*eye(n-m) (NON SYMMETRIC?)
lmiterm([1 2 1 P],1,1); % LMI #1: P
lmiterm([1 2 2 epsilon],.5*1,-eye(n-m),'s'); % LMI #1: -epsilon*eye(n-m) (NON SYMMETRIC?)
lmiterm([1 3 1 0],zeros(m,n-m)); % LMI #1: zeros(m,n-m)
lmiterm([1 3 2 0],zeros(m,n-m)); % LMI #1: zeros(m,n-m)
lmiterm([1 3 3 epsilon],.5*1,-eye(m),'s'); % LMI #1: -epsilon*eye(m) (NON SYMMETRIC?)
lmiterm([-2 1 1 P],1,1); % LMI #2: P
yingzhu=getlmis;
[tmin,xfeas]=feasp(yingzhu)
P=dec2mat(yingzhu,xfeas,P)
epsilon=dec2mat(yingzhu,xfeas,epsilon)