A.<X,Y> = PolynomialRing(QQ, order='lex')
f1 = X*Y
f2 = (X-Y)*X
J = A.ideal([f1,f2])
G=J.groebner_basis()
G

var('x,y,z,u,v')
parametric_plot3d([u*v,u^2,1],(u,-2,2),(v,-1,1))