added NRmatrix I/O