pub fn hessenberg(matrix: &mut Array) -> Result<Map, Box<EvalAltResult>>
Calculates the QR decomposition of a matrix
let matrix = eye(5); let h_results = hessenberg(matrix); assert_eq(h_results, #{"h": eye(5), "q": eye(5)});