/**/ use QQ[x[1..8],y[1..8]];
/**/ HPPI8 := [x[1]^I*x[I+2]*y[2]^(I+1) -y[1]^I*y[I+2]*x[2]^(I+1) | I in 1..6];
/**/ BL := toric(ideal(HPPI8), [x[1],y[2]]);
/**/ len(gens(BL));
340
/**/ use QQ[x[1..9],y[1..9]];
/**/ HPPI9 := [x[1]^I*x[I+2]*y[2]^(I+1) -y[1]^I*y[I+2]*x[2]^(I+1) | I in 1..7];
/**/ BL := toric(ideal(HPPI9), [x[1],y[2]]);
/**/ len(gens(BL));
798
/**/ use R ::= QQ[x,y,z,w];
/**/ toric(ideal(x*z-y^2, x*w-y*z));
ideal(-y^2 +x*z, -y*z +x*w, z^2 -y*w)
/**/ toric(ideal(x*z-y^2, x*w-y*z), [y]);
ideal(-y^2 +x*z, -y*z +x*w, z^2 -y*w)
/**/ use R ::= QQ[x,y,z];
/**/ toric([[1,3,2],[3,4,8]]);
ideal(-x^16 +y^2*z^5)
/**/ toric(mat([[1,3,2],[3,4,8]]));
ideal(-x^16 +y^2*z^5)
/**/ toric(mat([[1,3,2],[3,-4,8]]));
ideal(-x^32 +y^2*z^13)
|