Shared2017-12-31-011954.sagewsOpen in CoCalc
A.<X,Y> = PolynomialRing(QQ, order='lex')
f1 = X*Y
f2 = (X-Y)*X
J = A.ideal([f1,f2])
G=J.groebner_basis()
G
[X^2, X*Y]
var('x,y,z,u,v')
parametric_plot3d([u*v,u^2,1],(u,-2,2),(v,-1,1))
(x, y, z, u, v)
3D rendering not yet implemented