Function solve

pub fn solve<M: MatrixTrait + LinearAlgebra<M>>(
    A: &M,
    b: &M,
    sk: SolveKind,
) -> M