rr[j_, k_, l_, p_, q_] := rr15[j, k, l, p, q]
 
rr15[f_, sx_, a_, m_Integer, n_Integer] := Module[{na, w, zx}, 
     If[n < 1, Return[no$$iterations$$error]]; na = 0; zx = rat[sx, m]; 
      Label[o1]; If[na == n, Goto[o2]]; zx = r15f[f, zx, a, m]; na = na + 1; 
      Goto[o1]; Label[o2]; Return[zx]; ]

rat[x_, k_Integer] := Rationalize[x, 1/10^(k + 2)]
 
r15f[f_, sx_, a_, n_Integer] := r15fa[f, rat[sx, n], a, n]
 
r15fa[f_, sx_, a_, n_Integer] := (AccuracyGoal -> n; PrecisionGoal -> n; 
     WorkingPrecision -> n + 7; Module[{u, s, zx, o, nu, du, x1a, x2a, x3a, 
       x4a, x5a, x6a, x7a, x8a, x9a, x10a, x11a, x12a, x13a, x14a, x15a, fx1, 
       fx2, fx3, fx4, fx5, fx6, fx7, fx8, fx9, fx10, fx11, fx12, fx13, fx14, 
       fx15, fe, w}, {x1a, x2a, x3a, x4a, x5a, x6a, x7a, x8a, x9a, x10a, 
         x11a, x12a, x13a, x14a, x15a} = sx; 
       s = {x1 -> x1a, x2 -> x2a, x3 -> x3a, x4 -> x4a, x5 -> x5a, x6 -> x6a, 
         x7 -> x7a, x8 -> x8a, x9 -> x9a, x10 -> x10a, x11 -> x11a, 
         x12 -> x12a, x13 -> x13a, x14 -> x14a, x15 -> x15a}; 
       {fx1, fx2, fx3, fx4, fx5, fx6, fx7, fx8, fx9, fx10, fx11, fx12, fx13, 
         fx14, fx15, fe} = N[{D[f, x1], D[f, x2], D[f, x3], D[f, x4], 
           D[f, x5], D[f, x6], D[f, x7], D[f, x8], D[f, x9], D[f, x10], 
           D[f, x11], D[f, x12], D[f, x13], D[f, x14], D[f, x15], f} /. s, 
         n]; w = {fx1, fx2, fx3, fx4, fx5, fx6, fx7, fx8, fx9, fx10, fx11, 
         fx12, fx13, fx14, fx15, fe}; w = Rationalize[w, 1/10^(n + 2)]; 
       {fx1, fx2, fx3, fx4, fx5, fx6, fx7, fx8, fx9, fx10, fx11, fx12, fx13, 
         fx14, fx15, fe} = w; nu = -(a*fe); du = fx1^2 + fx2^2 + fx3^2 + 
         fx4^2 + fx5^2 + fx6^2 + fx7^2 + fx8^2 + fx9^2 + fx10^2 + fx11^2 + 
         fx12^2 + fx13^2 + fx14^2 + fx15^2; u = nu/du; 
       zx = N[sx + u*{fx1, fx2, fx3, fx4, fx5, fx6, fx7, fx8, fx9, fx10, 
            fx11, fx12, fx13, fx14, fx15}, n]; Return[zx]; Null; ])

rrs[z1_,z2_,z3_] := rrb[f,z1,9/5,z2,z3]

rrc[z1_, z2_, z3_, z4_] := rrb[f, z1, z2, z3, z4]
 
rra[v1_, v2_, v3_, v4_] := rrb[v1,v2,9/5,v3,v4]
 
rrb[f_, sx_, a_, m_Integer, n_Integer] := Module[{na, w, zx}, 
     If[n < 1, Return[no$$iterations$$error]]; na = 0; zx = rat[sx, m]; 
      Label[o1]; If[na == n, Goto[o2]]; zx = rfb[f, zx, a, m]; na = na + 1; 
      Goto[o1]; Label[o2]; Return[zx]; ]
 
rfb[f_, sx_, a_, n_Integer] := rf1b[f, rat[sx, n], a, n]
 
