This is my codes for the GN benchmark.

function [A]=B128(Zout)

A=zeros(128,128);

Label=[ones(32,1);2*ones(32,1);3*ones(32,1);4*ones(32,1)];

Node=1:128;

C=[Node',Label];

pin=(16-Zout)/31;

pout=Zout/96;

R=nan*ones(8500,2);

t=1;

for i=1:4

SN=Node(C(:,2)==i);

DN=Node(C(:,2)~=i);

for j=1:32

SNp=SN(SN~=SN(j));

S=[DN,SNp];

for m=1:127

T=[SN(j),S(m)];

a=Chong(T,R);

if a==1

test=rand(1);

if m