function [A, b]=elimpivot(A, b)
% ein einzelner Eliminationsschritt, mit Pivotisierung
% durch Zeilenvertauschen
[vmax imax]=max(abs(A(:,1)'))
if vmax<eps
    error('Pivot zu klein')
end
% und jetzt wird vertauscht
helprow=A(1,:);
A(1,:)=A(imax,:);
A(imax,:)=helprow;
help=b(1,1);
b(1,1)=b(imax,1);
b(imax,1)=help;
% und jetzt tun wir so, als waere nichts passiert:
% 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


