

% a)
disp('Teilaufgabe a)')
A = [2 11 14; 4 2 2; 2 26 11];
b = [5;4;-13];
N=3;
index=1:N;
pivot=1;
[A,b,index] = gaussElim(A,b,N,index,pivot);


% b)

disp('Teilaufgabe b)')
N=2;
index=1:N;
A=ones(2,2);
b=[1;2];
pivot=0;
for l=6:20
    A(1,1)=10^(-l);
    [x,index] = gaussElim(A,b,N,index,pivot);
end

index=1:N;
pivot=1;
for l=6:20
    A(1,1)=10^(-l);
    [A,b,index] = gaussElim(A,b,N,index,pivot);
end


% c)
disp('Teilaufgabe c)')

N = 20;
index=1:N;
for j = 1:N
    for k=1:N
        H(j,k)=1/(j+k-1);
    end
end

x0=ones(N,1);
b=H*x0;
pivot=0;
[A,b,index] = gaussElim(H,b,N,index,pivot);
pivot=1;
[A,b,index] = gaussElim(H,b,N,index,pivot);