这个lingo程序到底错哪了！？！？！？

Question

原题：
!对某厂一，二，三产品下一年各季度的合同预订数如下
各季度订货数为：
一产品1500,1000,2000,1200
二产品1500,1500,1200,1500
三产品1000,2000,1500,2500
该三种产品1季度初无库存
要求在4季度末各库存150件
已知该厂每季度生产工时为15000小时
生产一，二，三产品每件分别需要2、4、3小时
因更换工艺装备，产品一在2季度无法生产
规定当产品不能按期交货时，产品每件每迟交
产品一一个季度赔偿20元，产品二、三赔10元
又生产出来产品不在本季度交货的，每件每季度的库存费用为5元
问该厂应如何安排生产，使总的赔偿加库存的费用为最小。;
-----------------------------
model:

sets:
season/1..4/:;
product/1..3/:r,a;!a是生产用时,r赔偿;
link(season,product):x,p,s;!p销量,s是余量;
endsets

min=@sum(season(i):
@sum(season(j):
@if(s(i,j) #ge# 0, 5*s(i,j), -r(j)*s(i,j))));

@for(season(i):@sum(product(j):a(j)*x(i,j))<=15000);
x21=0;
@for(product(j):s(4,j)=150);
@for(link(i,j):s(i,j)=@if(i #eq# 1,x(i,j)-p(i,j),x(i,j)-p(i,j)+s(i-1,j)));
@for(link(i,j):@free(s);@gin(x));

data:
p=
1500 1500 1000
1000 1500 2000
2000 1200 1500
1200 1500 2500;
r=20 10 10;
a=2 4 3;
enddata

end

Answer 1

model:
sets:
season/1..4/:;
product/1..3/:r,a;!a是生产用时,r赔偿;
link(season,product):x,p,s;!p销量,s是余量;
endsets

min=@sum(season(i):
@sum(product(j): !将season(j)改成product(j);
@if(s(i,j) #ge# 0, 5*s(i,j), -r(j)*s(i,j))));

@for(season(i):@sum(product(j):a(j)*x(i,j))<=15000);
x21=0;
@for(product(j):s(4,j)=150);
@for(link(i,j):s(i,j)=@if(i #eq# 1,x(i,j)-p(i,j),x(i,j)-p(i,j)+s(i-1,j)));
@for(link(i,j):@free(s);@gin(x));

data:
p=
1500 1500 1000
1000 1500 2000
2000 1200 1500
1200 1500 2500;
r=20 10 10;
a=2 4 3;
enddata

end