pub fn solve(A: &Matrix, b: &Matrix, sk: SolveKind) -> Matrix