(*                            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];
newe5=Coefficient[workterm, e5Cmt]; workterm=Expand[workterm-newe5 e5Cmt];
newf5=Coefficient[workterm, f5Cmt]; workterm=Expand[workterm-newf5 f5Cmt];

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];
newk7=Coefficient[workterm, k7Cmt]; workterm=Expand[workterm-newk7 k7Cmt];
newl7=Coefficient[workterm, l7Cmt]; workterm=Expand[workterm-newl7 l7Cmt];
newm7=Coefficient[workterm, m7Cmt]; workterm=Expand[workterm-newm7 m7Cmt];
newn7=Coefficient[workterm, n7Cmt]; workterm=Expand[workterm-newn7 n7Cmt];
newo7=Coefficient[workterm, o7Cmt]; workterm=Expand[workterm-newo7 o7Cmt];
newp7=Coefficient[workterm, p7Cmt]; workterm=Expand[workterm-newp7 p7Cmt];
newq7=Coefficient[workterm, q7Cmt]; workterm=Expand[workterm-newq7 q7Cmt];
newr7=Coefficient[workterm, r7Cmt]; workterm=Expand[workterm-newr7 r7Cmt];

(* 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",newe5]; Save["Coeff",newf5];

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];
Save["Coeff",newk7]; Save["Coeff",newl7];
Save["Coeff",newm7]; Save["Coeff",newn7];
Save["Coeff",newo7]; Save["Coeff",newp7];
Save["Coeff",newq7]; Save["Coeff",newr7];

(* 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]; 
newe5p=Together[newe5]; Save["Coeffp",newe5p]; 
newf5p=Together[newf5]; Save["Coeffp",newf5p]; 

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]; 
newk7p=Together[newk7]; Save["Coeffp",newk7p]; 
newl7p=Together[newl7]; Save["Coeffp",newl7p]; 
newm7p=Together[newm7]; Save["Coeffp",newm7p]; 
newn7p=Together[newn7]; Save["Coeffp",newn7p]; 
newo7p=Together[newo7]; Save["Coeffp",newo7p]; 
newp7p=Together[newp7]; Save["Coeffp",newp7p]; 
newq7p=Together[newq7]; Save["Coeffp",newq7p]; 
newr7p=Together[newr7]; Save["Coeffp",newr7p]; 

