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]]; zy = r42f[f, zx, a, m, st43d]; na = na + 1; Goto[o1]; Label[o2]; Return[zx]; ] rat[y_, k_Integer] := Rationalize[y, 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_, sy_, a_, n_Integer, st43c_] := r42fa[f, rat[sy, n], a, n, st43c] r42fa[f_, sy_, a_, n_Integer, st43b_] := (AccuracyGoal -> n; PrecisionGoal -> n; WorkingPrecision -> n + 7; Module[{u, zy, nu, du, fe, w43, v42}, w43 = N[sby[st43b, sy], n]; w43 = Rationalize[w43, 1/10^(n + 2)]; {v42, fe} = w43; nu = -(a*fe); du = sqsum[v42]; u = nu/du; zy = N[sy + u*v42, n]; Return[zy]; Null;]) sby[phi_, sya_] := Module[{yaa, r}, yaa = sb[sya]; r = phi /. yaa; 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} = w; a = {y1 -> j1, y2 -> j2, y3 -> j3, y4 -> j4, y5 -> j5, y6 -> j6, y7 -> j7, y8 -> j8, y9 -> j9, y10 -> j10, y11 -> j11, y12 -> j12, y13 -> j13, y14 -> j14, y15 -> j15, y16 -> j16, y17 -> j17, y18 -> j18, y19 -> j19, y20 -> j20, y21 -> j21}; Return[a]; ] sqsum[jj_] := jj . jj rrs[z1_, z2_, z3_] := rrb2[f, z1, 9/5, z2, z3] rrb2[f_, sy_, a_, m_Integer, n_Integer] := Module[{na, w, zy, st22d}, If[n < 1, Return[no$$iterations$$error]]; st22d = rf1ba[f, n]; na = 0; zy = rat[sy, m]; Label[o1]; If[na == n, Goto[o2]]; zy = rfbb[f, zy, a, m, st22d]; na = na + 1; Goto[o1]; Label[o2]; Return[zy]; ] rfbb[f_, sy_, a_, n_Integer, st22d_] := rf1bb[f, rat[sy, n], a, n, st22d] rf1bb[f_, sy_, a_, n_Integer, st22c_] := (AccuracyGoal -> n; PrecisionGoal -> n; WorkingPrecision -> n + 7; Module[{u, zy, du, b, fe, fe2, st21, st22}, b = a; st22 = N[sby[st22c, sy], n]; st22 = Rationalize[st22, 1/10^(n + 2)]; {st21, fe} = st22; du = sqsum[st21]; u = -(fe/du); Goto[o2]; Label[o1]; b = (2*b)/3; Label[o2]; zy = N[sy + b*u*st21, n]; fe2 = N[sby[f, zy], n]; If[fe2 < fe, Goto[o3]]; If[b < 6^(-n), Return[gonetoosmall]]; Goto[o1]; Label[o3]; Return[zy]; 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_, sy_, a_, m_Integer, n_Integer] := Module[{na, w, zy}, If[n < 1, Return[no$$iterations$$error]]; na = 0; zy = rat[sy, m]; Label[o1]; If[na == n, Goto[o2]]; zy = FR1c[f, zy, a, m]; na = na + 1; Goto[o1]; Label[o2]; Return[zy]; ] FR1c[f_, sy_, a_, n_Integer] := FR1ac[f, rat[sy, n], a, n] FR1ac[f_, sy_, a_, n_Integer] := (AccuracyGoal -> n; PrecisionGoal -> n; WorkingPrecision -> n + 7; Module[{u, zy, du, b, fe, fe2, st21, st22, se}, b = a; st22 = {FR1[sy, n], N[fsb[sy], n]}; st22 = Rationalize[st22, 1/10^(n + 2)]; {st21, fe} = st22; du = sqsum[st21]; u = -(fe/du); Goto[o2]; Label[o1]; b = (6*b)/7; Label[o2]; zy = N[sy + b*u*st21, n]; fe2 = N[sby[f, zy], n]; If[fe2 < fe, Goto[o3]]; If[b < 5^(-n), Goto[o4]]; Goto[o1]; Label[o4]; zy = rrc[zy, a, n, 1]; Label[o3]; Return[zy]; Null; ]) FR1[u1_, u2_] := FRT1B[ff, u1, u2] FRT1B[w1_, w2_, w3_] := Module[{zy, o}, o = FRT1A[w1, w2, w3]; zy = yy21 /. o; Return[zy]; ] FRT1A[ll_, sy_, n_Integer] := Module[{ry, o, y1a, y2a, y3a, y4a, y5a, y6a, y7a, y8a, y9a, y10a, y11a, y12a, y13a, y14a, y15a, y16a, y17a, y18a, y19a, y20a, y21a, e7, e8, e9, e10, e11, e12, e13, e14, e15, e16, e17, e18, e19, e20, e21, e22, e23, e24, e25, e26, e27}, ry = Rationalize[sy, 1/10^(n + 2)]; {y1a, y2a, y3a, y4a, y5a, y6a, y7a, y8a, y9a, y10a, y11a, y12a, y13a, y14a, y15a, y16a, y17a, y18a, y19a, y20a, y21a} = ry; {e7, e8, e9, e10, e11, e12, e13, e14, e15, e16, e17, e18, e19, e20, e21, e22, e23, e24, e25, e26, e27} = 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}, {y1, y1a}, {y2, y2a}, {y3, y3a}, {y4, y4a}, {y5, y5a}, {y6, y6a}, {y7, y7a}, {y8, y8a}, {y9, y9a}, {y10, y10a}, {y11, y11a}, {y12, y12a}, {y13, y13a}, {y14, y14a}, {y15, y15a}, {y16, y16a}, {y17, y17a}, {y18, y18a}, {y19, y19a}, {y20, y20a}, {y21, y21a}, {AccuracyGoal -> n, MaxIterations -> 1, WorkingPrecision -> n + 7}]; o = N[o, n]; Return[o]; ] yy21 = {y1, y2, y3, y4, y5, y6, y7, y8, y9, y10, y11, y12, y13, y14, y15, y16, y17, y18, y19, y20, y21} fsb[z_] := sby[f, z]