If Cut/Copy and Paste fails, then click here for download.
fmcTitle("tp086"):
# 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:
a24a := 0.4:
a24b := 4:
a24 := a24a:
a_1_1 := -16:
a_1_2 := 2:
a_1_3 := 0:
a_1_4 := 1:
a_1_5 := 0:
a_2_1 := 0:
a_2_2 := -2:
a_2_3 := 0:
a_2_4 := a24:
a_2_5 := 2:
a_3_1 := -3.5:
a_3_2 := 0:
a_3_3 := 2:
a_3_4 := 0:
a_3_5 := 0:
a_4_1 := 0:
a_4_2 := -2:
a_4_3 := 0:
a_4_4 := -4:
a_4_5 := -1:
a_5_1 := 0:
a_5_2 := -9:
a_5_3 := -2:
a_5_4 := 1:
a_5_5 := -2.8:
a_6_1 := 2:
a_6_2 := 0:
a_6_3 := -4:
a_6_4 := 0:
a_6_5 := 0:
a_7_1 := -1:
a_7_2 := -1:
a_7_3 := -1:
a_7_4 := -1:
a_7_5 := -1:
a_8_1 := -1:
a_8_2 := -2:
a_8_3 := -3:
a_8_4 := -2:
a_8_5 := -1:
a_9_1 := 1:
a_9_2 := 2:
a_9_3 := 3:
a_9_4 := 4:
a_9_5 := 5:
a_10_1 := 1:
a_10_2 := 1:
a_10_3 := 1:
a_10_4 := 1:
a_10_5 := 1:
b_1 := -40:
b_2 := -2:
b_3 := -0.25:
b_4 := -4:
b_5 := -4:
b_6 := -1:
b_7 := -40:
b_8 := -60:
b_9 := 5:
b_10 := 1:
c_1_1 := 30:
c_1_2 := -20:
c_1_3 := -10:
c_1_4 := 32:
c_1_5 := -10:
c_2_1 := -20:
c_2_2 := 39:
c_2_3 := -6:
c_2_4 := -31:
c_2_5 := 32:
c_3_1 := -10:
c_3_2 := -6:
c_3_3 := 10:
c_3_4 := -6:
c_3_5 := -10:
c_4_1 := 32:
c_4_2 := -31:
c_4_3 := -6:
c_4_4 := 39:
c_4_5 := -20:
c_5_1 := -10:
c_5_2 := 32:
c_5_3 := -10:
c_5_4 := -20:
c_5_5 := 30:
d_1 := 4:
d_2 := 8:
d_3 := 10:
d_4 := 6:
d_5 := 2:
e_1 := -15:
e_2 := -27:
e_3 := -36:
e_4 := -18:
e_5 := -12:
sumax_1_0 := 0:
sumax_2_0 := 0:
sumax_3_0 := 0:
sumax_4_0 := 0:
sumax_5_0 := 0:
sumax_6_0 := 0:
sumax_7_0 := 0:
sumax_8_0 := 0:
sumax_9_0 := 0:
sumax_10_0 := 0:
sumex_0 := 0:
sumcxxi_1_0 := 0:
sumcxxi_2_0 := 0:
sumcxxi_3_0 := 0:
sumcxxi_4_0 := 0:
sumcxxi_5_0 := 0:
sumcxx_0 := 0:
sumdx_0 := 0:
fmcInitialValue(x_1, 0):
fmcStrongLowerBound(x_1, 0):
fmcInitialValue(x_2, 0):
fmcStrongLowerBound(x_2, 0):
fmcInitialValue(x_3, 0):
fmcStrongLowerBound(x_3, 0):
fmcInitialValue(x_4, 0):
fmcStrongLowerBound(x_4, 0):
fmcInitialValue(x_5, 1):
fmcStrongLowerBound(x_5, 0):
sumax_1_1 := sumax_1_0 + a_1_1*x_1:
sumax_1_2 := sumax_1_1 + a_1_2*x_2:
sumax_1_3 := sumax_1_2 + a_1_3*x_3:
sumax_1_4 := sumax_1_3 + a_1_4*x_4:
sumax_1_5 := sumax_1_4 + a_1_5*x_5:
sumax_2_1 := sumax_2_0 + a_2_1*x_1:
sumax_2_2 := sumax_2_1 + a_2_2*x_2:
sumax_2_3 := sumax_2_2 + a_2_3*x_3:
sumax_2_4 := sumax_2_3 + a_2_4*x_4:
sumax_2_5 := sumax_2_4 + a_2_5*x_5:
sumax_3_1 := sumax_3_0 + a_3_1*x_1:
sumax_3_2 := sumax_3_1 + a_3_2*x_2:
sumax_3_3 := sumax_3_2 + a_3_3*x_3:
sumax_3_4 := sumax_3_3 + a_3_4*x_4:
sumax_3_5 := sumax_3_4 + a_3_5*x_5:
sumax_4_1 := sumax_4_0 + a_4_1*x_1:
sumax_4_2 := sumax_4_1 + a_4_2*x_2:
sumax_4_3 := sumax_4_2 + a_4_3*x_3:
sumax_4_4 := sumax_4_3 + a_4_4*x_4:
sumax_4_5 := sumax_4_4 + a_4_5*x_5:
sumax_5_1 := sumax_5_0 + a_5_1*x_1:
sumax_5_2 := sumax_5_1 + a_5_2*x_2:
sumax_5_3 := sumax_5_2 + a_5_3*x_3:
sumax_5_4 := sumax_5_3 + a_5_4*x_4:
sumax_5_5 := sumax_5_4 + a_5_5*x_5:
sumax_6_1 := sumax_6_0 + a_6_1*x_1:
sumax_6_2 := sumax_6_1 + a_6_2*x_2:
sumax_6_3 := sumax_6_2 + a_6_3*x_3:
sumax_6_4 := sumax_6_3 + a_6_4*x_4:
sumax_6_5 := sumax_6_4 + a_6_5*x_5:
sumax_7_1 := sumax_7_0 + a_7_1*x_1:
sumax_7_2 := sumax_7_1 + a_7_2*x_2:
sumax_7_3 := sumax_7_2 + a_7_3*x_3:
sumax_7_4 := sumax_7_3 + a_7_4*x_4:
sumax_7_5 := sumax_7_4 + a_7_5*x_5:
sumax_8_1 := sumax_8_0 + a_8_1*x_1:
sumax_8_2 := sumax_8_1 + a_8_2*x_2:
sumax_8_3 := sumax_8_2 + a_8_3*x_3:
sumax_8_4 := sumax_8_3 + a_8_4*x_4:
sumax_8_5 := sumax_8_4 + a_8_5*x_5:
sumax_9_1 := sumax_9_0 + a_9_1*x_1:
sumax_9_2 := sumax_9_1 + a_9_2*x_2:
sumax_9_3 := sumax_9_2 + a_9_3*x_3:
sumax_9_4 := sumax_9_3 + a_9_4*x_4:
sumax_9_5 := sumax_9_4 + a_9_5*x_5:
sumax_10_1 := sumax_10_0 + a_10_1*x_1:
sumax_10_2 := sumax_10_1 + a_10_2*x_2:
sumax_10_3 := sumax_10_2 + a_10_3*x_3:
sumax_10_4 := sumax_10_3 + a_10_4*x_4:
sumax_10_5 := sumax_10_4 + a_10_5*x_5:
ineq_1 := sumax_1_5 - b_1:
ineq_2 := sumax_2_5 - b_2:
ineq_3 := sumax_3_5 - b_3:
ineq_4 := sumax_4_5 - b_4:
ineq_5 := sumax_5_5 - b_5:
ineq_6 := sumax_6_5 - b_6:
ineq_7 := sumax_7_5 - b_7:
ineq_8 := sumax_8_5 - b_8:
ineq_9 := sumax_9_5 - b_9:
ineq_10 := sumax_10_5 - b_10:
sumex_1 := sumex_0 + e_1*x_1:
sumex_2 := sumex_1 + e_2*x_2:
sumex_3 := sumex_2 + e_3*x_3:
sumex_4 := sumex_3 + e_4*x_4:
sumex_5 := sumex_4 + e_5*x_5:
sumcxxi_1_1 := sumcxxi_1_0 + c_1_1*x_1*x_1:
sumcxxi_1_2 := sumcxxi_1_1 + c_1_2*x_1*x_2:
sumcxxi_1_3 := sumcxxi_1_2 + c_1_3*x_1*x_3:
sumcxxi_1_4 := sumcxxi_1_3 + c_1_4*x_1*x_4:
sumcxxi_1_5 := sumcxxi_1_4 + c_1_5*x_1*x_5:
sumcxxi_2_1 := sumcxxi_2_0 + c_2_1*x_2*x_1:
sumcxxi_2_2 := sumcxxi_2_1 + c_2_2*x_2*x_2:
sumcxxi_2_3 := sumcxxi_2_2 + c_2_3*x_2*x_3:
sumcxxi_2_4 := sumcxxi_2_3 + c_2_4*x_2*x_4:
sumcxxi_2_5 := sumcxxi_2_4 + c_2_5*x_2*x_5:
sumcxxi_3_1 := sumcxxi_3_0 + c_3_1*x_3*x_1:
sumcxxi_3_2 := sumcxxi_3_1 + c_3_2*x_3*x_2:
sumcxxi_3_3 := sumcxxi_3_2 + c_3_3*x_3*x_3:
sumcxxi_3_4 := sumcxxi_3_3 + c_3_4*x_3*x_4:
sumcxxi_3_5 := sumcxxi_3_4 + c_3_5*x_3*x_5:
sumcxxi_4_1 := sumcxxi_4_0 + c_4_1*x_4*x_1:
sumcxxi_4_2 := sumcxxi_4_1 + c_4_2*x_4*x_2:
sumcxxi_4_3 := sumcxxi_4_2 + c_4_3*x_4*x_3:
sumcxxi_4_4 := sumcxxi_4_3 + c_4_4*x_4*x_4:
sumcxxi_4_5 := sumcxxi_4_4 + c_4_5*x_4*x_5:
sumcxxi_5_1 := sumcxxi_5_0 + c_5_1*x_5*x_1:
sumcxxi_5_2 := sumcxxi_5_1 + c_5_2*x_5*x_2:
sumcxxi_5_3 := sumcxxi_5_2 + c_5_3*x_5*x_3:
sumcxxi_5_4 := sumcxxi_5_3 + c_5_4*x_5*x_4:
sumcxxi_5_5 := sumcxxi_5_4 + c_5_5*x_5*x_5:
sumcxx_1 := sumcxx_0 + sumcxxi_1_5:
sumcxx_2 := sumcxx_1 + sumcxxi_2_5:
sumcxx_3 := sumcxx_2 + sumcxxi_3_5:
sumcxx_4 := sumcxx_3 + sumcxxi_4_5:
sumcxx_5 := sumcxx_4 + sumcxxi_5_5:
sumdx_1 := sumdx_0 + d_1*x_1^3:
sumdx_2 := sumdx_1 + d_2*x_2^3:
sumdx_3 := sumdx_2 + d_3*x_3^3:
sumdx_4 := sumdx_3 + d_4*x_4^3:
sumdx_5 := sumdx_4 + d_5*x_5^3:
mf := sumex_5 + sumcxx_5 + sumdx_5:
myminfun := mf:
myabsdevnod0_0 := myownabs(myminfun + 32.34867896572271):
myreldevnod0_0 := myownabs(myminfun/32.34867896572271 + 1):
myabsdevnod0_1 := myownabs(x_1 - 0.3):
myreldevnod0_1 := myownabs(x_1/0.3 - 1):
myabsdevnod0_2 := myownabs(x_2 - 0.3334676065346071):
myreldevnod0_2 := myownabs(x_2/0.3334676065346071 - 1):
myabsdevnod0_3 := myownabs(x_3 - 0.4):
myreldevnod0_3 := myownabs(x_3/0.4 - 1):
myabsdevnod0_4 := myownabs(x_4 - 0.4283101047816988):
myreldevnod0_4 := myownabs(x_4/0.4283101047816988 - 1):
myabsdevnod0_5 := myownabs(x_5 - 0.2239648735607981):
myreldevnod0_5 := myownabs(x_5/0.2239648735607981 - 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, ineq_1):
fmcInequality(ci_2, ineq_2):
fmcInequality(ci_3, ineq_3):
fmcInequality(ci_4, ineq_4):
fmcInequality(ci_5, ineq_5):
fmcInequality(ci_6, ineq_6):
fmcInequality(ci_7, ineq_7):
fmcInequality(ci_8, ineq_8):
fmcInequality(ci_9, ineq_9):
fmcInequality(ci_10, ineq_10):
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