1 | initial version |

Not sure what `x^y^2`

was supposed to be in your polynomial `f`

.
Anyway, try something like this.

```
sage: R.<x,y> = PolynomialRing(ZZ, order='deglex')
sage: f = x^5 * y^2 + x^3 * y^2 - y + 1
sage: F = (x * y^2 - x, x - y^3)
sage: f.reduce(F)
x^5 + x^3 - y + 1
```

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.