1 | initial version |

Looks like a bug. Here is an alternative method, for polynomial equations with a finite number of solutions:

```
sage: R.<x,y> = PolynomialRing(QQ)
sage: I = R.ideal([x*x*x-y*y-10.5, 3.0*x*y+y-4.6])
sage: I.variety(RR)
[{y: 0.601783026716651, x: 2.21465035058553}]
```

You might also like `AA`

(the field of real algebraic numbers) instead of `RR`

, for exact computations.

2 | No.2 Revision |

Looks like a ~~bug. ~~bug in sympy. Here is an alternative method, for polynomial equations with a finite number of solutions:

```
sage: R.<x,y> = PolynomialRing(QQ)
sage: I = R.ideal([x*x*x-y*y-10.5, 3.0*x*y+y-4.6])
sage: I.variety(RR)
[{y: 0.601783026716651, x: 2.21465035058553}]
```

You might also like `AA`

(the field of real algebraic numbers) instead of `RR`

, for exact computations.

Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.