clear;clc;

format shortG

options = optimset('Display','off');

H = [1.32706 -0.26629 -0.10485 2.37964
    -2.14291 1.11803 0.89604 0.92137
    -0.27878 -0.86637 -0.43404 0.18012
    -0.35438 0.09057 0.22529 -0.09379
    0.42988 -1.50226 2.11735 -0.66289
    0.34362 0.63772 -0.49262 -0.80495
    -0.00402 -0.09919 -1.29207 0.04081
    -0.60345 -0.65027 -0.68906 -0.44298
    -0.04899 -0.38018 -0.83804 -0.83195
    1.33196 1.91824 0.61199 -0.68537]

q = zeros(10,4);
w = zeros(10,1);

for j=1:10
   
    c=[-H(j,4),0,0,0];
   
    a=[H(1,4),-H(1,1:3);
        H(2,4),-H(2,1:3);
        H(3,4),-H(3,1:3);
        H(4,4),-H(4,1:3);
        H(5,4),-H(5,1:3);
        H(6,4),-H(6,1:3);
        H(7,4),-H(7,1:3);
        H(8,4),-H(8,1:3);
        H(9,4),-H(9,1:3);
        H(10,4),-H(10,1:3);
        H(j,4),0,0,0];
   
    b=[0;0;0;0;0;0;0;0;0;0;1];
    aeq=[0,H(j,1:3)];
    beq=[1];
    lb=[0,0,0,0];
    ub=[1,1,1,1];
   
    [xx,fval,exitflag,output] = linprog(c,a,b,aeq,beq,lb,ub,[],options);

    q(j,:) = xx';
    w(j,:) = exitflag';
   
end

q
w