Using lexicographic order compute a Gröbner basis for $$I=\langle f_1=xy^2-xy+y,f_2=xy-z^2,f_3=x-yz^4\rangle\subset \Bbb R[x,y,z]$$
I was strictly using these notes to compute a Gröbner basis.
When the computation gets messy I used a computer software at least to know the answer, but the software gave me
$$\{z^8-z^6-2 z^4-z^2, y-z^6+2 z^4, x-z^6+z^4+z^2\}$$
as a Gröbner basis, which I wouldn't even imagine, because the elements in $I$ are not included. I was expecting the answer to be like $\langle f_1,f_2,f_3,\ldots\rangle$.
The Gröbner basis depends on the monomial order you choose. If you use lexicographical order you do get (in M2; commented out -- is the answer):
But if you use graded reverse lexicographical order (default for M2):
In either case the ideals are the same, what needs to be true for a Gröbner basis (that it is a generating set and) that the initial terms generates the initial ideal for that monomial order.
This means that your original generators could be part of the list but that the initial terms aren't needed for generating the initial ideal.