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.

Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.