x=sdpvar(2,1); A=[0 x(1) 0;x(1) x(2) 0;0 0 1+x(1)]; F = set(A>0); solvesdp(F,x(1),sdpsettings('solver','penbmi','savedebug',1));