use mathru::algebra::linear::matrix::{General, HessenbergDecomposition, UpperHessenberg};
use mathru::matrix;
fn main() {
let a: General<f64> = matrix![ 1.0, 5.0, 3.0;
1.0, 0.0, -7.0;
3.0, 8.0, 9.0];
let (_q, _h): (General<f64>, UpperHessenberg<f64>) = a.dec_hessenberg().qh();
}