Function identity 

pub fn identity<T: Element>(n: usize) -> FerrayResult<Array<T, Ix2>>

Create a 2-D identity matrix of size n x n.

Analogous to numpy.identity().