clear all;
% Gauss-sches Eliminationsverfahren naiv, ohne Pivotisierung
n=5;
A=rand(n,n)    % eine Zufallsmatrix
Asave=A;       % Kopie wird gespeichert 
b=rand(n,1)    % ein Zufallsvektor
bsave=b;       % Kopie wird gespeichert
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)]=...
          elimnaiv(A(k:end,k:end),b(k:end,1))
end
x=rwe(A,b)               % Rueckwaertseinsetzen 
rest=Asave*x-bsave       % Ergebnistest
x=Asave\bsave            % was macht MATLAB? 
restMATLAB=Asave*x-bsave
