model:
sets:
dmu/1..6/:m;
ii/1..4/:s;
oi/1..6/:a;
iv(ii,dmu):x;
ov(oi,dmu):y;
endsets
data:
x=100 98 96 105 70 95
30 43 28 30 18 26
25 30 24 25 32 37
85 78 75 81 91 82;
y=8 5 7 6 9 7
9 8 7 5 7 6
7 6 5 8 6 6
9 7 8 6 8 6
5 2 6 5 7 4
8 6 7 8 9 5;
enddata
min=v-e-(s1+s2+s3+s4+a1+a2+a3+a4+a5+a6);
@for(dmu(n):
@for(ii(i):@sum(dmu(t):m(t)*x(i,t)+s(i))=x(i,2));
@for(oi(j):@sum(dmu(t):m(t)*y(j,t)-a(j))=v*y(j,2)));
@sum(dmu(n):m(n))=1;
End