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!
