autoDn:= function(G)
local a,b,aims,bims,maps,autos,abims;
a:=G.1; b:=G.2;
aims:=Filtered(Elements(G), i -> Order(a) = Order(i));
bims:=Filtered(Elements(G), i -> Order(b) = Order(i));
abims:= Cartesian(aims,bims);
maps:= List(abims, i -> GroupHomomorphismByImages(G,G,[a,b],i));
maps:= Filtered(maps, i -> i <> fail);
autos:= Filtered(maps, IsInjective);
return autos;
end;