Revision history [back]
 | 1 | initial version |
Try this:
sage: g = f.polynomial(y)
sage: g.change_ring(g.base_ring().fraction_field()).factor()
(9) * (y^2 + 16*x^5) * (y^4 + 16*x^2*y^3 + 15*x^3*y^2 + x^6)
I'm not sure whether that's any better than your original
sage: f.factor()
(-8) * (-x^5 + y^2) * (x^6 - 2*x^3*y^2 - x^2*y^3 + y^4)
since as John already indicated, there is no fundamental difference.
