If G is a permutation group, I am aware of sage commands like G.subgroups(), G.normal_subgroups() and G.maximal_normal_subgroups() which returns the set of all subgroups , set of all normal subgroups and set of all maximal normal subgroups respectively.
Is there any such command for set of all maximal subgroups? Please help.
P.S. If there is no such command, can you help me to write a function which does it in an EFFICIENT way?
