做规划用lingo

sets:
 cangku /ck1 .. ck6/:a;
 kehu /kh1 ..kh8 /:d;
tp(cangku,kehu):c,x;
endsets

DATA:
 A=60,55,51,43,41,52;
 D=35,37,22,32,41,32,43,38;
 c=6,2,6,7,4,2,5,9
 4,9,5,3,8,5,8,2
 5,2,1,9,7,4,3,3
 7,6,7,3,9,2,7,1
 2,3,9,5,7,2,6,5
 5,5,2,2,8,1,4,3;
ENDDATA
min=@sum(tp(i,j):c(i,j)*x(i,j));
@for(cangku(i):@sum(kehu(j):x(i,j))<=a(i));
@for(kehu(j):@sum(cangku(i):x(i,j))=d(j));