pub fn subst_bw<A: Copy + Field>(u: &Mat<A>, b: &Col<A>) -> Col<A>
Solve U x = b where U is an upper triangular matrix.