麻烦哪位大神可以帮忙看看这个lingo代码哪出问题了吗

麻烦哪位大神可以帮忙看看这个代码哪出问题了吗

model:
!集合定义;
sets:
!p:每个回收点的回收量,f:回收点到处理中心的单位运输成本;
consumer/c1…c16/:p,f;
!a:处理中心固定成本,b:单位运营成本,c:容量限制,z:处理中心是否被选;
warehouse/w1…w16/:a,b,c,z;
!f:回收点到处理中心的单位运输成本,v:从回收点运往处理中心的运量;
links(consumer,warehouse):v;
endsets
!给出已知数据;
data:
p=42 30 15 30 30 30 45 36 15 15 30 36 60 36 9 9;
f=0 3 4 5 14 15 20 6.6 10 8.7 10 6.3 9.5 15 13 17
3 0 6 8 20 21 27 8.9 16 12 17 4.1 12 18 17 26
4 6 0 7.3 11 12 18 7.1 11 9.1 6.7 8.2 10 16 13 14
5 8 7.3 0 18 21 38 9.4 8.9 18 14 5.6 19 21 20 23
14 20 11 18 0 3 17 10 12 20 11 18 21 34 22 25
15 21 12 21 3 0 18 6.1 7.2 21 13 19 22 32 20 23
20 27 18 38 17 18 0 19 21 20 15 25 20 25 10 10
6.6 8.9 7.1 9.4 10 6.1 19 0 2.7 17 8.6 11 18 23 15 18
10 16 11 8.9 12 7.2 24 2.7 0 20 12 15 20 25 19 21
8.7 12 9.1 18 20 21 20 17 20 0 12 16 1.1 6.8 12 14
10 17 6.7 14 11 13 15 8.6 12 12 0 15 12 17 6.2 9.6
6.3 4.1 8.2 5.6 18 19 25 11 15 16 15 0 12 18 18 21
9.5 12 10 19 21 22 20 18 20 1.1 12 12 0 6 12 13
15 18 16 21 34 32 25 23 25 6.8 17 18 6 0 13 13
13 17 13 20 22 20 10 15 19 12 6.2 18 12 13 0 4.5
17 26 14 23 25 23 10 18 21 14 9.6 21 13 13 4.5 0;
a=120000 120000 120000 120000 120000 120000;
b=30 30 30 30 30 30;
c=85 85 85 85 85 85;
enddata
!目标函数;
[email protected](link1(i,j):v(i,j)*f(i,j))[email protected](link2(i,j):v(i,j)*b(j)[email protected](warhouse(j):a(j)*z(j));!!!;
@for(warehouse(j):
@sum(consumer(i):v(i,j)=p(i));
@for(consumer(i):
@sum(warehouse(j):v(i,j)<=c(j));
@sum(warehouse(j)????(j))=6;
@bin(z);!!!;
end

一直改不了麻烦哪位大神可以帮忙看看这个lingo代码哪出问题了吗