A = [-1 2;-3 -4];
X = sdpvar(2,2);
I = eye(2,2);

C = [X >= I, A'*X+X*A <= -I];

solvesdp(C);

X_feasible_2 = double(X)