If Cut/Copy and Paste fails, then click here for download.
fmcTitle("tp070"):
# 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:
e1 := -1:
e2 := 1:
d1 := 7.658:
d2 := 7.685:
e := e1:
d := d1:
yobs_1 := 0.00189:
yobs_2 := 0.1038:
yobs_3 := 0.268:
yobs_4 := 0.506:
yobs_5 := 0.577:
yobs_6 := 0.604:
yobs_7 := 0.725:
yobs_8 := 0.898:
yobs_9 := 0.947:
yobs_10 := 0.845:
yobs_11 := 0.702:
yobs_12 := 0.528:
yobs_13 := 0.385:
yobs_14 := 0.257:
yobs_15 := 0.159:
yobs_16 := 0.0869:
yobs_17 := 0.0453:
yobs_18 := 0.01509:
yobs_19 := 0.00189:
c_1 := 0.1:
c_2 := 1:
c_3 := c_2 + 1:
c_4 := c_3 + 1:
c_5 := c_4 + 1:
c_6 := c_5 + 1:
c_7 := c_6 + 1:
c_8 := c_7 + 1:
c_9 := c_8 + 1:
c_10 := c_9 + 1:
c_11 := c_10 + 1:
c_12 := c_11 + 1:
c_13 := c_12 + 1:
c_14 := c_13 + 1:
c_15 := c_14 + 1:
c_16 := c_15 + 1:
c_17 := c_16 + 1:
c_18 := c_17 + 1:
c_19 := c_18 + 1:
s_0 := 0:
fmcInitialValue(x_1, 2):
fmcStrongLowerBound(x_1, 0.00001):
fmcStrongUpperBound(x_1, 100):
fmcInitialValue(x_2, 4):
fmcStrongLowerBound(x_2, 0.00001):
fmcStrongUpperBound(x_2, 100):
fmcInitialValue(x_3, 0.04):
fmcStrongLowerBound(x_3, 0.00001):
fmcStrongUpperBound(x_3, 1):
fmcInitialValue(x_4, 2):
fmcStrongLowerBound(x_4, 0.00001):
fmcStrongUpperBound(x_4, 100):
b := x_3 + (1 - x_3)*x_4:
ycal_1 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_1/d)^(x_2 - 1) * exp(x_2 -
b*c_1*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_1/7.658)^(x_1 - 1)
* exp(x_1 - b*c_1*x_1/(7.658*x_4)):
ycal_2 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_2/d)^(x_2 - 1) * exp(x_2 -
b*c_2*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_2/7.658)^(x_1 - 1)
* exp(x_1 - b*c_2*x_1/(7.658*x_4)):
ycal_3 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_3/d)^(x_2 - 1) * exp(x_2 -
b*c_3*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_3/7.658)^(x_1 - 1)
* exp(x_1 - b*c_3*x_1/(7.658*x_4)):
ycal_4 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_4/d)^(x_2 - 1) * exp(x_2 -
b*c_4*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_4/7.658)^(x_1 - 1)
* exp(x_1 - b*c_4*x_1/(7.658*x_4)):
ycal_5 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_5/d)^(x_2 - 1) * exp(x_2 -
b*c_5*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_5/7.658)^(x_1 - 1)
* exp(x_1 - b*c_5*x_1/(7.658*x_4)):
ycal_6 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_6/d)^(x_2 - 1) * exp(x_2 -
b*c_6*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_6/7.658)^(x_1 - 1)
* exp(x_1 - b*c_6*x_1/(7.658*x_4)):
ycal_7 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_7/d)^(x_2 - 1) * exp(x_2 -
b*c_7*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_7/7.658)^(x_1 - 1)
* exp(x_1 - b*c_7*x_1/(7.658*x_4)):
ycal_8 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_8/d)^(x_2 - 1) * exp(x_2 -
b*c_8*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_8/7.658)^(x_1 - 1)
* exp(x_1 - b*c_8*x_1/(7.658*x_4)):
ycal_9 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_9/d)^(x_2 - 1) * exp(x_2 -
b*c_9*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_9/7.658)^(x_1 - 1)
* exp(x_1 - b*c_9*x_1/(7.658*x_4)):
ycal_10 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_10/d)^(x_2 - 1) * exp(x_2 -
b*c_10*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_10/7.658)^(x_1 -
1) * exp(x_1 - b*c_10*x_1/(7.658*x_4)):
ycal_11 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_11/d)^(x_2 - 1) * exp(x_2 -
b*c_11*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_11/7.658)^(x_1 -
1) * exp(x_1 - b*c_11*x_1/(7.658*x_4)):
ycal_12 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_12/d)^(x_2 - 1) * exp(x_2 -
b*c_12*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_12/7.658)^(x_1 -
1) * exp(x_1 - b*c_12*x_1/(7.658*x_4)):
ycal_13 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_13/d)^(x_2 - 1) * exp(x_2 -
b*c_13*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_13/7.658)^(x_1 -
1) * exp(x_1 - b*c_13*x_1/(7.658*x_4)):
ycal_14 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_14/d)^(x_2 - 1) * exp(x_2 -
b*c_14*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_14/7.658)^(x_1 -
1) * exp(x_1 - b*c_14*x_1/(7.658*x_4)):
ycal_15 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_15/d)^(x_2 - 1) * exp(x_2 -
b*c_15*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_15/7.658)^(x_1 -
1) * exp(x_1 - b*c_15*x_1/(7.658*x_4)):
ycal_16 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_16/d)^(x_2 - 1) * exp(x_2 -
b*c_16*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_16/7.658)^(x_1 -
1) * exp(x_1 - b*c_16*x_1/(7.658*x_4)):
ycal_17 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_17/d)^(x_2 - 1) * exp(x_2 -
b*c_17*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_17/7.658)^(x_1 -
1) * exp(x_1 - b*c_17*x_1/(7.658*x_4)):
ycal_18 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_18/d)^(x_2 - 1) * exp(x_2 -
b*c_18*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_18/7.658)^(x_1 -
1) * exp(x_1 - b*c_18*x_1/(7.658*x_4)):
ycal_19 := (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832)^(1/2) * (c_19/d)^(x_2 - 1) * exp(x_2 -
b*c_19*x_2/7.658) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832)^(1/2) * (c_19/7.658)^(x_1 -
1) * exp(x_1 - b*c_19*x_1/(7.658*x_4)):
s_1 := s_0 + (yobs_1 - ycal_1)^2:
s_2 := s_1 + (yobs_2 - ycal_2)^2:
s_3 := s_2 + (yobs_3 - ycal_3)^2:
s_4 := s_3 + (yobs_4 - ycal_4)^2:
s_5 := s_4 + (yobs_5 - ycal_5)^2:
s_6 := s_5 + (yobs_6 - ycal_6)^2:
s_7 := s_6 + (yobs_7 - ycal_7)^2:
s_8 := s_7 + (yobs_8 - ycal_8)^2:
s_9 := s_8 + (yobs_9 - ycal_9)^2:
s_10 := s_9 + (yobs_10 - ycal_10)^2:
s_11 := s_10 + (yobs_11 - ycal_11)^2:
s_12 := s_11 + (yobs_12 - ycal_12)^2:
s_13 := s_12 + (yobs_13 - ycal_13)^2:
s_14 := s_13 + (yobs_14 - ycal_14)^2:
s_15 := s_14 + (yobs_15 - ycal_15)^2:
s_16 := s_15 + (yobs_16 - ycal_16)^2:
s_17 := s_16 + (yobs_17 - ycal_17)^2:
s_18 := s_17 + (yobs_18 - ycal_18)^2:
s_19 := s_18 + (yobs_19 - ycal_19)^2:
mf := s_19:
myminfun := mf:
myabsdevnod0_0 := myownabs(myminfun - 0.007498463574427645):
myreldevnod0_0 := myownabs(myminfun/0.007498463574427645 - 1):
myabsdevnod0_1 := myownabs(x_1 - 12.27697912756719):
myreldevnod0_1 := myownabs(x_1/12.27697912756719 - 1):
myabsdevnod0_2 := myownabs(x_2 - 4.631748162745852):
myreldevnod0_2 := myownabs(x_2/4.631748162745852 - 1):
myabsdevnod0_3 := myownabs(x_3 - 0.3128646302166193):
myreldevnod0_3 := myownabs(x_3/0.3128646302166193 - 1):
myabsdevnod0_4 := myownabs(x_4 - 2.029282825337289):
myreldevnod0_4 := myownabs(x_4/2.029282825337289 - 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 := myabsdevnod0_4:
myreldevnod1_4 := myreldevnod0_4:
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, b):
fmcEscortFloat(x_1, x_1):
fmcEscortFloat(x_2, x_2):
fmcEscortFloat(x_3, x_3):
fmcEscortFloat(x_4, x_4):
fmcEscortFloat('myminfun', myminfun):
fmcEscortFloat('zmyabsdevmax', zmyabsdevmax):
fmcEscortFloat('zmyreldevmax', zmyreldevmax):
fmcControlMinimum(1, 0):
# End Model