mcp500_1;
makeQ;
[phi, X, y] = pdmaxcut( Q);
