% Teste die verschiedenen Loeser mit grossen Zufallsmatrizen

A = (2*rand(200)-1)./rand(200);
b = rand(200,1);

x_octave           = A\b;
x_einfach          = gauss_einfach(A,b);
x_zeilenpivot      = gauss_zeilenpivot(A,b);
x_verhaeltnispivot = gauss_verhaeltnispivot(A,b);
x_vollpivot        = gauss_vollpivot(A,b);

% [x_octave x_einfach x_zeilenpivot x_verhaeltnispivot x_vollpivot]

% Jetzt berechnen wir die Residuen

r_octave           = A*x_octave - b;
r_einfach          = A*x_einfach - b;
r_zeilenpivot      = A*x_zeilenpivot - b;
r_verhaeltnispivot = A*x_verhaeltnispivot - b;
r_vollpivot        = A*x_vollpivot - b;

% [r_octave r_einfach r_zeilenpivot r_verhaeltnispivot r_vollpivot]

% Und hier die Norm der Residuen:

[norm(r_octave)
 norm(r_einfach)
 norm(r_zeilenpivot)
 norm(r_verhaeltnispivot)
 norm(r_vollpivot)]
