If Cut/Copy and Paste fails, then click here for download.
fmcTitle("tp070")$
/* Source version 1 */
/* Simple APM canonicalizer version 1.3 */
/* FMCMAC backend version 1.2 */
myownabs(x) := fmc_m_abs(x)$
e1 : -1$
e2 : 1$
d1 : 7.658b0$
d2 : 7.685b0$
e : e1$
d : d1$
yobs_1 : 0.00189b0$
yobs_2 : 0.1038b0$
yobs_3 : 0.268b0$
yobs_4 : 0.506b0$
yobs_5 : 0.577b0$
yobs_6 : 0.604b0$
yobs_7 : 0.725b0$
yobs_8 : 0.898b0$
yobs_9 : 0.947b0$
yobs_10 : 0.845b0$
yobs_11 : 0.702b0$
yobs_12 : 0.528b0$
yobs_13 : 0.385b0$
yobs_14 : 0.257b0$
yobs_15 : 0.159b0$
yobs_16 : 0.0869b0$
yobs_17 : 0.0453b0$
yobs_18 : 0.01509b0$
yobs_19 : 0.00189b0$
c_1 : 0.1b0$
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.00001b0)$
fmcStrongUpperBound(x_1, 100)$
fmcInitialValue(x_2, 4)$
fmcStrongLowerBound(x_2, 0.00001b0)$
fmcStrongUpperBound(x_2, 100)$
fmcInitialValue(x_3, 0.04b0)$
fmcStrongLowerBound(x_3, 0.00001b0)$
fmcStrongUpperBound(x_3, 1)$
fmcInitialValue(x_4, 2)$
fmcStrongLowerBound(x_4, 0.00001b0)$
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.2832b0)^(1/2) * (c_1/d)^(x_2 - 1) * exp(x_2 -
b*c_1*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_1/7.658b0)^(x_1 - 1)
* exp(x_1 - b*c_1*x_1/(7.658b0*x_4))$
ycal_2 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_2/d)^(x_2 - 1) * exp(x_2 -
b*c_2*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_2/7.658b0)^(x_1 - 1)
* exp(x_1 - b*c_2*x_1/(7.658b0*x_4))$
ycal_3 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_3/d)^(x_2 - 1) * exp(x_2 -
b*c_3*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_3/7.658b0)^(x_1 - 1)
* exp(x_1 - b*c_3*x_1/(7.658b0*x_4))$
ycal_4 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_4/d)^(x_2 - 1) * exp(x_2 -
b*c_4*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_4/7.658b0)^(x_1 - 1)
* exp(x_1 - b*c_4*x_1/(7.658b0*x_4))$
ycal_5 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_5/d)^(x_2 - 1) * exp(x_2 -
b*c_5*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_5/7.658b0)^(x_1 - 1)
* exp(x_1 - b*c_5*x_1/(7.658b0*x_4))$
ycal_6 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_6/d)^(x_2 - 1) * exp(x_2 -
b*c_6*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_6/7.658b0)^(x_1 - 1)
* exp(x_1 - b*c_6*x_1/(7.658b0*x_4))$
ycal_7 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_7/d)^(x_2 - 1) * exp(x_2 -
b*c_7*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_7/7.658b0)^(x_1 - 1)
* exp(x_1 - b*c_7*x_1/(7.658b0*x_4))$
ycal_8 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_8/d)^(x_2 - 1) * exp(x_2 -
b*c_8*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_8/7.658b0)^(x_1 - 1)
* exp(x_1 - b*c_8*x_1/(7.658b0*x_4))$
ycal_9 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_9/d)^(x_2 - 1) * exp(x_2 -
b*c_9*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_9/7.658b0)^(x_1 - 1)
* exp(x_1 - b*c_9*x_1/(7.658b0*x_4))$
ycal_10 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_10/d)^(x_2 - 1) * exp(x_2 -
b*c_10*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_10/7.658b0)^(x_1 -
1) * exp(x_1 - b*c_10*x_1/(7.658b0*x_4))$
ycal_11 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_11/d)^(x_2 - 1) * exp(x_2 -
b*c_11*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_11/7.658b0)^(x_1 -
1) * exp(x_1 - b*c_11*x_1/(7.658b0*x_4))$
ycal_12 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_12/d)^(x_2 - 1) * exp(x_2 -
b*c_12*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_12/7.658b0)^(x_1 -
1) * exp(x_1 - b*c_12*x_1/(7.658b0*x_4))$
ycal_13 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_13/d)^(x_2 - 1) * exp(x_2 -
b*c_13*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_13/7.658b0)^(x_1 -
1) * exp(x_1 - b*c_13*x_1/(7.658b0*x_4))$
ycal_14 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_14/d)^(x_2 - 1) * exp(x_2 -
b*c_14*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_14/7.658b0)^(x_1 -
1) * exp(x_1 - b*c_14*x_1/(7.658b0*x_4))$
ycal_15 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_15/d)^(x_2 - 1) * exp(x_2 -
b*c_15*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_15/7.658b0)^(x_1 -
1) * exp(x_1 - b*c_15*x_1/(7.658b0*x_4))$
ycal_16 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_16/d)^(x_2 - 1) * exp(x_2 -
b*c_16*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_16/7.658b0)^(x_1 -
1) * exp(x_1 - b*c_16*x_1/(7.658b0*x_4))$
ycal_17 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_17/d)^(x_2 - 1) * exp(x_2 -
b*c_17*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_17/7.658b0)^(x_1 -
1) * exp(x_1 - b*c_17*x_1/(7.658b0*x_4))$
ycal_18 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_18/d)^(x_2 - 1) * exp(x_2 -
b*c_18*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_18/7.658b0)^(x_1 -
1) * exp(x_1 - b*c_18*x_1/(7.658b0*x_4))$
ycal_19 : (1 + 1/(12*x_2))^e * x_3*b^x_2 * (x_2/6.2832b0)^(1/2) * (c_19/d)^(x_2 - 1) * exp(x_2 -
b*c_19*x_2/7.658b0) + (1 + 1/(12*x_1))^e * (1 - x_3) * (b/x_4)^x_1 * (x_1/6.2832b0)^(1/2) * (c_19/7.658b0)^(x_1 -
1) * exp(x_1 - b*c_19*x_1/(7.658b0*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.007498463574427645b0)$
myreldevnod0_0 : myownabs(myminfun/0.007498463574427645b0 - 1)$
myabsdevnod0_1 : myownabs(x_1 - 12.27697912756719b0)$
myreldevnod0_1 : myownabs(x_1/12.27697912756719b0 - 1)$
myabsdevnod0_2 : myownabs(x_2 - 4.631748162745852b0)$
myreldevnod0_2 : myownabs(x_2/4.631748162745852b0 - 1)$
myabsdevnod0_3 : myownabs(x_3 - 0.3128646302166193b0)$
myreldevnod0_3 : myownabs(x_3/0.3128646302166193b0 - 1)$
myabsdevnod0_4 : myownabs(x_4 - 2.029282825337289b0)$
myreldevnod0_4 : myownabs(x_4/2.029282825337289b0 - 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 */