# Model tp085 # Source version 1 # Simple APM canonicalizer version 1.3 # AMPL backend version 1.1 param coefy5a := .004324; param coefy5b := .00423; param coefy5 := coefy5a; param a_2 := 17.505; param a_3 := 11.275; param a_4 := 214.228; param a_5 := 7.458; param a_6 := .961; param a_7 := 1.612; param a_8 := .146; param a_9 := 107.99; param a_10 := 922.693; param a_11 := 926.832; param a_12 := 18.766; param a_13 := 1072.163; param a_14 := 8961.448; param a_15 := .063; param a_16 := 71084.33; param a_17 := 2802713; param b_2 := 1053.6667; param b_3 := 35.03; param b_4 := 665.585; param b_5 := 584.463; param b_6 := 265.916; param b_7 := 7.046; param b_8 := .222; param b_9 := 273.366; param b_10 := 1286.105; param b_11 := 1444.046; param b_12 := 537.141; param b_13 := 3247.039; param b_14 := 26844.086; param b_15 := .386; param b_16 := 140000; param b_17 := 12146108; param c_10 := 123/7523; var x_1 >= 704.4148, <= 906.3855; let x_1 := 900; var x_2 >= 68.6, <= 288.88; let x_2 := 80; var x_3 >= 0, <= 134.75; let x_3 := 115; var x_4 >= 193, <= 287.0966; let x_4 := 267; var x_5 >= 25, <= 84.1988; let x_5 := 27; var y_1 = x_2 + x_3 + 41.6; var c_1 = .024*x_4 - 4.62; var y_2 = 12.5/c_1 + 12; var c_2 = .0003535*x_1^2 + .5311*x_1 + .08705*y_2*x_1; var c_3 = .052*x_1 + 78 + .002377*y_2*x_1; var y_3 = c_2/c_3; var y_4 = 19*y_3; var c_4 = .04782*(x_1 - y_3) + .1956*(x_1 - y_3)^2/x_2 + .6376*y_4 + 1.594*y_3; var c_5 = 100*x_2; var c_6 = x_1 - y_3 - y_4; var c_7 = .95 - c_4/c_5; var y_5 = c_6*c_7; var y_6 = x_1 - y_5 - y_4 - y_3; var c_8 = (y_5 + y_4)*.995; var y_7 = c_8/y_1; var y_8 = c_8/3798; var c_9 = y_7 - .0663*y_7/y_8 - .3153; var y_9 = 96.82/c_9 + .321*y_1; var y_10 = 1.29*y_5 + 1.258*y_4 + 2.29*y_3 + 1.71*y_6; var y_11 = 1.71*x_1 - .452*y_4 + .58*y_3; var c_11 = (1.75*y_2)*(.995*x_1); var c_12 = .995*y_10 + 1998; var y_12 = c_10*x_1 + c_11/c_12; var y_13 = c_12 - 1.75*y_2; var y_14 = 3623 + 64.4*x_2 + 58.4*x_3 + 146312/(y_9 + x_5); var c_13 = .995*y_10 + 60.8*x_2 + 48*x_4 - .1121*y_14 - 5095; var y_15 = y_13/c_13; var y_16 = 148000 - 331000*y_15 + 40*y_13 - 61*y_15*y_13; var c_14 = 2324*y_10 - 28740000*y_2; var y_17 = 14130000 - 1328*y_10 - 531*y_11 + c_14/c_12; var c_15 = y_13/y_15 - y_13/.52; var c_16 = 1.104 - .72*y_15; var c_17 = y_9 + x_5; var myminfun = (-5.843e-7)*y_17 + 1.17e-4*y_14 + 2.358e-5*y_13 + 1.502e-6*y_16 + .0321*y_12 + coefy5*y_5 + 1.0e-4*c_15/c_16 + 37.48*y_2/c_12 + .1365; var myabsdevnod0_0 = abs(myminfun + 1.905155258534784); var myreldevnod0_0 = abs(myminfun/1.905155258534784 + 1); var myabsdevnod0_1 = abs(x_1 - 705.1745370700908); var myreldevnod0_1 = abs(x_1/705.1745370700908 - 1); var myabsdevnod0_2 = abs(x_2 - 68.6); var myreldevnod0_2 = abs(x_2/68.6 - 1); var myabsdevnod0_3 = abs(x_3 - 102.9); var myreldevnod0_3 = abs(x_3/102.9 - 1); var myabsdevnod0_4 = abs(x_4 - 282.3249315936603); var myreldevnod0_4 = abs(x_4/282.3249315936603 - 1); var myabsdevnod0_5 = abs(x_5 - 37.58411642580555); var myreldevnod0_5 = abs(x_5/37.58411642580555 - 1); var myabsdevnod1_0 = (1/2)*(abs(myabsdevnod0_0 - myabsdevnod0_1) + (myabsdevnod0_0 + myabsdevnod0_1)); var myreldevnod1_0 = (1/2)*(abs(myreldevnod0_0 - myreldevnod0_1) + (myreldevnod0_0 + myreldevnod0_1)); var myabsdevnod1_2 = (1/2)*(abs(myabsdevnod0_2 - myabsdevnod0_3) + (myabsdevnod0_2 + myabsdevnod0_3)); var myreldevnod1_2 = (1/2)*(abs(myreldevnod0_2 - myreldevnod0_3) + (myreldevnod0_2 + myreldevnod0_3)); var myabsdevnod1_4 = (1/2)*(abs(myabsdevnod0_4 - myabsdevnod0_5) + (myabsdevnod0_4 + myabsdevnod0_5)); var myreldevnod1_4 = (1/2)*(abs(myreldevnod0_4 - myreldevnod0_5) + (myreldevnod0_4 + myreldevnod0_5)); var myabsdevnod2_0 = (1/2)*(abs(myabsdevnod1_0 - myabsdevnod1_2) + (myabsdevnod1_0 + myabsdevnod1_2)); var myreldevnod2_0 = (1/2)*(abs(myreldevnod1_0 - myreldevnod1_2) + (myreldevnod1_0 + myreldevnod1_2)); var myabsdevnod2_4 = myabsdevnod1_4; var myreldevnod2_4 = myreldevnod1_4; var myabsdevnod3_0 = (1/2)*(abs(myabsdevnod2_0 - myabsdevnod2_4) + (myabsdevnod2_0 + myabsdevnod2_4)); var myreldevnod3_0 = (1/2)*(abs(myreldevnod2_0 - myreldevnod2_4) + (myreldevnod2_0 + myreldevnod2_4)); var zmyabsdevmax = myabsdevnod3_0; var zmyreldevmax = myreldevnod3_0; minimize obj: myminfun; subject to ci_1: 1.5*x_2 - x_3 >= 0; subject to ci_2: y_1 - 213.1 >= 0; subject to ci_3: 405.23 - y_1 >= 0; subject to ci_4: y_2 - a_2 >= 0; subject to ci_5: y_3 - a_3 >= 0; subject to ci_6: y_4 - a_4 >= 0; subject to ci_7: y_5 - a_5 >= 0; subject to ci_8: y_6 - a_6 >= 0; subject to ci_9: y_7 - a_7 >= 0; subject to ci_10: y_8 - a_8 >= 0; subject to ci_11: y_9 - a_9 >= 0; subject to ci_12: y_10 - a_10 >= 0; subject to ci_13: y_11 - a_11 >= 0; subject to ci_14: y_12 - a_12 >= 0; subject to ci_15: y_13 - a_13 >= 0; subject to ci_16: y_14 - a_14 >= 0; subject to ci_17: y_15 - a_15 >= 0; subject to ci_18: y_16 - a_16 >= 0; subject to ci_19: y_17 - a_17 >= 0; subject to ci_20: b_2 - y_2 >= 0; subject to ci_21: b_3 - y_3 >= 0; subject to ci_22: b_4 - y_4 >= 0; subject to ci_23: b_5 - y_5 >= 0; subject to ci_24: b_6 - y_6 >= 0; subject to ci_25: b_7 - y_7 >= 0; subject to ci_26: b_8 - y_8 >= 0; subject to ci_27: b_9 - y_9 >= 0; subject to ci_28: b_10 - y_10 >= 0; subject to ci_29: b_11 - y_11 >= 0; subject to ci_30: b_12 - y_12 >= 0; subject to ci_31: b_13 - y_13 >= 0; subject to ci_32: b_14 - y_14 >= 0; subject to ci_33: b_15 - y_15 >= 0; subject to ci_34: b_16 - y_16 >= 0; subject to ci_35: b_17 - y_17 >= 0; subject to ci_36: y_4 - (28/72)*y_5 >= 0; subject to ci_37: 21 - 3496*y_2/c_12 >= 0; subject to ci_38: 62212/c_17 - 110.6 - y_1 >= 0; solve; display x_1; display x_2; display x_3; display x_4; display x_5; display myminfun; display zmyabsdevmax; display zmyreldevmax; # End Model