xrrceb[a3_, b5_] := Module[{r}, Clear[s1, s1a, s1b, s1c, sy0, s2, t1, t2, u]; t1 = dat[]; {s1a, s1b, s1c} = a3; s2 = b5; u = $xrrceb$; s1 = {sy0, s1b, s1c}; Save["start.data.gen", memo, u, s1, s2, t1]; r = xrrc[a3, b5]; Clear[OUT]; OUT = r; t2 = dat[]; Save["out.xrrceb", ps, memo, OUT, s1, s2, t1, t2]; Return[r]; ] dat[] := Module[{r, t1, t2, t3, t4, t5, t6}, {t1, t2, t3, t4, t5, t6} = Date[]; r = {t3, t2, t1, t4, t5}; Return[r]; ] xrrc[rrc3_, std5_] := Module[{v1, md, acc, st1, st2, jp, rep, x, xa, ct}, {st1, st2, jp, rep, st3} = std5; {v1, md, acc} = rrc3; x = v1; ct = 0; Label[o1]; x = rrc[x, md, acc, st1]; ct = ct + 1; xa = x; x = rrc[x, md, acc, st2]; x = (1 + jp)*x - jp*xa; If[ct < rep, Goto[o1]]; x = rrc[x, md, acc, st3]; Return[x]; ] rrc[z1_, z2_, z3_, z4_] := rrb2[f, z1, z2, z3, z4] 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]; ] rf1ba[f_, n_Integer] := (AccuracyGoal -> n; PrecisionGoal -> n; WorkingPrecision -> n + 7; Module[{st22a}, st22a = {{D[f, y1], D[f, y2], D[f, y3], D[f, y4], D[f, y5], D[f, y6], D[f, y7], D[f, y8], D[f, y9], D[f, y10], D[f, y11], D[f, y12], D[f, y13], D[f, y14], D[f, y15], D[f, y16], D[f, y17], D[f, y18], D[f, y19], D[f, y20], D[f, y21]}, f}; Return[st22a]; ]) rat[y_, k_Integer] := Rationalize[y, 1/10^(k + 2)] 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; ]) sby[phi_, sya_] := Module[{yaa, r}, yaa = sb[sya]; r = phi /. yaa; Return[r]; Null; ] sb[wq_] := s21b[wq] s21b[w_] := Module[{j1, j2, j3, a, j4, j5, j6, j7, j8, j9, j10, j11, j12, j13, j14, j15, j16, j17, j18, j19, j20, j21}, {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 OUT = OUT xrrce[rrc3_, std5_] := Module[{r, t1, t2, s, p, q}, t1 = dat[]; r = xrrc[rrc3, std5]; Clear[STARTCLOCK, ENDCLOCK, OUT, SY00, std0a, rrc0a]; OUT = r; rrc0a = rrc3; std0a = std5; t2 = dat[]; {s, p, q} = rrc3; SY00 = s; STARTCLOCK = t1; ENDCLOCK = t2; Save["out.xrrce", STARTCLOCK, ps, memo, OUT, SY00, std0a, rrc0a, ENDCLOCK]; ]