solvesdp([constraints, x0 == [3;1]],objective);
double(u{1})