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

data:
  X =
  1177  187   794
  1137 2157  5473
   970 1565  7128
   620  454  2094
  2055 1974  7903
   492  536  2171
   850  874  6030
   326  385  3115
  1248.146 1728  4462
  1222 1576 11192
   785.660 1917  1835
   492  455  1631
  1403 1379  4322
   511  420  2736;

  Y =
    2  26  42  122
  110 241 519  626
   98 176 213  881
   24  41 117  274
   37 253 330  744
   28  52 224  333
  170 107 108  644
  240  48 360  321
  184 180 579  532
   33 142 456  710
   62 286 133  129
   36  53 247  270
   67 175 331  729
   71  55 237  304;
  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;

  -0.365 * v(1) + 1.00  * v(2) <= 0;
   0.095 * v(1) - 1.00  * v(2) <= 0;
  -0.037 * v(1) + 1.00  * v(3) <= 0;
   0.009 * v(1) - 1.00  * v(3) <= 0;
  -0.037 * v(2) + 0.095 * v(3) <= 0;
   0.009 * v(2) - 0.365 * v(3) <= 0;
   
  -19.6  * u(1) + 1.1   * u(2) <= 0;
   14.68 * u(1) - 4.04  * u(2) <= 0;
  -0.412 * u(1) + 1.1   * u(3) <= 0;
   0.104 * u(1) - 4.04  * u(3) <= 0;
  -2.0   * u(1) + 1.1   * u(4) <= 0;
   1.34  * u(1) - 4.04  * u(4) <= 0;
  -0.412 * u(2) + 14.68 * u(3) <= 0;
   0.104 * u(2) - 19.6  * u(3) <= 0;
  -2.0   * u(2) + 14.68 * u(4) <= 0;
   1.34  * u(2) - 19.6  * u(4) <= 0;
  -2.0   * u(3) + 0.104 * u(4) <= 0;
   1.34  * u(3) - 0.412 * u(4) <= 0;