Maxima's facilities for manipulating equations and
solutions look rather weak. It would be nice if there was
a way for SOLVE to return a ROOTSOF object if it is not
able to solve an equation in term of radicals. One should
be able to manipulate ROOTSOF objects: for instance,
LSUM and LPRODUCT should be able to compute sums or
products over all roots of a polynomial.