rf1b[f_, sx_, a_, n_Integer] := (AccuracyGoal -> n; PrecisionGoal -> n; 
     WorkingPrecision -> n + 7; Module[{u, s, zx, du, b, fx1, fx2, fx3, fx4, 
       fx5, fx6, fx7, fx8, fx9, fx10, fx11, fx12, fx13, fx14, fx15, fe, fe2, 
       w, st16}, b = a; st16 = N[sbx[{D[f, x1], D[f, x2], D[f, x3], D[f, x4], 
           D[f, x5], D[f, x6], D[f, x7], D[f, x8], D[f, x9], D[f, x10], 
           D[f, x11], D[f, x12], D[f, x13], D[f, x14], D[f, x15], f}, sx], 
         n]; st16 = Rationalize[st16, 1/10^(n + 2)]; 
       {fx1, fx2, fx3, fx4, fx5, fx6, fx7, fx8, fx9, fx10, fx11, fx12, fx13, 
         fx14, fx15, fe} = st16; du = fx1^2 + fx2^2 + fx3^2 + fx4^2 + fx5^2 + 
         fx6^2 + fx7^2 + fx8^2 + fx9^2 + fx10^2 + fx11^2 + fx12^2 + fx13^2 + 
         fx14^2 + fx15^2; u = -(fe/du); Goto[o2]; Label[o1]; b = (2*b)/3; 
       Label[o2]; zx = N[sx + b*u*{fx1, fx2, fx3, fx4, fx5, fx6, fx7, fx8, 
            fx9, fx10, fx11, fx12, fx13, fx14, fx15}, n]; 
       fe2 = N[sbx[f, zx], n]; If[fe2 < fe, Goto[o3]]; 
       If[b < 6^(-n), Return[gonetoosmall]]; Goto[o1]; Label[o3]; Return[zx]; 
       Null; ])
 
sbx[phi_, sxa_] := Module[{xaa, r}, xaa = sb[sxa]; r = phi /. xaa; Return[r]; 
      Null; ]
 
sb[wq_] := s15b[wq]
 
s15b[w_] := Module[{j1, j2, j3, a, j4, j5, j6, j7, j8, j9, j10, j11, j12, 
      j13, j14, j15}, {j1, j2, j3, j4, j5, j6, j7, j8, j9, j10, j11, j12, 
        j13, j14, j15} = w; a = {x1 -> j1, x2 -> j2, x3 -> j3, x4 -> j4, 
        x5 -> j5, x6 -> j6, x7 -> j7, x8 -> j8, x9 -> j9, x10 -> j10, 
        x11 -> j11, x12 -> j12, x13 -> j13, x14 -> j14, x15 -> j15}; 
      Return[a]; ]
 
FRTB[w1_, w2_, w3_] := Module[{zx, o}, o = FRTA[w1, w2, w3]; zx = xx15 /. o; 
      Return[zx]; ]
 
FRTA[ll_, sx_, n_Integer] := Module[{rx, o, x1a, x2a, x3a, x4a, x5a, x6a, 
      x7a, x8a, x9a, x10a, x11a, x12a, x13a, x14a, x15a}, 
     rx = Rationalize[sx, 1/10^(n + 2)]; {x1a, x2a, x3a, x4a, x5a, x6a, x7a, 
        x8a, x9a, x10a, x11a, x12a, x13a, x14a, x15a} = rx; 
      o = FindRoot[ll == zz15, {x1, x1a}, {x2, x2a}, {x3, x3a}, {x4, x4a}, 
        {x5, x5a}, {x6, x6a}, {x7, x7a}, {x8, x8a}, {x9, x9a}, {x10, x10a}, 
        {x11, x11a}, {x12, x12a}, {x13, x13a}, {x14, x14a}, {x15, x15a}, 
        {AccuracyGoal -> n, WorkingPrecision -> n + 7}]; o = N[o, n]; 
      Return[o]; ]
 
zz15 = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}
 
xx15 = {x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, x15}
 
FRT[w1_, w2_, w3_] := Module[{za, zb, zc}, za = FRTA[w1, w2, w3]; 
      zb = FRTB[w1, w2, w3]; zc = {zb, za}; Return[zc]; ]
 
SFRB[w1_, w2_, w3_, v_] := Module[{zx, o}, o = SFRA[w1, w2, w3, v]; 
      zx = xx15 /. o; Return[zx]; ]
 
