function ans = gauss(n, a2)
    for i=1:n-1
        main_elem = a2(i, i);
        for j=i:n+1
            a2(i, j) /= main_elem;
        end

        for j=i+1:n
            curr_elem = a2(j, i);
            for k=i:n+1
                a2(j, k) = (a2(j, k) - a2(i, k)*(curr_elem));
            end
        end
    end
end