L=[2 -1 0 0 -1
-1 3 -1 -1 0
0 -1 2 -1 0
0 -1 -1 3 -1
-1 0 0 -1 2];
[phi, X, y] = maxcut( L);
disp(['upper bound on optimal value is  ',num2str(sum(y)/4)]);
[vals,inds]=sort(4*X(:,1));
disp(['unsorted and sorted values of column 1, and index positions  ']);
             [4*X(:,1) vals inds]
