(*                         Ren-Cang Li, June 1, 1996
                           na.rcli@na-net.ornl.gov   *)

Print["Compute term1"];
term1=Expand[c[1][x,y]];
(* Save["b1rec",term1]; *)

Print["Compute term2"];
y=Expand[Normal[Series[y,{t,0,6}]]];
x=Normal[Series[x,{t,0,6}]];
term2=Expand[Normal[Series[c[2][x,y],{t,0,7}]]];
(* Save["b2rec",term2];  *)

Print["Compute term3"];
y=Expand[Normal[Series[y,{t,0,5}]]];
x=Normal[Series[x,{t,0,5}]];
tmp=Expand[c[3][xx,yy]]; LL=Length[tmp]; term3=0;
Do[{tmp1=tmp[[i]]; term3=term3+termaction[tmp1]},{i,LL}];
term3=Expand[term3];
(* Save["b3rec",term3];  *)

Print["Compute term4"];
y=Expand[Normal[Series[y,{t,0,4}]]];
x=Normal[Series[x,{t,0,4}]];
term4=Expand[termaction[c[4][xx,yy]]];
(* Save["b4rec",term4];  *)

Print["Compute term5"];
y=Expand[Normal[Series[y,{t,0,3}]]];
x=Normal[Series[x,{t,0,3}]];
tmp=Expand[c[5][xx,yy]]; LL=Length[tmp]; term5=0;
Do[{tmp1=tmp[[i]]; term5=term5+termaction[tmp1]},{i,LL}];
term5=Expand[term5];
(* Save["b5rec",term5];  *)

Print["Compute term6"];
y=Expand[Normal[Series[y,{t,0,2}]]];
x=Normal[Series[x,{t,0,2}]];
tmp=Expand[c[6][xx,yy]]; LL=Length[tmp]; term6=0;
Do[{tmp1=tmp[[i]]; term6=term6+termaction[tmp1]},{i,LL}];
term6=Expand[term6];
(* Save["b6rec",term6];  *)

Print["Compute term7"];
y=Expand[Normal[Series[y,{t,0,1}]]];
x=Normal[Series[x,{t,0,1}]];
tmp=Expand[c[7][xx,yy]]; LL=Length[tmp]; term7=0;
Do[{tmp1=tmp[[i]]; term7=term7+termaction[tmp1]},{i,LL}];
term7=Expand[term7];
(* Save["b7rec",term7];  *)
