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

data:
	X =
	 8.939  64.3
	 8.625  99.0
	10.813  99.6
	10.638  96.0
	 6.240  96.2
	 4.719  79.9;
	Y =
	25.2  223
	28.2  287
	29.4  317
	26.4  291
	27.2  295
	25.5  222;
	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;

	-v(1) + 0.3*v(2) <= 0;
	 v(1) -   4*v(2) <= 0;
	-u(1) + 0.3*u(2) <= 0;
	 u(1) -   4*u(2) <= 0;