pub fn inv(matrix: &ArrayView2<'_, f64>) -> Result<Array2<f64>, String>Expand description
Real matrix inversion via Gauss-Jordan elimination with partial pivoting. Works for any square NxN matrix (not just 2x2); returns an honest error for non-square or singular/near-singular input.