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