function x=gausspivot(A, b)
% Gauss-Elimination mit Zeilenpivotisierung, fuer A*x=b
[m, n]=size(A);
if m~=n 
   error('Matrix ist unsymmetrisch');
end
for k=1:n-1    % Schleife ueber Eliminationsschritte
    % das wirkt auf die Restmatrizen
    eliminationsschritt=k % nur zur Ausgabe, damit man sieht, wo man ist
    % und hier der Standard-Eliminationsschritt
    [A(k:end,k:end), b(k:end,1)]=elimpivot(A(k:end,k:end),b(k:end,1))
end
x=rwe(A,b)               % Rueckwaertseinsetzen 

