PRump( G, p )
PRump
returns the p-rump of an ag group G for a prime p.
The p-rump of a group G is the normal closure under G of the subgroup generated by the commutators and p.th powers of the generators of G.
gap> PRump( s4, 2 ); Subgroup( s4, [ b, c, d ] ) gap> PRump( s4, 3 ); s4
GAP 3.4.4