查看: 2626|回复: 4

请教一lingo程序，谢谢！

[复制链接]

字体大小: 正常放大

yangwj78

2 主题	4 听众	87 积分

升级 86.32%

该用户从未签到

电梯直达

1^#

发表于 2009-6-30 09:34 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

由于刚接触lingo，很多东西都不明白，编的程序总是不对，请大家指正，谢谢！
model:
sets:
dmu/a1..a39/:ee;
out/o1..o4/;
do(dmu,out):f,w;
endsets
data:
n=4;
f=41.374 0.76 0.01214 10.7516
7.5057 0.692 0.00089 0.05207
64.977 2.85 0.01989 54.8142
160.75 1.252 0.02476 53.4763
39.785 5.563 0.02646 5.57088
408.5  1.298 0.09064 17.1816
391.58 19.82 0.01033 0.20892
188.83 29.45 0.00955 0.76588
250.24 4.361 0.01199 0.05094
5.3832 0.058 0.00084 0.30294
115.85 5.526 0.01165 0.0156
8.173  0.345 0.00093 0.04326
92.117 9.302 0.00211 0.20876
83.279 5.697 0.01185 0.03052
5.7716 0.076 0.00156 4.93978
2018.5 39.17 0.04799 0.05951
6.2366 0.179 0.00086 0.00539
2.8194 0.132 0.00042 5.47648
42.478 10.07 0.03094 5.20861
168.74 65.45 0.03044 0.09468
164.89 4.455 0.00601 0.25325
322.08 9.937 0.03895 0.24721
18.543 1.284 0.01 0.08281
3.064  0.26 0.00113 2.90088
21.045 0.93 0.08005 0.73361
43.218 4.642 0.02656 1.09635
35.178 10.17 0.0581 0.0869
11.856 0.793 0.00219 0.16704
7.3769 0.264 0.00224 0.14108
18.473 4.966 0.00324 0.33516
17.068 2.088 0.00218 0.01928
5.4787 0.196 0.00081 0.02805
3.0761 0.435 0.00047 0.01578
19.818 0.757 0.00199 0.01822
5.1734 0.343 0.00085 0.02332
109.35 11.87 0.01363 0.8445
32.229 0.812 0.21538 4.7541
59.355 33.63 0.03417 1.57874
133.09 3.28 0.0013 0.02362;
enddata
max=1/@sum(dmu:ee);
@for(dmu(i)
ee=@sum(out(j)|j #le# n:f(i,j)*w(i,j);
@sum(do:w)=1;
@for(dmu(k):
@sum(out(j)|j #le# n:f(i,j)*w(i,j))>=1
)
);
end

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

madio

3万主题	1311 听众	5万积分

TA的每日心情

	奋斗 2024-7-1 22:21

签到天数: 2014 天

[LV.Master]伴坛终老

自我介绍: 数学中国站长

群组: 数学建模培训课堂1

群组: 数学中国美赛辅助报名

群组: Matlab讨论组

群组: 2013认证赛A题讨论群组

群组: 2013认证赛C题讨论群组

2^#

发表于 2009-6-30 23:09 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

给你修改了一下，没有错误了，但是我这里的是测试版没有运行出结果，你去正式版上运行一下看看！

model:
sets:
dmu/a1..a39/:ee;
out/o1..o4/;
do(dmu,out):f,w;
endsets
data:
n=4;
f=41.374 0.76 0.01214 10.7516
7.5057 0.692 0.00089 0.05207
64.977 2.85 0.01989 54.8142
160.75 1.252 0.02476 53.4763
39.785 5.563 0.02646 5.57088
408.5  1.298 0.09064 17.1816
391.58 19.82 0.01033 0.20892
188.83 29.45 0.00955 0.76588
250.24 4.361 0.01199 0.05094
5.3832 0.058 0.00084 0.30294
115.85 5.526 0.01165 0.0156
8.173  0.345 0.00093 0.04326
92.117 9.302 0.00211 0.20876
83.279 5.697 0.01185 0.03052
5.7716 0.076 0.00156 4.93978
2018.5 39.17 0.04799 0.05951
6.2366 0.179 0.00086 0.00539
2.8194 0.132 0.00042 5.47648
42.478 10.07 0.03094 5.20861
168.74 65.45 0.03044 0.09468
164.89 4.455 0.00601 0.25325
322.08 9.937 0.03895 0.24721
18.543 1.284 0.01 0.08281
3.064  0.26 0.00113 2.90088
21.045 0.93 0.08005 0.73361
43.218 4.642 0.02656 1.09635
35.178 10.17 0.0581 0.0869
11.856 0.793 0.00219 0.16704
7.3769 0.264 0.00224 0.14108
18.473 4.966 0.00324 0.33516
17.068 2.088 0.00218 0.01928
5.4787 0.196 0.00081 0.02805
3.0761 0.435 0.00047 0.01578
19.818 0.757 0.00199 0.01822
5.1734 0.343 0.00085 0.02332
109.35 11.87 0.01363 0.8445
32.229 0.812 0.21538 4.7541
59.355 33.63 0.03417 1.57874
133.09 3.28 0.0013 0.02362;
enddata
max=1/@sum(dmu:ee);
@for(dmu(i):
ee=@sum(out(j)|j #le# n:f(i,j)*w(i,j));
@sum(do:w)=1;
@for(dmu(k):
      @sum(out(j)|j #le# n:f(i,j)*w(i,j))>=1
)
);
end

我也说一句

发表

使用道具举报

yangwj78

2 主题	4 听众	87 积分

升级 86.32%

该用户从未签到

3^#

发表于 2009-7-1 08:28 |只看该作者 |招呼Ta 关注Ta

谢谢兄弟，我运行了下，结果不对（每个DMU的EE都等于1）。扩展模型如下：
根据扩展模型（部分），有3个地方还需要修改：（1）目标函数有误：我的本意是目标函数为每个DMU的EE的倒数取最大值,而不是39项之和的倒数取最大值（此问题我一直解决不了）；（2）权约束有误：模型的本意是对每个DMU，W的4项和为1（这个问题我已解决；（3）循环时DMU没有按照想象的进行（此问题我可以解决）。请兄弟再帮忙看看。
MODEL:
[_1] MAX= 1 / ( EE_A1 + EE_A2 + EE_A3 + EE_A4 + EE_A5 +
   EE_A6 + EE_A7 + EE_A8 + EE_A9 + EE_A10 + EE_A11 + EE_A12
   + EE_A13 + EE_A14 + EE_A15 + EE_A16 + EE_A17 + EE_A18 +
   EE_A19 + EE_A20 + EE_A21 + EE_A22 + EE_A23 + EE_A24 +
   EE_A25 + EE_A26 + EE_A27 + EE_A28 + EE_A29 + EE_A30 +
   EE_A31 + EE_A32 + EE_A33 + EE_A34 + EE_A35 + EE_A36 +
   EE_A37 + EE_A38 + EE_A39 ) ;
[_2] - 41.374 * W_A1_O1 - 0.76 * W_A1_O2 - 0.01214 * W_A1_O3 - 10.7516 *
W_A1_O4 + EE_A1 = 0 ;
[_3] W_A1_O1 + W_A1_O2 + W_A1_O3 + W_A1_O4 + W_A2_O1 + W_A2_O2 + W_A2_O3
+ W_A2_O4 + W_A3_O1 + W_A3_O2 + W_A3_O3 + W_A3_O4 + W_A4_O1 + W_A4_O2 +
W_A4_O3 + W_A4_O4 + W_A5_O1 + W_A5_O2 + W_A5_O3 + W_A5_O4 + W_A6_O1 +
W_A6_O2 + W_A6_O3 + W_A6_O4 + W_A7_O1 + W_A7_O2 + W_A7_O3 + W_A7_O4 +
W_A8_O1 + W_A8_O2 + W_A8_O3 + W_A8_O4 + W_A9_O1 + W_A9_O2 + W_A9_O3 +
W_A9_O4 + W_A10_O1 + W_A10_O2 + W_A10_O3 + W_A10_O4 + W_A11_O1 +
W_A11_O2 + W_A11_O3 + W_A11_O4 + W_A12_O1 + W_A12_O2 + W_A12_O3 +
W_A12_O4 + W_A13_O1 + W_A13_O2 + W_A13_O3 + W_A13_O4 + W_A14_O1 +
W_A14_O2 + W_A14_O3 + W_A14_O4 + W_A15_O1 + W_A15_O2 + W_A15_O3 +
W_A15_O4 + W_A16_O1 + W_A16_O2 + W_A16_O3 + W_A16_O4 + W_A17_O1 +
W_A17_O2 + W_A17_O3 + W_A17_O4 + W_A18_O1 + W_A18_O2 + W_A18_O3 +
W_A18_O4 + W_A19_O1 + W_A19_O2 + W_A19_O3 + W_A19_O4 + W_A20_O1 +
W_A20_O2 + W_A20_O3 + W_A20_O4 + W_A21_O1 + W_A21_O2 + W_A21_O3 +
W_A21_O4 + W_A22_O1 + W_A22_O2 + W_A22_O3 + W_A22_O4 + W_A23_O1 +
W_A23_O2 + W_A23_O3 + W_A23_O4 + W_A24_O1 + W_A24_O2 + W_A24_O3 +
W_A24_O4 + W_A25_O1 + W_A25_O2 + W_A25_O3 + W_A25_O4 + W_A26_O1 +
W_A26_O2 + W_A26_O3 + W_A26_O4 + W_A27_O1 + W_A27_O2 + W_A27_O3 +
W_A27_O4 + W_A28_O1 + W_A28_O2 + W_A28_O3 + W_A28_O4 + W_A29_O1 +
W_A29_O2 + W_A29_O3 + W_A29_O4 + W_A30_O1 + W_A30_O2 + W_A30_O3 +
W_A30_O4 + W_A31_O1 + W_A31_O2 + W_A31_O3 + W_A31_O4 + W_A32_O1 +
W_A32_O2 + W_A32_O3 + W_A32_O4 + W_A33_O1 + W_A33_O2 + W_A33_O3 +
W_A33_O4 + W_A34_O1 + W_A34_O2 + W_A34_O3 + W_A34_O4 + W_A35_O1 +
W_A35_O2 + W_A35_O3 + W_A35_O4 + W_A36_O1 + W_A36_O2 + W_A36_O3 +
W_A36_O4 + W_A37_O1 + W_A37_O2 + W_A37_O3 + W_A37_O4 + W_A38_O1 +
W_A38_O2 + W_A38_O3 + W_A38_O4 + W_A39_O1 + W_A39_O2 + W_A39_O3 +
W_A39_O4 = 1 ;
[_4] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_5] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_6] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_7] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_8] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_9] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_10] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_11] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_12] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_13] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_14] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_15] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_16] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_17] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_18] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_19] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_20] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_21] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_22] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_23] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_24] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_25] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_26] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_27] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_28] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_29] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_30] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_31] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_32] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_33] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_34] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_35] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_36] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_37] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_38] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_39] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_40] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_41] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_42] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_43] - 7.5057 * W_A2_O1 - 0.692 * W_A2_O2 - 0.00089 * W_A2_O3 - 0.05207
* W_A2_O4 + EE_A2 = 0 ;
[_44] W_A1_O1 + W_A1_O2 + W_A1_O3 + W_A1_O4 + W_A2_O1 + W_A2_O2 +
W_A2_O3 + W_A2_O4 + W_A3_O1 + W_A3_O2 + W_A3_O3 + W_A3_O4 + W_A4_O1 +
W_A4_O2 + W_A4_O3 + W_A4_O4 + W_A5_O1 + W_A5_O2 + W_A5_O3 + W_A5_O4 +
W_A6_O1 + W_A6_O2 + W_A6_O3 + W_A6_O4 + W_A7_O1 + W_A7_O2 + W_A7_O3 +
W_A7_O4 + W_A8_O1 + W_A8_O2 + W_A8_O3 + W_A8_O4 + W_A9_O1 + W_A9_O2 +
W_A9_O3 + W_A9_O4 + W_A10_O1 + W_A10_O2 + W_A10_O3 + W_A10_O4 + W_A11_O1
+ W_A11_O2 + W_A11_O3 + W_A11_O4 + W_A12_O1 + W_A12_O2 + W_A12_O3 +
W_A12_O4 + W_A13_O1 + W_A13_O2 + W_A13_O3 + W_A13_O4 + W_A14_O1 +
W_A14_O2 + W_A14_O3 + W_A14_O4 + W_A15_O1 + W_A15_O2 + W_A15_O3 +
W_A15_O4 + W_A16_O1 + W_A16_O2 + W_A16_O3 + W_A16_O4 + W_A17_O1 +
W_A17_O2 + W_A17_O3 + W_A17_O4 + W_A18_O1 + W_A18_O2 + W_A18_O3 +
W_A18_O4 + W_A19_O1 + W_A19_O2 + W_A19_O3 + W_A19_O4 + W_A20_O1 +
W_A20_O2 + W_A20_O3 + W_A20_O4 + W_A21_O1 + W_A21_O2 + W_A21_O3 +
W_A21_O4 + W_A22_O1 + W_A22_O2 + W_A22_O3 + W_A22_O4 + W_A23_O1 +
W_A23_O2 + W_A23_O3 + W_A23_O4 + W_A24_O1 + W_A24_O2 + W_A24_O3 +
W_A24_O4 + W_A25_O1 + W_A25_O2 + W_A25_O3 + W_A25_O4 + W_A26_O1 +
W_A26_O2 + W_A26_O3 + W_A26_O4 + W_A27_O1 + W_A27_O2 + W_A27_O3 +
W_A27_O4 + W_A28_O1 + W_A28_O2 + W_A28_O3 + W_A28_O4 + W_A29_O1 +
W_A29_O2 + W_A29_O3 + W_A29_O4 + W_A30_O1 + W_A30_O2 + W_A30_O3 +
W_A30_O4 + W_A31_O1 + W_A31_O2 + W_A31_O3 + W_A31_O4 + W_A32_O1 +
W_A32_O2 + W_A32_O3 + W_A32_O4 + W_A33_O1 + W_A33_O2 + W_A33_O3 +
W_A33_O4 + W_A34_O1 + W_A34_O2 + W_A34_O3 + W_A34_O4 + W_A35_O1 +
W_A35_O2 + W_A35_O3 + W_A35_O4 + W_A36_O1 + W_A36_O2 + W_A36_O3 +
W_A36_O4 + W_A37_O1 + W_A37_O2 + W_A37_O3 + W_A37_O4 + W_A38_O1 +
W_A38_O2 + W_A38_O3 + W_A38_O4 + W_A39_O1 + W_A39_O2 + W_A39_O3 +
W_A39_O4 = 1 ;
[_45] 7.5057 * W_A2_O1 + 0.692 * W_A2_O2 + 0.00089 * W_A2_O3 + 0.05207 *
W_A2_O4 >= 1 ;
[_46] 7.5057 * W_A2_O1 + 0.692 * W_A2_O2 + 0.00089 * W_A2_O3 + 0.05207 *
W_A2_O4 >= 1 ;
[_47] 7.5057 * W_A2_O1 + 0.692 * W_A2_O2 + 0.00089 * W_A2_O3 + 0.05207 *
W_A2_O4 >= 1 ;

