arithmetic.clj: more work on matrix diagonalisation.