G := Sym(11); involutions := [ [ G!(1, 8)(2, 10)(3, 5)(6, 9)(7, 11), G!(1, 3)(2, 9)(4, 6)(5, 10)(7, 11), G!(1, 4)(2, 11)(5, 10)(6, 9)(7, 8) ] , [ G!(1, 9)(2, 8)(3, 11)(4, 7)(5, 10), G!(1, 11)(2, 3)(4, 10)(5, 8)(6, 9), G!(1, 10)(2, 4)(3, 6)(7, 11)(8, 9) ] , [ G!(1, 11)(2, 10)(4, 9)(5, 7)(6, 8), G!(1, 10)(2, 9)(3, 7)(5, 6)(8, 11), G!(1, 9)(2, 11)(3, 10)(4, 5)(6, 7) ] , [ G!(1, 3)(2, 9)(4, 5)(6, 7)(8, 11), G!(1, 4)(2, 11)(3, 10)(5, 7)(6, 8), G!(1, 8)(2, 10)(3, 7)(4, 9)(5, 6) ] , [ G!(1, 11)(2, 8)(3, 10)(4, 7)(5, 6), G!(1, 10)(2, 3)(4, 9)(5, 8)(6, 7), G!(1, 9)(2, 4)(3, 6)(5, 7)(8, 11) ] , [ G!(1, 3)(2, 7)(5, 11)(6, 9)(8, 10), G!(1, 4)(2, 5)(3, 9)(6, 10)(7, 11), G!(1, 8)(2, 6)(4, 11)(5, 10)(7, 9) ] ]; geometries := [* *]; for i in [1..#involutions] do G1 := sub; G2 := sub; G3 := sub; Append(~geometries,CosetGeometry(G,{G1,G2,G3})); end for;