slau

function X = Gauss(A, B)
    r = rows(A)
    c = columns(A)
    Ag = [A B]

    if (rank(A) == rank(Ag))
        if (rank(A) == c)
            X = "Единственное решение";
            for k = 1:r
                if (Ag(k,k) == 0)

                    X = "Деление на ноль"
                else

                    Ag(k,:) = Ag(k,:) / Ag(k,k);
                    for i = k+1:c
                        Ag(i,:) = Ag(i,:)  - Ag(k,:)*Ag(i,k);
                    end
                end
            end
            X = Ag(end,end);
            for k = r - 1:-1:1
               X = [Ag(k, c + 1) - sum(Ag(k, c - (length(X) - 1) : c) * X); X];
            end
        else
            X = "Неопределена";
        end
    else
        X = "Нет решений";
    end
end