matrix_fact
