Hmmm, This works: R.<x> = PolynomialRing(QQ,"x") a = x^2 +3*x b = x^3 - x^2 - x+1 c = a/b c.partial_fraction_decomposition() (0, [-1/2/(x + 1), (3/2*x + 1/2)/(x^2 - 2*x + 1)])
So, it likely has to do with the specific a and b since in the original, a/b works (with a remainder).
Note that this interface is not consistent with the partial fraction decomposition for symbolic expressions which is just f.partial_fraction(). I'd still like to be able to convert between the two.
Thanks, Tim. --- Tim Lahey PhD Candidate, Systems Design Engineering University of Waterloo http://www.linkedin.com/in/timlahey
smime.p7s
Description: S/MIME cryptographic signature