It is io difficult to do this.

obj… R =e= sum((i,j,k),a(j,k)* b(j)*c(j,k)*x(i,j,k)) + sum((i,j,k),a(j,k)* b(j)*c(j,k)*z(i,j,k)) + sum((i,j,k), x(i,j,k)*h(i,j)*m(i,j)) + sum((i,j,k), z(i,j,k)*h(i,j)*m(i,j));
con1(i,k)… sum((j),a(j,k)* b(j)

*x(i,j,k)) =g= d(i,k);*

con2(i,k)… sum((j),a(j,k)b(j)*z(i,j,k)) =l= d(i,k);

con2(i,k)… sum((j),a(j,k)

From con2 (The number of crops in tonnes (left) is less than the demand of factory), I found a lot of problems because the objeective is to minimize cost. It seems con 2 will not enable. So, I would like to merge to one constraint or create sum loop, like if con1 works then con2 not work and when con2 works then con1 not work. Could you recommend how to solve this coding? x(i,j,k) and z(i,j,k) are binary variable that I use for differentiate these 2 constrants. Please let me know if you have any comments or suggestions.

Thanks in advance !