memo2 = {this,is,the,corrected,programming,of,2002d03}

memo3 = {I,found,that,a,more,complicated,program,form,existed,allowing,
     for,a,modified,steps,procedure,but,it,need,not,be,used,if,one,
     just,uses,1,for,the,first,input,argument,of,the,function,called,
     $$eeseq4$$,and,so,I,saved,that,structure,in,case,of,occasional,use}

memo1 = {portion, of, this, file, below, and, including, this, advisory, is, 
     OK, for, input, to, a, session, of, MATHEMATICA}
 
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]; ]