Revision history [back]

One point from infinity variety

Given an ideal, I want to compute a non-zero point in the variety of that ideal. This can be done with .variety() if the ideal has 0-dimension. How can I get a non-zero point in the variety if the ideal is not 0-dimensional?

Example:

R.<x,y> = PolynomialRing(QQ)
I = R.ideal([x^2-y^2])
I.variety()


yields an error. How can I get a point like (1,1) which lies in the variety?

Thanks!