%%%% lyapunov_eps.m

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