pub fn inv_u(u: Matrix) -> Matrix
Expand description

Inverse of upper triangular matrix

Examples

#[macro_use]
extern crate peroxide;
use peroxide::fuga::*;

fn main() {
    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));
}