pub fn lstsq_ndarray<T>(a: &Array2<T>, b: &Array1<T>) -> LapackResult<Array1<T>>where T: Field + Clone + Zeroable + Real,
Solves the least squares problem min ||Ax - b||.