use cratevec;
use crate;
/// TODO: I cann't figure out why `DMatrix::<f64>::zeros(nrows,ncols)` doesn't work, so init a matrix in this way.
/// TODO: I cann't figure out why `DVector::<f64>::zeros(nrows)` doesn't work, so init a vector in this way.
/// TODO: I cann't figure out why `DVector::<f64>::identity(nrows,ncols)` doesn't work, so init a vector in this way.