1 | initial version |

I don't think solve does what you hope it does. It seems from your code that you want to select elements from a groebner basis of bounded degree (I guess n is the length of G?) and want to find the solutions to those equations.

If the equations you select still define a 0-dimensional variety and you're interested in solutions that take values in the ground field, you might be able to do something like

```
J=ideal([f for f in I.groebner_basis() if f.degree() lt num])
J.variety()
```

Putting your selected polynomials in a dictionary is probably not a good idea. You'd be better off putting them in a list, i.e.

```
P=[]
...
P.append(G[i])
...
```

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.