clear
F=-[100 400 200 200 100 0 200 800 0 0 100 100 100 100 600 267 153 99 451 30 100 33 33 34 800];
A=[]; B=[];
Aeq=[ones(1,5) zeros(1,20); zeros(1,5) ones(1,5) zeros(1,15); zeros(1,10) ones(1,5) zeros(1,10);
zeros(1,15) ones(1,5) zeros(1,5); zeros(1,20) ones(1,5);
zeros(1,0) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4);
zeros(1,1) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,3);
zeros(1,2) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,2);
zeros(1,3) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,1);
zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,0)];
Beq = ones(1,10);
[x1,min_fval1]=bintprog(F,A,B,Aeq,Beq);
clear
F=-[100 400 200 200 100 0 200 800 0 0 100 100 100 100 600 267 153 99 451 30 zeros(1,5)];
A=[zeros(1,0) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4);
zeros(1,1) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,3);
zeros(1,2) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,2);
zeros(1,3) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,1);
zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,0)];
B=ones(1,5);
Aeq=[ones(1,5) zeros(1,20); zeros(1,5) ones(1,5) zeros(1,15);
zeros(1,10) ones(1,5) zeros(1,10); zeros(1,15) ones(1,5) zeros(1,5); zeros(1,20) ones(1,5)];
Beq=[ones(1,4) 0];
[x2,min_fval2]=bintprog(F,A,B,Aeq,Beq);
clear
F=-[100 400 200 200 100 0 200 800 0 0 100 100 100 100 600 267 153 99 451 30 zeros(1,5)];
A=[zeros(1,5) ones(1,5) zeros(1,15); zeros(1,15) ones(1,5) zeros(1,5)];
B=2*ones(1,2);
Aeq=[ones(1,5) zeros(1,20); zeros(1,10) ones(1,5) zeros(1,10); zeros(1,20) ones(1,5);
zeros(1,0) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4);
zeros(1,1) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,3);
zeros(1,2) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,2);
zeros(1,3) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,1);
zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,0)];
Beq=[ones(1,2) 0 ones(1,5)];
[x3,min_fval3]=bintprog(F,A,B,Aeq,Beq);
clear
F=-[100 700 200 0 0 0 200 800 0 0 100 100 100 100 600 871 0 99 0 30 0 33 33 34 900];
A=[]; B=[];
Aeq = [ones(1,3) zeros(1,22); zeros(1,5) ones(1,5) zeros(1,15);
zeros(1,10) ones(1,5) zeros(1,10); zeros(1,15) 1 0 1 0 1 zeros(1,5);
zeros(1,21) ones(1,4); zeros(1,0) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,9);
zeros(1,1) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,9) 1 zeros(1,3);
zeros(1,2) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,2);
zeros(1,8) 1 zeros(1,4) 1 zeros(1,9) 1 zeros(1,1);
zeros(1,9) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,0)];
Beq=ones(1,10);
[x4,min_fval4]=bintprog(F,A,B,Aeq,Beq);
clear
F=-[100 700 200 0 0 0 200 800 0 0 100 100 100 100 600 871 0 99 0 30...
0 33 33 34 900 250 250 0 250 250 111 1 0 333 555];
A=[]; B=[];
Aeq=[ones(1,3) zeros(1,32); zeros(1,5) ones(1,5) zeros(1,25);
zeros(1,10) 1 0 ones(1,3) zeros(1,20); zeros(1,15) 1 0 1 0 1 zeros(1,15);
zeros(1,21) ones(1,4) zeros(1,10); zeros(1,25) ones(1,2) 0 ones(1,2) zeros(1,5);
zeros(1,30) ones(1,2) 0 ones(1,2);
zeros(1,0) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,9) 1 zeros(1,4) 1 zeros(1,4);
zeros(1,1) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,9) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,3);
zeros(1,2) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,12);
zeros(1,8) 1 zeros(1,4) 1 zeros(1,9) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,1);
zeros(1,9) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,0)];
Beq=[ones(1,10) 2*ones(1,2)];
[x5,min_fval5]=bintprog(F,A,B,Aeq,Beq);
clear
F=-[100 400 200 200 100 0 200 800 0 0 100 100 100 100 600 267 153 99 451 30 100 33 33 34 800];
A=[ones(1,5) zeros(1,20); zeros(1,5) ones(1,5) zeros(1,15);
zeros(1,10) ones(1,5) zeros(1,10); zeros(1,15) ones(1,5) zeros(1,5);
zeros(1,20) ones(1,5);...
zeros(1,0) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4);
zeros(1,1) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,3);
zeros(1,2) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,2);
zeros(1,3) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,1);
zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,0)];
B=ones(1,10) ;
Aeq=ones(1,25);
Beq=3;
[x6,min_fval6]=bintprog(F,A,B,Aeq,Beq);