If Cut/Copy and Paste fails, then click here for download.
fmcTitle("tp087r2"): # Source version 7 # Simple APM canonicalizer version 1.3 # FMCMAP backend version 1.2 myownabs := proc(x::algebraic) return(fmc_m_abs(x)): end proc: bv1 := 1.48477: bv2 := 1.48577: bv := bv1: ubx5hsv1 := 1000: ubx5hsv2 := 10000: ubx5hsv := ubx5hsv1: a := 131.078: b := bv: c := 0.90798: arg := 1.47588: d := cos(arg): e := sin(arg): ivx5hs := 198.175: lbx5hs := -1000: lbx5r := -18: lbx5 := lbx5r: ubx5hs := ubx5hsv: ubx5r := -10.7: ubx5 := ubx5r: ivx5 := (1/2)*(ivx5hs + ubx5 - abs(ivx5hs - ubx5)): sign1_0 := 1: sign1_2 := -1: sign2_0 := 1: sign2_3 := -1: fmcInitialValue(x_1, 390): fmcStrongLowerBound(x_1, 0): fmcStrongUpperBound(x_1, 400): fmcInitialValue(x_2, 1000): fmcStrongLowerBound(x_2, 0): fmcStrongUpperBound(x_2, 1000): fmcInitialValue(x_3, 419.5): fmcStrongLowerBound(x_3, 340): fmcStrongUpperBound(x_3, 420): fmcInitialValue(x_4, 340.5): fmcStrongLowerBound(x_4, 340): fmcStrongUpperBound(x_4, 420): fmcInitialValue(x_5, ivx5): fmcStrongLowerBound(x_5, lbx5): fmcStrongUpperBound(x_5, ubx5): fmcInitialValue(x_6, 0.5): fmcStrongLowerBound(x_6, 0): fmcStrongUpperBound(x_6, 0.5236): fmcInitialValue(x_7, 0): fmcStrongLowerBound(x_7, 0): fmcStrongUpperBound(x_7, 300): fmcInitialValue(x_8, 90): fmcStrongLowerBound(x_8, 0): fmcStrongUpperBound(x_8, 300): fmcInitialValue(x_9, 0): fmcStrongLowerBound(x_9, 0): fmcStrongUpperBound(x_9, 900): fmcInitialValue(x_10, 900): fmcStrongLowerBound(x_10, 0): fmcStrongUpperBound(x_10, 900): fmcInitialValue(x_11, 0): fmcStrongLowerBound(x_11, 0): fmcStrongUpperBound(x_11, 800): fmcInitialValue(x_12, 800): fmcStrongLowerBound(x_12, 0): fmcStrongUpperBound(x_12, 800): fmcInitialValue(x_13, 1): fmcStrongLowerBound(x_13, -1): fmcStrongUpperBound(x_13, 1): fmcInitialValue(x_14, 1): fmcStrongLowerBound(x_14, -1): fmcStrongUpperBound(x_14, 1): fmcInitialValue(x_15, 1): fmcStrongLowerBound(x_15, -1): fmcStrongUpperBound(x_15, 1): cons_1 := 300 - x_1 - x_3*x_4*cos(b - x_6)/a + c*d*x_3^2/a: cons_2 := (-1)*x_2 - x_3*x_4*cos(b + x_6)/a + c*d*x_4^2/a: cons_3 := (-1)*x_5 - x_3*x_4*sin(b + x_6)/a + c*e*x_4^2/a: cons_4 := 200 - x_3*x_4*sin(b - x_6)/a + c*e*x_3^2/a: comp1_1 := x_1 - 300: comp2_1 := x_2 - 100: comp2_2 := x_2 - 200: comp1nnxa_1 := x_7: comp1nnxb_1 := x_8: comp2nnxa_1 := x_9: comp2nnxb_1 := x_10: comp2nnxa_2 := x_11: comp2nnxb_2 := x_12: cons_5 := comp1_1 + comp1nnxa_1 - comp1nnxb_1: cons_6 := comp2_1 + comp2nnxa_1 - comp2nnxb_1: cons_7 := comp2_2 + comp2nnxa_2 - comp2nnxb_2: comp1abs_1 := comp1nnxa_1 + comp1nnxb_1: comp2abs_1 := comp2nnxa_1 + comp2nnxb_1: comp2abs_2 := comp2nnxa_2 + comp2nnxb_2: sign1_1 := x_13: sign2_1 := x_14: sign2_2 := x_15: cons_8 := sign1_1*comp1abs_1 - comp1_1: cons_9 := sign2_1*comp2abs_1 - comp2_1: cons_10 := sign2_2*comp2abs_2 - comp2_2: step1_1 := (1 + sign1_0)*(1 - sign1_1)/4: step1_2 := (1 + sign1_1)*(1 - sign1_2)/4: step2_1 := (1 + sign2_0)*(1 - sign2_1)/4: step2_2 := (1 + sign2_1)*(1 - sign2_2)/4: step2_3 := (1 + sign2_2)*(1 - sign2_3)/4: f1_1 := 30*x_1: f1_2 := 31*x_1: f2_1 := 28*x_2: f2_2 := 29*x_2: f2_3 := 30*x_2: f1comb_1 := f1_1*step1_1: f1comb_2 := f1comb_1 + comp1nnxa_1*comp1nnxb_1 + f1_2*step1_2: f2comb_1 := f2_1*step2_1: f2comb_2 := f2comb_1 + comp2nnxa_1*comp2nnxb_1 + f2_2*step2_2: f2comb_3 := f2comb_2 + comp2nnxa_2*comp2nnxb_2 + f2_3*step2_3: mf := f1comb_2 + f2comb_3: myminfun := mf: myabsdevnod0_0 := myownabs(myminfun - 8853.5398913296): myreldevnod0_0 := myownabs(myminfun/8853.5398913296 - 1): myabsdevnod0_1 := myownabs(x_1 - 201.78466304432): myreldevnod0_1 := myownabs(x_1/201.78466304432 - 1): myabsdevnod0_2 := myownabs(x_2 - 100): myreldevnod0_2 := myownabs(x_2/100 - 1): myabsdevnod0_3 := myownabs(x_3 - 383.0709952729388): myreldevnod0_3 := myownabs(x_3/383.0709952729388 - 1): myabsdevnod0_4 := myownabs(x_4 - 420): myreldevnod0_4 := myownabs(x_4/420 - 1): myabsdevnod0_5 := myownabs(x_5 + 10.90760562936625): myreldevnod0_5 := myownabs(x_5/10.90760562936625 + 1): myabsdevnod0_6 := myownabs(x_6 - 0.07314814840034025): myreldevnod0_6 := myownabs(x_6/0.07314814840034025 - 1): myabsdevnod0_7 := myownabs(x_7 - 98.21533695567999): myreldevnod0_7 := myownabs(x_7/98.21533695567999 - 1): myabsdevnod0_8 := myownabs(x_8): myreldevnod0_8 := myownabs(x_8): myabsdevnod0_9 := myownabs(x_9): myreldevnod0_9 := myownabs(x_9): myabsdevnod0_10 := myownabs(x_10): myreldevnod0_10 := myownabs(x_10): myabsdevnod0_11 := myownabs(x_11 - 100): myreldevnod0_11 := myownabs(x_11/100 - 1): myabsdevnod0_12 := myownabs(x_12): myreldevnod0_12 := myownabs(x_12): myabsdevnod0_13 := myownabs(x_13 + 1): myreldevnod0_13 := myownabs(x_13 + 1): myabsdevnod0_14 := myownabs(x_14 + 1): myreldevnod0_14 := myownabs(x_14 + 1): myabsdevnod0_15 := myownabs(x_15 + 1): myreldevnod0_15 := myownabs(x_15 + 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)): myabsdevnod1_6 := (1/2)*(myownabs(myabsdevnod0_6 - myabsdevnod0_7) + (myabsdevnod0_6 + myabsdevnod0_7)): myreldevnod1_6 := (1/2)*(myownabs(myreldevnod0_6 - myreldevnod0_7) + (myreldevnod0_6 + myreldevnod0_7)): myabsdevnod1_8 := (1/2)*(myownabs(myabsdevnod0_8 - myabsdevnod0_9) + (myabsdevnod0_8 + myabsdevnod0_9)): myreldevnod1_8 := (1/2)*(myownabs(myreldevnod0_8 - myreldevnod0_9) + (myreldevnod0_8 + myreldevnod0_9)): myabsdevnod1_10 := (1/2)*(myownabs(myabsdevnod0_10 - myabsdevnod0_11) + (myabsdevnod0_10 + myabsdevnod0_11)): myreldevnod1_10 := (1/2)*(myownabs(myreldevnod0_10 - myreldevnod0_11) + (myreldevnod0_10 + myreldevnod0_11)): myabsdevnod1_12 := (1/2)*(myownabs(myabsdevnod0_12 - myabsdevnod0_13) + (myabsdevnod0_12 + myabsdevnod0_13)): myreldevnod1_12 := (1/2)*(myownabs(myreldevnod0_12 - myreldevnod0_13) + (myreldevnod0_12 + myreldevnod0_13)): myabsdevnod1_14 := (1/2)*(myownabs(myabsdevnod0_14 - myabsdevnod0_15) + (myabsdevnod0_14 + myabsdevnod0_15)): myreldevnod1_14 := (1/2)*(myownabs(myreldevnod0_14 - myreldevnod0_15) + (myreldevnod0_14 + myreldevnod0_15)): 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 := (1/2)*(myownabs(myabsdevnod1_4 - myabsdevnod1_6) + (myabsdevnod1_4 + myabsdevnod1_6)): myreldevnod2_4 := (1/2)*(myownabs(myreldevnod1_4 - myreldevnod1_6) + (myreldevnod1_4 + myreldevnod1_6)): myabsdevnod2_8 := (1/2)*(myownabs(myabsdevnod1_8 - myabsdevnod1_10) + (myabsdevnod1_8 + myabsdevnod1_10)): myreldevnod2_8 := (1/2)*(myownabs(myreldevnod1_8 - myreldevnod1_10) + (myreldevnod1_8 + myreldevnod1_10)): myabsdevnod2_12 := (1/2)*(myownabs(myabsdevnod1_12 - myabsdevnod1_14) + (myabsdevnod1_12 + myabsdevnod1_14)): myreldevnod2_12 := (1/2)*(myownabs(myreldevnod1_12 - myreldevnod1_14) + (myreldevnod1_12 + myreldevnod1_14)): 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)): myabsdevnod3_8 := (1/2)*(myownabs(myabsdevnod2_8 - myabsdevnod2_12) + (myabsdevnod2_8 + myabsdevnod2_12)): myreldevnod3_8 := (1/2)*(myownabs(myreldevnod2_8 - myreldevnod2_12) + (myreldevnod2_8 + myreldevnod2_12)): myabsdevnod4_0 := (1/2)*(myownabs(myabsdevnod3_0 - myabsdevnod3_8) + (myabsdevnod3_0 + myabsdevnod3_8)): myreldevnod4_0 := (1/2)*(myownabs(myreldevnod3_0 - myreldevnod3_8) + (myreldevnod3_0 + myreldevnod3_8)): zmyabsdevmax := myabsdevnod4_0: zmyreldevmax := myreldevnod4_0: obj := myminfun: fmcMinimum(obj): fmcEquality(ce_1, cons_1): fmcEquality(ce_2, cons_2): fmcEquality(ce_3, cons_3): fmcEquality(ce_4, cons_4): fmcEquality(ce_5, cons_5): fmcEquality(ce_6, cons_6): fmcEquality(ce_7, cons_7): fmcEquality(ce_8, cons_8): fmcEquality(ce_9, cons_9): fmcEquality(ce_10, cons_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(x_6, x_6): fmcEscortFloat(x_7, x_7): fmcEscortFloat(x_8, x_8): fmcEscortFloat(x_9, x_9): fmcEscortFloat(x_10, x_10): fmcEscortFloat(x_11, x_11): fmcEscortFloat(x_12, x_12): fmcEscortFloat(x_13, x_13): fmcEscortFloat(x_14, x_14): fmcEscortFloat(x_15, x_15): fmcEscortFloat('myminfun', myminfun): fmcEscortFloat('zmyabsdevmax', zmyabsdevmax): fmcEscortFloat('zmyreldevmax', zmyreldevmax): fmcControlMinimum(1, 0): # End Model