哪位帮忙看看这个lingo程序的毛病，自己建的模型，求赐教！

pangpangpang3

sets:
Layer/1..3/;
subedge/1..3/:v;
bian/edge1..edge16/:C;
sink/1,2,3/:T;
path/N1,N2,N3/;
!邻接边;
!adjacentpair(bian,bian)/edge1,edge4 edge1,edge5 edge2,edge6 edge2,edge7
edge2,edge8 edge3,edge9 edge3,edge10 edge5,edge11 edge6,edge11 edge8,edge12
edge9,edge12 edge11,edge13 edge11,edge14 edge12,edge15 edge12,edge16/:AD;
adjacentpair(bian,bian):AD;
!各信宿的路径;
PathXSink(sink,path):Pathsink;
!各边是否属于该路径;
PathXSinkXEdge(sink,Path,bian):x;
!各路径上不同层的流;
FlowXPathXSink(path,sink,Layer):y;
!该信宿是否能解码改成信息;
ZXLayer(sink,Layer):z;
BXLayer(Layer):b;
!虚拟子链路消耗的路径容量;
UseXPathXLayer(sink,path,Layer,bian,subedge):w;
!虚拟链路的码率;
GXSubedge(bian,subedge):u;
endsets
data:
AD=0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 1 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 0 0 0 0 0 0 0 0 1 0 0 0 0 0
   0 0 0 0 0 0 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 0 0 0 0 0 0 0 0 1 0 0 0 0
   0 0 0 0 0 0 0 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 0 0 0 0 0 0 0 0 1 1 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 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 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
  
Pathsink=1 1 0
         1 1 1
         1 1 0;
x=1 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0
  0 1 0 0 0 1 0 0 0 0 1 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 1 0 0 0 0 0 1 0 0 1 0 0
  0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0
  0 0 1 0 0 0 0 0 1 0 0 1 0 0 1 0
  0 1 0 0 0 0 0 1 0 0 0 1 0 0 0 1
  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 0 0 0 0 0 0;
C= 1 3 3 1
   1 1 3 1
   2 2 1 2
   1 1 2 2;  
enddata
[obj]max=@sum(ZXLayer(i,l):z(i,l)*b(l));

@for(sink(i):
     @for(layer(j)|j#LE#L-1:z(i,j)>=z(i,j+1));
     @for(layer(j):
              @sum(path(k):y(k,i,j))=b(l)*z(i,l);
               @for(path(d):
                    @for(bian(e):
                      @sum(subedge(r)|(r#GE#j)#and#(z(i,j)#EQ#1):w(i,p,j,e,r))=x(i,p,e)*y(k,i,j);
                      @for(bian(f):
                          @for(subedge(r)|AD(e,f)#eq#1:w(i,p,j,f,r)<=w(i,p,j,f,r)););););););
@for(sink(i):
@for(bian(e):
  @sum(subedge(r):u(e,r))<=c(e);     
  @for(layer(j):
    @for(subedge(r):
        @sum(subedge(h)|h#LE#r:@sum(path(Nt):w(i,Nt,j,e,h)))<=u(e,h);););););   
               
@for(ZXLayer:@bin(z));

@for(FlowXPathXSink(Nt,i,j):@gin(y));
            
@for(GXSubedge(e,r):@gin(u));

@for(UseXPathXLayer(i3,d,j3,e,r):@gin(w));
           
@for(BXLayer(j4):@gin(b));
         ! @bnd(1,b,3);!);
end

pangpangpang3

不用看了，已经查出来错误了，结贴！