dev:=function(g,bb)
return Union(List(bb,x->Orbit(g,x,OnSets)));
end;

g:=Group((4,5,6)(7,8,9)(10,11,12)(13,14,15),(1,4)(2,5)(3,6)(8,13)(9,10)(11,15));

bb2:=[[1, 2, 3, 4, 5], [1, 2, 3, 7, 11], [1, 2, 4, 7, 11], [1, 2, 4, 7, 14],
[1, 2, 4, 8, 10], [1, 2, 4, 8, 11], [1, 2, 4, 9, 12], [1, 2, 4, 9, 15],
[1, 2, 4, 10, 13], [1, 2, 4, 12, 14], [1, 2, 4, 13, 15], [1, 2, 7, 8, 9],
[1, 2, 7, 8, 14], [1, 2, 7, 12, 15], [1, 2, 10, 11, 12], [1, 2, 10, 11, 15],
[1, 2, 13, 14, 15], [1, 4, 7, 8, 12], [1, 4, 7, 9, 11], [1, 4, 8, 9, 14],
[1, 4, 8, 12, 13], [1, 4, 8, 13, 14], [1, 4, 9, 10, 11], [1, 4, 11, 12, 15],
[1, 4, 11, 14, 15], [1, 7, 8, 10, 11], [1, 7, 8, 10, 13], [1, 7, 8, 11, 15],
[1, 7, 8, 14, 15], [1, 7, 10, 11, 14], [1, 7, 10, 13, 15], [1, 7, 11, 12, 14],
[1, 7, 12, 13, 14], [1, 10, 11, 13, 15], [7, 8, 9, 10, 13], [7, 8, 11, 12, 13]];

bb3:=[[1, 2, 3, 4, 8], [1, 2, 3, 7, 10], [1, 2, 4, 5, 9], [1, 2, 4, 5, 14], [1, 2, 4, 7, 8], 
[1, 2, 4, 7, 11], [1, 2, 4, 7, 15], [1, 2, 4, 9, 12], [1, 2, 4, 9, 13], [1, 2, 4, 10, 12], 
[1, 2, 4, 10, 15], [1, 2, 4, 11, 13], [1, 2, 4, 14, 15], [1, 2, 7, 8, 11], [1, 2, 7, 8, 14], 
[1, 2, 7, 11, 12], [1, 2, 7, 12, 13], [1, 2, 7, 12, 15], [1, 2, 7, 14, 15], [1, 2, 10, 11, 14], 
[1, 2, 10, 14, 15], [1, 4, 7, 8, 9], [1, 4, 7, 8, 10], [1, 4, 7, 9, 15], [1, 4, 7, 12, 14], 
[1, 4, 8, 9, 14], [1, 4, 8, 11, 12], [1, 4, 8, 11, 13], [1, 4, 8, 12, 13], [1, 4, 8, 12, 14], 
[1, 4, 8, 14, 15], [1, 4, 9, 10, 14], [1, 4, 9, 11, 12], [1, 4, 9, 11, 15], [1, 4, 11, 14, 15], 
[1, 7, 8, 10, 13], [1, 7, 8, 10, 14], [1, 7, 8, 11, 12], [1, 7, 8, 11, 15], [1, 7, 8, 12, 15], 
[1, 7, 10, 12, 14], [1, 7, 10, 13, 15], [1, 7, 11, 12, 14], [1, 7, 11, 13, 14], [1, 7, 13, 14, 15], 
[1, 10, 11, 12, 13], [7, 8, 9, 13, 14], [7, 8, 10, 11, 15], [7, 8, 11, 12, 13]];


