MathGroup Archive 2004

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Extracting Coefficients and Powers


Bruce W. Colletti wrote:
> I have a sum of terms, each of the form "a * x^r" or "a * 1/x^r", where 
> a and (positive) r are real numbers (if a is missing, it's presumably 1 
> or -1 and if there is no x-term, r is understood to be 0). 
> 
> In my application, the complete expression is given by Expand[p], where 
> p is a product of terms having form a t^b + c t^d (p is actually a 
> factorial moment generating function).
> 
> I want to extract all coefficients into a list (including 1 or -1 
> coefficients), and exponents into another (to include 0 for any constant 
> term).  For instance,
> 
> x^3.2 + 5 - 3 / x^5.2 - x^5
> 
> yields the coefficients' list {1, 5, -3, -1} and exponents' list {3.2, 
> 0, -5.2,5}.
> 
> How would I build these lists?  Must Cases[ ] be used or are there 
> built-in functions?
> 
> Thanks.
> 
> Bruce
> 

If you are willing to work with rational exponents you can do as below 
to get exponent/coefficient list pairs. These can be further processed 
in an abvious way to get the separate lists you seek.

In[15]:= InputForm[ipoly = Internal`NestedTermsList[Rationalize[x^3.2 +
   5 - 3 / x^5.2 - x^5], x]]
Out[15]//InputForm= {{5, -1}, {16/5, 1}, {0, 5}, {-26/5, -3}}

In[16]:= InputForm[poly = Internal`FromNestedTermsList[ipoly, x]]
Out[16]//InputForm= 5 - 3/x^(26/5) + x^(16/5) - x^5


Daniel Lichtblau
Wolfram Research


  • Prev by Date: Re: Uniform design
  • Next by Date: Re: Trouble : simple stuffs for beginner
  • Previous by thread: Extracting Coefficients and Powers
  • Next by thread: Re: Extracting Coefficients and Powers