SFRA[ll_, sx_, n_Integer, kk_Integer] := 
    Module[{rx, o, k, x1a, x2a, x3a, x4a, x5a, x6a, x7a, x8a, x9a, x10a, 
      x11a, x12a, x13a, x14a, x15a}, k = kk/(1 + kk); 
      rx = Rationalize[sx, 1/10^(n + 2)]; {x1a, x2a, x3a, x4a, x5a, x6a, x7a, 
        x8a, x9a, x10a, x11a, x12a, x13a, x14a, x15a} = rx; 
      o = FindRoot[ll == zz15, {x1, k*x1a, x1a/k}, {x2, k*x2a, x2a/k}, 
        {x3, k*x3a, x3a/k}, {x4, k*x4a, x4a/k}, {x5, k*x5a, x5a/k}, 
        {x6, k*x6a, x6a/k}, {x7, k*x7a, x7a/k}, {x8, k*x8a, x8a/k}, 
        {x9, k*x9a, x9a/k}, {x10, k*x10a, x10a/k}, {x11, k*x11a, x11a/k}, 
        {x12, k*x12a, x12a/k}, {x13, k*x13a, x13a/k}, {x14, k*x14a, x14a/k}, 
        {x15, k*x15a, x15a/k}, {AccuracyGoal -> n, WorkingPrecision -> 
          n + 7}]; o = N[o, n]; Return[o]; Null; ]
 
SFR[w1_, w2_, w3_, v_] := Module[{za, zb, zc}, za = SFRA[w1, w2, w3, v]; 
      zb = SFRB[w1, w2, w3, v]; zc = {zb, za}; Return[zc]; ]
 
fsb[z_] := sbx[f, z]
 
ffsb[v_] := sbx[ff, v]
 
ff0 = eqtsxx15

