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

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

Print["Compute term3"];
y=Expand[Normal[Series[y,{t,0,5}]]];
x=Normal[Series[x,{t,0,5}]];
tmp=Expand[nc3[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 term5"];
y=Expand[Normal[Series[y,{t,0,3}]]];
x=Normal[Series[x,{t,0,3}]];
tmp=Expand[nc5[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 term7"];
y=Expand[Normal[Series[y,{t,0,1}]]];
x=Normal[Series[x,{t,0,1}]];
tmp=Expand[nc7[xx,yy]]; LL=Length[tmp]; term7=0;
Do[{tmp1=tmp[[i]]; term7=term7+termaction[tmp1]},{i,LL}];
term7=Expand[term7];
(* Save["b7rec",term7];  *)
