function [A] = LU(A)

% function [A] =LU(A)
% LU Zerlegung der Matrix A
%
% Input:  A   Matrix
%        
% Output: A Matrix
%



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