
% Optimization Toolbox.
%% The Paley graph
% The graph we use is the Paley graph with n=17. The vertices of this graph
% correspond to the integers mod 17 with an edge between two vertices
% exactly when their difference is a quadratic residue mod 17. The |paley|
% function creates such a graph. 
g = graph;
paley(g,17);
%% Maximum clique
% We use |omega| to find the clique number and a maximum clique.
[w,S] = omega(g);
disp(['The clique number is ', int2str(w)])
notS = setdiff(1:17,S);
p = partition({S,not