Computing all power integral bases in orders of totally real cyclic sextic number fields

An algorithm is given for determining all power integral bases in orders of totally real cyclic sextic number fields. The orders considered are in most cases the maximal orders of the fields. The corresponding index form equation is reduced to a relative Thue equation of degree 3 over the quadratic subfield and to some inhomogeneous Thue equations of degree 3 over the rationals. At the end of the paper, numerical examples are given.