model:
sets:
sj/1..3/:a,c;
si/1..10/:y;
sji(sj,si):x;
endsets
data:
y=29118.40,32092.00,34946.30,38564.30,37783.60,39610.10,40141.90,42484.20,24408.19,26971.44;
x=1316.08,1518.17,1780.42,2163.80,2695.82,3206.58,3786.52,4436.18,4854.00,5666.20
330.00,590.90,740.00,890.40,1006.70,1199.20,1464.80,1730.20,1961.10,2302.40
1052.88,1733.88,2559.38,3330.14,4297.60,5264.87,6685.20,8116.56,9521.57,11858.20;
enddata
min=@sum(sj(j):c(j)*x(j,i));
@for(si(i):
@sum(sj(j):a(j)*x(j,i))-0.2*@sum(sj(j):c(j)*x(j,i))<=y(i));
@for(si(i):
@sum(sj(j):a(j)*x(j,i))+0.2*@sum(sj(j):c(j)*x(j,i))>=y(i));
end