# On function quo_rem

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.