/**/ use P ::= QQ[x,y];
/**/ Points := mat([[1, 2], [3, 4], [5, 6]]);
/**/ I := IdealOfPoints(P, Points);
/**/ I;
ideal(x -y +1, y^3 -12*y^2 +44*y -48)
/**/ K := NewFractionField(NewPolyRing(QQ, "a"));
/**/ use K;
/**/ Points := mat([[1,2,0], [3,4,a], [5,1,6]]);
/**/ use P ::= K[x,y,z], Lex;
/**/ I := IdealOfPoints(P, Points);
/**/ indent(I);
ideal(
z^3 +(-a -6)*z^2 +(6*a)*z,
y +((-a -12)/(6*a^2 -36*a))*z^2 +((a^2 +72)/(6*a^2 -36*a))*z -2,
x +((2*a -6)/(3*a^2 -18*a))*z^2 +((-2*a^2 +36)/(3*a^2 -18*a))*z -1
)
|