function B=gausslu(B);
% Matrix b wird obere und untere Dreiecksmatrix zerlegt, so daß
% B=LU
% U ist auf oberes Dreieck von B gespeichert
% L wird unterhalb der Diagonalen gespeichert 
% Eins-Diagonale von L wird weggelassen
% Roman Unger

% Groesse von B ermitteln
[m,n]=size(B);

for k=1:n-1
    B(k+1:n,k)=B(k+1:n,k)/B(k,k);
    for  j=k+1:n
       for  i=k+1:n
	  B(i,j)=B(i,j)-B(i,k)*B(k,j);
       end;
    end;
end;


