Remove Matrix::det() and use Matrix::det_via_lu() as the only method to calculate the determinant

[cm-jones]

Matrix::det_via_lu() is more readable and less computationally complex.