add function to compute orbit size under action of conjugation by permutation-x-transposition
add stability figure and combined stability regions figure