use nalgebra::{DMatrix};
#[macro_use]
extern crate ndarray;
//use ndarray::Array;
use ndarray::{Array,ArrayD, Dim, Ix, Ix2, IxDyn};


fn main() {
    let xdata = [0.05,  0.18,0.31, 0.42, 0.5];
    let ydata = [0.12, 0.22,  0.35, 0.38, 0.49];
    let dmaty = DMatrix::from_iterator(5,1,ydata.iter().cloned());
    let dmat = DMatrix::from_iterator(5,1,xdata.iter().cloned());
    let qr = dmat.qr();
    let q = qr.q();
    let r = qr.r();
    println!("{}  {}",q,r);
    let rinv = r.try_inverse().unwrap();
    println!("{}",rinv.dot(&dmaty.transpose()));
}