FUNCTION: polysolve - find roots of a set of polynomials using Grobner bases

CALLING SEQUENCE:
   polysolve(f,x,optns)

PARAMETERS:
   f     - list of functions whose zeroes are to be found
   x     - list of unknown variables,
   optns - (optional) string or set of strings for the options,

SYNOPSIS:
- polysolve determines the roots of the polynomials using the
  Groebner Basis method.

- This method gives, in principle, all the roots of the polynomials
  (over the complex numbers), but in some cases, it can require a
  huge amount of computational work.

- If the string `verbose` is included in the options, the Groebner bases
  will be displayed.

EXAMPLES:
> polysolve([x+y,x*y-y+1], [x,y], verbose);

The Groebner bases are:
[[x+y, y^2+y-1]]

             [{y = RootOf(_Z^2+_Z-1), x = -RootOf(_Z^2+_Z-1)}],


SEE ALSO: groebner, grobner[gsolve]
