New answers tagged

0

If this is for convenience, then you could write a wrapped for JordanDecomposition, e.g. something like this: Clear[jd] jd[m_] := Module[{s, j}, {s, j} = JordanDecomposition[m]; {s, Inverse[s], j} ] Then for example: a = {{27, 48, 81}, {-6, 0, 0}, {1, 0, 3}}; jd[a] (* Out: { (* the similarity matrix s *) {{3, 18, 2}, {-3, -9, -(1/4)}, {1, 2, 0}}, ...



Top 50 recent answers are included