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

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

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

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

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

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

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

Save["Coeff",newa6]; Save["Coeff",newb6]; 
Save["Coeff",newc6]; Save["Coeff",newd6];
Save["Coeff",newe6]; 

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

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

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

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

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

