function f=ld(x)
A=[1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 32
15 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 10 1 0 0 0 0 2 5
17 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5 1 0 1 0 0 1 9 0 0 0 0 0 0 0
1 1 0 0 0 0 1 13 1 0 0 0 0 2 8 2 0 0 0 0 2 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0
12 0 0 1 0 0 0 0 0 0 0 0 0 0 7 2 0 0 0 0 2 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0
15 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 18 0 0 0 0 0 0 0
11 1 0 1 0 0 0 3 1 0 0 0 0 0 4 2 0 0 1 0 1 10 0 0 0 0 0 0 0 0 0 0 0 0 0 0
1 0 1 0 0 0 0 0 0 0 0 0 1 0 3 2 0 0 0 0 3 20 1 0 0 0 0 1 2 0 0 0 0 0 0 0
1 0 0 0 1 0 0 0 0 0 0 0 0 0 4 1 0 0 0 0 0 2 1 0 1 0 0 0 24 0 0 0 0 0 0 0
3 0 0 0 1 0 0 0 0 0 0 0 0 0 4 0 0 0 0 0 1 25 0 0 0 0 0 1 0 0 0 0 0 0 0 0
1 0 0 1 0 0 0 1 1 1 0 0 0 0 2 1 0 0 1 1 0 10 3 0 0 0 0 3 9 0 0 0 0 0 0 0
5 0 0 0 0 0 0 1 1 0 1 0 0 0 0 1 0 0 0 0 0 5 2 0 1 1 0 1 16 0 0 0 0 0 0 0
1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 16 3 1 0 0 0 2 9
1 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 7 3 0 1 0 0 2 18 0 0 0 0 1 0 0
14 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 10 1 1 1 0 0 0 7 0 0 0 0 1 0 0
2 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5 1 1 1 0 0 1 22
];
for j=1:16
for i=1:36
B(j,i)=x(i)*A(j,i);
end
end
f=sum(B(:));
f=-f;
Aeq=[ones(1,6) zeros(1,30)
zeros(1,6) ones(1,6) zeros(1,24)
zeros(1,12) ones(1,6) zeros(1,18)
zeros(1,18) ones(1,6) zeros(1,12)
zeros(1,24) ones(1,6) zeros(1,6)
zeros(1,30) ones(1,6)];
beq=[1
1
1
1
1
1];
x0=[0.999999999994889 0.000000000001012 0.000000000001003 0.000000000001022 0.000000000001052 0.000000000001022 0.000000000008805 0.999999999883916 0.000000000080863 0.000000000008805 0.000000000008805 0.000000000008805 0.000000000001417 0.000000000001502 0.999999999990798 0.000000000003448 0.000000000001417 0.000000000001417 0.00000000000095 0.000000000000942 0.000000000001011 0.999999999995115 0.00000000000104 0.000000000000942 0.000000000000833 0.000000000000797 0.000000000000797 0.000000000000846 0.999999999995906 0.000000000000821 0.000000000001329 0.000000000001295 0.000000000001278 0.000000000001312 0.000000000001365 0.999999999993421
];
lb=zeros(1,36);
options = optimset('LargeScale','off')
[x(1,:),fval]=fmincon(@ld,x0,[],[],Aeq,beq,lb,[],[],options)