Not good for teaching (and I just got a complain for that). Of course I can do map(lambda x: x.radical_expression(), M.eigenvalues()) or something similar. Or begin with Matrix(SR, ...).
But is there a default setting to change output type of AlgebraicNumber? -- Jori Mäntysalo