If Cut/Copy and Paste fails, then click here for download.
Model tp087r1 ! Source version 7 ! Simple APM canonicalizer version 1.3 ! APM backend version 1.1 Parameters 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)) midget = 2.0e-38 c2 = 2 c4 = 4 c8 = 8 c16 = 16 ct10 = -2186927100809997/c16 ct11 = 9112324500347/c4 ct12 = -121498919997/c8 ct13 = 101249775/c2 ct14 = -1350001/c16 ct15 = 225/c4 ct20 = -2999100089997/c16 ct21 = 37493500189/c4 ct22 = -1499879997/c8 ct23 = 3749925/c2 ct24 = -150001/c16 ct25 = 75/c4 ct30 = -191985600359997/c16 ct31 = 299987000067 ct32 = -23999519997/c8 ct33 = 14999925 ct34 = -600001/c16 ct35 = 75/c2 sign1_0 = 1 sign1_3 = -1 sign2_0 = 1 sign2_5 = -1 End Parameters Variables x_1 = 390, >= 0, <= 400 x_2 = 1000, >= 0, <= 1000 x_3 = 419.5, >= 340, <= 420 x_4 = 340.5, >= 340, <= 420 x_5 = ivx5, >= lbx5, <= ubx5 x_6 = 0.5, >= 0, <= 0.5236 obj End Variables Intermediates 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 - 299 comp1_2 = x_1 - 301 comp2_1 = x_2 - 99 comp2_2 = x_2 - 101 comp2_3 = x_2 - 199 comp2_4 = x_2 - 201 sign1x_1 = comp1_1/sqrt(comp1_1^2 + midget) sign1x_2 = comp1_2/sqrt(comp1_2^2 + midget) sign2x_1 = comp2_1/sqrt(comp2_1^2 + midget) sign2x_2 = comp2_2/sqrt(comp2_2^2 + midget) sign2x_3 = comp2_3/sqrt(comp2_3^2 + midget) sign2x_4 = comp2_4/sqrt(comp2_4^2 + midget) sign1_1 = sign1x_1 sign1_2 = sign1x_2 sign2_1 = sign2x_1 sign2_2 = sign2x_2 sign2_3 = sign2x_3 sign2_4 = sign2x_4 step1_1 = (1 + sign1_0)*(1 - sign1_1)/4 step1_2 = (1 + sign1_1)*(1 - sign1_2)/4 step1_3 = (1 + sign1_2)*(1 - sign1_3)/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 step2_4 = (1 + sign2_3)*(1 - sign2_4)/4 step2_5 = (1 + sign2_4)*(1 - sign2_5)/4 f1_1 = 30*x_1 f1_2 = (((((ct15)*x_1 + ct14)*x_1 + ct13)*x_1 + ct12)*x_1 + ct11)*x_1 + ct10 f1_3 = 31*x_1 f2_1 = 28*x_2 f2_2 = (((((ct25)*x_2 + ct24)*x_2 + ct23)*x_2 + ct22)*x_2 + ct21)*x_2 + ct20 f2_3 = 29*x_2 f2_4 = (((((ct35)*x_2 + ct34)*x_2 + ct33)*x_2 + ct32)*x_2 + ct31)*x_2 + ct30 f2_5 = 30*x_2 f1comb_1 = f1_1*step1_1 f1comb_2 = f1comb_1 + f1_2*step1_2 f1comb_3 = f1comb_2 + f1_3*step1_3 f2comb_1 = f2_1*step2_1 f2comb_2 = f2comb_1 + f2_2*step2_2 f2comb_3 = f2comb_2 + f2_3*step2_3 f2comb_4 = f2comb_3 + f2_4*step2_4 f2comb_5 = f2comb_4 + f2_5*step2_5 mf = f1comb_3 + f2comb_5 myminfun = mf myabsdevnod0_0 = abs(myminfun - 8854.997011124169) myreldevnod0_0 = abs(myminfun/8854.997011124169 - 1) myabsdevnod0_1 = abs(x_1 - 202.7039339773582) myreldevnod0_1 = abs(x_1/202.7039339773582 - 1) myabsdevnod0_2 = abs(x_2 - 99.06588958062389) myreldevnod0_2 = abs(x_2/99.06588958062389 - 1) myabsdevnod0_3 = abs(x_3 - 383.1404626009829) myreldevnod0_3 = abs(x_3/383.1404626009829 - 1) myabsdevnod0_4 = abs(x_4 - 420) myreldevnod0_4 = abs(x_4/420 - 1) myabsdevnod0_5 = abs(x_5 + 11.11782763722497) myreldevnod0_5 = abs(x_5/11.11782763722497 + 1) myabsdevnod0_6 = abs(x_6 - 0.07238952863601766) myreldevnod0_6 = abs(x_6/0.07238952863601766 - 1) myabsdevnod1_0 = (1/2)*(abs(myabsdevnod0_0 - myabsdevnod0_1) + (myabsdevnod0_0 + myabsdevnod0_1)) myreldevnod1_0 = (1/2)*(abs(myreldevnod0_0 - myreldevnod0_1) + (myreldevnod0_0 + myreldevnod0_1)) myabsdevnod1_2 = (1/2)*(abs(myabsdevnod0_2 - myabsdevnod0_3) + (myabsdevnod0_2 + myabsdevnod0_3)) myreldevnod1_2 = (1/2)*(abs(myreldevnod0_2 - myreldevnod0_3) + (myreldevnod0_2 + myreldevnod0_3)) myabsdevnod1_4 = (1/2)*(abs(myabsdevnod0_4 - myabsdevnod0_5) + (myabsdevnod0_4 + myabsdevnod0_5)) myreldevnod1_4 = (1/2)*(abs(myreldevnod0_4 - myreldevnod0_5) + (myreldevnod0_4 + myreldevnod0_5)) myabsdevnod1_6 = myabsdevnod0_6 myreldevnod1_6 = myreldevnod0_6 myabsdevnod2_0 = (1/2)*(abs(myabsdevnod1_0 - myabsdevnod1_2) + (myabsdevnod1_0 + myabsdevnod1_2)) myreldevnod2_0 = (1/2)*(abs(myreldevnod1_0 - myreldevnod1_2) + (myreldevnod1_0 + myreldevnod1_2)) myabsdevnod2_4 = (1/2)*(abs(myabsdevnod1_4 - myabsdevnod1_6) + (myabsdevnod1_4 + myabsdevnod1_6)) myreldevnod2_4 = (1/2)*(abs(myreldevnod1_4 - myreldevnod1_6) + (myreldevnod1_4 + myreldevnod1_6)) myabsdevnod3_0 = (1/2)*(abs(myabsdevnod2_0 - myabsdevnod2_4) + (myabsdevnod2_0 + myabsdevnod2_4)) myreldevnod3_0 = (1/2)*(abs(myreldevnod2_0 - myreldevnod2_4) + (myreldevnod2_0 + myreldevnod2_4)) zmyabsdevmax = myabsdevnod3_0 zmyreldevmax = myreldevnod3_0 End Intermediates Equations obj = myminfun cons_1 = 0 cons_2 = 0 cons_3 = 0 cons_4 = 0 End Equations End Model