Hi, I write the following code:
R.<x,y,z> = PolynomialRing(GF(29), order='degrevlex');
f = x*y + z;
g = x + y;
f.quo_rem(g)
and the result is
(x, -x^2 + z).
Is it a really correct result of the function quo_rem?
Because of the 'degrevlex' (x > y > z and y^2 < x*y < x^2), I expect something like this
(y, -y^2 + z).
Thank you very much in advance for your help.
