x_o := MaximalOrder(X^6+X^5+X^4+X^3+X^2+X+1); x_u:=ExceptionalUnits(x_o); ExceptionalUnitOrbit(x_u[2]);