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

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

newa1=Coefficient[workterm, X[1]]; workterm=Expand[workterm-newa1  X[1]];

newa3=Coefficient[workterm, X[3]]; workterm=Expand[workterm-newa3  X[3]];

newb4=Coefficient[workterm,b4Cmt]; workterm=Expand[workterm-newb4 b4Cmt];

newa5=Coefficient[workterm, X[5]]; workterm=Expand[workterm-newa5  X[5]];
newb5=Coefficient[workterm,b5Cmt]; workterm=Expand[workterm-newb5 b5Cmt];

newb6=Coefficient[workterm,b6Cmt]; workterm=Expand[workterm-newb6 b6Cmt];
newc6=Coefficient[workterm,c6Cmt]; workterm=Expand[workterm-newc6 c6Cmt];

newa7=Coefficient[workterm, X[7]]; workterm=Expand[workterm-newa7  X[7]];
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];

newb8=Coefficient[workterm,b8Cmt]; workterm=Expand[workterm-newb8 b8Cmt];
newc8=Coefficient[workterm,c8Cmt]; workterm=Expand[workterm-newc8 c8Cmt];
newd8=Coefficient[workterm,d8Cmt]; workterm=Expand[workterm-newd8 d8Cmt];
newe8=Coefficient[workterm,e8Cmt]; workterm=Expand[workterm-newe8 e8Cmt];
newf8=Coefficient[workterm,f8Cmt]; workterm=Expand[workterm-newf8 f8Cmt];

newa9=Coefficient[workterm, X[9]]; workterm=Expand[workterm-newa9  X[9]];
newb9=Coefficient[workterm,b9Cmt]; workterm=Expand[workterm-newb9 b9Cmt];
newc9=Coefficient[workterm,c9Cmt]; workterm=Expand[workterm-newc9 c9Cmt];
newd9=Coefficient[workterm,d9Cmt]; workterm=Expand[workterm-newd9 d9Cmt];
newe9=Coefficient[workterm,e9Cmt]; workterm=Expand[workterm-newe9 e9Cmt];
newf9=Coefficient[workterm,f9Cmt]; workterm=Expand[workterm-newf9 f9Cmt];
newg9=Coefficient[workterm,g9Cmt]; workterm=Expand[workterm-newg9 g9Cmt];
newh9=Coefficient[workterm,h9Cmt]; workterm=Expand[workterm-newh9 h9Cmt];

newb10=Coefficient[workterm,b10Cmt]; workterm=Expand[workterm-newb10 b10Cmt];
newc10=Coefficient[workterm,c10Cmt]; workterm=Expand[workterm-newc10 c10Cmt];
newd10=Coefficient[workterm,d10Cmt]; workterm=Expand[workterm-newd10 d10Cmt];
newe10=Coefficient[workterm,e10Cmt]; workterm=Expand[workterm-newe10 e10Cmt];
newf10=Coefficient[workterm,f10Cmt]; workterm=Expand[workterm-newf10 f10Cmt];
newg10=Coefficient[workterm,g10Cmt]; workterm=Expand[workterm-newg10 g10Cmt];
newh10=Coefficient[workterm,h10Cmt]; workterm=Expand[workterm-newh10 h10Cmt];
newi10=Coefficient[workterm,i10Cmt]; workterm=Expand[workterm-newi10 i10Cmt];
newj10=Coefficient[workterm,j10Cmt]; workterm=Expand[workterm-newj10 j10Cmt];
newk10=Coefficient[workterm,k10Cmt]; workterm=Expand[workterm-newk10 k10Cmt];
newl10=Coefficient[workterm,l10Cmt]; workterm=Expand[workterm-newl10 l10Cmt];

newa11=Coefficient[workterm, X[11]]; workterm=Expand[workterm-newa11  X[11]];
newb11=Coefficient[workterm,b11Cmt]; workterm=Expand[workterm-newb11 b11Cmt];
newc11=Coefficient[workterm,c11Cmt]; workterm=Expand[workterm-newc11 c11Cmt];
newd11=Coefficient[workterm,d11Cmt]; workterm=Expand[workterm-newd11 d11Cmt];
newe11=Coefficient[workterm,e11Cmt]; workterm=Expand[workterm-newe11 e11Cmt];
newf11=Coefficient[workterm,f11Cmt]; workterm=Expand[workterm-newf11 f11Cmt];
newg11=Coefficient[workterm,g11Cmt]; workterm=Expand[workterm-newg11 g11Cmt];
newh11=Coefficient[workterm,h11Cmt]; workterm=Expand[workterm-newh11 h11Cmt];
newi11=Coefficient[workterm,i11Cmt]; workterm=Expand[workterm-newi11 i11Cmt];
newj11=Coefficient[workterm,j11Cmt]; workterm=Expand[workterm-newj11 j11Cmt];
newk11=Coefficient[workterm,k11Cmt]; workterm=Expand[workterm-newk11 k11Cmt];
newl11=Coefficient[workterm,l11Cmt]; workterm=Expand[workterm-newl11 l11Cmt];
newm11=Coefficient[workterm,m11Cmt]; workterm=Expand[workterm-newm11 m11Cmt];
newn11=Coefficient[workterm,n11Cmt]; workterm=Expand[workterm-newn11 n11Cmt];
newo11=Coefficient[workterm,o11Cmt]; workterm=Expand[workterm-newo11 o11Cmt];
newp11=Coefficient[workterm,p11Cmt]; workterm=Expand[workterm-newp11 p11Cmt];
newq11=Coefficient[workterm,q11Cmt]; workterm=Expand[workterm-newq11 q11Cmt];
newr11=Coefficient[workterm,r11Cmt]; workterm=Expand[workterm-newr11 r11Cmt];

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

