function [L, U] = lu0(A)
	n = size(A);
	L = zeros(n);
	U = A;
	for j = 1:n
		L(j,j) = 1;
		L(j+1:n, j) = U(j+1:n, j)./U(j,j);
%		U(j+1:n, j) = zeros(n-j, 1);
		U(j+1:n, j:n) = U(j+1:n, j:n) - L(j+1:n, j) * U(j, j:n);
	end
%	L
%	A
%	for i = 1:n
%		M(i,1:i-1) = L(i,1:i-1);
%		M(i,i:n) = A(i,i:n);
%	end