c=[1;8;4;1;2;7;6;5;3;5;4;8;4;1;6;3];
a=[];
b=[];
aeq=[1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0;
    0,0,0,0,1,1,1,1,0,0,0,0,0,0,0,0;
    0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,0;
    0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1;
    1,0,0,0,1,0,0,0,1,0,0,0,1,0,0,0;
    0,1,0,0,0,1,0,0,0,1,0,0,0,1,0,0;
    0,0,1,0,0,0,1,0,0,0,1,0,0,0,1,0;
    0,0,0,1,0,0,0,1,0,0,0,1,0,0,0,1];
beq=ones(8,1);
vlb=zeros(16,1);
vub=ones(16,1);
[x,y]=linprog(c,a,b,aeq,beq,vlb,vub)
x =
    0.0000
    0.0000
    0.0000
    1.0000
    1.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    1.0000
    0.0000
    0.0000
    1.0000
    0.0000
    0.0000
y =
8.0000
c=[16;15;19;22;17;21;19;18;24;22;18;17;17;19;22;16];
a=[];
b=[];
aeq=[1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0;
    0,0,0,0,1,1,1,1,0,0,0,0,0,0,0,0;
    0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,0;
    0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1;
    1,0,0,0,1,0,0,0,1,0,0,0,1,0,0,0;
    0,1,0,0,0,1,0,0,0,1,0,0,0,1,0,0;
    0,0,1,0,0,0,1,0,0,0,1,0,0,0,1,
                        0,0/0,1/16/17/22/19/aeq/数学/18/1,0/  
                        
                          0,0/0,1/16/17/22/19/aeq/数学/18/1,0/  
                        
                        
                        
                        
                    -->