gap> G := JohnsonScheme(5,2);; gap> IsPPolynomialScheme(G); true gap> AllPPolynomialOrdering(G); [ [ 1, 2 ], [ 2, 1 ] ] gap> PrintArray(IntersectionArray(G, 1)); [ [ *, 1, 4 ], [ 0, 3, 2 ], [ 6, 2, * ] ]