rr[v1_, v2_, v3_, v4_] := rr42[f, v1, v2, v3, v4] rr42[f_, sx_, a_, m_Integer, n_Integer] := Module[{na, w, zx, st43d}, If[n < 1, Return[no$$iterations$$error]]; na = 0; zx = rat[sx, m]; st43d = rf1ba[f, m]; Label[o1]; If[na == n, Goto[o2]]; zx = r42f[f, zx, a, m, st43d]; na = na + 1; Goto[o1]; Label[o2]; Return[zx]; ] rat[x_, k_Integer] := Rationalize[x, 1/10^(k + 2)] rf1ba[f_, n_Integer] := (AccuracyGoal -> n; PrecisionGoal -> n; WorkingPrecision -> n + 7; Module[{st43a}, st43a = {{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], D[f, x16], D[f, x17], D[f, x18], D[f, x19], D[f, x20], D[f, x21], D[f, x22], D[f, x23], D[f, x24], D[f, x25], D[f, x26], D[f, x27], D[f, x28], D[f, x29], D[f, x30], D[f, x31], D[f, x32], D[f, x33], D[f, x34], D[f, x35], D[f, x36], D[f, x37], D[f, x38], D[f, x39], D[f, x40], D[f, x41]}, D[f, x42], f}; Return[st43a]; ]) r42f[f_, sx_, a_, n_Integer, st43c_] := r42fa[f, rat[sx, n], a, n, st43c] r42fa[f_, sx_, a_, n_Integer, st43b_] := (AccuracyGoal -> n; PrecisionGoal -> n; WorkingPrecision -> n + 7; Module[{u, zx, nu, du, fe, w43, v42}, w43 = N[sbx[st43b, sx], n]; w43 = Rationalize[w43, 1/10^(n + 2)]; {v42, fe} = w43; nu = -(a*fe); du = sqsum[v42]; u = nu/du; zx = N[sx + u*v42, n]; Return[zx]; Null;]) sbx[phi_, sxa_] := Module[{xaa, r}, xaa = sb[sxa]; r = phi /. xaa; Return[r]; Null; ] sb[wq_] := s42b[wq] s42b[w_] := Module[{j1, j2, j3, a, j4, j5, j6, j7, j8, j9, j10, j11, j12, j13, j14, j15, j16, j17, j18, j19, j20, j21, j22, j23, j24, j25, j26, j27, j28, j29, j30, j31, j32, j33, j34, j35, j36, j37, j38, j39, j40, j41, j42}, {j1, j2, j3, j4, j5, j6, j7, j8, j9, j10, j11, j12, j13, j14, j15, j16, j17, j18, j19, j20, j21, j22, j23, j24, j25, j26, j27, j28, j29, j30, j31, j32, j33, j34, j35, j36, j37, j38, j39, j40, j41, j42} = 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, x16 -> j16, x17 -> j17, x18 -> j18, x19 -> j19, x20 -> j20, x21 -> j21, x22->j22, x23->j23, x24->j24, x25->j25, x26->j26, x27->j27, x28->j28, x29-> j29, x30->j30, x31->j31, x32->j32, x33->j33, x34->j34, x35->j35, x36->j36, x37->j37, x38->j38, x39->j39, x40->j40, x41->j41, x42->j42}; Return[a]; ] sqsum[jj_] := jj . jj rrs[z1_, z2_, z3_] := rrb2[f, z1, 9/5, z2, z3] rrb2[f_, sx_, a_, m_Integer, n_Integer] := Module[{na, w, zx, st43d}, If[n < 1, Return[no$$iterations$$error]]; st43d = rf1ba[f, n]; na = 0; zx = rat[sx, m]; Label[o1]; If[na == n, Goto[o2]]; zx = rfbb[f, zx, a, m, st43d]; na = na + 1; Goto[o1]; Label[o2]; Return[zx]; ] rfbb[f_, sx_, a_, n_Integer, st43d_] := rf1bb[f, rat[sx, n], a, n, st43d] rf1bb[f_, sx_, a_, n_Integer, st43c_] := (AccuracyGoal -> n; PrecisionGoal -> n; WorkingPrecision -> n + 7; Module[{u, zx, du, b, fe, fe2, st42, st43}, b = a; st43 = N[sbx[st43c, sx], n]; st43 = Rationalize[st43, 1/10^(n + 2)]; {st42, fe} = st43; du = sqsum[st42]; u = -(fe/du); Goto[o2]; Label[o1]; b = (2*b)/3; Label[o2]; zx = N[sx + b*u*st42,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; ]) rre[v1_, v2_, v3_, v4_] := Module[{r}, Clear[OUT, procedure, t1, t2]; t1 = dat[]; r = rr[v1, v2, v3, v4]; OUT = r; t2 = dat[]; procedure = {v2, v3, v4}; Save["out.rre", ps, memo, OUT, procedure, t1, t2]; Return[r]; ] OUT = OUT dat[] := Module[{r, t1, t2, t3, t4, t5, t6}, {t1, t2, t3, t4, t5, t6} = Date[]; r = {t3, t2, t1, t4, t5}; Return[r]; ] rrce[v1_, v2_, v3_, v4_] := Module[{r, s, t1, t2}, t1 = dat[]; r = rrc[v1, v2, v3, v4]; Clear[OUT]; OUT = r; t2 = dat[]; s = {v2, v3, v4}; Save["out.rrce", ps, memo, OUT, s, t1, t2]; Return[r]; ] rrc[z1_, z2_, z3_, z4_] := rrb2[f, z1, z2, z3, z4] rrfrce[v1_, v2_, v3_, v4_] := Module[{r, s, t1, t2}, t1 = dat[]; r = rrfrc[v1, v2, v3, v4]; Clear[OUT]; OUT = r; t2 = dat[]; s = {v2, v3, v4}; Save["out.rrfrce", ps, memo, OUT, s, t1, t2]; Return[r]; ] rrfrc[z1_, z2_, z3_, z4_] := rrfrb[f, z1, z2, z3, z4] rrfrb[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 = FR1c[f, zx, a, m]; na = na + 1; Goto[o1]; Label[o2]; Return[zx]; ] FR1c[f_, sx_, a_, n_Integer] := FR1ac[f, rat[sx, n], a, n] FR1ac[f_, sx_, a_, n_Integer] := (AccuracyGoal -> n; PrecisionGoal -> n; WorkingPrecision -> n + 7; Module[{u, zx, du, b, fe, fe2, st21, st22, se}, b = a; st43 = {FR1[sx, n], N[fsb[sx], n]}; st43 = Rationalize[st43, 1/10^(n + 2)]; {st42, fe} = st43; du = sqsum[st42]; u = -(fe/du); Goto[o2]; Label[o1]; b = (6*b)/7; Label[o2]; zx = N[sx + b*u*st42, n]; fe2 = N[sbx[f, zx], n]; If[fe2 < fe, Goto[o3]]; If[b < 5^(-n), Goto[o4]]; Goto[o1]; Label[o4]; zx = rrc[zx, a, n, 1]; Label[o3]; Return[zx]; Null; ]) FR1[u1_, u2_] := FRT1B[ff, u1, u2] FRT1B[w1_, w2_, w3_] := Module[{zx, o}, o = FRT1A[w1, w2, w3]; zx = xx42 /. o; Return[zx]; ] FRT1A[ll_, sx_, n_Integer] := Module[{rx, o, x1a, x2a, x3a, x4a, x5a, x6a, x7a, x8a, x9a, x10a, x11a, x12a, x13a, x14a, x15a, x16a, x17a, x18a, x19a, x20a, x21a, x22a, x23a, x24a, x25a, x26a, x27a, x28a, x29a, x30a, x31a, x32a, x33a, x34a, x35a, x36a, x37a, x38a, x39a, x40a, x41a, x42a, e7, e8, e9, e10, e11, e12, e13, e14, e15, e16, e17, e18, e19, e20, e21, e22, e23, e24, e25, e26, e27, e28, e29, e30, e31, e32, e33, e34, e35, e36, e37, e38, e39, e40, e41, e42, e43, e44, e45, e46, e47, e48}, rx = Rationalize[sx, 1/10^(n + 2)]; {x1a, x2a, x3a, x4a, x5a, x6a, x7a, x8a, x9a, x10a, x11a, x12a, x13a, x14a, x15a, x16a, x17a, x18a, x19a, x20a, x21a, x22a, x23a, x24a, x25a, x26a, x27a, x28a, x29a, x30a, x31a, x32a, x33a, x34a, x35a, x36a, x37a, x38a, x39a, x40a, x41a, x42a} = rx; {e7, e8, e9, e10, e11, e12, e13, e14, e15, e16, e17, e18, e19, e20, e21, e22, e23, e24, e25, e26, e27 e28, e29, e30, e31, e32, e33, e34, e35, e36, e37, e38, e39, e40, e41, e42, e43, e44, e45, e46, e47, e48} = ll; o = FindRoot[{e7 == 0, e8 == 0, e9 == 0, e10 == 0, e11 == 0, e12 == 0, e13 == 0, e14 == 0, e15 == 0, e16 == 0, e17 == 0, e18 == 0, e19 == 0, e20 == 0, e21 == 0, e22 == 0, e23 == 0, e24 == 0, e25 == 0, e26 == 0, e27 == 0, e28 == 0, e29 == 0, e30 == 0, e31 == 0, e32 == 0, e33 == 0, e34 == 0, e35 == 0, e36 == 0, e37 == 0, e38 == 0, e39 == 0, e40 == 0, e41 == 0, e42 == 0, e43 == 0, e44 == 0, e45 == 0, e46 == 0, e47 == 0, e48 == 0}, {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}, {x16, x16a}, {x17, x17a}, {x18, x18a}, {x19, x19a}, {x20, x20a}, {x21, x21a}, {x22, x22a}, {x23, x23a}, {x24, x24a}, {x25, x25a}, {x26, x26a}, {x27, x27a}, {x28, x28a}, {x29, x29a}, {x30, x30a}, {x31, x31a}, {x32, x32a}, {x33, x33a}, {x34, x34a}, {x35, x35a}, {x36, x36a}, {x37, x37a}, {x38, x38a}, {x39, x39a}, {x40, x40a}, {x41, x41a}, {x42, x42a}, {AccuracyGoal -> n, MaxIterations -> 1, WorkingPrecision -> n + 7}]; o = N[o, n]; Return[o]; ] xx42 = {x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, x15, x16, x17, x18, x19, x20, x21, x22, x23, x24, x25, x26, x27, x28, x29, x30, x31, x32, x33, x34, x35, x36, x37, x38, x39, x40, x41, x42} fsb[z_] := sbx[f, z]