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

t=1;
term=Expand[term1+term3+term5+term7+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]];

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

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

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

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",newa5]; Save["Coeff",newb5]; 

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

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

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

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

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

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