Function read_matrix 

fn read_matrix(
    reader: &mut Reader<'_>,
    rows: usize,
    columns: usize,
) -> Result<DMatrix<Complex>>