groupact:= function(a,b)
local f,g,H,h,w,i;
g:= Group(a,b);
H:= Elements(Group(a,b));
h:= Length(H);
w:= [];
i:=1;
  repeat
     if a=a^H[i] and b=b^H[i] then
       Append(w, [H[i]]);
     fi;
     i:= i+1;
  until i = h+1;
return w;
end;