pub fn subst_fw<A: Copy + Field>(l: &Mat<A>, b: &Col<A>) -> Col<A>
Solve L x = b where L is a lower triangular matrix.