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

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

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

newa2=Coefficient[workterm, a2Cmt]; workterm=Expand[workterm-newa2 a2Cmt];

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

newa4=Coefficient[workterm, a4Cmt]; workterm=Expand[workterm-newa4 a4Cmt];
newb4=Coefficient[workterm, b4Cmt]; workterm=Expand[workterm-newb4 b4Cmt];
newc4=Coefficient[workterm, c4Cmt]; workterm=Expand[workterm-newc4 c4Cmt];

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];

newa6=Coefficient[workterm, a6Cmt]; workterm=Expand[workterm-newa6 a6Cmt];
newb6=Coefficient[workterm, b6Cmt]; workterm=Expand[workterm-newb6 b6Cmt];
newc6=Coefficient[workterm, c6Cmt]; workterm=Expand[workterm-newc6 c6Cmt];
newd6=Coefficient[workterm, d6Cmt]; workterm=Expand[workterm-newd6 d6Cmt];
newe6=Coefficient[workterm, e6Cmt]; workterm=Expand[workterm-newe6 e6Cmt];
newf6=Coefficient[workterm, f6Cmt]; workterm=Expand[workterm-newf6 f6Cmt];
newg6=Coefficient[workterm, g6Cmt]; workterm=Expand[workterm-newg6 g6Cmt];
newh6=Coefficient[workterm, h6Cmt]; workterm=Expand[workterm-newh6 h6Cmt];
newi6=Coefficient[workterm, i6Cmt]; workterm=Expand[workterm-newi6 i6Cmt];

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",newa2]; 

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

Save["Coeff",newa4]; Save["Coeff",newb4]; 
Save["Coeff",newc4];

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

Save["Coeff",newa6]; Save["Coeff",newb6]; 
Save["Coeff",newc6]; Save["Coeff",newd6];
Save["Coeff",newe6]; Save["Coeff",newf6];
Save["Coeff",newg6]; Save["Coeff",newh6];
Save["Coeff",newi6];

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]; 

newa2p=Together[newa2]; Save["Coeffp",newa2p]; 

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

newa4p=Together[newa4]; Save["Coeffp",newa4p]; 
newb4p=Together[newb4]; Save["Coeffp",newb4p]; 
newc4p=Together[newc4]; Save["Coeffp",newc4p]; 

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]; 

newa6p=Together[newa6]; Save["Coeffp",newa6p]; 
newb6p=Together[newb6]; Save["Coeffp",newb6p]; 
newc6p=Together[newc6]; Save["Coeffp",newc6p]; 
newd6p=Together[newd6]; Save["Coeffp",newd6p]; 
newe6p=Together[newe6]; Save["Coeffp",newe6p]; 
newf6p=Together[newf6]; Save["Coeffp",newf6p]; 
newg6p=Together[newg6]; Save["Coeffp",newg6p]; 
newh6p=Together[newh6]; Save["Coeffp",newh6p]; 
newi6p=Together[newi6]; Save["Coeffp",newi6p]; 

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]; 

