function [A, b]=elimnaiv(A, b)
% ein einzelner naiver Eliminationsschritt, S. 30 des Buches
if abs(A(1,1))<eps
    error('Diagonalelement zu klein')
end
% Dies wird der Vektor m des Skripts
m=A(:,1)/A(1,1);  % renormierte erste Spalte
m(1,1)=0;         % vorne eine Null rein
% Aneu=(E-m*e_1')*A=A-m*A(1,:)
A=A-m*A(1,:);     % dann A transformieren
% Hier wird zuviel gemacht, 
% weil die erste Zeile und Spalte von A mittransformiert werden
b=b-m*b(1,1);     % und b transformieren
