C = [52 59 43.5 47.8 50 74 80 68 180 200;
        18 180 10 8 4 140 35 47 7 4;
        3 5.9 11 6 14.95 9.35 4.53 18 0.817 0.5;
        0.2 8.25 0.05 0.01 0 9 0.342 0 0 0;
        0.1 22 0.001 0.001 0 20 26.38 4 0 0;
        -5.2 -2.5 0 1.09 -2 29 50.7 5.24 100 100;
        1.806 2.5 5.084 3.142 5.055 2.6 1.7 9.5 1.5 2;
        2.65 2.5 2.71 2.851 2.98 2.67 2.51 2.87 1.03 0.9]
   d = [58.491;96.779;2.004;3.324;7.738;7.934;3.771;2.688]
A=ones(1,10);
A(1,10)=0;A(1,9)=0;
b=1
   Aeq = [1 1 1 1 1 1 1 1 1 1]
   beq =1
   lb = zeros(1,10)
   ub = ones(1,10)
   [x,resnorm,residual,exitfag,output] = lsqlin(C,d,A,b,Aeq,beq,lb,ub)