
memo1 = {portion, of, this, file, below, and, including, this, advisory, is, 
     OK, for, input, to, a, session, of, MATHEMATICA}
 
eeseqe[m1_, m2_, m3_, m4_, m5_] := Module[{r, s}, 
     Clear[OUT, PSB00, STD00, TIMESTART, TIMEFINISH, EPS00, ACC00]; 
      TIMESTART = Date[]; r = eeseq[m1, m2, m3, m4, m5]; OUT = r; 
      PSB00 = m2; EPS00 = m3; STD00 = m4; ACC00 = m5; TIMEFINISH = Date[]; 
      Save["out.eeseqe", memo, cmemo, PSB00, EPS00, STD00, ACC00, 
       TIMESTART, TIMEFINISH, OUT]; Return[r]; ]
  
eeseq[systart_, psb0_, epsdata_, stepsdata_, accy_Integer] := 
    Module[{psba, sya, e3a, e4a, e5a, b, k1a, k2a, k3a, ja, sta, st2a, acc}, 
     psba = psb0; acc = accy; sya = systart; {e3a, e4a, e5a} = epsdata; 
      {b, k1a, k2a, k3a, ja, sta, st2a} = stepsdata; 
      Return[eeseq4[psba, e3a, e4a, e5a, b, k1a, k2a, k3a, ja, sta, st2a, 
        sya, acc]]; ]
 
eeseq4[psb_, e30_, e40_, e50_, b_, k1_, k2_, k3_, j_, st_Integer, 
     st2_Integer, sy0_, acc_Integer] := 
    Module[{ffav, ct, ct2, syva, e3v, e4v, e5v, r}, 
     ffav = ff0 /. psb; ct = 0; ct2 = 0; syva = sy0; e3v = e30; e4v = e40; 
      e5v = e50; r = {e30, e40, e50, sy0}; Label[o1]; 
      {e3v, e4v, e5v} = eestep4[b, e3v, e4v, e5v, k1, k2, k3, j]; 
      syva = systep[syva, e3v, e4v, e5v, ffav, acc]; ct = ct + 1; 
      ct2 = ct2 + 1; If[ct2 < st2, Goto[o1]]; ct2 = 0; 
      r = r + u^ct*{e3v, e4v, e5v, syva}; If[ct < st, Goto[o1]]; Return[r]; ]
 
eestep4[b_, c1_, c2_, c3_, j1_, j2_, j3_, m_] := 
    b*eestep3[c1/b, c2/b, c3/b, j1, j2, j3, m]
 
eestep3[a1_, a2_, a3_, k1_, k2_, k3_, j_] := 
    Module[{b1, b2, b3, s1, s2, s3, r}, b1 = 1/a1; b2 = 1/a2; b3 = 1/a3; 
      s1 = k1*b1^j; s1 = IntegerPart[s1]; s2 = k2*b2^j; s2 = IntegerPart[s2]; 
      s3 = k3*b3^j; s3 = IntegerPart[s3]; b1 = b1 + s1; b2 = b2 + s2; 
      b3 = b3 + s3; r = {1/b1, 1/b2, 1/b3}; Return[r]; ]
 
systep[insy_, e3a_, e4a_, e5a_, ffa_, acc_] := 
    Module[{syv}, ff = ffa /. {e3 -> e3a, e4 -> e4a, e5 -> e5a}; 
      f = sqsum[ff]; syv = rrs[insy, acc - 3, 2]; syv = N[syv, acc - 5]; 
      syv = FRTB[ff, syv, acc]; syv = N[syv, acc - 3]; Return[syv]; ]
 
steptest[epsdata_, stepsdata_] := Module[{bv, d1, d2, d3, k1, k2, k3, n, r, q1, q2}, 
     {d1, d2, d3} = epsdata; {bv, k1, k2, k3, n, q1, q2} = stepsdata; 
      r = eestep4[bv, d1, d2, d3, k1, k2, k3, n]; Return[r]; ]
 