我也说一句

发表

使用道具举报

yangwj78

2 主题	4 听众	87 积分

升级 86.32%

该用户从未签到

4^#

发表于 2009-7-1 14:44 |只看该作者 |招呼Ta 关注Ta

本帖最后由 yangwj78 于 2009-7-1 14:46 编辑

我改写后的模型和扩展模型如下，错误之处在于目标函数不对，如有时间的话，请指正：
model:
sets:
dmu/a1..a39/;
out/o1..o4/;
do(dmu,out):f,w;
endsets
data:
f=41.374 0.76 0.01214 10.7516
7.5057 0.692 0.00089 0.05207
64.977 2.85 0.01989 54.8142
160.75 1.252 0.02476 53.4763
39.785 5.563 0.02646 5.57088
408.5  1.298 0.09064 17.1816
391.58 19.82 0.01033 0.20892
188.83 29.45 0.00955 0.76588
250.24 4.361 0.01199 0.05094
5.3832 0.058 0.00084 0.30294
115.85 5.526 0.01165 0.0156
8.173  0.345 0.00093 0.04326
92.117 9.302 0.00211 0.20876
83.279 5.697 0.01185 0.03052
5.7716 0.076 0.00156 4.93978
2018.5 39.17 0.04799 0.05951
6.2366 0.179 0.00086 0.00539
2.8194 0.132 0.00042 5.47648
42.478 10.07 0.03094 5.20861
168.74 65.45 0.03044 0.09468
164.89 4.455 0.00601 0.25325
322.08 9.937 0.03895 0.24721
18.543 1.284 0.01 0.08281
3.064  0.26 0.00113 2.90088
21.045 0.93 0.08005 0.73361
43.218 4.642 0.02656 1.09635
35.178 10.17 0.0581 0.0869
11.856 0.793 0.00219 0.16704
7.3769 0.264 0.00224 0.14108
18.473 4.966 0.00324 0.33516
17.068 2.088 0.00218 0.01928
5.4787 0.196 0.00081 0.02805
3.0761 0.435 0.00047 0.01578
19.818 0.757 0.00199 0.01822
5.1734 0.343 0.00085 0.02332
109.35 11.87 0.01363 0.8445
32.229 0.812 0.21538 4.7541
59.355 33.63 0.03417 1.57874
133.09 3.28 0.0013 0.02362;
enddata
@for(dmu(m)|m #le# 39:
max=1/(@sum(out(i):f(m,i)*w(m,i)));
@sum(out(i):w(m,i))=1;
@for(dmu(k):
@sum(out(j):f(k,j)*w(m,j))>=1;
  )
);
end