Save["Coeff",newb4]; 

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

Save["Coeff",newb6]; Save["Coeff",newc6]; 

Save["Coeff",newa7]; Save["Coeff",newb7]; 
Save["Coeff",newc7]; Save["Coeff",newd7]; 

Save["Coeff",newb8]; Save["Coeff",newc8]; 
Save["Coeff",newd8]; Save["Coeff",newe8]; 
Save["Coeff",newf8]; 

Save["Coeff",newa9]; Save["Coeff",newb9]; 
Save["Coeff",newc9]; Save["Coeff",newd9]; 
Save["Coeff",newe9]; Save["Coeff",newf9]; 
Save["Coeff",newg9]; Save["Coeff",newh9];

Save["Coeff",newb10]; Save["Coeff",newc10]; 
Save["Coeff",newd10]; Save["Coeff",newe10]; 
Save["Coeff",newf10]; Save["Coeff",newg10]; 
Save["Coeff",newh10]; Save["Coeff",newi10]; 
Save["Coeff",newj10]; Save["Coeff",newk10]; 
Save["Coeff",newl10]; 

Save["Coeff",newa11]; Save["Coeff",newb11]; 
Save["Coeff",newc11]; Save["Coeff",newd11]; 
Save["Coeff",newe11]; Save["Coeff",newf11]; 
Save["Coeff",newg11]; Save["Coeff",newh11];
Save["Coeff",newi11]; Save["Coeff",newj11]; 
Save["Coeff",newk11]; Save["Coeff",newl11]; 
Save["Coeff",newm11]; Save["Coeff",newn11]; 
Save["Coeff",newo11]; Save["Coeff",newp11];
Save["Coeff",newq11]; Save["Coeff",newr11]; 

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

newa1p=Together[newa1]; Save["Coeffp",newa1p]; 
newa3p=Together[newa3]; Save["Coeffp",newa3p]; 

newb4p=Together[newb4]; Save["Coeffp",newb4p]; 

newa5p=Together[newa5]; Save["Coeffp",newa5p]; 
newb5p=Together[newb5]; Save["Coeffp",newb5p]; 

newb6p=Together[newb6]; Save["Coeffp",newb6p]; 
newc6p=Together[newc6]; Save["Coeffp",newc6p]; 

newa7p=Together[newa7]; Save["Coeffp",newa7p]; 
newb7p=Together[newb7]; Save["Coeffp",newb7p]; 
newc7p=Together[newc7]; Save["Coeffp",newc7p]; 
newd7p=Together[newd7]; Save["Coeffp",newd7p]; 

newb8p=Together[newb8]; Save["Coeffp",newb8p]; 
newc8p=Together[newc8]; Save["Coeffp",newc8p]; 
newd8p=Together[newd8]; Save["Coeffp",newd8p]; 
newe8p=Together[newe8]; Save["Coeffp",newe8p]; 
newf8p=Together[newf8]; Save["Coeffp",newf8p]; 

newa9p=Together[newa9]; Save["Coeffp",newa9p]; 
newb9p=Together[newb9]; Save["Coeffp",newb9p]; 
newc9p=Together[newc9]; Save["Coeffp",newc9p]; 
newd9p=Together[newd9]; Save["Coeffp",newd9p]; 
newe9p=Together[newe9]; Save["Coeffp",newe9p]; 
newf9p=Together[newf9]; Save["Coeffp",newf9p]; 
newg9p=Together[newg9]; Save["Coeffp",newg9p]; 
newh9p=Together[newh9]; Save["Coeffp",newh9p];

newb10p=Together[newb10]; Save["Coeffp",newb10p]; 
newc10p=Together[newc10]; Save["Coeffp",newc10p]; 
newd10p=Together[newd10]; Save["Coeffp",newd10p]; 
newe10p=Together[newe10]; Save["Coeffp",newe10p]; 
newf10p=Together[newf10]; Save["Coeffp",newf10p]; 
newg10p=Together[newg10]; Save["Coeffp",newg10p]; 
newh10p=Together[newh10]; Save["Coeffp",newh10p];
newi10p=Together[newi10]; Save["Coeffp",newi10p]; 
newj10p=Together[newj10]; Save["Coeffp",newj10p]; 
newk10p=Together[newk10]; Save["Coeffp",newk10p]; 
newl10p=Together[newl10]; Save["Coeffp",newl10p]; 

newa11p=Together[newa11]; Save["Coeffp",newa11p]; 
newb11p=Together[newb11]; Save["Coeffp",newb11p]; 
newc11p=Together[newc11]; Save["Coeffp",newc11p]; 
newd11p=Together[newd11]; Save["Coeffp",newd11p]; 
newe11p=Together[newe11]; Save["Coeffp",newe11p]; 
newf11p=Together[newf11]; Save["Coeffp",newf11p]; 
newg11p=Together[newg11]; Save["Coeffp",newg11p]; 
newh11p=Together[newh11]; Save["Coeffp",newh11p];
newi11p=Together[newi11]; Save["Coeffp",newi11p]; 
newj11p=Together[newj11]; Save["Coeffp",newj11p]; 
newk11p=Together[newk11]; Save["Coeffp",newk11p]; 
newl11p=Together[newl11]; Save["Coeffp",newl11p]; 
newm11p=Together[newm11]; Save["Coeffp",newm11p]; 
newn11p=Together[newn11]; Save["Coeffp",newn11p]; 
newo11p=Together[newo11]; Save["Coeffp",newo11p]; 
newp11p=Together[newp11]; Save["Coeffp",newp11p];
newq11p=Together[newq11]; Save["Coeffp",newq11p]; 
newr11p=Together[newr11]; Save["Coeffp",newr11p];
