model:
sets:
dmu/1..10/:S,T,P;
II/1..4/:s1;
OI/1..2/:s2;
inv(II,dmu):X;
outv(OI,dmu):Y;
endsets

data:
P=1,0,0,0,0,0,0,0,0,0;
X= 36,36,36,36,35,33,30,30,31,31,
16353,16353,16353,16353,16793,16718,15542,15542,15757,15757,
4688,5604,5906,6036,6203,4528,6071,7024,6886,6151,
0.362,0.319,0.286,0.279,0.271,0.31,0.441,0.525,0.623,0.573;

Y= 23391,21570,30185,28209,25612,22036,20238,18006,17092,15539,
261169,254275,253798,243464,223476,207337,225843,221021,215482,196067;
enddata

max=@sum(dmu:P*T);
@for(dmu(j):S(j)=@sum(II(i):w(i)*X(i,j)));
@for(dmu(j):T(j)=@sum(OI(i):u(i)*Y(i,j)));
S(j)>=T(j);
@sum(dmu:P*S)=1;
end