fmcmap/tp085.map


If Cut/Copy and Paste fails, then click here for download.


fmcTitle("tp085"):
# Source version 1
# Simple APM canonicalizer version 1.3
# FMCMAP back end version 1.1
asin := proc(x::algebraic) local y: y := arcsin(x): RETURN(y): end:
atan := proc(x::algebraic) local y: y := arctan(x): RETURN(y): end:
myownabs := proc(x::algebraic) local y: y := fmc_m_abs(x): RETURN(y): end:
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

Stephan K.H. Seidl