Function identity 

pub fn identity<T>() -> Matrix2d<T>
where
    T: Float,

Create an identity matrix.