pub fn inv_u(u: Matrix) -> Matrix
Inverse of upper triangular matrix
extern crate peroxide; use peroxide::*; let u = matrix(c!(2,2,0,1), 2, 2, Row); assert_eq!(inv_u(u), matrix(c!(0.5,-1,0,1), 2, 2, Row));