Poly would not work because the expression is not a polynomial! Take a look at the apart_list() function, which gives a structured format of the partial fraction decomposition.
Aaron Meurer On Thu, Jan 9, 2014 at 10:12 AM, Hans-Martin Tröger <[email protected]> wrote: > Hi, > i have given a Polynominal with > > a=((1.0*x + 1.0)*(x + 9.0)*(x + 81.0)*(x + 729.0))/((3.0*x + 1.0)*(3.0*x + > 9.0)*(3.0*x + 81.0)*(3.0*x + 729.0)) > > with > > the code > >>>> a=((1.0*x + 1.0)*(x + 9.0)*(x + 81.0)*(x + 729.0))/((3.0*x + 1.0)*(3.0*x >>>> + 9.0)*(3.0*x + 81.0)*(3.0*x + 729.0)) >>>> temp = sp.apart(a) >>>> print temp > 0.0123456790123457 + 0.728907627865961/(3.0*x + 1.0) + > 4.37544642857143/(1.0*x + 243.0) + 1.584375/(1.0*x + 27.0) + > 0.546180555555556/(1.0*x + 3.0) >>>> > > , i get the solution of the partial fraction decomposition. > > All I want now is to extract the poles and the coefficients. I tryed for 3 > hours to transform it into a real polynominal expression with > sympy.Poly([a],x), but all went wrong. > How can I do this an can I than extract the coefficients with > sympy.Poly.coeff ??? > > Thanks for help > Hans-Martin > > Also other solutions are very welcome > > -- > You received this message because you are subscribed to the Google Groups > "sympy" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to [email protected]. > To post to this group, send email to [email protected]. > Visit this group at http://groups.google.com/group/sympy. > For more options, visit https://groups.google.com/groups/opt_out. -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To post to this group, send email to [email protected]. Visit this group at http://groups.google.com/group/sympy. For more options, visit https://groups.google.com/groups/opt_out.
