/* Compute the class group for some metacyclic groups and compares the result
   with [CR, Theorem 50.25].
*/

pqr := [ [5,4,2], [7,3,2], [7,6,3], [11,5,4], [13,3,3],
         [13,4,8], [13,6,4],  [11,10,2], [13, 12, 2]
       ];

for v in pqr do
    if not TestMetacyclic(v[1], v[2], v[3])  then
        print "Katastrophe für pqr = ", v;
    else 
        print "pqr = ", v, " ist OK!!!!";
    end if;
    print "\n************************************************************\n";
end for;