pub fn det(matrix: &ArrayView2<'_, f64>) -> Result<f64, String>
Real matrix determinant via LU decomposition with partial pivoting. Works for any square NxN matrix (not just 2x2).