lu_nopivoting.m
33bba04159d6master

lu_nopivoting.m
View Options

			function [L, U] = lu_nopivoting(A)
			[Ni, Nj] = size(A);
			assert(Ni == Nj, "The input must be diagonal");

			N = Nj;
			assert(N > 0, "The input must non empty");

			L = eye(N); % if zeros doesn't give the same result
			U = A; % if zeros doesn't give the same result

			for k=1:(N-1)
			% computing LU
			for i=(k+1):N
			L(i,k) = U(i,k) / U(k,k);
			U(i,:) = U(i,:) - L(i,k) * U(k,:);
			end
			end
			end