A = [0,1,0,0,0,0,0,0,1,0,0,0,0,0,0;
0,0,1,0,1,0,1,0,0,0,0,0,0,0,0;
0,1,0,0,0,1,0,1,0,0,0,0,0,0,0;
0,0,1,0,0,0,0,0,0,0,0,1,0,0,0;
1,0,0,0,0,0,0,0,0,1,0,0,0,0,0;
0,0,0,0,0,0,0,0,0,1,1,0,0,0,0;
0,0,0,0,0,0,0,0,0,1,1,0,0,0,0;
0,0,0,1,0,0,0,0,0,0,1,0,0,0,0;
0,0,0,0,1,1,0,0,0,1,0,0,0,0,0;
0,0,0,0,0,0,0,0,0,0,0,0,1,0,0;
0,0,0,0,0,0,0,0,0,0,0,0,0,0,1;
0,0,0,0,0,0,1,1,0,0,1,0,0,0,0;
0,0,0,0,0,0,0,0,1,0,0,0,0,1,0;
0,0,0,0,0,0,0,0,0,1,1,0,1,0,1;
0,0,0,0,0,0,0,0,0,0,0,1,0,1,0];

q = 0.15;

[G,p] = Gmatrix(A,q);
p