FUNCTION: uf_roots - compute real roots for a unit form

CALLING SEQUENCE:
uf_roots(q)
uf_roots(q,m)

PARAMETERS:
q - matrix of the unit form
m - optional parameter, maximal value for root entries. Default value is 6.

SYNOPSIS:

• The function computes all real roots (with entries up to m) performing sequences of reflections on base vectors.

• The function returns a list of roots as a vector, in which the entries are roots.

EXAMPLES:

> q:=uf_init([2,[1,2,-2]]);

> uf_roots(q);

> uf_roots(q,2);

> r:=uf_roots(q,8);

> r[vectdim(r)];