model:
sets:
inputs/ip1..ip3/;
badoutputs/bo1..bo2/;
provinces/p1..p5/;
links1(inputs,provinces):xt;
links2(badoutputs,provinces):bt;
goodoutputs(provinces):yt;
weighting(provinces):zt;
endsets
max=y_t;
@sum(provinces(K):zt(K)*yt(K))>=y_t;
@for(provinces(K):@for(badoutputs(J):@sum(provinces(K):zt(K)*bt(J,K))=bt(J,K)));
@for(provinces(K):@for(inputs(N):@sum(provinces(K):zt(K)*xt(N,K))<=xt(N,K)));
@for(provinces(K):zt(K)>=0);
data:
xt=337.33       240.41       475.86       309.73       189.76     

566.50       455.98       2555.43       1154.10       856.60     

2291.00       1680.00       4918.00       4395.00       1871.00     
bt=30.97       21.99       63.46       47.45       37.99     

6087.17       4254.68        13769.51        12185.54        6108.83      
yt= 211.91     
  156.13     
  325.10     
  170.34     
  127.04     
；

enddata
end