In Mathematica 9, the function PrimitiveRoot[n] does not always return the SMALLEST primitive root of n for composite n that have a primitive root.
How does Mathematica select the particular primitive root that it does return?
