那位大侠知道这个lingo 程序错在哪

ljwabc115

model:
title
sets:
   supplynumber/1/:w;
   distributionnumber/1..5/:z,f,ml;
   clientnumber/1..10/:h;
   productcategory/1..2/;
   link21(supplynumber,distributionnumber):c;
   link22(distributionnumber,clientnumber):d;
   link23(productcategory,distributionnumber):m;
   link24(productcategory,supplynumber):a;
   link31(productcategory,supplynumber,distributionnumber):x;
   link32(productcategory,distributionnumber,clientnumber):y;
endsets
data:
   e=1;
   c=30 35 20 26 39;
   d=7  4  6  7  5  3  1  2  2  3
     5  2  1  8  6  2  3  1  7  9
     2  6  5  1  1  8  4  3  4  6
     3  9  3  4  5  4  7  7  3  4
     2  2  2  2  7  4  3  9  5  8;
   m=35  25  30  50  60
     25  40  20  20  40;
   f=3000  3500  3100  2500  3200;
   a=300
     300;
   ml=700  800  700  600  700;
enddata
   min=@sum(link31(r,k,i):c(k,i)*e*x(r,k,i)+
       @sum(link32(r,i,j):d(i,j)*e*y(r,i,j)+
       @sum(link31(r,k,i):m(r,i)*x(r,k,i))+
       @sum(distributionnumber(i):z(i)*f(i));
       @for(supplynumber(k):
       @sum(distributionnumber(i):x(1,k,i))<=a(1,k));
       @for(supplynumber(k):
       @sum(distributionnumber(i):x(2,k,i))<=a(2,k));
       @for(clientnumber(j):
       @sum(distributionnumber(i):y(1,i,j))>=d(1,j));
       @for(clientnumber(j):
       @sum(distributionnumber(i):y(2,i,j))>=d(2,j));
       @for(distributionnumber(i):
       @sum(productcategory(r):
       @sum(clientnumber(j):y(r,i,j)))<=z(i)*ml(i));
       @for(distributionnumber(i):
       @sum(productcategory(r):
       @sum(supplynumber(k):y(r,i,j)))<=z(i)*ml(i));
       @for(productcategory(r):
       @for(distributionnumber(i):
       @sum(supplynumber(k):x(r,k,i)=@sum(clientnumber(j):y(r,i,l)))));
       @for(<distributionnumber(i):@bin(z(i)));
end