我也说一句

发表

使用道具举报

yangwj78

2 主题	4 听众	87 积分

升级 86.32%

该用户从未签到

5^#

发表于 2009-7-1 14:45 |只看该作者 |招呼Ta 关注Ta

部分扩展模型：

MODEL:
[_1559] MAX= 1 / ( ( 133.09 * W_A39_O1 + 3.28 * W_A39_O2 +
0.0013 * W_A39_O3 + 0.02362 * W_A39_O4 ) ) ;
[_2] W_A1_O1 + W_A1_O2 + W_A1_O3 + W_A1_O4 = 1 ;
[_3] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
W_A1_O4 >= 1 ;
[_4] 7.5057 * W_A1_O1 + 0.692 * W_A1_O2 + 0.00089 * W_A1_O3 + 0.05207 *
W_A1_O4 >= 1 ;
[_5] 64.977 * W_A1_O1 + 2.85 * W_A1_O2 + 0.01989 * W_A1_O3 + 54.8142 *
W_A1_O4 >= 1 ;
[_6] 160.75 * W_A1_O1 + 1.252 * W_A1_O2 + 0.02476 * W_A1_O3 + 53.4763 *
W_A1_O4 >= 1 ;
[_7] 39.785 * W_A1_O1 + 5.563 * W_A1_O2 + 0.02646 * W_A1_O3 + 5.57088 *
W_A1_O4 >= 1 ;
[_8] 408.5 * W_A1_O1 + 1.298 * W_A1_O2 + 0.09064 * W_A1_O3 + 17.1816 *
W_A1_O4 >= 1 ;
[_9] 391.58 * W_A1_O1 + 19.82 * W_A1_O2 + 0.01033 * W_A1_O3 + 0.20892 *
W_A1_O4 >= 1 ;
[_10] 188.83 * W_A1_O1 + 29.45 * W_A1_O2 + 0.00955 * W_A1_O3 + 0.76588 *
W_A1_O4 >= 1 ;
[_11] 250.24 * W_A1_O1 + 4.361 * W_A1_O2 + 0.01199 * W_A1_O3 + 0.05094 *
W_A1_O4 >= 1 ;
[_12] 5.3832 * W_A1_O1 + 0.058 * W_A1_O2 + 0.00084 * W_A1_O3 + 0.30294 *
W_A1_O4 >= 1 ;
[_13] 115.85 * W_A1_O1 + 5.526 * W_A1_O2 + 0.01165 * W_A1_O3 + 0.0156 *
W_A1_O4 >= 1 ;
[_14] 8.173 * W_A1_O1 + 0.345 * W_A1_O2 + 0.0009300000000000001 *
W_A1_O3 + 0.04326 * W_A1_O4 >= 1 ;
[_15] 92.117 * W_A1_O1 + 9.302 * W_A1_O2 + 0.00211 * W_A1_O3 + 0.20876 *
W_A1_O4 >= 1 ;
[_16] 83.279 * W_A1_O1 + 5.697 * W_A1_O2 + 0.01185 * W_A1_O3 + 0.03052 *
W_A1_O4 >= 1 ;
[_17] 5.7716 * W_A1_O1 + 0.076 * W_A1_O2 + 0.00156 * W_A1_O3 + 4.93978 *
W_A1_O4 >= 1 ;
[_18] 2018.5 * W_A1_O1 + 39.17 * W_A1_O2 + 0.04799 * W_A1_O3 + 0.05951 *
W_A1_O4 >= 1 ;
[_19] 6.2366 * W_A1_O1 + 0.179 * W_A1_O2 + 0.00086 * W_A1_O3 + 0.00539 *
W_A1_O4 >= 1 ;
[_20] 2.8194 * W_A1_O1 + 0.132 * W_A1_O2 + 0.00042 * W_A1_O3 + 5.47648 *
W_A1_O4 >= 1 ;
[_21] 42.478 * W_A1_O1 + 10.07 * W_A1_O2 + 0.03094 * W_A1_O3 + 5.20861 *
W_A1_O4 >= 1 ;
[_22] 168.74 * W_A1_O1 + 65.45 * W_A1_O2 + 0.03044 * W_A1_O3 + 0.09468 *
W_A1_O4 >= 1 ;
[_23] 164.89 * W_A1_O1 + 4.455 * W_A1_O2 + 0.00601 * W_A1_O3 + 0.25325 *
W_A1_O4 >= 1 ;
[_24] 322.08 * W_A1_O1 + 9.936999999999999 * W_A1_O2 + 0.03895 * W_A1_O3
+ 0.24721 * W_A1_O4 >= 1 ;
[_25] 18.543 * W_A1_O1 + 1.284 * W_A1_O2 + 0.01 * W_A1_O3 + 0.08281 *
W_A1_O4 >= 1 ;
[_26] 3.064 * W_A1_O1 + 0.26 * W_A1_O2 + 0.00113 * W_A1_O3 + 2.90088 *
W_A1_O4 >= 1 ;
[_27] 21.045 * W_A1_O1 + 0.9300000000000001 * W_A1_O2 + 0.08005 *
W_A1_O3 + 0.73361 * W_A1_O4 >= 1 ;
[_28] 43.218 * W_A1_O1 + 4.642 * W_A1_O2 + 0.02656 * W_A1_O3 + 1.09635 *
W_A1_O4 >= 1 ;
[_29] 35.178 * W_A1_O1 + 10.17 * W_A1_O2 + 0.0581 * W_A1_O3 +
0.08690000000000001 * W_A1_O4 >= 1 ;
[_30] 11.856 * W_A1_O1 + 0.793 * W_A1_O2 + 0.00219 * W_A1_O3 + 0.16704 *
W_A1_O4 >= 1 ;
[_31] 7.3769 * W_A1_O1 + 0.264 * W_A1_O2 + 0.00224 * W_A1_O3 + 0.14108 *
W_A1_O4 >= 1 ;
[_32] 18.473 * W_A1_O1 + 4.966 * W_A1_O2 + 0.00324 * W_A1_O3 + 0.33516 *
W_A1_O4 >= 1 ;
[_33] 17.068 * W_A1_O1 + 2.088 * W_A1_O2 + 0.00218 * W_A1_O3 + 0.01928 *
W_A1_O4 >= 1 ;
[_34] 5.4787 * W_A1_O1 + 0.196 * W_A1_O2 + 0.00081 * W_A1_O3 + 0.02805 *
W_A1_O4 >= 1 ;
[_35] 3.0761 * W_A1_O1 + 0.435 * W_A1_O2 + 0.00047 * W_A1_O3 + 0.01578 *
W_A1_O4 >= 1 ;
[_36] 19.818 * W_A1_O1 + 0.757 * W_A1_O2 + 0.00199 * W_A1_O3 + 0.01822 *
W_A1_O4 >= 1 ;
[_37] 5.1734 * W_A1_O1 + 0.343 * W_A1_O2 + 0.00085 * W_A1_O3 + 0.02332 *
W_A1_O4 >= 1 ;
[_38] 109.35 * W_A1_O1 + 11.87 * W_A1_O2 + 0.01363 * W_A1_O3 + 0.8445 *
W_A1_O4 >= 1 ;
[_39] 32.229 * W_A1_O1 + 0.8120000000000001 * W_A1_O2 + 0.21538 *
W_A1_O3 + 4.7541 * W_A1_O4 >= 1 ;
[_40] 59.355 * W_A1_O1 + 33.63 * W_A1_O2 + 0.03417 * W_A1_O3 + 1.57874 *
W_A1_O4 >= 1 ;
[_41] 133.09 * W_A1_O1 + 3.28 * W_A1_O2 + 0.0013 * W_A1_O3 + 0.02362 *
W_A1_O4 >= 1 ;
[_42] 1 / ( ( 7.5057 * W_A2_O1 + 0.692 * W_A2_O2 + 0.00089 *
W_A2_O3 + 0.05207 * W_A2_O4 ) ) ;
[_43] W_A2_O1 + W_A2_O2 + W_A2_O3 + W_A2_O4 = 1 ;
[_44] 41.374 * W_A2_O1 + 0.76 * W_A2_O2 + 0.01214 * W_A2_O3 + 10.7516 *
W_A2_O4 >= 1 ;
[_45] 7.5057 * W_A2_O1 + 0.692 * W_A2_O2 + 0.00089 * W_A2_O3 + 0.05207 *
W_A2_O4 >= 1 ;