function count_hilbert(n)
    A = hilb(n)
    x = ones(n)
    b = zeros(n)
    b = A * x

    gauss = A \ b
    inverse = inv(A) * b

    gaussError = norm(x - gauss) / norm(x)
    inverseError = norm(x - inverse) / norm(x)
    println("n = ", n, ", rank(A) = ", rank(A), ", cond(A) = ", cond(A), ", gaussError = ", gaussError, ", inverseError = ", inverseError)
    println("B: ", b)
    println("X: ", gauss)
end