bb5:=[[1, 2, 3, 4, 5], [1, 2, 3, 4, 8], [1, 2, 3, 7, 8], [1, 2, 3, 7, 11], [1, 2, 4, 5, 9], 
[1, 2, 4, 5, 14], [1, 2, 4, 7, 9], [1, 2, 4, 7, 12], [1, 2, 4, 7, 13], [1, 2, 4, 7, 14], 
[1, 2, 4, 7, 15], [1, 2, 4, 8, 11], [1, 2, 4, 8, 12], [1, 2, 4, 8, 15], [1, 2, 4, 9, 11], 
[1, 2, 4, 9, 14], [1, 2, 4, 9, 15], [1, 2, 4, 10, 11], [1, 2, 4, 10, 12], [1, 2, 4, 10, 13], 
[1, 2, 4, 10, 14], [1, 2, 4, 11, 13], [1, 2, 4, 12, 15], [1, 2, 4, 13, 15], [1, 2, 7, 8, 9], 
[1, 2, 7, 8, 11], [1, 2, 7, 8, 12], [1, 2, 7, 10, 13], [1, 2, 7, 10, 14], [1, 2, 7, 10, 15], 
[1, 2, 7, 11, 14], [1, 2, 7, 12, 13], [1, 2, 7, 12, 14], [1, 2, 7, 13, 15], [1, 2, 10, 11, 12], 
[1, 2, 10, 11, 13], [1, 2, 10, 14, 15], [1, 2, 13, 14, 15], [1, 4, 7, 8, 10], [1, 4, 7, 8, 13], 
[1, 4, 7, 8, 15], [1, 4, 7, 9, 15], [1, 4, 7, 11, 12], [1, 4, 7, 11, 15], [1, 4, 7, 12, 14], 
[1, 4, 8, 9, 10], [1, 4, 8, 9, 11], [1, 4, 8, 9, 14], [1, 4, 8, 9, 15], [1, 4, 8, 10, 12], 
[1, 4, 8, 10, 14], [1, 4, 8, 11, 12], [1, 4, 8, 12, 13], [1, 4, 8, 13, 14], [1, 4, 8, 14, 15], 
[1, 4, 9, 10, 12], [1, 4, 9, 11, 14], [1, 4, 11, 12, 14], [1, 4, 11, 12, 15], [1, 4, 11, 14, 15], 
[1, 7, 8, 9, 13], [1, 7, 8, 10, 11], [1, 7, 8, 10, 12], [1, 7, 8, 10, 13], [1, 7, 8, 11, 12], 
[1, 7, 8, 11, 15], [1, 7, 8, 12, 14], [1, 7, 8, 14, 15], [1, 7, 10, 12, 13], [1, 7, 10, 12, 14], 
[1, 7, 10, 12, 15], [1, 7, 10, 13, 15], [1, 7, 10, 14, 15], [1, 7, 11, 12, 13], [1, 7, 11, 13, 14], 
[1, 7, 11, 13, 15], [1, 7, 11, 14, 15], [1, 7, 13, 14, 15], [1, 10, 11, 12, 13], [1, 10, 11, 14, 15], 
[7, 8, 9, 10, 15], [7, 8, 9, 13, 14], [7, 8, 10, 11, 15], [7, 8, 11, 12, 13]];

d2:=dev(g,bb2);

d3:=dev(g,bb3);

d5:=dev(g,bb5);

g1:=Group((1,13,12)(2,10,8)(3,7,4)(5,14,11)(6,9,15),(1,4,8,14,3)(2,5,12,6,13)(7,9,10,11,15));

bb4:=[[1, 2, 3, 4, 6], [1, 2, 3, 4, 7], [1, 2, 3, 4, 12], [1, 2, 3, 4, 14], [1, 2, 3, 7, 8], 
[1, 2, 4, 6, 8], [1, 2, 4, 7, 10], [1, 2, 4, 8, 10], [1, 2, 4, 8, 13], [1, 2, 4, 10, 14], 
[1, 2, 4, 11, 12], [1, 2, 4, 11, 15], [1, 2, 4, 14, 15], [1, 2, 5, 7, 12], [1, 2, 5, 8, 12], 
[1, 2, 5, 8, 14], [1, 2, 5, 8, 15], [1, 2, 5, 9, 15], [1, 2, 5, 12, 13], [1, 2, 6, 7, 12], 
[1, 2, 6, 8, 13], [1, 2, 6, 10, 12], [1, 2, 6, 12, 13], [1, 2, 6, 12, 15], [1, 2, 7, 8, 15]];

d4:=dev(g1,bb4);