function returnValue=main
% Setze zunaechst die Konstanten:
n = 9;
max_fehler = zeros(1, 11);
condition = zeros(1, 11);

for n = n:20
    H = hilb(n);
    b = zeros(1, n);
    for i = 1:length(H(1:n))
        b(i) = sum(H(i:n));
    end

    % Einfluss auf diese Methode nicht gegeben ?
    x = H\b';

    m_diff = 0;
    for i = 1:length(x)
        if m_diff < abs(x(i)) -1
            m_diff = abs(x(i)) -1;
        end
    end

    max_fehler(n) = m_diff;
    condition(n) = log(cond(H));
end
%max_fehler = max_fehler(10:length(max_fehler));
%condition = condition(10:length(condition));
plot(max_fehler, condition)
plot(condition, max_fehler)
end