fmcmac/tp085.mac


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


fmcTitle("tp085")$
/* Source version 1 */
/* Simple APM canonicalizer version 1.3 */
/* FMCMAC backend version 1.2 */
myownabs(x) := fmc_m_abs(x)$
coefy5a : .004324b0$
coefy5b : .00423b0$
coefy5 : coefy5a$
a_2 : 17.505b0$
a_3 : 11.275b0$
a_4 : 214.228b0$
a_5 : 7.458b0$
a_6 : .961b0$
a_7 : 1.612b0$
a_8 : .146b0$
a_9 : 107.99b0$
a_10 : 922.693b0$
a_11 : 926.832b0$
a_12 : 18.766b0$
a_13 : 1072.163b0$
a_14 : 8961.448b0$
a_15 : .063b0$
a_16 : 71084.33b0$
a_17 : 2802713$
b_2 : 1053.6667b0$
b_3 : 35.03b0$
b_4 : 665.585b0$
b_5 : 584.463b0$
b_6 : 265.916b0$
b_7 : 7.046b0$
b_8 : .222b0$
b_9 : 273.366b0$
b_10 : 1286.105b0$
b_11 : 1444.046b0$
b_12 : 537.141b0$
b_13 : 3247.039b0$
b_14 : 26844.086b0$
b_15 : .386b0$
b_16 : 140000$
b_17 : 12146108$
c_10 : 123/7523$
fmcInitialValue(x_1, 900)$
fmcStrongLowerBound(x_1, 704.4148b0)$
fmcStrongUpperBound(x_1, 906.3855b0)$
fmcInitialValue(x_2, 80)$
fmcStrongLowerBound(x_2, 68.6b0)$
fmcStrongUpperBound(x_2, 288.88b0)$
fmcInitialValue(x_3, 115)$
fmcStrongLowerBound(x_3, 0)$
fmcStrongUpperBound(x_3, 134.75b0)$
fmcInitialValue(x_4, 267)$
fmcStrongLowerBound(x_4, 193)$
fmcStrongUpperBound(x_4, 287.0966b0)$
fmcInitialValue(x_5, 27)$
fmcStrongLowerBound(x_5, 25)$
fmcStrongUpperBound(x_5, 84.1988b0)$
y_1 : x_2 + x_3 + 41.6b0$
c_1 : .024b0*x_4 - 4.62b0$
y_2 : 12.5b0/c_1 + 12$
c_2 : .0003535b0*x_1^2 + .5311b0*x_1 + .08705b0*y_2*x_1$
c_3 : .052b0*x_1 + 78 + .002377b0*y_2*x_1$
y_3 : c_2/c_3$
y_4 : 19*y_3$
c_4 : .04782b0*(x_1 - y_3) + .1956b0*(x_1 - y_3)^2/x_2 + .6376b0*y_4 + 1.594b0*y_3$
c_5 : 100*x_2$
c_6 : x_1 - y_3 - y_4$
c_7 : .95b0 - 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)*.995b0$
y_7 : c_8/y_1$
y_8 : c_8/3798$
c_9 : y_7 - .0663b0*y_7/y_8 - .3153b0$
y_9 : 96.82b0/c_9 + .321b0*y_1$
y_10 : 1.29b0*y_5 + 1.258b0*y_4 + 2.29b0*y_3 + 1.71b0*y_6$
y_11 : 1.71b0*x_1 - .452b0*y_4 + .58b0*y_3$
c_11 : (1.75b0*y_2)*(.995b0*x_1)$
c_12 : .995b0*y_10 + 1998$
y_12 : c_10*x_1 + c_11/c_12$
y_13 : c_12 - 1.75b0*y_2$
y_14 : 3623 + 64.4b0*x_2 + 58.4b0*x_3 + 146312/(y_9 + x_5)$
c_13 : .995b0*y_10 + 60.8b0*x_2 + 48*x_4 - .1121b0*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/.52b0$
c_16 : 1.104b0 - .72b0*y_15$
c_17 : y_9 + x_5$
myminfun : (-5.843b-7)*y_17 + 1.17b-4*y_14 + 2.358b-5*y_13 + 1.502b-6*y_16 + .0321b0*y_12 + coefy5*y_5 +
1.0b-4*c_15/c_16 + 37.48b0*y_2/c_12 + .1365b0$
myabsdevnod0_0 : myownabs(myminfun + 1.905155258534784b0)$
myreldevnod0_0 : myownabs(myminfun/1.905155258534784b0 + 1)$
myabsdevnod0_1 : myownabs(x_1 - 705.1745370700908b0)$
myreldevnod0_1 : myownabs(x_1/705.1745370700908b0 - 1)$
myabsdevnod0_2 : myownabs(x_2 - 68.6b0)$
myreldevnod0_2 : myownabs(x_2/68.6b0 - 1)$
myabsdevnod0_3 : myownabs(x_3 - 102.9b0)$
myreldevnod0_3 : myownabs(x_3/102.9b0 - 1)$
myabsdevnod0_4 : myownabs(x_4 - 282.3249315936603b0)$
myreldevnod0_4 : myownabs(x_4/282.3249315936603b0 - 1)$
myabsdevnod0_5 : myownabs(x_5 - 37.58411642580555b0)$
myreldevnod0_5 : myownabs(x_5/37.58411642580555b0 - 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.5b0*x_2 - x_3)$
fmcInequality(ci_2, y_1 - 213.1b0)$
fmcInequality(ci_3, 405.23b0 - 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.6b0 - 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