pub fn characteristic_polynomial(m: &[Vec<Rational>]) -> Option<Poly>
det(x·I − M) by Faddeev–LeVerrier; None for a matrix that is not square.
det(x·I − M)
None