请教一lingo程序

yangwj78

最近写论文，需要用lingo编一个求解最大值的线性规划程序，但是总是出错，请大家帮忙看看。我QQ421200956，欢迎赐教！
model:
sets:
dmu/a1..a30/:s,w;
factor/fac1..fac5/;
do(dmu,factor):f;
endsets
data:
f=
592 126 3913 254 3227
658 141 2445 249 1644
917 294 9151 523 9506
637 208 6626 241 6316
397 97 3050 129 4504
1303 355 9106 490 8963
563 170 3751 209 3008
807 253 5975 382 4053
1389 240 4874 692 4912
1731 579 8118 1152 7318
1076 314 5222 646 5016
439 185 4575 282 3647
560 157 2579 332 2007
331 133 2028 125 1743
1578 548 9012 1198 9089
993 381 7244 529 6618
863 274 6046 570 5105
488 187 4894 240 3452
2791 549 8376 1408 6344
407 108 2438 169 4152
102 13 407 31 327
339 116 3277 108 1713
817 253 6750 463 4606
202 73 4315 108 3347
333 90 3364 277 1887
694 143 3022 172 2497
315 94 2687 122 2772
132 19 739 30 589
95 26 827 29 1041
336 68 3280 130 1790 ;
enddata
@for(dmu(k):
    max=s(k);
    (1-s(k))*f(k,5)-@sum(factor(j)|j # eq # 5:f(i,j)*w(i))=0;
    (1+s(k))*f(k,4)-@sum(factor(j)|j &#35; eq &#35; 4:f(i,j)*w(i))<=0;
    f(k,1)-@sum(factor(j)|j &#35; eq &#35; 1:f(i,j)*w(i))>=0;
    f(k,2)-@sum(factor(j)|j &#35; eq &#35; 2:f(i,j)*w(i))>=0;
    f(k,3)-@sum(factor(j)|j &#35; eq &#35; 3:f(i,j)*w(i))>=0;
    @sum(w(i))=1);
end