The jordan command finds the Jordan form of a matrix.
| J=P−1AP |
jordan(A), in Maple mode, only returns the matrix J.
Input not in Maple mode:
| jordan([[4,1,1],[1,4,1],[1,1,4]]) |
|
| jordan([[4,1,1],[1,4,1],[1,1,4]],P) |
|
| P |
|
If A is symmetric and has eigenvalues with multiple orders, the matrix P returned by jordan(A) will contain orthogonal eigenvectors (not always of norm equal to 1); that is, PTP will be a diagonal matrix where the diagonal is the square norm of the eigenvectors.
| jordan([[4,1,1],[1,4,1],[1,1,4]]) |
|