eqtsxx15 = {2*(-((3 - 6*x1 + 4*x1^2)*(e4 - 4*(-1 + e4)*x1 + 4*(-1 + e4)*x1^2)*
         (b*e5*(-1 + x2)*(-1 + x3) + 2*x2*x7 - x3*((-2 + x2)*x6 + x2*x7))) + 
       x1*(6 - 3*e4 + 2*(-9 + 7*e4)*x1 - 20*(-1 + e4)*x1^2 + 
         8*(-1 + e4)*x1^3)*(b*e5 - x3*(-2 + b*e5 + 2*x6) + 
         x2*(2 - b*e5 - 2*x7 + x3*(-2 + b*e5 + x6 + x7)))), 
     e3*(6 - 12*x1 + 8*x1^2)*(1 - (-1 + e4)*(-1 + 2*x1)^3)*
       (b*e5*(-1 + x3) + (-2 + x3)*(-e5 + (-1 + e5)*x3)*x4 + 
        x3*(2 - e5 + (-1 + e5)*x3)*(-1 + x5 + x6)) - 
      (b*e5*(-1 + x3) + x3*(-2 + e5 + x3 - e5*x3)*x6 + 
        (-2 + x3)*(-e5 + (-1 + e5)*x3)*x7)*
       (((1 - 2*x1)*x1*(6 - 3*e4 + 2*(-9 + 7*e4)*x1 - 20*(-1 + e4)*x1^2 + 
           8*(-1 + e4)*x1^3)*(3*b*e5 + x3*(2 - 3*b*e5 + 2*x6) + 
           x2*(2 - 3*b*e5 + x3*(-2 + 3*b*e5 - x6 - x7) + 2*x7)))/
         (2*(e5 + (-1 + e5)*x2*(-1 + x3) + x3 - e5*x3)) + 
        ((6 - 12*x1 + 8*x1^2)*((-1 + x1 + (-1 + e4)*(-1 + 2*x1)^3)*
            (b*e5*(-1 + x2)*(-1 + x3) - x2*(-2 + x3)*x4 + (-2 + x2)*x3*
              (-1 + x5 + x6)) + x1*(-((-2 + x2)*x3*x5) + x2*(-2 + x3)*
              (-1 + x4 + x7))))/(2 + 2*(-1 + e5)*(-1 + x2)*(-1 + x3))), 
     e3*(6 - 12*x1 + 8*x1^2)*(1 - (-1 + e4)*(-1 + 2*x1)^3)*
       ((-2 + x2)*(-e5 + (-1 + e5)*x2)*x5 + x2*(2 - e5 + (-1 + e5)*x2)*
         (-1 + x4 + x7)) - (b*e5*(-1 + x2) + (-2 + x2)*(-e5 + (-1 + e5)*x2)*
         x6 + x2*(-2 + e5 + x2 - e5*x2)*x7)*
       (((1 - 2*x1)*x1*(6 - 3*e4 + 2*(-9 + 7*e4)*x1 - 20*(-1 + e4)*x1^2 + 
           8*(-1 + e4)*x1^3)*(3*b*e5 + x3*(2 - 3*b*e5 + 2*x6) + 
           x2*(2 - 3*b*e5 + x3*(-2 + 3*b*e5 - x6 - x7) + 2*x7)))/
         (2*(e5 + (-1 + e5)*x2*(-1 + x3) + x3 - e5*x3)) + 
        ((6 - 12*x1 + 8*x1^2)*(x1*(b*e5*(-1 + x2)*(-1 + x3) - 
             x2*(-2 + x3)*x4 + (-2 + x2)*x3*(-1 + x5 + x6)) + 
           (-1 + x1 + (-1 + e4)*(-1 + 2*x1)^3)*(-((-2 + x2)*x3*x5) + 
             x2*(-2 + x3)*(-1 + x4 + x7))))/(2 + 2*(-1 + e5)*(-1 + x2)*
           (-1 + x3))), e3*(6 - 12*x1 + 8*x1^2)*
       (1 - (-1 + e4)*(-1 + 2*x1)^3) + 
      ((1 - 2*x1)*x1*(6 - 3*e4 + 2*(-9 + 7*e4)*x1 - 20*(-1 + e4)*x1^2 + 
         8*(-1 + e4)*x1^3)*(3*b*e5 + x3*(2 - 3*b*e5 + 2*x6) + 
         x2*(2 - 3*b*e5 + x3*(-2 + 3*b*e5 - x6 - x7) + 2*x7)))/
       (2*(e5 + (-1 + e5)*x2*(-1 + x3) + x3 - e5*x3)) + 
      ((6 - 12*x1 + 8*x1^2)*((-1 + x1 + (-1 + e4)*(-1 + 2*x1)^3)*
          (b*e5*(-1 + x2)*(-1 + x3) - x2*(-2 + x3)*x4 + (-2 + x2)*x3*
            (-1 + x5 + x6)) + x1*(-((-2 + x2)*x3*x5) + x2*(-2 + x3)*
            (-1 + x4 + x7))))/(2 + 2*(-1 + e5)*(-1 + x2)*(-1 + x3)), 
     e3*(6 - 12*x1 + 8*x1^2)*(1 - (-1 + e4)*(-1 + 2*x1)^3)*
       (e3*(-2 + x2)*(e5 + (-1 + e5)*x2*(-1 + x3) + x3 - e5*x3) + 
        (1 - x3)*(b*e5*(-1 + x2) + x2*(-2 + e5 + x2 - e5*x2)*x4 - 
          (-2 + x2)*(-e5 + (-1 + e5)*x2)*(-1 + x5 + x6))) - 
      (1 - x3)*(b*e5*(-1 + x2) + (-2 + x2)*(-e5 + (-1 + e5)*x2)*x6 + 
        x2*(-2 + e5 + x2 - e5*x2)*x7)*
       (((1 - 2*x1)*x1*(6 - 3*e4 + 2*(-9 + 7*e4)*x1 - 20*(-1 + e4)*x1^2 + 
           8*(-1 + e4)*x1^3)*(3*b*e5 + x3*(2 - 3*b*e5 + 2*x6) + 
           x2*(2 - 3*b*e5 + x3*(-2 + 3*b*e5 - x6 - x7) + 2*x7)))/
         (2*(e5 + (-1 + e5)*x2*(-1 + x3) + x3 - e5*x3)) + 
        ((6 - 12*x1 + 8*x1^2)*((-1 + x1 + (-1 + e4)*(-1 + 2*x1)^3)*
            (b*e5*(-1 + x2)*(-1 + x3) - x2*(-2 + x3)*x4 + (-2 + x2)*x3*
              (-1 + x5 + x6)) + x1*(-((-2 + x2)*x3*x5) + x2*(-2 + x3)*
              (-1 + x4 + x7))))/(2 + 2*(-1 + e5)*(-1 + x2)*(-1 + x3))), 
     e3*(6 - 12*x1 + 8*x1^2)*(1 - (-1 + e4)*(-1 + 2*x1)^3)*
       (e3*(-2 + x3)*(e5 + (-1 + e5)*x2*(-1 + x3) + x3 - e5*x3) + 
        (1 - x2)*(x3*(-2 + e5 + x3 - e5*x3)*x5 - (-2 + x3)*
           (-e5 + (-1 + e5)*x3)*(-1 + x4 + x7))) - 
      (1 - x2)*(b*e5*(-1 + x3) + x3*(-2 + e5 + x3 - e5*x3)*x6 + 
        (-2 + x3)*(-e5 + (-1 + e5)*x3)*x7)*
       (((1 - 2*x1)*x1*(6 - 3*e4 + 2*(-9 + 7*e4)*x1 - 20*(-1 + e4)*x1^2 + 
           8*(-1 + e4)*x1^3)*(3*b*e5 + x3*(2 - 3*b*e5 + 2*x6) + 
           x2*(2 - 3*b*e5 + x3*(-2 + 3*b*e5 - x6 - x7) + 2*x7)))/
         (2*(e5 + (-1 + e5)*x2*(-1 + x3) + x3 - e5*x3)) + 
        ((6 - 12*x1 + 8*x1^2)*(x1*(b*e5*(-1 + x2)*(-1 + x3) - 
             x2*(-2 + x3)*x4 + (-2 + x2)*x3*(-1 + x5 + x6)) + 
           (-1 + x1 + (-1 + e4)*(-1 + 2*x1)^3)*(-((-2 + x2)*x3*x5) + 
             x2*(-2 + x3)*(-1 + x4 + x7))))/(2 + 2*(-1 + e5)*(-1 + x2)*
           (-1 + x3))), e3*(6 - 12*x1 + 8*x1^2)*
       (1 - (-1 + e4)*(-1 + 2*x1)^3) + 
      ((1 - 2*x1)*x1*(6 - 3*e4 + 2*(-9 + 7*e4)*x1 - 20*(-1 + e4)*x1^2 + 
         8*(-1 + e4)*x1^3)*(3*b*e5 + x3*(2 - 3*b*e5 + 2*x6) + 
         x2*(2 - 3*b*e5 + x3*(-2 + 3*b*e5 - x6 - x7) + 2*x7)))/
       (2*(e5 + (-1 + e5)*x2*(-1 + x3) + x3 - e5*x3)) + 
      ((6 - 12*x1 + 8*x1^2)*(x1*(b*e5*(-1 + x2)*(-1 + x3) - x2*(-2 + x3)*x4 + 
           (-2 + x2)*x3*(-1 + x5 + x6)) + (-1 + x1 + (-1 + e4)*(-1 + 2*x1)^3)*
          (-((-2 + x2)*x3*x5) + x2*(-2 + x3)*(-1 + x4 + x7))))/
       (2 + 2*(-1 + e5)*(-1 + x2)*(-1 + x3)), 0, 0, 0, 0, 0, 0, 0, 0}
  
sqsum[kk_] := Module[{r, f1, f2, f3, f4, f5, f6, f7, f8, f9, f10, f11, f12, 
      f13, f14, f15}, {f1, f2, f3, f4, f5, f6, f7, f8, f9, f10, f11, f12, 
        f13, f14, f15} = kk; r = f1*f1 + f2^2 + f3*f3 + f4^2 + f5*f5 +
     f6^2 + f7^2 + f8*f8 + f9*f9 + f10*f10 + f11^2 + f12^2 + f13^2 + 
      f14^2 +  f15^2; Return[r]; ]

s3sum[kk_] := Module[{r, f1, f2, f3}, {f1, f2, f3} = kk; r = f1 + f2 + f3; 
      Return[r]; ]
   
shapval = {1/3 + (-2*b1 + b2 + b3)/6, 1/3 + (b1 - 2*b2 + b3)/6, 
     1/3 + (b1 + b2 - 2*b3)/6}
 
nucleo = {1/3, 1/3, 1/3}
 
memo = {nucleolus,value,vector,is,correct,ONLY,when,b1b2b3,quantities,are,small,enough}

end = lastline