This can be done in a few steps in probably any computer algebra package. You take the generators of your original ideal I, and bi-homogenize them, as described in the question. Then saturate with respect to the two hyperplanes at infinity, which are defined by the equation x0y0.
For example, the diagonal in mathbbA3timesmathbbA3 is defined by x1−y1, x2−y2, and x3−y3. If I wanted to use this to compute the ideal of the diagonal in mathbbP3timesmathbbP3, I would use the following commands in Macaulay2:
r = QQ[x0,x1,x2,x3,y0,y1,y2,y3]
i = ideal(x1*y0-y1*x0, x2*y0-y2*x0, x3*y0-y3*x0)
saturate(i, x0*y0)
The code in Singular would be:
ring r = 0, (x0,x1,x2,x3,y0,y1,y2,y3), dp;
ideal i = x1*y0-y1*x0, x2*y0-y2*x0, x3*y0-y3*x0;
LIB "elim.lib";
sat(i, x0*y0);
No comments:
Post a Comment