%%%% lyapunov2.m

clear; format compact
% --- ステップ 1 ---------
A = [ 0  1
      2 -3 ];
% --- ステップ 2 ---------
P = sdpvar(2,2,'sy');
% --- ステップ 3 ---------
LMI = [];                   % ---> set に代わる新しい記述形式 []
M = P*A + A'*P;
LMI = [LMI, M < 0];         % ---> set に代わる新しい記述形式 []
LMI = [LMI, P > eye(2)];    % ---> set に代わる新しい記述形式 []
% --- ステップ 4 ---------
sol = solvesdp(LMI)
% ------------------------
if sol.problem ~= 1 
  % --- ステップ 5 -------
  P_feas = double(P)  
  sign(eig(double(P))')
  sign(eig(double(-M))')
end
