solution of system of quadratic eqns.



> Perhaps solve(), or algsys() can be made to recognize that the problem
> requires a call to_poly_solve() and compute accordingly.

Alternative:

load(grobner);

myalgsys(e,v) := (
  e : map(lambda([s], rhs(s)-lhs(s)),e),
  algsys(poly_reduced_grobner(e, block([listconstvars : false], listofvars(e))),v));

Well maybe this should check for polynomial equations before calling  poly_reduced_grobner. This code  a bit redundant: algys
doesn't know that the system is triangular.  I'm not sure that the author of the grobner package is interested in maintaining it.
Maybe it is bug free :)

> By the way, did you mean Maxima 5.27? 

Yes.

--Barton