%%%% lyapunov2.m

clear; format compact
% --- ƒXƒeƒbƒv 1 ---------
A = [ 0  1
      2 -3 ];
% --- ƒXƒeƒbƒv 2 ---------
P = sdpvar(2,2,'sy');
% --- ƒXƒeƒbƒv 3 ---------
LMI = set([]);
M = P*A + A'*P;
LMI = LMI + set(M < 0);
LMI = LMI + set(P > eye(2));
% --- ƒXƒeƒbƒv 4 ---------
sol = solvesdp(LMI)
% ------------------------
if sol.problem ~= 1 
  % --- ƒXƒeƒbƒv 5 -------
  P_feas = double(P)  
  sign(eig(double(P))')
  sign(eig(double(-M))')
end
