sets:
Dmus/1..19/;
Inputs/1..3/: v;
Outputs/1/: u;
DxI(Dmus,Inputs): X;
DxO(Dmus,Outputs): Y;
endsets

data:
X =
179.0 3969.10 1005.5
 96.0  896.32  194.6
 94.0  803.17  200.0  
141.0 1629.90  564.1
147.0 3264.40  617.4
261.0 1427.50  347.1
448.0 6365.30 1680.4
223.0 4911.00 1319.4
279.0 7779.30 1785.6
214.0 2307.50  690.0  
367.0 4334.70 1268.9
 97.0 1870.70  416.5
    143.0 2730.60  779.8
 95.0   540.61  305.0  
241.0  9509.20 1892.1
147.0  1561.10  718.5
306.0  1699.00  306.3
182.0   881.16  182.3
500.0  2548.90  650.1;

    Y =
3531.4
 841.3
 690.5
1498.1
3347.7
1263.2
5734.0  
4107.7
6964.2
1961.0
4207.2
1752.7
2667.6
 457.2
9529.9
1360.6
1689.7
 750.3
2439.5;
enddata

k=1;

max = @sum(Outputs(r): Y(k,r)*u(r));

@for(Dmus(j): 
   @sum(Outputs(r): Y(j,r)*u(r)) 
   - @sum(Inputs(i): X(j,i)*v(i)) <= 0; 
);

@sum(Inputs(i): X(k,i)*v(i)) = 1;

@for(Inputs(i):  v(i) >= 0.000001; );
@for(Outputs(r): u(r) >= 0.000001; );