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