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

t=1;
term=Expand[term1+term3+term5+term7];
workterm=term;

newa1=Coefficient[workterm,     X]; workterm=Expand[workterm-newa1     X];
newb1=Coefficient[workterm,     Y]; workterm=Expand[workterm-newb1     Y];

newa3=Coefficient[workterm, a3Cmt]; workterm=Expand[workterm-newa3 a3Cmt];
newb3=Coefficient[workterm, b3Cmt]; workterm=Expand[workterm-newb3 b3Cmt];

newa5=Coefficient[workterm, a5Cmt]; workterm=Expand[workterm-newa5 a5Cmt];
newb5=Coefficient[workterm, b5Cmt]; workterm=Expand[workterm-newb5 b5Cmt];
newc5=Coefficient[workterm, c5Cmt]; workterm=Expand[workterm-newc5 c5Cmt];
newd5=Coefficient[workterm, d5Cmt]; workterm=Expand[workterm-newd5 d5Cmt];

newa7=Coefficient[workterm, a7Cmt]; workterm=Expand[workterm-newa7 a7Cmt];
newb7=Coefficient[workterm, b7Cmt]; workterm=Expand[workterm-newb7 b7Cmt];
newc7=Coefficient[workterm, c7Cmt]; workterm=Expand[workterm-newc7 c7Cmt];
newd7=Coefficient[workterm, d7Cmt]; workterm=Expand[workterm-newd7 d7Cmt];
newe7=Coefficient[workterm, e7Cmt]; workterm=Expand[workterm-newe7 e7Cmt];
newf7=Coefficient[workterm, f7Cmt]; workterm=Expand[workterm-newf7 f7Cmt];
newg7=Coefficient[workterm, g7Cmt]; workterm=Expand[workterm-newg7 g7Cmt];
newh7=Coefficient[workterm, h7Cmt]; workterm=Expand[workterm-newh7 h7Cmt];
newi7=Coefficient[workterm, i7Cmt]; workterm=Expand[workterm-newi7 i7Cmt];
newj7=Coefficient[workterm, j7Cmt]; workterm=Expand[workterm-newj7 j7Cmt];

(* Print["workterm is"]; Print[workterm]; *)
If[workterm==0,
   Print["Recursion formulas are correctly computed."],
   Print["Something is wrong. The computed formulas can not be used."]];

(* Dump Coefficients *)
Print["Dump Coefficients"];

Save["Coeff",newa1]; Save["Coeff",newb1]; 

Save["Coeff",newa3]; Save["Coeff",newb3]; 

Save["Coeff",newa5]; Save["Coeff",newb5]; 
Save["Coeff",newc5]; Save["Coeff",newd5];

Save["Coeff",newa7]; Save["Coeff",newb7]; 
Save["Coeff",newc7]; Save["Coeff",newd7];
Save["Coeff",newe7]; Save["Coeff",newf7];
Save["Coeff",newg7]; Save["Coeff",newh7];
Save["Coeff",newi7]; Save["Coeff",newj7];

(* Another Way to dump them *)
Print["Another Way to dump them"];

newa1p=Together[newa1]; Save["Coeffp",newa1p]; 
newb1p=Together[newb1]; Save["Coeffp",newb1p]; 

newa3p=Together[newa3]; Save["Coeffp",newa3p]; 
newb3p=Together[newb3]; Save["Coeffp",newb3p]; 

newa5p=Together[newa5]; Save["Coeffp",newa5p]; 
newb5p=Together[newb5]; Save["Coeffp",newb5p]; 
newc5p=Together[newc5]; Save["Coeffp",newc5p]; 
newd5p=Together[newd5]; Save["Coeffp",newd5p]; 

newa7p=Together[newa7]; Save["Coeffp",newa7p]; 
newb7p=Together[newb7]; Save["Coeffp",newb7p]; 
newc7p=Together[newc7]; Save["Coeffp",newc7p]; 
newd7p=Together[newd7]; Save["Coeffp",newd7p]; 
newe7p=Together[newe7]; Save["Coeffp",newe7p]; 
newf7p=Together[newf7]; Save["Coeffp",newf7p]; 
newg7p=Together[newg7]; Save["Coeffp",newg7p]; 
newh7p=Together[newh7]; Save["Coeffp",newh7p]; 
newi7p=Together[newi7]; Save["Coeffp",newi7p]; 
newj7p=Together[newj7]; Save["Coeffp",newj7p]; 

