/**/ 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
)
