working through factored_matrix.py