groebner basis[{x^3 - p, y-( x^2 + x)}, {x,y,p}]