# Revision history [back]

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.