lingo中属性下标越界

MODEL:
sets:
gd/1..12/:a,b;
ch/1..13/:x,y;
links(gd,ch):;

endsets
data:
s=0;
a=0.50 1.41 3.00 3.37 3.40 4.72 4.72 5.43 7.57 8.38 8.98 9.50;
b=2.00 3.50 1.50 3.51 5.50 2.00 6.24 4.10 2.01 4.50 3.41 0.80;
enddata
f1=@if((a(i)-x(j)#LE#x(j+1)-b(i))#AND#(b(i)-y(k)#LE#y(k+1)-b(i)),@IF((a(i)-x(j)#LE#0.05)#AND#(b(i)-y(k)#LE#0.05),1,0),@if((a(i)-x(j)#LE#x(j+1)-a(i))#AND#(b(i)-y(k)#GT#y(k+1)-b(i)),@IF((a(i)-x(j)#LE#0.05)#AND#(y(k+1)-b(i)#LE#0.05),1,0),@if((a(i)-x(j)#GT#x(j+1)-a(i))#AND#(b(i)-y(k)#LE#y(k+1)-b(i)),@IF((x(j+1)-a(i)#LE#0.05)#AND#(b(i)-y(k)#LE#0.05),1,0),@if((a(i)-x(j)#GT#x(j+1)-a(i))#AND#(b(i)-y(k)#GT#y(k+1)-b(i)),@IF((x(j+1)-a(i)#LE#0.05)#AND#(y(k+1)-b(i)#LE#0.05),1,0),0))));
max=@sum(links(gd,ch):f1);
@for(links(i,j):0<=a(i)-x(j));
@for(links(i,j):a(i)-x(j)<=1);
@for(links(i,k):0<=b(i)-y(k));
@for(links(i,k):b(i)-y(k)<=1);
@for(ch(j):x(j+1)-x(j)=1);
@for(ch(k):y(k+1)-y(k)=1);
end

勘探部门在某地区找矿。初步勘探时期已零散地在若干位置上钻井，取得了地质资料。进入系统勘探时期后，要在一个区域内按纵横等距的网格点来布置井位，进行“撒网式”全面钻探。由于钻一口井的费用很高，如果新设计的井位与原有井位重合（或相当接近），便可利用旧井的地质资料，不必打这口新井。因此，应该尽量利用旧井，少打新井，以节约钻探费用。比如钻一口新井的费用为500万元，利用旧井资料的费用为10万元，则利用一口旧井就节约费用490万元。
设平面上有n个点pi，其坐标为(ai,bi),i=1,2,…,n，表示已有的n个井位。新布置的井位是一个正方形网格n的所有结点（所谓“正方形网格”是指每个格子都是正方形的网格；结点是指纵线和横线的交叉点）。假定每个格子的边长（井位的纵横间距）都是1单位（比如100米）。整个网格是可以在平面上任意移动的。若一个已知点pi与某个网格结点xi的距离不超过给定误差ε（=0.05单位），则认为pi处的旧井资料可以利用，不必在结点xi处打新井。

第一句有大错 根本看不出你是什么意思 没法改 其它的在这里
MODEL:
sets:
gd/1..12/:a,b;
ch/1..13/:x,y;
links(gd,ch):f1;

endsets
data:
s=0;
a=0.50 1.41 3.00 3.37 3.40 4.72 4.72 5.43 7.57 8.38 8.98 9.50;
b=2.00 3.50 1.50 3.51 5.50 2.00 6.24 4.10 2.01 4.50 3.41 0.80;
enddata
!@for(links(i,j):f1(i,j)=@if((a(i)-x(j)#LE#x(j+1)-b(i))#AND#(b(i)-y(k)#LE#y(k+1)-b(i)),@IF((a(i)-x(j)#LE#0.05)#AND#(b(i)-y(k)#LE#0.05),1,0),@if((a(i)-x(j)#LE#x(j+1)-a(i))#AND#(b(i)-y(k)#GT#y(k+1)-b(i)),@IF((a(i)-x(j)#LE#0.05)#AND#(y(k+1)-b(i)#LE#0.05),1,0),@if((a(i)-x(j)#GT#x(j+1)-a(i))#AND#(b(i)-y(k)#LE#y(k+1)-b(i)),@IF((x(j+1)-a(i)#LE#0.05)#AND#(b(i)-y(k)#LE#0.05),1,0),@if((a(i)-x(j)#GT#x(j+1)-a(i))#AND#(b(i)-y(k)#GT#y(k+1)-b(i)),@IF((x(j+1)-a(i)#LE#0.05)#AND#(y(k+1)-b(i)#LE#0.05),1,0),0)))));
max=@sum(links:f1);
@for(links(i,j):0<=a(i)-x(j));
@for(links(i,j):a(i)-x(j)<=1);
@for(links(i,k):0<=b(i)-y(k));
@for(links(i,k):b(i)-y(k)<=1);
@for(ch(j)|j#lt#13:x(j+1)-x(j)=1);
@for(ch(k)|k#lt#13:y(k+1)-y(k)=1);
end