菜鸟求解

丢丢1122

求大神们帮忙检查到底哪出错了，lingo提示错误出现在@for(shijian(i)|i#gt#1: w(i)=w(i-1)+7*x11*30/t(1,1)+7*x21*30/t(2,l)+7*x31*30/t(3,l)+16*x12*30/t(1,2)+16*x22*30/t(2,2)+16*x32*30/t(3,2)-0.4*y11*30/m(1,1)-0.4*y21*30/m(2,1)-0.4*y31*30/m(3,1)-0.4*y41*30/m(4,1)-0.004*yl2*30/m(1,2)-0.004*y22*30/m(2,2)-0.004*y32*30/m(3,2)-0.004*y42*30/m(4,2));
这行，实在找不错哪里不对了


model：
sets:

chuanxing/l,2/:;

chexing/1,2/:;

vendors/1.. 3/:;

customers/1..4/:;

link1(vendors,chuanxing):t,e;

link2(customers,chexing):m,f;

shijian/1..12/:w;

endsets

data:

e=175,352,

462,928,

259,0;

f=3.45,0.052,

4.19,0.063,

6.8,0.1,

12.63,0.19;

t=15,15,

  40,40,

  20,20;

m=1,1,

  1,1,

  2,2,

  3,3;

w=100,,,,,,,,,,,;

enddata

min=490*xll + 1120*xl2+420*x21+960*x22+385*x31+880*x32+xll*30/t(l,l)*e(l,l)+xl2*30/t(l,2)*e(l,2)+x21*30/t(2,l)*e(2,l)+x22*30/t(2,2)*e(2,2)+x31*30/t(3,l)*e(3,l)+x32*30/t(3,2)*e(3,2)+yll*30/m(l,l)*f(l,l)+yl2*30/m(l,2)*f(l,2)+y21*30/m(2,l)*f(2,l)+y22*30/m(2,2)*f(2,2)+y31*30/m(3,l)*f(3,l)+y32*30/m(3,2)*f(3,2)+y41*30/m(4,1 )*f(4,1 )+y42*30/m(4,2)*f(4,2)+0.006+0.000336*(x11*30/t( 1,1 )+x21*30/t(2,l )+x31*30/t(3,l))+0.000768*(xl2*30/t(l,2)+x22*30/t(2,2)+x32*30/t(3,2))-0.0000192*(yll*30/m(l,l)+y21*30/m(2,l)+y31*30/m(3,l)+y41*30/m(4,l))-0.000000192*(yl2*30/m(l,2)+y22*30/m(2,2)+y32*30/m(3,2)+y42*30/m(4,2));

! gongyingyueshu;

7*xll*30/t(l,l)+7*x21*30/t(2,l)+7*x31*30/t(3,l)+16*xl2*30/t(l,2)+16*x22*30/t(2,2)+16*x32*30/t(3,2)<=370.67;

0.4*yll*30/m(l,l)+0.004*yl2*30/m(l,2)>=173.75;

0.4*y21*30/m(2,l)+0.004*y22*30/m(2,2)>=57.92;

0.4*y31*30/m(3,l)+0.004*y32*30/m(3,2)>=115.83;

0.4*y41*30/m(4,1 )+0.004*y42*30/m(4,2)>=23.17;

! gangkouzhuangchenengliyueshu;

yll*30/m(l,l)+y21*30/m(2,l)+y31*30/m(3,l)+y41*30/m(4,l)<=1000;

! kucunnengliyueshu;

7*x11*30/t(l,l)+7*x21*30/t(2,l)+7*x31*30/t(3,l)+16*xl2*30/t(l,2)+16*x22+30/t(2,2)+16*x32*30/t(3,2)-0.4*yll*30/m(l,l)-0.4*y21*30/m(2,l)-0.4*y31*30/m(3,1)-0.4*y41*30/m(4,1)-0.004*yl2*30/m(l,2)-0.004*y22*30/m(2,2)-0.004*y32*30/m(3,2)-0.004*y42*30/m(4,2)<=475;

! anquankucunyueshu;

@for(shijian(i)|i#gt#1: w(i)=w(i-1)+7*x11*30/t(1,1)+7*x21*30/t(2,l)+7*x31*30/t(3,l)+16*x12*30/t(1,2)+16*x22*30/t(2,2)+16*x32*30/t(3,2)-0.4*y11*30/m(1,1)-0.4*y21*30/m(2,1)-0.4*y31*30/m(3,1)-0.4*y41*30/m(4,1)-0.004*yl2*30/m(1,2)-0.004*y22*30/m(2,2)-0.004*y32*30/m(3,2)-0.004*y42*30/m(4,2));

@for(shijian(i)|i#ge#1: w(i)>=125);

@for(shijian(i)|i#ge#1: w(i)<=600);

! zhengshuyueshu;

@gin(xll);

@gin(xl2);

@gin(x21);

@gin(x22);

@gin(x31);

@gin(x32);

@gin(yll);

@gin(yl2);

@gin(y21);

@gin(y22);

@gin(y31);

@gin(y32);

@gin(y41);

@gin(y42);

x31>=l;

x32=0;

x22=@if(x21#le#l#and#x21#ge#0,x22#gt#0,0);

x21=@if(x22#le#l#and#x22#ge#0,x21#gt#0,0);

end