pub fn matvec(a: &[Vec<Rational>], x: &[Rational]) -> Vec<Rational>
A·x over the rationals, computed exactly.