rrfrce[v1_, v2_, v3_, v4_] := Module[{r, s, t1, t2}, t1 = Date[]; r = rrfrc[v1, v2, v3, v4]; Clear[OUT]; OUT = r; t2 = Date[]; 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; ])