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,
16353163, 53163, 53163, 53167, 93167, 18155, 42155, 42157, 57157, 57,
4688560 4590 6603 6620 3452 8607 1702 4688 6615 1,
0.362, 0.319, 0.286, 0.279, 0.271, 0.31, 0.441, 0.525, 0.623, 0.573;

Y=23391215, 70301, 85282, 09256, 12220, 36202, 38180, 06170, 92155, 39
261169254 275253 798243 464223 476207 337225 843221 021215 482196 067;
Enddata

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